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.

Seasonal drought stress in the Amazon: Reconciling models and observations

The Amazon Basin is crucial to global circulatory and carbon patterns due to the large areal extent and large flux magnitude. Biogeophysical models have had difficulty reproducing the annual cycle of net ecosystem exchange (NEE) of carbon in some regions of the Amazon, generally simulating uptake during the wet season and efflux during seasonal drought. In reality, the opposite occurs. Observational and modeling studies have identified several mechanisms that explain the observed annual cycle, including: (1) deep soil columns that can store large water amount, (2) the ability of deep roots to access moisture at depth when near-surface soil dries during annual drought, (3) movement of water in the soil via hydraulic redistribution, allowing for more efficient uptake of water during the wet season, and moistening of near-surface soil during the annual drought, and (4) photosynthetic response to elevated light levels as cloudiness decreases during the dry season. We incorporate these mechanisms into the third version of the Simple Biosphere model (SiB3) both singly and collectively, and confront the results with observations. For the forest to maintain function through seasonal drought, there must be sufficient water storage in the soil to sustain transpiration through the dry season in addition to the ability of the roots to access the stored water. We find that individually, none of these mechanisms by themselves produces a simulation of the annual cycle of NEE that matches the observed. When these mechanisms are combined into the model, NEE follows the general trend of the observations, showing efflux during the wet season and uptake during seasonal drought.

Seismic activity triggered by water wells in the Parana Basin, Brazil

Triggered seismicity is commonly associated with deep water reservoirs or injection wells where water is injected at high pressure into the reservoir rock. However, earth tremors related solely to the opening of groundwater wells are extremely rare. Here we present a clear case of seismicity induced by pore-pressure changes following the drilling of water wells that exploit a confined aquifer in the intracratonic Parana Basin of southeastern Brazil. Since 2004, shallow seismic activity, with magnitudes up to 2.9 and intensities V MM, has been observed near deep wells (120-200 m) that were drilled in early 2003 near the town of Bebedouro. The wells were drilled for irrigation purposes, cross a sandstone layer about 60-80 m thick and extract water from a confined aquifer in fractured zones between basalt flow layers. Seismic activity, mainly event swarms, has occurred yearly since 2004, mostly during the rainy season when the wells are not pumped. During the dry season when the wells are pumped almost continuously, the activity is very low. A seismographic network, installed in March 2005, has located more than 2000 microearthquakes. The events are less than 1 km deep (mostly within the 0.5 km thick basalt layer) and cover an area roughly 1.5 km x 5 km across. The seismicity generally starts in a small area and expands to larger distances with an equivalent hydraulic diffusivity ranging from 0.06 to 0.6 m(2)/s. Geophysical and geothermal logging of several wells in the area showed that water from the shallow sandstone aquifer enters the well at the top and usually forms waterfalls. The waterfalls flow down the sides of the wells and feed the confined, fractured aquifer in the basalt layer at the bottom. Two seismic areas are observed: the main area surrounds several wells that are pumped continuously during the dry season, and a second area near another well (about 10 km from the first area) that is not used for irrigation and not pumped regularly. The main area displays cyclic annual activity, but the second area does not. We explain the earthquake swarms as being triggered by pore pressure diffusion in the fractured basalt layer due to additional pressure from the newly connected surface aquifer. This reaches critically prestressed areas up to a few kilometers away from the wells. During periods of continuous pumping, the reduction of pore pressure in the confined aquifer stops the seismic activity. Our study suggests that this kind of activity may be more common than previously thought and implies that many other cases of small tremors associated with the drilling of water wells may have gone unnoticed.

THE ENDOPARASITE PILOSTYLES ULEI (APODANTHACEAE - CUCURBITALES) INFLUENCES WOOD STRUCTURE IN THREE HOST SPECIES OF MIMOSA

Pilostyles species (Apodanthaceae) are endoparasites in stems of the plant family Fabaceae. The body comprises masses of parenchyma in the host bark and cortex, with sinkers, comprising groups of twisted tracheal elements surrounded by parenchyma that enter the secondary xylem of the host plant. Here we report for the first time the effects of Pilostyles parasitism on host secondary xylem. We obtained healthy and parasitized stems from Mimosa foliolosa, M. maguirei and M. setosa and compared vessel element length, fiber length, vessel diameter and vessel frequency, measured through digital imaging. Also, tree height and girth were compared between healthy and parasitized M. setosa. When parasitized, plant size, vessel diameter, vessel element length and fiber length are all less than in healthy plants. Also, vessel frequency is greater and vessels are narrower in parasitized stems. These responses to parasitism are similar to those observed in stressed plants. Thus, hosts respond to the parasite by changing its wood micromorphology in favour of increased hydraulic safety.

Evaluation of a bounded high order upwind scheme for 3D incompressible free surface flow computations

fit the context of normalized variable formulation (NVF) of Leonard and total variation diminishing (TVD) constraints of Harten. this paper presents an extension of it previous work by the authors for solving unsteady incompressible flow problems. The main contributions of the paper are threefold. First, it presents the results of the development and implementation of a bounded high order upwind adaptative QUICKEST scheme in the 3D robust code (Freeflow), for the numerical solution of the full incompressible Navier-Stokes equations. Second, it reports numerical simulation results for 1D hock tube problem, 2D impinging jet and 2D/3D broken clam flows. Furthermore, these results are compared with existing analytical and experimental data. And third, it presents the application of the numerical method for solving 3D free surface flow problems. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved,