Deterministic sensitivity analysis for a model for flow in porous media

A deterministic method for sensitivity analysis is developed and applied to a mathematical model for the simulation of flow in porous media. The method is based on the singular value decomposition (SVD) of the Jacobian matrix of the model. It is a local approach to sensitivity analysis providing a hierarchical classification of the directions in both the input space and of those in the output space reflecting the degree of sensitiveness of the latter to the former. Its low computational cost, in comparison with that of statistical approaches, allows the study of the variability of the results of the sensitivity analysis due to the variations of the input parameters of the model, and thus it can provide a quality criterion for the validity of more classical probabilistic global approaches. For the example treated here, however, this variability is weak, and deterministic and statistical methods yield similar sensitivity results.     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.     

OSL signal resetting in young deposits determined with a pulsed photon‐stimulated luminescence (PPSL) unit

Optical Stimulated Luminescence (OSL) is a technique that can be used for dating geological materials deposited within the last half‐million years, including sediments transported by air, water or gravity, as well as rocks heated at high temperatures. Recently, several studies have shown that OSL can also provide information on sediment transport. The pulsed photon‐stimulated luminescence (PPSL) unit (also known as a portable OSL reader) developed by the Scottish Universities Environmental Research Centre is an instrument designed to read luminescence signals from bulk (untreated) sediment samples comprising poly‐mineral and poly‐grain fractions. In this contribution, we evaluate the potential of the PPSL unit to assess the degree of OSL signal resetting in 27 young deposits (<2 ka) transported by different geomorphic agents in volcanic, coastal and fluvial depositional settings located in Mexico. Our results are in agreement with previous findings that used the Risø TL/OSL reader, confirming that sediments transported by debrisflows contain the highest inherited luminesce signals. Infrared stimulation (IRSL) values in volcanic ash, lavas, and sand beach and dune deposits exhibit low scatter. However, with blue stimulation (BLSL) these samples reveal a large degree of scattering, attributed to charge transfer in the case of the coastal deposits and to the low sensitivity of quartz in the case of volcanic material. The luminescence signals of fluvial sediments exhibit a highly scattered distribution in both IRSL and BLSL. We conclude that the use of a PPSL unit is a simple approach to assess the degree of OSL signal resetting in deposits sourced from different geological environments. This research contributes to previous studies that have investigated new applications of the PPSL unit to assist in OSL dating of geological materials.     

Dissolved nitrous oxide (N2O) dynamics in agricultural field drains and headwater streams in an intensive arable catchment

Indirect nitrous oxide (N₂O) emissions produced by nitrogen (N) leaching into surface water and groundwater bodies are poorly understood in comparison to direct N₂O emissions from soils. In this study, dissolved N₂O concentrations were measured weekly in both lowland headwater streams and subsurface agricultural field drain discharges over a 2‐year period (2013–2015) in an intensive arable catchment, Norfolk, UK. All field drain and stream water samples were found to have dissolved N₂O concentrations higher than the water–air equilibrium concentration, illustrating that all sites were acting as a net source of N₂O emissions to the atmosphere. Soil texture was found to significantly influence field drain N₂O dynamics, with mean concentrations from drains in clay loam soils (5.3 μg N L^?1) being greater than drains in sandy loam soils (4.0 μg N L^?1). Soil texture also impacted upon the relationships between field drain N₂O concentrations and other water quality parameters (pH, flow rate, and nitrate (NO₃) and nitrite (NO₂) concentrations), highlighting possible differences in N₂O production mechanisms in different soil types. Catchment antecedent moisture conditions influenced the storm event mobilisation of N₂O in both field drains and streams, with the greatest concentration increases recorded during precipitation events preceded by prolonged wet conditions. N₂O concentrations also varied seasonally, with the lowest mean concentrations typically occurring during the summer months (JJA). Nitrogen fertiliser application rates and different soil inversion regimes were found to have no effect on dissolved N₂O concentrations, whereas higher N₂O concentrations recorded in field drains under a winter cover crop compared to fallow fields revealed cover crops are an ineffective greenhouse gas emission mitigation strategy. Overall, this study highlights the complex interactions governing the dynamics of dissolved N₂O concentrations in field drains and headwater streams in a lowland intensive agricultural catchment.     

Evaluation of Hydraulic Properties of Aquitards Using Earthquake‐Triggered Groundwater Variation

The hydraulic properties of aquitards are not easily obtained because monitoring wells are usually installed in aquifers for groundwater resources management. Earthquake‐induced crust stress (strain) triggers groundwater level variations over a short period of time in a large area. These groundwater anomalies can be used to investigate aquifer systems. This study uses a poroelastic model to fit the postseismic variations of groundwater level triggered by the Chi‐Chi earthquake to evaluate the hydraulic properties of aquitards in the Jhoushuei River alluvial fan (JRAF), Taiwan. Six of the adopted eight wells with depths of 70 to 130 m showed good agreement with the recovery theory. The mean hydraulic conductivities (K) of the aquifers for the eight wells are 1.62 × 10^?4 to 9.06 × 10^?4 m/s, and the thicknesses are 18.8 to 46.1 m. The thicknesses of the aquitards are 11.3 to 42.0 m. Under the isotropic assumption for K, the estimated values of K for the aquitards are 3.0 × 10^?8 to 2.1 × 10^?6 m/s, corresponding to a silty medium. The results match the values obtained for the geological material of the drilling core and those reported in previous studies. The estimated values were combined with those given in previous studies to determine the distribution of K in the first two aquitards in the JRAF. The distribution patterns of the aquitards reflect the sedimentary environments and fit the geological material. The proposed technique can be used to evaluate the K value of aquitards using inverse methods. The inversion results can be used in hydrogeological analyses, contaminant modeling, and subsidence evaluation.     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Bioaccumulation et elimination du 65Zn par Gammarus aequicauda martimov

Uptake of ⁶⁵Zn by Gammarus aequicauda results in a concentration factor of approximately 50 and it is characterised by a maximal accumulation level attained as early as the third day of the experiment. After ingestion by the gammarids of twenty contaminated meals which are distributed over a 45-day period, a transfer of ⁶⁵Zn between the food and the consumer can be observed but there is no sign of the biomagnification phenomenon. The retention rate of the radionuclide is only about 1·5% and the ⁶⁵Zn concentration in the gammarids remains markedly inferior to that of the ingested food.The elimination of the ⁶⁵Zn fixed directly from water takes place according to an exponential model which corresponds to the existence of three biological half-lives of the radionuclide, $\text{Tb}{}_1\text{? 0·5 day, Tb}{}_2\text{? 3·5 days and Tb}{}_3\text{? days}$. The excretion of ⁶⁵Zn accumulated from food is a simple exponential phenomenon which corresponds to only one biological half-life of approximately 17 days.     

Severe soil erosion during a 3‐day exceptional rainfall event: combining modelling and field data for a fallow cereal field

Exceptional rainfall events cause significant losses of soil, although few studies have addressed the validation of model predictions at field scale during severe erosive episodes. In this study, we evaluate the predictive ability of the enhanced Soil Erosion and Redistribution Tool (SERT‐2014) model for mapping and quantifying soil erosion during the exceptional rainfall event (~235 mm) that affected the Central Spanish Pyrenees in October 2012. The capacity of the simulation model is evaluated in a fallow cereal field (1.9 ha) at a high spatial scale (1 × 1 m). Validation was performed with field‐quantified rates of soil loss in the rills and ephemeral gullies and also with a detailed map of soil redistribution. The SERT‐2014 model was run for the six rainfall sub‐events that made up the exceptional event, simulating the different hydrological responses of soils with maximum runoff depths ranging between 40 and 1017 mm. Predicted average and maximum soil erosion was 11 and 117 Mg ha^?1 event^?1, respectively. Total soil loss and sediment yield to the La Reina gully amounted to 16.3 and 9.0 Mg event^?1. These rates are in agreement with field estimations of soil loss of 20.0 Mg event^?1. Most soil loss (86%) occurred during the first sub‐event. Although soil accumulation was overestimated in the first sub‐event because of the large amount of detached soil, the enhanced SERT‐2014 model successfully predicted the different spatial patterns and values of soil redistribution for each sub‐event. Further research should focus on stream transport capacity. Copyright © 2014 John Wiley & Sons, Ltd.