3D geoelectrical monitoring to follow hydric soil behaviour during rainfall

RESUME Durant les dernières années, les méthodes électriques ont souvent été utilisées pour l'investigation des structures de subsurface. L'imagerie électrique (Electrical Resistivity Tomography, ERT) est une technique de prospection non-invasive et spatialement intégrée. La méthode ERT a subi des améliorations significatives avec le développement de nouveaux algorithmes d'inversion et le perfectionnement des techniques d'acquisition. La technologie multicanale et les ordinateurs de dernière génération permettent la collecte et le traitement de données en quelques heures. Les domaines d'application sont nombreux et divers: géologie et hydrogéologie, génie civil et géotechnique, archéologie et études environnementales. En particulier, les méthodes électriques sont souvent employées dans l'étude hydrologique de la zone vadose. Le but de ce travail est le développement d'un système de monitorage 3D automatique, non- invasif, fiable, peu coûteux, basé sur une technique multicanale et approprié pour suivre les variations de résistivité électrique dans le sous-sol lors d'événements pluvieux. En raison des limitations techniques et afin d'éviter toute perturbation physique dans la subsurface, ce dispositif de mesure emploie une installation non-conventionnelle, où toutes les électrodes de courant sont placées au bord de la zone d'étude. Le dispositif le plus approprié pour suivre les variations verticales et latérales de la résistivité électrique à partir d'une installation permanente a été choisi à l'aide de modélisations numériques. Les résultats démontrent que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle et plus apte à détecter les variations latérales et verticales de la résistivité électrique, et cela malgré la configuration non-conventionnelle des électrodes. Pour tester l'efficacité du système proposé, des données de terrain ont été collectées sur un site d'étude expérimental. La technique de monitorage utilisée permet de suivre le processus d'infiltration 3D pendant des événements pluvieux. Une bonne corrélation est observée entre les résultats de modélisation numérique et les données de terrain, confirmant par ailleurs que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle. La nouvelle technique de monitorage 3D de résistivité électrique permet de caractériser les zones d'écoulement préférentiel et de caractériser le rôle de la lithologie et de la pédologie de manière quantitative dans les processus hydrologiques responsables d'écoulement de crue. ABSTRACT During the last years, electrical methods were often used for the investigation of subsurface structures. Electrical resistivity tomography (ERT) has been reported to be a useful non-invasive and spatially integrative prospecting technique. The ERT method provides significant improvements, with the developments of new inversion algorithms, and the increasing efficiency of data collection techniques. Multichannel technology and powerful computers allow collecting and processing resistivity data within few hours. Application domains are numerous and varied: geology and hydrogeology, civil engineering and geotechnics, archaeology and environmental studies. In particular, electrical methods are commonly used in hydrological studies of the vadose zone. The aim of this study was to develop a multichannel, automatic, non-invasive, reliable and inexpensive 3D monitoring system designed to follow electrical resistivity variations in soil during rainfall. Because of technical limitations and in order to not disturb the subsurface, the proposed measurement device uses a non-conventional electrode set-up, where all the current electrodes are located near the edges of the survey grid. Using numerical modelling, the most appropriate arrays were selected to detect vertical and lateral variations of the electrical resistivity in the framework of a permanent surveying installation system. The results show that a pole-dipole array has a better resolution than a pole-pole array and can successfully follow vertical and lateral resistivity variations despite the non-conventional electrode configuration used. Field data are then collected at a test site to assess the efficiency of the proposed monitoring technique. The system allows following the 3D infiltration processes during a rainfall event. A good correlation between the results of numerical modelling and field data results can be observed since the field pole-dipole data give a better resolution image than the pole-pole data. The new device and technique makes it possible to better characterize the zones of preferential flow and to quantify the role of lithology and pedology in flood- generating hydrological processes.

Getting the most sulfate from soil: Regulation of sulfate uptake transporters in Arabidopsis.

Sulfur (S) is an essential macronutrient for all living organisms. Plants require large amounts of sulfate for growth and development, and this serves as a major entry point of sulfate into the food web. Plants acquire S in its ionic form from the soil; they have evolved tightly controlled mechanisms for the regulation of sulfate uptake in response to its external and internal availability. In the model plant Arabidopsis thaliana, the first key step in sulfate uptake is presumed to be carried out exclusively by only two high-affinity sulfate transporters: SULTR1;1 and SULTR1;2. A better understanding of the mode of regulation for these two transporters is crucial because they constitute the first determinative step in balancing sulfate in respect to its supply and demand. Here, we review the recent progress achieved in our comprehension of (i) mechanisms that regulate these two high-affinity sulfate transporters at the transcriptional and post-transcriptional levels, and (ii) their structure-function relationship. Such progress is important to enable biotechnological and agronomic strategies aimed at enhancing sulfate uptake and improving crop yield in S-deficient soils.

Fungi, bacteria and soil pH: the oxalate-carbonate pathway as a model for metabolic interaction

The oxalatecarbonate pathway involves the oxidation of calcium oxalate to low-magnesium calcite and represents a potential long-term terrestrial sink for atmospheric CO2. In this pathway, bacterial oxalate degradation is associated with a strong local alkalinization and subsequent carbonate precipitation. In order to test whether this process occurs in soil, the role of bacteria, fungi and calcium oxalate amendments was studied using microcosms. In a model system with sterile soil amended with laboratory cultures of oxalotrophic bacteria and fungi, the addition of calcium oxalate induced a distinct pH shift and led to the final precipitation of calcite. However, the simultaneous presence of bacteria and fungi was essential to drive this pH shift. Growth of both oxalotrophic bacteria and fungi was confirmed by qPCR on the frc (oxalotrophic bacteria) and 16S rRNA genes, and the quantification of ergosterol (active fungal biomass) respectively. The experiment was replicated in microcosms with non-sterilized soil. In this case, the bacterial and fungal contribution to oxalate degradation was evaluated by treatments with specific biocides (cycloheximide and bronopol). Results showed that the autochthonous microflora oxidized calcium oxalate and induced a significant soil alkalinization. Moreover, data confirmed the results from the model soil showing that bacteria are essentially responsible for the pH shift, but require the presence of fungi for their oxalotrophic activity. The combined results highlight that the interaction between bacteria and fungi is essential to drive metabolic processes in complex environments such as soil.

Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE

PECUBE is a three-dimensional thermal-kinematic code capable of solving the heat production-diffusion-advection equation under a temporally varying surface boundary condition. It was initially developed to assess the effects of time-varying surface topography (relief) on low-temperature thermochronological datasets. Thermochronometric ages are predicted by tracking the time-temperature histories of rock-particles ending up at the surface and by combining these with various age-prediction models. In the decade since its inception, the PECUBE code has been under continuous development as its use became wider and addressed different tectonic-geomorphic problems. This paper describes several major recent improvements in the code, including its integration with an inverse-modeling package based on the Neighborhood Algorithm, the incorporation of fault-controlled kinematics, several different ways to address topographic and drainage change through time, the ability to predict subsurface (tunnel or borehole) data, prediction of detrital thermochronology data and a method to compare these with observations, and the coupling with landscape-evolution (or surface-process) models. Each new development is described together with one or several applications, so that the reader and potential user can clearly assess and make use of the capabilities of PECUBE. We end with describing some developments that are currently underway or should take place in the foreseeable future. (C) 2012 Elsevier B.V. All rights reserved.

Derivation of soil-specific streaming potential electrical parameters from hydrodynamic characteristics of partially saturated soils

Water movement in unsaturated soils gives rise to measurable electrical potential differences that are related to the flow direction and volumetric fluxes, as well as to the soil properties themselves. Laboratory and field data suggest that these so-called streaming potentials may be several orders of magnitudes larger than theoretical predictions that only consider the influence of the relative permeability and electrical conductivity on the self potential (SP) data. Recent work has improved predictions somewhat by considering how the volumetric excess charge in the pore space scales with the inverse of water saturation. We present a new theoretical approach that uses the flux-averaged excess charge, not the volumetric excess charge, to predict streaming potentials. We present relationships for how this effective excess charge varies with water saturation for typical soil properties using either the water retention or the relative permeability function. We find large differences between soil types and the predictions based on the relative permeability function display the best agreement with field data. The new relationships better explain laboratory data than previous work and allow us to predict the recorded magnitudes of the streaming potentials following a rainfall event in sandy loam, whereas previous models predict values that are three orders of magnitude too small. We suggest that the strong signals in unsaturated media can be used to gain information about fluxes (including very small ones related to film flow), but also to constrain the relative permeability function, the water retention curve, and the relative electrical conductivity function.

Burkholderia sartisoli sp. nov., isolated from a polycyclic aromatic hydrocarbon-contaminated soil.

A Gram-negative, rod-shaped, aerobic bacterium, designated strain RP007(T), was isolated from a polycyclic aromatic hydrocarbon-contaminated soil in New Zealand. Two additional strains were recovered from a compost heap in Belgium (LMG 18808) and from the rhizosphere of maize in the Netherlands (LMG 24204). The three strains had virtually identical 16S rRNA gene sequences and whole-cell protein profiles, and they were identified as members of the genus Burkholderia, with Burkholderia phenazinium as their closest relative. Strain RP007(T) had a DNA G+C content of 63.5 mol% and could be distinguished from B. phenazinium based on a range of biochemical characteristics. Strain RP007(T) showed levels of DNA-DNA relatedness towards the type strain of B. phenazinium and those of other recognized Burkholderia species of less than 30 %. The results of 16S rRNA gene sequence analysis, DNA-DNA hybridization experiments and physiological and biochemical tests allowed the differentiation of strain RP007(T) from all recognized species of the genus Burkholderia. Strains RP007(T), LMG 18808 and LMG 24204 are therefore considered to represent a single novel species of the genus Burkholderia, for which the name Burkholderia sartisoli sp. nov. is proposed. The type strain is RP007(T) (=LMG 24000(T) =CCUG 53604(T) =ICMP 13529(T)).

Association of hemolytic activity of Pseudomonas entomophila, a versatile soil bacterium, with cyclic lipopeptide production.

Pseudomonas entomophila is an entomopathogenic bacterium that is able to infect and kill Drosophila melanogaster upon ingestion. Its genome sequence suggests that it is a versatile soil bacterium closely related to Pseudomonas putida. The GacS/GacA two-component system plays a key role in P. entomophila pathogenicity, controlling many putative virulence factors and AprA, a secreted protease important to escape the fly immune response. P. entomophila secretes a strong diffusible hemolytic activity. Here, we showed that this activity is linked to the production of a new cyclic lipopeptide containing 14 amino acids and a 3-C(10)OH fatty acid that we called entolysin. Three nonribosomal peptide synthetases (EtlA, EtlB, EtlC) were identified as responsible for entolysin biosynthesis. Two additional components (EtlR, MacAB) are necessary for its production and secretion. The P. entomophila GacS/GacA two-component system regulates entolysin production, and we demonstrated that its functioning requires two small RNAs and two RsmA-like proteins. Finally, entolysin is required for swarming motility, as described for other lipopeptides, but it does not participate in the virulence of P. entomophila for Drosophila. While investigating the physiological role of entolysin, we also uncovered new phenotypes associated with P. entomophila, including strong biocontrol abilities.

Radiocarbon dating of modern groundwater: the role of the unsaturated zone

Biological and physical processes occurring in soils may lead to significant isotopic changes between the isotopic compositions of atmospheric CO2 and of soil CO2. Also, during water and gas transport from the soil surface to the water table, isotopic changes likely occur due to numerous physical processes such as gas production and diffusion, water advection, and gas-water-rock interactions. In most cases, these changes are not included in the correction models developed for groundwater dating, whereas they can significantly impact the calculation of the 14C age. We explore the role of these processes using: i) experimental data from two aquifer sites (Fontainebleau sands and Astian sands, France), ii) a distributed model to simulate the 14C activities of soil CO2, and iii) numerical simulations in order to highlight the role of the physical processes.¦The 13C content in soil CO2 showed seasonal variations and highlighted the competition between CO2 production and CO2 diffusion. Their respective contributions played a significant role in defining the isotopic composition of CO2 at the water table. On both study sites, variations of the 14C activity in soil CO2 reflect the competition between the fluxes of root derived-CO2 and organic matter derived-CO2. Since the nuclear weapon tests in the fifties and sixties, soil CO2 became significantly depleted in 14C compared to modern atmospheric CO2. Models that take into account this 14C depletion in soil CO2 for dating modern groundwater would lead to apparent younger 14C ages than models that only consider the 14C activity in atmospheric CO2. Moreover, since 2000-2005, the inverse effect is observed as soil CO2 is enriched in 14C compared to atmospheric CO2.¦Therefore, we conclude that the isotopic composition of CO2 at the water table have to be taken into account for the dating of modern groundwater. This requires a systematic sampling of soil CO2 and the measurement of its 13C and 14C contents. We used this information in a numerical simulation to calculate the evolution of isotopic composition of CO2 from the soil surface to the water table. This simulation integrated physical processes in the unsaturated zone (e.g. CO2 production and diffusion, water advection, etc.) and gas-water-rock interactions.

High-resolution numerical modelling of flow-vegetation interactions

In this paper, we present and apply a new three-dimensional model for the prediction of canopy-flow and turbulence dynamics in open-channel flow. The approach uses a dynamic immersed boundary technique that is coupled in a sequentially staggered manner to a large eddy simulation. Two different biomechanical models are developed depending on whether the vegetation is dominated by bending or tensile forces. For bending plants, a model structured on the Euler-Bernoulli beam equation has been developed, whilst for tensile plants, an N-pendula model has been developed. Validation against flume data shows good agreement and demonstrates that for a given stem density, the models are able to simulate the extraction of energy from the mean flow at the stem-scale which leads to the drag discontinuity and associated mixing layer.

Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Natural Arabidopsis brx loss-of-function alleles confer root adaptation to acidic soil.

Soil acidification is a major agricultural problem that negatively affects crop yield. Root systems counteract detrimental passive proton influx from acidic soil through increased proton pumping into the apoplast, which is presumably also required for cell elongation and stimulated by auxin. Here, we found an unexpected impact of extracellular pH on auxin activity and cell proliferation rate in the root meristem of two Arabidopsis mutants with impaired auxin perception, axr3 and brx. Surprisingly, neutral to slightly alkaline media rescued their severely reduced root (meristem) growth by stimulating auxin signaling, independent of auxin uptake. The finding that proton pumps are hyperactive in brx roots could explain this phenomenon and is consistent with more robust growth and increased fitness of brx mutants on overly acidic media or soil. Interestingly, the original brx allele was isolated from a natural stock center accession collected from acidic soil. Our discovery of a novel brx allele in accessions recently collected from another acidic sampling site demonstrates the existence of independently maintained brx loss-of-function alleles in nature and supports the notion that they are advantageous in acidic soil pH conditions, a finding that might be exploited for crop breeding.

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.

Interplay between wheat cultivars, biocontrol pseudomonads, and soil.

There is a significant potential to improve the plant-beneficial effects of root-colonizing pseudomonads by breeding wheat genotypes with a greater capacity to sustain interactions with these bacteria. However, the interaction between pseudomonads and crop plants at the cultivar level, as well as the conditions which favor the accumulation of beneficial microorganisms in the wheat rhizosphere, is largely unknown. Therefore, we characterized the three Swiss winter wheat (Triticum aestivum) cultivars Arina, Zinal, and Cimetta for their ability to accumulate naturally occurring plant-beneficial pseudomonads in the rhizosphere. Cultivar performance was measured also by the ability to select for specific genotypes of 2,4-diacetylphloroglucinol (DAPG) producers in two different soils. Cultivar-specific differences were found; however, these were strongly influenced by the soil type. Denaturing gradient gel electrophoresis (DGGE) analysis of fragments of the DAPG biosynthetic gene phlD amplified from natural Pseudomonas rhizosphere populations revealed that phlD diversity substantially varied between the two soils and that there was a cultivar-specific accumulation of certain phlD genotypes in one soil but not in the other. Furthermore, the three cultivars were tested for their ability to benefit from Pseudomonas inoculants. Interestingly, Arina, which was best protected against Pythium ultimum infection by inoculation with Pseudomonas fluorescens biocontrol strain CHA0, was the cultivar which profited the least from the bacterial inoculant in terms of plant growth promotion in the absence of the pathogen. Knowledge gained of the interactions between wheat cultivars, beneficial pseudomonads, and soil types allows us to optimize cultivar-soil combinations for the promotion of growth through beneficial pseudomonads. Additionally, this information can be implemented by breeders into a new and unique breeding strategy for low-input and organic conditions.

Analyzing complex rock slope deformation at Stampa, western Norway, by integrating geomorphology, kinematics and numerical modeling

The unstable rock slope, Stampa, above the village of Flåm, Norway, shows signs of both active and postglacial gravitational deformation over an area of 11 km2. Detailed structural field mapping, annual differential Global Navigation Satellite System (GNSS) surveys, as well as geomorphic analysis of high-resolution digital elevation models based on airborne and terrestrial laser scanning indicate that slope deformation is complex and spatially variable. Numerical modeling was used to investigate the influence of former rockslide activity and to better understand the failure mechanism. Field observations, kinematic analysis and numerical modeling indicate a strong structural control of the unstable area. Based on the integration of the above analyses, we propose that the failure mechanism is dominated by (1) a toppling component, (2) subsiding bilinear wedge failure and (3) planar sliding along the foliation at the toe of the unstable slope. Using differential GNSS, 18 points were measured annually over a period of up to 6 years. Two of these points have an average yearly movement of around 10 mm/year. They are located at the frontal cliff on almost completely detached blocks with volumes smaller than 300,000 m3. Large fractures indicate deep-seated gravitational deformation of volumes reaching several 100 million m3, but the movement rates in these areas are below 2 mm/year. Two different lobes of prehistoric rock slope failures were dated with terrestrial cosmogenic nuclides. While the northern lobe gave an average age of 4,300 years BP, the southern one resulted in two different ages (2,400 and 12,000 years BP), which represent most likely multiple rockfall events. This reflects the currently observable deformation style with unstable blocks in the northern part in between Joasete and Furekamben and no distinct blocks but a high rockfall activity around Ramnanosi in the south. With a relative susceptibility analysis it is concluded that small collapses of blocks along the frontal cliff will be more frequent. Larger collapses of free-standing blocks along the cliff with volumes > 100,000 m3, thus large enough to reach the fjord, cannot be ruled out. A larger collapse involving several million m3 is presently considered of very low likelihood.