Heated Optical Fiber for Distributed Soil-Moisture Measurements: A Lysimeter Experiment

An Actively Heated Fiber Optics (AHFO) method to estimate soil moisture is tested and the analysis technique improved on. The measurements were performed in a lysimeter uniformly packed with loam soil with variable water content profiles. In the first meter of the soil profi le, 30 m of fiber optic cable were installed in a 12 loops coil. The metal sheath armoring the fiber cable was used as an electrical resistance heater to generate a heat pulse, and the soil response was monitored with a Distributed Temperature Sensing (DTS) system. We study the cooling following three continuous heat pulses of 120 s at 36 W m(-1) by means of long-time approximation of radial heat conduction. The soil volumetric water contents were then inferred from the estimated thermal conductivities through a specifically calibrated model relating thermal conductivity and volumetric water content. To use the pre-asymptotic data we employed a time correction that allowed the volumetric water content to be estimated with a precision of 0.01-0.035 (m(3) m(-3)). A comparison of the AHFO measurements with soil-moisture measurements obtained with calibrated capacitance-based probes gave good agreement for wetter soils [discrepancy between the two methods was less than 0.04 (m(3) m(-3))]. In the shallow drier soils, the AHFO method underestimated the volumetric water content due to the longertime required for the temperature increment to become asymptotic in less thermally conductive media [discrepancy between the two methods was larger than 0.1 (m(3) m(-3))]. The present work suggests that future applications of the AHFO method should include longer heat pulses, that longer heating and cooling events are analyzed, and, temperature increments ideally be measured with higher frequency.

On the diurnal soil water content dynamics during evaporation using dielectric methods

The water content dynamics in the upper soil surface during evaporation is a key element in land-atmosphere exchanges. Previous experimental studies have suggested that the soil water content increases at the depth of 5 to 15 cm below the soil surface during evapo- ration, while the layer in the immediate vicinity of the soil surface is drying. In this study, the dynamics of water content profiles exposed to solar radiative forcing was monitored at a high temporal resolution using dielectric methods both in the presence and absence of evaporation. A 4-d comparison of reported moisture content in coarse sand in covered and uncovered buckets using a commercial dielectric-based probe (70 MHz ECH2O-5TE, Decagon Devices, Pullman, WA) and the standard 1-GHz time domain reflectometry method. Both sensors reported a positive correlation between temperature and water content in the 5- to 10-cm depth, most pronounced in the morning during heating and in the afternoon during cooling. Such positive correlation might have a physical origin induced by evaporation at the surface and redistribution due to liquid water fluxes resulting from the temperature- gradient dynamics within the sand profile at those depths. Our experimental data suggest that the combined effect of surface evaporation and temperature-gradient dynamics should be considered to analyze experimental soil water profiles. Additional effects related to the frequency of operation and to protocols for temperature compensation of the dielectric sensors may also affect the probes' response during large temperature changes.

Testing the influence of topography and material properties on catchment-scale soil moisture patterns using remotely sensed vegetation patterns in a humid temperate catchment, northern Britain

In order to evaluate the relationship between the apparent complexity of hillslope soil moisture and the emergent patterns of catchment hydrological behaviour and water quality, we need fine-resolution catchment-wide data on soil moisture characteristics. This study proposes a methodology whereby vegetation patterns obtained from high-resolution orthorectified aerial photographs are used as an indicator of soil moisture characteristics. This enables us to examine a set of hypotheses regarding what drives the spatial patterns of soil moisture at the catchment scale (material properties or topography). We find that the pattern of Juncus effusus vegetation is controlled largely by topography and mediated by the catchment's material properties. Characterizing topography using the topographic index adds value to the soil moisture predictions relative to slope or upslope contributing area (UCA). However, these predictions depart from the observed soil moisture patterns at very steep slopes or low UCAs. Copyright (c) 2012 John Wiley & Sons, Ltd.

Changes in reproductive investment with altitude in an alpine plant

Aims: In perennial species, the allocation of resources to reproduction results in a reduction of allocation to vegetative growth and, therefore, impacts future reproductive success. As a consequence, variation in this trade-off is among the most important driving forces in the life-history evolution of perennial plants and can lead to locally adapted genotypes. In addition to genetic variation, phenotypic plasticity might also contribute to local adaptation of plants to local conditions by mediating changes in reproductive allocation. Knowledge on the importance of genetic and environmental effects on the trade-off between reproduction and vegetative growth is therefore essential to understand how plants may respond to environmental changes. Methods: We conducted a transplant experiment along an altitudinal gradient from 425 m to 1921 m in the front range of the Western Alps of Switzerland to assess the influence of both altitudinal origin of populations and altitude of growing site on growth, reproductive investment and local adaptation in Poa alpina. Important findings: In our study, the investment in reproduction increased with plant size. Plant growth and the relative importance of reproductive investment decreased in populations originating from higher altitudes compared to populations originating from lower altitudes. The changes in reproductive investment were mainly explained by differences in plant size. In contrast to genetic effects, phenotypic plasticity of all traits measured was low and not related to altitude. As a result, the population from the lowest altitude of origin performed best at all sites. Our results indicate that in P. alpina genetic differences in growth and reproductive investment are related to local conditions affecting growth, i.e. interspecific competition and soil moisture content.

Identification of opsA, a gene involved in solute stress mitigation and survival in soil, in the polycyclic aromatic hydrocarbon-degrading bacterium Novosphingobium sp. strain LH128.

The aim of this study was to identify genes involved in solute and matric stress mitigation in the polycyclic aromatic hydrocarbon (PAH)-degrading Novosphingobium sp. strain LH128. The genes were identified using plasposon mutagenesis and by selection of mutants that showed impaired growth in a medium containing 450 mM NaCl as a solute stress or 10% (wt/vol) polyethylene glycol (PEG) 6000 as a matric stress. Eleven and 14 mutants showed growth impairment when exposed to solute and matric stresses, respectively. The disrupted sequences were mapped on a draft genome sequence of strain LH128, and the corresponding gene functions were predicted. None of them were shared between solute and matric stress-impacted mutants. One NaCl-affected mutant (i.e., NA7E1) with a disruption in a gene encoding a putative outer membrane protein (OpsA) was susceptible to lower NaCl concentrations than the other mutants. The growth of NA7E1 was impacted by other ions and nonionic solutes and by sodium dodecyl sulfate (SDS), suggesting that opsA is involved in osmotic stress mitigation and/or outer membrane stability in strain LH128. NA7E1 was also the only mutant that showed reduced growth and less-efficient phenanthrene degradation in soil compared to the wild type. Moreover, the survival of NA7E1 in soil decreased significantly when the moisture content was decreased but was unaffected when soluble solutes from sandy soil were removed by washing. opsA appears to be important for the survival of strain LH128 in soil, especially in the case of reduced moisture content, probably by mitigating the effects of solute stress and retaining membrane stability.

Oxygen-sensing reporter strain of Pseudomonas fluorescens for monitoring the distribution of low-oxygen habitats in soil.

The root-colonizing bacterium Pseudomonas fluorescens CHA0 was used to construct an oxygen-responsive biosensor. An anaerobically inducible promoter of Pseudomonas aeruginosa, which depends on the FNR (fumarate and nitrate reductase regulation)-like transcriptional regulator ANR (anaerobic regulation of arginine deiminase and nitrate reductase pathways), was fused to the structural lacZ gene of Escherichia coli. By inserting the reporter fusion into the chromosomal attTn7 site of P. fluorescens CHA0 by using a mini-Tn7 transposon, the reporter strain, CHA900, was obtained. Grown in glutamate-yeast extract medium in an oxystat at defined oxygen levels, the biosensor CHA900 responded to a decrease in oxygen concentration from 210 x 10(2) Pa to 2 x 10(2) Pa of O(2) by a nearly 100-fold increase in beta-galactosidase activity. Half-maximal induction of the reporter occurred at about 5 x 10(2) Pa. This dose response closely resembles that found for E. coli promoters which are activated by the FNR protein. In a carbon-free buffer or in bulk soil, the biosensor CHA900 still responded to a decrease in oxygen concentration, although here induction was about 10 times lower and the low oxygen response was gradually lost within 3 days. Introduced into a barley-soil microcosm, the biosensor could report decreasing oxygen concentrations in the rhizosphere for a 6-day period. When the water content in the microcosm was raised from 60% to 85% of field capacity, expression of the reporter gene was elevated about twofold above a basal level after 2 days of incubation, suggesting that a water content of 85% caused mild anoxia. Increased compaction of the soil was shown to have a faster and more dramatic effect on the expression of the oxygen reporter than soil water content alone, indicating that factors other than the water-filled pore space influenced the oxygen status of the soil. These experiments illustrate the utility of the biosensor for detecting low oxygen concentrations in the rhizosphere and other soil habitats.

Speciation of naturally-accumulated uranium in an organic-rich soil of an alpine region (Switzerland)

Very high concentrations of uranium (up to 4000 ppm) were found in a natural soil in the Dischma valley, an alpine region in the Grisons canton in Switzerland. The goal of this study was to examine the redox state and the nature of uranium binding in the soil matrix in order to understand the accumulation mechanism. Pore water profiles collected from Dischma soil revealed the establishment of anoxic conditions with increasing soil depth. A combination of chemical extraction methods and spectroscopy was applied to characterize the redox state and binding environment of uranium in the soil. Bicarbonate extraction under anoxic conditions released most of the uranium indicating that uranium occurs predominantly in the hexavalent form. Surprisingly, the uranium redox state did not vary greatly as a function of depth. X-ray absorption near edge spectroscopy (XANES), confirmed that uranium was present as a mixture of U(VI) and U(IV) with U(VI) dominating. Sequential extractions of soil samples showed that the dissolution of solid organic matter resulted in the simultaneous release of the majority of the soil uranium content (>95%). Extended X-ray absorption fine structure (EXAFS) spectroscopy also revealed that soil-associated uranium in the soil matrix was mainly octahedrally coordinated, with an average of 1.7 axial (at about 1.76 Å) and 4.6 to 5.3 equatorial oxygen atoms (at about 2.36 Å) indicating the dominance of a uranyl-like (UO22+) structure presumably mixed with some U(IV). An additional EXAFS signal (at about 3.2 Å) identified in some spectra suggested that uranium was also bound (via an oxygen atom) to a light element such as carbon, phosphorus or silicon. Gamma spectrometric measurements of soil profiles failed to identify uranium long-life daughter products in the soil which is an indication that uranium originates elsewhere and was transported to its current location by water. Finally, it was found that the release of uranium from the soil was significantly promoted at very low pH values (pH 2) and increased with increasing pH values (between pH 5 and 9).

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data

Time-lapse geophysical monitoring and inversion are valuable tools in hydrogeology for monitoring changes in the subsurface due to natural and forced (tracer) dynamics. However, the resulting models may suffer from insufficient resolution, which leads to underestimated variability and poor mass recovery. Structural joint inversion using cross-gradient constraints can provide higher-resolution models compared with individual inversions and we present the first application to time-lapse data. The results from a synthetic and field vadose zone water tracer injection experiment show that joint 3-D time-lapse inversion of crosshole electrical resistance tomography (ERT) and ground penetrating radar (GPR) traveltime data significantly improve the imaged characteristics of the point injected plume, such as lateral spreading and center of mass, as well as the overall consistency between models. The joint inversion method appears to work well for cases when one hydrological state variable (in this case moisture content) controls the time-lapse response of both geophysical methods. Citation: Doetsch, J., N. Linde, and A. Binley (2010), Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data, Geophys. Res. Lett., 37, L24404, doi: 10.1029/2010GL045482.

Improving the prediction of plant species distribution and community composition by adding edaphic to topo-climatic variables

Questions Soil properties have been widely shown to influence plant growth and distribution. However, the degree to which edaphic variables can improve models based on topo-climatic variables is still unclear. In this study, we tested the roles of seven edaphic variables, namely (1) pH; (2) the content of nitrogen and of (3) phosphorus; (4) silt; (5) sand; (6) clay and (7) carbon-to-nitrogen ratio, as predictors of species distribution models in an edaphically heterogeneous landscape. We also tested how the respective influence of these variables in the models is linked to different ecological and functional species characteristics. Location The Western Alps, Switzerland. Methods With four different modelling techniques, we built models for 115 plant species using topo-climatic variables alone and then topo-climatic variables plus each of the seven edaphic variables, one at a time. We evaluated the contribution of each edaphic variable by assessing the change in predictive power of the model. In a second step, we evaluated the importance of the two edaphic variables that yielded the largest increase in predictive power in one final set of models for each species. Third, we explored the change in predictive power and the importance of variables across plant functional groups. Finally, we assessed the influence of the edaphic predictors on the prediction of community composition by stacking the models for all species and comparing the predicted communities with the observed community. Results Among the set of edaphic variables studied, pH and nitrogen content showed the highest contributions to improvement of the predictive power of the models, as well as the predictions of community composition. When considering all topo-climatic and edaphic variables together, pH was the second most important variable after degree-days. The changes in model results caused by edaphic predictors were dependent on species characteristics. The predictions for the species that have a low specific leaf area, and acidophilic preferences, tolerating low soil pH and high humus content, showed the largest improvement by the addition of pH and nitrogen in the model. Conclusions pH was an important predictor variable for explaining species distribution and community composition of the mountain plants considered in our study. pH allowed more precise predictions for acidophilic species. This variable should not be neglected in the construction of species distribution models in areas with contrasting edaphic conditions.

On the use of spatially discrete data to compute energy and mass balance

In many practical applications the state of field soils is monitored by recording the evolution of temperature and soil moisture at discrete depths. We theoretically investigate the systematic errors that arise when mass and energy balances are computed directly from these measurements. We show that, even with no measurement or model errors, large residuals might result when finite difference approximations are used to compute fluxes and storage term. To calculate the limits set by the use of spatially discrete measurements on the accuracy of balance closure, we derive an analytical solution to estimate the residual on the basis of the two key parameters: the penetration depth and the distance between the measurements. When the thickness of the control layer for which the balance is computed is comparable to the penetration depth of the forcing (which depends on the thermal diffusivity and on the forcing period) large residuals arise. The residual is also very sensitive to the distance between the measurements, which requires accurately controlling the position of the sensors in field experiments. We also demonstrate that, for the same experimental setup, mass residuals are sensitively larger than the energy residuals due to the nonlinearity of the moisture transport equation. Our analysis suggests that a careful assessment of the systematic mass error introduced by the use of spatially discrete data is required before using fluxes and residuals computed directly from field measurements.

Estimating vadose zone hydraulic properties using ground penetrating radar: The impact of prior information

A number of geophysical methods, such as ground-penetrating radar (GPR), have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, the stochastic inversion of such data within a coupled geophysical-hydrological framework may allow for the effective estimation of vadose zone hydraulic parameters and their corresponding uncertainties. A critical issue in stochastic inversion is choosing prior parameter probability distributions from which potential model configurations are drawn and tested against observed data. A well chosen prior should reflect as honestly as possible the initial state of knowledge regarding the parameters and be neither overly specific nor too conservative. In a Bayesian context, combining the prior with available data yields a posterior state of knowledge about the parameters, which can then be used statistically for predictions and risk assessment. Here we investigate the influence of prior information regarding the van Genuchten-Mualem (VGM) parameters, which describe vadose zone hydraulic properties, on the stochastic inversion of crosshole GPR data collected under steady state, natural-loading conditions. We do this using a Bayesian Markov chain Monte Carlo (MCMC) inversion approach, considering first noninformative uniform prior distributions and then more informative priors derived from soil property databases. For the informative priors, we further explore the effect of including information regarding parameter correlation. Analysis of both synthetic and field data indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when we combine these data with a realistic, informative prior.

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.