The domestication of water: water management in the ancient world and its prehistoric origins in the Jordan Valley

The ancient civilizations were dependent upon sophisticated systems of water management. The hydraulic engineering works found in ancient Angkor (ninth to thirteenth century AD), the Aztec city of Tenochtitlan (thirteenth to fifteenth century AD), Byzantine Constantinople (fourth to sixth century AD) and Nabatean Petra (sixth century BC to AD 106) are particularly striking because each of these is in localities of the world that are once again facing a water crisis. Without water management, such ancient cities would never have emerged, nor would the urban communities and towns from which they developed. Indeed, the ‘domestication’ of water marked a key turning point in the cultural trajectory of each region of the world where state societies developed. This is illustrated by examining the prehistory of water management in the Jordan Valley, identifying the later Neolithic (approx. 8300–6500 years ago) as a key period when significant investment in water management occurred, laying the foundation for the development of the first urban communities of the Early Bronze Age.

A model of the three-dimensional evolution of Arctic melt ponds on first-year and multiyear sea ice

During winter the ocean surface in polar regions freezes over to form sea ice. In the summer the upper layers of sea ice and snow melts producing meltwater that accumulates in Arctic melt ponds on the surface of sea ice. An accurate estimate of the fraction of the sea ice surface covered in melt ponds is essential for a realistic estimate of the albedo for global climate models. We present a melt-pond–sea-ice model that simulates the three-dimensional evolution of melt ponds on an Arctic sea ice surface. The advancements of this model compared to previous models are the inclusion of snow topography; meltwater transport rates are calculated from hydraulic gradients and ice permeability; and the incorporation of a detailed one-dimensional, thermodynamic radiative balance. Results of model runs simulating first-year and multiyear sea ice are presented. Model results show good agreement with observations, with duration of pond coverage, pond area, and ice ablation comparing well for both the first-year ice and multiyear ice cases. We investigate the sensitivity of the melt pond cover to changes in ice topography, snow topography, and vertical ice permeability. Snow was found to have an important impact mainly at the start of the melt season, whereas initial ice topography strongly controlled pond size and pond fraction throughout the melt season. A reduction in ice permeability allowed surface flooding of relatively flat, first-year ice but had little impact on the pond coverage of rougher, multiyear ice. We discuss our results, including model shortcomings and areas of experimental uncertainty.

Groundwater residence time distributions in peatlands: implications for peat decomposition and accumulation

Peat soils consist of poorly decomposed plant detritus, preserved by low decay rates, and deep peat deposits are globally significant stores in the carbon cycle. High water tables and low soil temperatures are commonly held to be the primary reasons for low peat decay rates. However, recent studies suggest a thermodynamic limit to peat decay, whereby the slow turnover of peat soil pore water may lead to high concentrations of phenols and dissolved inorganic carbon. In sufficient concentrations, these chemicals may slow or even halt microbial respiration, providing a negative feedback to peat decay. We document the analysis of a simple, one-dimensional theoretical model of peatland pore water residence time distributions (RTDs). The model suggests that broader, thicker peatlands may be more resilient to rapid decay caused by climate change because of slow pore water turnover in deep layers. Even shallow peat deposits may also be resilient to rapid decay if rainfall rates are low. However, the model suggests that even thick peatlands may be vulnerable to rapid decay under prolonged high rainfall rates, which may act to flush pore water with fresh rainwater. We also used the model to illustrate a particular limitation of the diplotelmic (i.e., acrotelm and catotelm) model of peatland structure. Model peatlands of contrasting hydraulic structure exhibited identical water tables but contrasting RTDs. These scenarios would be treated identically by diplotelmic models, although the thermodynamic limit suggests contrasting decay regimes. We therefore conclude that the diplotelmic model be discarded in favor of model schemes that consider continuous variation in peat properties and processes.

On the Nusselt number for frazil ice growth—a correction to ”Frazil evolution in channels“ by Lars Hammar and Hung-Tao Shen

The growth (melt) rate of frazil ice is governed by heat transfer away from (towards) the ice crystal, which can be represented by the Nusselt number. We discuss choices for the Nusselt number and turbulent length scale appropriate for frazil ice and note an inaccuracy in the study ”Frazil evolution in channels“ by Lars Hammar and Hung-Tao Shen, which has also led to potentially significant errors in several other papers. We correct this error and suggest an appropriate strategy for determining the Nusselt number applicable to frazil ice growth and melting.

Computer simulations of pyroclastic flows from dome collapse

We present a method of simulating both the avalanche and surge components of pyroclastic flows generated by lava collapsing from a growing Pelean dome. This is used to successfully model the pyroclastic flows generated on 12 May 1996 by the Soufriere Hills volcano, Montserrat. In simulating the avalanche component we use a simple 3-fold parameterisation of flow acceleration for which we choose values using an inverse method. The surge component is simulated by a 1D hydraulic balance of sedimentation of clasts and entrainment of air away from the avalanche source. We show how multiple simulations based on uncertainty of the starting conditions and parameters, specifically location and size (mass flux), could be used to map hazard zones.

The influence of ocean flow on newly forming sea ice

The heat and mass balance of the Arctic Ocean is very sensitive to the growth and decay of sea ice and the interaction between the heat and salt fields in the oceanic boundary layer. The hydraulic roughness of sea ice controls the detailed nature of turbulent fluxes in the boundary layer and hence is an important ingredient in model parameterizations. We describe a novel mechanism for the generation of corrugations of the sea ice–ocean interface, present a mathematical analysis elucidating the mechanism, and present numerical calculations for geophysically relevant conditions. The mechanism relies on brine flows developing in the sea ice due to Bernoulli suction by flow of ocean past the interface. For oceanic shears at the ice interface of 0.2 s−1, we expect the corrugations to form with a wavelength dependent upon the permeability structure of the sea ice which is described herein. The mechanism should be particularly important during sea ice formation in wind-maintained coastal polynyas and in leads. This paper applies our earlier analyses of the fundamental instability to field conditions and extends it to take account of the anisotropic and heterogeneous permeability of sea ice.

The impact of uncertain precipitation data on insurance loss estimates using a flood catastrophe model

Catastrophe risk models used by the insurance industry are likely subject to significant uncertainty, but due to their proprietary nature and strict licensing conditions they are not available for experimentation. In addition, even if such experiments were conducted, these would not be repeatable by other researchers because commercial confidentiality issues prevent the details of proprietary catastrophe model structures from being described in public domain documents. However, such experimentation is urgently required to improve decision making in both insurance and reinsurance markets. In this paper we therefore construct our own catastrophe risk model for flooding in Dublin, Ireland, in order to assess the impact of typical precipitation data uncertainty on loss predictions. As we consider only a city region rather than a whole territory and have access to detailed data and computing resources typically unavailable to industry modellers, our model is significantly more detailed than most commercial products. The model consists of four components, a stochastic rainfall module, a hydrological and hydraulic flood hazard module, a vulnerability module, and a financial loss module. Using these we undertake a series of simulations to test the impact of driving the stochastic event generator with four different rainfall data sets: ground gauge data, gauge-corrected rainfall radar, meteorological reanalysis data (European Centre for Medium-Range Weather Forecasts Reanalysis-Interim; ERA-Interim) and a satellite rainfall product (The Climate Prediction Center morphing method; CMORPH). Catastrophe models are unusual because they use the upper three components of the modelling chain to generate a large synthetic database of unobserved and severe loss-driving events for which estimated losses are calculated. We find the loss estimates to be more sensitive to uncertainties propagated from the driving precipitation data sets than to other uncertainties in the hazard and vulnerability modules, suggesting that the range of uncertainty within catastrophe model structures may be greater than commonly believed.

Drought and resprouting plants

Many species have the ability to resprout vegetatively after a substantial loss of biomass induced by environmental stress, including drought. Many of the regions characterised by ecosystems where resprouting is common are projected to experience more frequent and intense drought during the 21st Century. However, in assessments of ecosystem response to drought disturbance there has been scant consideration of the resilience and post-drought recovery of resprouting species. Systematic differences in hydraulic and allocation traits suggest that resprouting species are more resilient to drought-stress than nonresprouting species. Evidence suggests that ecosystems dominated by resprouters recover from disturbance more quickly than ecosystems dominated by nonresprouters. The ability of resprouters to avoid mortality and withstand drought, coupled with their ability to recover rapidly, suggests that the impact of increased drought stress in ecosystems dominated by these species may be small. The strategy of resprouting needs to be modelled explicitly to improve estimates of future climate-change impacts on the carbon cycle, but this will require several important knowledge gaps to be filled before resprouting can be properly implemented.

New-on-line method for water isotope analysis of speleothem fluid inclusions using laser absorption spectroscopy (WS-CRDS)

A new online method to analyse water isotopes of speleothem fluid inclusions using a wavelength scanned cavity ring down spectroscopy (WS-CRDS) instrument is presented. This novel technique allows us simultaneously to measure hydrogen and oxygen isotopes for a released aliquot of water. To do so, we designed a new simple line that allows the online water extraction and isotope analysis of speleothem samples. The specificity of the method lies in the fact that fluid inclusions release is made on a standard water background, which mainly improves the δ D robustness. To saturate the line, a peristaltic pump continuously injects standard water into the line that is permanently heated to 140 °C and flushed with dry nitrogen gas. This permits instantaneous and complete vaporisation of the standard water, resulting in an artificial water background with well-known δ D and δ18O values. The speleothem sample is placed in a copper tube, attached to the line, and after system stabilisation it is crushed using a simple hydraulic device to liberate speleothem fluid inclusions water. The released water is carried by the nitrogen/standard water gas stream directly to a Picarro L1102-i for isotope determination. To test the accuracy and reproducibility of the line and to measure standard water during speleothem measurements, a syringe injection unit was added to the line. Peak evaluation is done similarly as in gas chromatography to obtain &delta D; and δ18O isotopic compositions of measured water aliquots. Precision is better than 1.5 ‰ for δ D and 0.4 ‰ for δ18O for water measurements for an extended range (−210 to 0 ‰ for δ D and −27 to 0 ‰ for δ18O) primarily dependent on the amount of water released from speleothem fluid inclusions and secondarily on the isotopic composition of the sample. The results show that WS-CRDS technology is suitable for speleothem fluid inclusion measurements and gives results that are comparable to the isotope ratio mass spectrometry (IRMS) technique.

Tracking the uncertainty in flood alerts driven by grand ensemble weather predictions

The incorporation of numerical weather predictions (NWP) into a flood warning system can increase forecast lead times from a few hours to a few days. A single NWP forecast from a single forecast centre, however, is insufficient as it involves considerable non-predictable uncertainties and can lead to a high number of false or missed warnings. Weather forecasts using multiple NWPs from various weather centres implemented on catchment hydrology can provide significantly improved early flood warning. The availability of global ensemble weather prediction systems through the ‘THORPEX Interactive Grand Global Ensemble’ (TIGGE) offers a new opportunity for the development of state-of-the-art early flood forecasting systems. This paper presents a case study using the TIGGE database for flood warning on a meso-scale catchment (4062 km2) located in the Midlands region of England. For the first time, a research attempt is made to set up a coupled atmospheric-hydrologic-hydraulic cascade system driven by the TIGGE ensemble forecasts. A probabilistic discharge and flood inundation forecast is provided as the end product to study the potential benefits of using the TIGGE database. The study shows that precipitation input uncertainties dominate and propagate through the cascade chain. The current NWPs fall short of representing the spatial precipitation variability on such a comparatively small catchment, which indicates need to improve NWPs resolution and/or disaggregating techniques to narrow down the spatial gap between meteorology and hydrology. The spread of discharge forecasts varies from centre to centre, but it is generally large and implies a significant level of uncertainties. Nevertheless, the results show the TIGGE database is a promising tool to forecast flood inundation, comparable with that driven by raingauge observation.

Efficient incorporation of channel cross-section geometry uncertainty into regional and global scale flood inundation models

This paper investigates the challenge of representing structural differences in river channel cross-section geometry for regional to global scale river hydraulic models and the effect this can have on simulations of wave dynamics. Classically, channel geometry is defined using data, yet at larger scales the necessary information and model structures do not exist to take this approach. We therefore propose a fundamentally different approach where the structural uncertainty in channel geometry is represented using a simple parameterization, which could then be estimated through calibration or data assimilation. This paper first outlines the development of a computationally efficient numerical scheme to represent generalised channel shapes using a single parameter, which is then validated using a simple straight channel test case and shown to predict wetted perimeter to within 2% for the channels tested. An application to the River Severn, UK is also presented, along with an analysis of model sensitivity to channel shape, depth and friction. The channel shape parameter was shown to improve model simulations of river level, particularly for more physically plausible channel roughness and depth parameter ranges. Calibrating channel Manning’s coefficient in a rectangular channel provided similar water level simulation accuracy in terms of Nash-Sutcliffe efficiency to a model where friction and shape or depth were calibrated. However, the calibrated Manning coefficient in the rectangular channel model was ~2/3 greater than the likely physically realistic value for this reach and this erroneously slowed wave propagation times through the reach by several hours. Therefore, for large scale models applied in data sparse areas, calibrating channel depth and/or shape may be preferable to assuming a rectangular geometry and calibrating friction alone.

Geotechnical systems that evolve with ecological processes

Geotechnical systems, such as landfills, mine tailings storage facilities (TSFs), slopes, and levees, are required to perform safely throughout their service life, which can span from decades for levees to “in perpetuity” for TSFs. The conventional design practice by geotechnical engineers for these systems utilizes the as-built material properties to predict its performance throughout the required service life. The implicit assumption in this design methodology is that the soil properties are stable through time. This is counter to long-term field observations of these systems, particularly where ecological processes such as plant, animal, biological, and geochemical activity are present. Plant roots can densify soil and/or increase hydraulic conductivity, burrowing animals can increase seepage, biological activity can strengthen soil, geochemical processes can increase stiffness, etc. The engineering soil properties naturally change as a stable ecological system is gradually established following initial construction, and these changes alter system performance. This paper presents an integrated perspective and new approach to this issue, considering ecological, geotechnical, and mining demands and constraints. A series of data sets and case histories are utilized to examine these issues and to propose a more integrated design approach, and consideration is given to future opportunities to manage engineered landscapes as ecological systems. We conclude that soil scientists and restoration ecologists must be engaged in initial project design and geotechnical engineers must be active in long-term management during the facility’s service life. For near-surface geotechnical structures in particular, this requires an interdisciplinary perspective and the embracing of soil as a living ecological system rather than an inert construction material.

Assessing the reliability of probabilistic flood inundation model predictions

An ability to quantify the reliability of probabilistic flood inundation predictions is a requirement not only for guiding model development but also for their successful application. Probabilistic flood inundation predictions are usually produced by choosing a method of weighting the model parameter space, but previous study suggests that this choice leads to clear differences in inundation probabilities. This study aims to address the evaluation of the reliability of these probabilistic predictions. However, a lack of an adequate number of observations of flood inundation for a catchment limits the application of conventional methods of evaluating predictive reliability. Consequently, attempts have been made to assess the reliability of probabilistic predictions using multiple observations from a single flood event. Here, a LISFLOOD-FP hydraulic model of an extreme (>1 in 1000 years) flood event in Cockermouth, UK, is constructed and calibrated using multiple performance measures from both peak flood wrack mark data and aerial photography captured post-peak. These measures are used in weighting the parameter space to produce multiple probabilistic predictions for the event. Two methods of assessing the reliability of these probabilistic predictions using limited observations are utilized; an existing method assessing the binary pattern of flooding, and a method developed in this paper to assess predictions of water surface elevation. This study finds that the water surface elevation method has both a better diagnostic and discriminatory ability, but this result is likely to be sensitive to the unknown uncertainties in the upstream boundary condition

Effect of temperature and strain rate on the compressive behaviour of supramolecular polyurethane

Supramolecular polyurethanes (SPUs) possess thermoresponsive and thermoreversible properties, and those characteristics are highly desirable in both bulk commodity and value-added applications such as adhesives, shape-memory materials, healable coatings and lightweight, impact-resistant structures (e.g. protection for mobile electronics). A better understanding of the mechanical properties, especially the rate and temperature sensitivity, of these materials are required to assess their suitability for different applications. In this paper, a newly developed SPU with tuneable thermal properties was studied, and the response of this SPU to compressive loading over strain rates from 10−3 to 104 s−1 was presented. Furthermore, the effect of temperature on the mechanical response was also demonstrated. The sample was tested using an Instron mechanical testing machine for quasi-static loading, a home-made hydraulic system for moderate rates and a traditional split Hopkinson pressure bars (SHPBs) for high strain rates. Results showed that the compression stress-strain behaviour was affected significantly by the thermoresponsive nature of SPU, but that, as expected for polymeric materials, the general trends of the temperature and the rate dependence mirror each other. However, this behaviour is more complicated than observed for many other polymeric materials, as a result of the richer range of transitions that influence the behaviour over the range of temperatures and strain rates tested.

Effect of enzymatic pretreatment and increasing the organic loading rate of lipid-rich wastewater treated in a hybrid UASB reactor

This study aimed at evaluating the effect of increasing organic loading rates and of enzyme pretreatment on the stability and efficiency of a hybrid upflow anaerobic sludge blanket reactor (UASBh) treating dairy effluent. The UASBh was submitted to the following average organic loading rates (OLR) 0.98 Kg.m(-3).d(-1), 4.58 Kg.m(-3).d(-1), 8.89 Kg.m(-3).d(-1) and 15.73 Kg.m(-3).d(-1), and with the higher value, the reactor was fed with effluent with and without an enzymatic pretreatment to hydrolyze fats. The hydraulic detention time was 24 h, and the temperature was 30 +/- 2 degrees C. The reactor was equipped with a superior foam bed and showed good efficiency and stability until an OLR of 8.89 Kg.m(-3).d(-1). The foam bed was efficient for solid retention and residual volatile acid concentration consumption. The enzymatic pretreatment did not contribute to the process stability, propitiating loss in both biomass and system efficiency. Specific methanogenic activity tests indicated the presence of inhibition after the sludge had been submitted to the pretreated effluent It was concluded that continuous exposure to the hydrolysis products or to the enzyme caused a dramatic drop in the efficiency and stability of the process, and the single exposure of the biomass to this condition did not inhibit methane formation. (C) 2011 Elsevier B.V. All rights reserved.