A field,laboratory, and literature review evaluation of the water retention curve of volcanic ash soils: How well do standard laboratory methods reflect field conditions?

Accurate determination of the water retention curve (WRC) of a soil is essential for the understanding and modelling of the subsurface hydrological, ecological, and biogeochemical processes. Volcanic ash soils with andic properties (Andosols) are recognized as important providers of ecological and hydrological services in mountainous regions worldwide due to their large fraction of small size particles (clay, silt, and organic matter) that gives them an outstanding water holding capacity. Previous comparative analyses of in situ (field) and standard laboratory methods for the determination of the WRC of Andosols showed contrasting results. Based on an extensive analysis of laboratory, experimental, and field measured WRCs of Andosols in combination with data extracted from the published literature we show that standard laboratory methods using small soil sample volumes (≤300 cm³) mimic the WRC of these soils only partially. The results obtained by the latter resemble only a small portion of the wet range of the Andosols' WRC (from saturation up to −5 kPa, or pF 1.7), but overestimate substantially their water content for higher matric potentials. This discrepancy occurs irrespective of site-specific land use and cover, soil properties, and applied method. The disagreement limits our capacity to infer correctly subsurface hydrological behaviour, as illustrated through the analysis of long-term soil moisture and matric potential data from an experimental site in the tropical Andes. These findings imply that results reported in past research should be used with caution and that future research should focus on determining laboratory methods that allow obtaining a correct characterization of the WRC of Andosols. For the latter, a set of recommendations and future directions to solve the identified methodological issues is proposed.     

Magnetic signature of hydrocarbon-contaminated soils and sediments at the former oil field H?nigsen, Germany

Magnetic properties of hydrocarbon (HC) containing soils and sediments from two sites (Site A and B) of the former oil-field H?nigsen were analyzed in order to determine whether magnetic methods can be employed to delineate HC contamination of soils and sediments. Magnetic parameters such as magnetic susceptibility and induced isothermal remanent magnetizations, as well as soil and sediment properties such as pH, iron content and water content, HC content and most probable number counts of iron-metabolizing microorganisms were determined. The magnetic concentration-dependent parameters for HC contaminated samples were 25 times higher in soils from Site A than in sediment samples from Site B. However, at Site B the magnetic susceptibility was still four times higher in comparison to lithologically similar non-contaminated sediment samples from a third Site C. Newly formed magnetite containing mainly single domain particles was responsible for the magnetic enhancement, whereas superparamagnetic grains represented only a minor component. Site A had an acidic pH compared to neutral pH at Site B, and a higher crystalline and bioavailable total iron content. Nevertheless, Site B samples contained significant numbers of both iron(II)-oxidizing and iron(III)-reducing microorganisms indicating that microbial iron cycling might have taken place at this site and potentially played a role for iron mineral transformation, including magnetite (trans)formation. The content of total non-polar hydrocarbons (TNPH) at Site A was one order of magnitude higher than at Site B. Only at Site A magnetic susceptibility correlated well with TNPH. Our results demonstrate that HC contaminated samples had an enhanced magnetite content compared to non-contaminated soils and sediments. Therefore, magnetic methods may provide a quick and cost-effective way to assess HC contamination in soils and sediments. However, more field sites and laboratory investigations are needed to reveal the complex nature of the processes involved.     

Does pore water velocity affect the reaction rates of adsorptive solute transport in soils? Demonstration with pore-scale modelling

Modelling adsorptive solute transport in soils needs a number of parameters to describe its reaction kinetics and the values of these parameters are usually determined from batch and displacement experiments. Some experimental results reveal that when describing the adsorption as first-order kinetics, its associated reaction rates are not constants but vary with pore water velocity. Explanation of this varies but an independent verification of each explanation is difficult because simultaneously measuring the spatiotemporal distributions of dissolved and adsorbed solutes in soils is formidable. Pore-scale modelling could play an important role to address this gap and has received increased attention over the past few years. This paper investigated the transport of adsorptive solute in a simple porous medium using pore-scale modelling. Fluid flow through the void space of the medium was assumed to be laminar and in saturated condition, and solute transport consisted of advection and molecular diffusion; the sorption and desorption occurring at the fluid–solid interface were modelled as linear first-order kinetics. Based on the simulated spatiotemporal distribution of dissolved and adsorbed solutes at pore scale, volumetric-average reaction kinetics at macroscopic scale and its associated reactive parameters were measured. Both homogeneous adsorption where the reaction rates at microscopic scale are constant, and heterogeneous adsorption where the reaction rates vary from site to site, were investigated. The results indicate that, in contrast to previously thought, the macroscopic reaction rates directly measured from the pore-scale simulations do not change with pore velocity under both homogeneous and heterogeneous adsorptions. In particular, we found that for the homogeneous adsorption, the macroscopic adsorption remains first-order kinetic and can be described by constant reaction rates, regardless of flow rate; whilst for the heterogeneous adsorption, the macroscopic adsorption kinetics continues not to be affected by flow rate but is no longer first-order kinetics that can be described by constant reaction rates. We discuss how these findings could help explain some contrary literature reports over the dependence of reaction rates on pore water velocity.     

Are spatial patterns of soil moisture at plot scales generalisable across catchments,climates, and other characteristics? A synthesis of synoptic soil moisture across the Mid-Atlantic

The spatial variation of soil moisture over very small areas (<100 m²) can have nonlinear impacts on cycling and flux rates resulting in bias if it is not considered, but measuring this variation is difficult over extensive temporal and spatial scales. Most studies examining spatial variation of soil moisture were conducted at hillslope (0.01 km²) to multi-catchment spatial scales (1000 km²). They found the greatest variation at mid wetness levels and the smallest variation at wet and dry wetness levels forming a concave down relationship. There is growing evidence that concave down relationships formed between spatial variation of soil moisture and average soil moisture are consistent across spatial scales spanning several orders of magnitude, but more research is needed at very small, plot scales (<100 m²). The goal of this study was to characterise spatial variation in shallow soil moisture at the plot scale by relating the mean of measurements collected in a plot to the standard deviation (SD). We combined data from a previous study with thousands of new soil moisture measurements from 212 plots in eight catchments distributed across the US Mid-Atlantic Region to (1) test for a generalisable mean–SD relationship at plot scales, (2) characterise how landcover, land use, season, and hillslope position contribute to differences in mean–SD relationships, and (3) use these generalised mean–SD relationships to quantify their impacts on catchment scale nitrification and denitrification potential. Our study found that 98% of all measurements formed a generalised mean–SD relationship like those observed at hillslope and catchment spatial scales. The remaining 2% of data comprised a mean–SD relationship with greater spatial variation that originated from two riparian plots reported in a previous study. Incorporating the generalised mean–SD relationship into estimates of nitrification and denitrification potential revealed strong bias that was even greater when incorporating mean–SD observations from the two riparian plots with significantly greater spatial variation.     

What factors drive seasonal variation of phytoplankton,protozoans and metazoans on leaves of Posidonia oceanica and in the water column along the coast of the Kerkennah Islands,Tunisia?

A hierarchical sampling design was used during two seasons (spring (May) and summer (August) 2006). Using this design, three regions of the Kerkennah Islands (Tunisia) were analyzed for the distribution of microalgal, protozoan and metazoan assemblages in two different habitats: (1) the water column; and (2) on Posidonia oceanica (L.) Delile (P. oceanica) leaves in shallow meadows. A total of 85 species were obtained. In particular, the diatom family Naviculacea consistently dominated (both numerically and in their diversity) the micro-algae in all regions for the two seasons of the study and in both habitats. In the Chergui region, which is the closest area to a source of impact, fast growing centric diatoms (such as Thalassionema, Rhizosolenia, Striatella, and Skeletonema) were identified as indicators of high organic matter and nutrient enrichment in water bodies. Protozoan and metazoan species abundance in the different regions indicate a non-random spatial and temporal distribution of the epiphytic organisms on leaves of P. oceanica that correlated with phytoplankton. The results also indicate that (1) the abundance of micro- and macroorganisms in the three regions were higher on P. oceanica leaves than in the water column for the two seasons; (2) environmental factors such as currents and tide influenced assemblages; and (3) the highest abundance was due to direct exposure to the polluted coast of Sfax and the effect of tidal asymmetries generating nutrient-rich inputs from the city.     

What is the source of baseflow in agriculturally fragmented catchments? Complex groundwater/surface‐water interactions in three tributary catchments of the Wabash River,Indiana, USA

Some conceptual models suggest that baseflow in agriculturally fragmented watersheds may contain little, if any, groundwater. This has critical implications for stream quality and ecosystem functioning. Here, we (a) identify the sources and flowpaths contributing to baseflow using ²²²Rn and ⁸⁷Sr/⁸⁶Sr and (b) quantify mean apparent ages of groundwater and baseflow using multiple isotopic tracers (CFC, SF₆, ³⁶Cl, and ³H) in 4 small (0.08 to 0.64 km²) tributary catchments to the Wabash River in Indiana, USA. ²²²Rn activities and ⁸⁷Sr/⁸⁶Sr ratios indicate that baseflow in 3 catchments is sourced primarily from groundwater; baseflow in the fourth is dominated by a source similar to agricultural run‐off. CFC‐12 data indicate that springs in 1 catchment are discharging significant proportions of water that recharged between 1974 (42 ± 2 years) and 1961 (55 ± 2 years). Those same springs have ³⁶Cl/Cl ratios between 1,381.08 ± 29.37 (×10^?15) and 1,530.64 ± 27.65 (×10^?15) indicating that a substantial proportion of the discharge likely recharged between 1975 (41 years) and 1950 (66 years). Groundwater samples collected from streambed mini‐piezometers in a separate catchment have CFC‐12 concentrations indicating that a large proportion of the recharge occurred between 1948 (68 ± 2 years) and 1950 (66 ± 2 years). Repeat sampling conducted in September 2015 after above‐average summer rainfall did not show significant decreases in mean apparent age. The relatively old ages observed in 3 of the catchments can be explained by geological complexities that are likely present in all 4 catchments, but overwhelmed by flow from the shallow phreatic aquifer in the fourth catchment.     

Do model parameters change under changing climate and land use in the upstream of the Lancang River Basin,China?

ABSTRACT

Hydrological model parameters may vary under the impacts of climate and land-use change. This study proposes a hydrological modelling framework based on the Soil and Water Assessment Tool (SWAT) to investigate the variability of model parameters in three different experiments and to assess the impacts of climate and/or land-use change on these parameters in the upstream of the Lancang River Basin, China. In Experiment 1 (E1), most parameters show clear temporal trends under changing climate and land use, implying that model parameters are strongly influenced by their combined effects. Experiments 2 (E2) and 3 (E3) investigate the separate impacts of land-use change and climate change, respectively. Due to the almost invisible changes in land use in E2, there is no change detected in the model parameters. Temporal trends are found in most parameters in E3 and over half of them show consistent trends with E1, which indicates that climate change has greater impacts on model parameter variability. The simulated extreme streamflow and sediment fluxes vary substantially with time-variant parameters, implying that the variations in model parameters do matter for hydrological prediction.     

Mixed layer depth (MLD) variability in the southern Bay of Biscay. Deepening of winter MLDs concurrent with generalized upper water warming trends?

Mixed layer depth (MLD) variability from seasonal to decadal time scales in the Bay of Biscay is studied in this work. A hydrographic time series running since 1991 in the study area, a climatology of the upper layer vertical structure based on the topology of this temperature profile time series and a one-dimensional water column model have been used for this purpose. The prevailing factors driving MLD variability have been determined with detail, and agreement with observations is achieved. Tests carried out to investigate climatological profile skill to reproduce the upper layer temporal evolution have demonstrated its ability to simulate variability at seasonal time scales and reproduce the most conspicuous events observed. This has enabled us to carry out a reconstruction of the MLD variability for the last 60 years in the study area. Favourable sequence of intense mixing events explains interannual differences and cases of extraordinary deepening of winter mixed layer. The negative phase of the Eastern Atlantic pattern seems to determine important interannual variability through intense episodes of cooling and mixing as in winter 2005 in the Bay of Biscay. Low-frequency variability is also observed. A very striking and unexpected shallower winter MLD during the 1970s and 1980s than those observed from 1995 has been found. Simulation results support this counter-intuitive outcome of shallower winter mixed layers concurrent with generalized upper water warming trends reported on several occasions for the area. The long-term trends in MLD seem related with decadal variability in the North Atlantic Oscillation, being in phase and opposition with other deepening-shallowing cycles found from subtropical-to-subpolar areas in the North Atlantic.     

Sublethal effects of a chlorinated and heated effluent on the physiology of the mussel,Mytilus edulis L.: A reduction in fitness for survival?

Coastal power stations entrain large volumes of cooling water, requiring biocidal treatment to prevent biological fouling. Discharged effluent is both heated and contaminated with residual traces of biocide and so it is necessary to quantify the impacts of this discharge. Cooling water from Heysham 2 nuclear power station, NW England, UK, is discharged to the intertidal area, via a culvert (to minimise erosion and maximise dilution and dispersion by directing the effluent into the receiving water at all states of the tide) within which the effluent is contained at low water. The culvert and surrounding coastal area support a population of blue mussels (Mytilus edulis). Mussel health was determined along a gradient of exposure, using three physiological indices: Scope for Growth, Gonad Mantle Index and Somatic Condition Index (K Factor). The Mussels within the culvert exhibited reduced physiological index values compared to an external site. A trend was identified down the length of the culvert, representing a gradient of exposure and indicating a potential negative effect on growth and reproductive output.     

Assessment of future groundwater recharge in semi‐arid regions under climate change scenarios (Serral‐Salinas aquifer,SE Spain). Could increased rainfall variability increase the recharge rate?

The projected impact of climate change on groundwater recharge is a challenge in hydrogeological research because substantial doubts still remain, particularly in arid and semi‐arid zones. We present a methodology to generate future groundwater recharge scenarios using available information about regional climate change projections developed in European Projects. It involves an analysis of regional climate model (RCM) simulations and a proposal for ensemble models to assess the impacts of climate change. Future rainfall and temperature series are generated by modifying the mean and standard deviation of the historical series in accordance with estimates of their change provoked by climate change. Future recharge series will be obtained by simulating these new series within a continuous balance model of the aquifer. The proposed method is applied to the Serral‐Salinas aquifer, located in a semi‐arid zone of south‐east Spain. The results show important differences depending on the RCM used. Differences are also observed between the series generated by imposing only the changes in means or also in standard deviations. An increase in rainfall variability, as expected under future scenarios, could increase recharge rates for a given mean rainfall because the number of extreme events increases. For some RCMs, the simulations predict total recharge increases over the historical values, even though climate change would produce a reduction in the mean rainfall and an increased mean temperature. A method based on a multi‐objective analysis is proposed to provide ensemble predictions that give more value to the information obtained from the best calibrated models. The ensemble of predictions estimates a reduction in mean annual recharge of 14% for scenario A2 and 58% for scenario A1B. Lower values of future recharge are obtained if only the change in the mean is imposed. Copyright © 2014 John Wiley & Sons, Ltd.