Measuring and modeling water-related soil-vegetation feedbacks in a fallow plot

Abstract. Land fallowing is one possible response to shortage of water for irrigation. Leaving the soil unseeded implies a change of the soil functioning that has an impact on the water cycle. The development of a soil crust in the open spaces between the patterns of grass weed affects the soil properties and the field-scale water balance. The objectives of this study are to test the potential of integrated non-invasive geophysical methods and ground-image analysis and to quantify the effect of the soil–vegetation interaction on the water balance of fallow land at the local- and plot scale. We measured repeatedly in space and time local soil saturation and vegetation cover over two small plots located in southern Sardinia, Italy, during a controlled irrigation experiment. One plot was left unseeded and the other was cultivated. The comparative analysis of ERT maps of soil moisture evidenced a considerably different hydrologic response to irrigation of the two plots. Local measurements of soil saturation and vegetation cover were repeated in space to evidence a positive feedback between weed growth and infiltration at the fallow plot. A simple bucket model captured the different soil moisture dynamics at the two plots during the infiltration experiment and was used to estimate the impact of the soil vegetation feedback on the yearly water balance at the fallow site

Monitoring and modelling of soil–plant interactions: the joint use of ERT, sap flow and eddy covariance data to characterize the volume of an orange tree root zone

Abstract. Mass and energy exchanges between soil, plants and atmosphere control a number of key environmental processes involving hydrology, biota and climate. The understanding of these exchanges also play a critical role for practical purposes e.g. in precision agriculture. In this paper we present a methodology based on coupling innovative data collection and models in order to obtain quantitative estimates of the key parameters of such complex flow system. In particular we propose the use of hydro-geophysical monitoring via "time-lapse" electrical resistivity tomography (ERT) in conjunction with measurements of plant transpiration via sap flow and evapotranspiration (ET) from eddy covariance (EC). This abundance of data is fed to spatially distributed soil models in order to characterize the distribution of active roots. We conducted experiments in an orange orchard in eastern Sicily (Italy), characterized by the typical Mediterranean semi-arid climate. The subsoil dynamics, particularly influenced by irrigation and root uptake, were characterized mainly by the ERT set-up, consisting of 48 buried electrodes on 4 instrumented micro-boreholes (about 1.2 m deep) placed at the corners of a square (with about 1.3 m long sides) surrounding the orange tree, plus 24 mini-electrodes on the surface spaced 0.1 m on a square grid. During the monitoring, we collected repeated ERT and time domain reflectometry (TDR) soil moisture measurements, soil water sampling, sap flow measurements from the orange tree and EC data. We conducted a laboratory calibration of the soil electrical properties as a function of moisture content and porewater electrical conductivity. Irrigation, precipitation, sap flow and ET data are available allowing for knowledge of the system's long-term forcing conditions on the system. This information was used to calibrate a 1-D Richards' equation model representing the dynamics of the volume monitored via 3-D ERT. Information on the soil hydraulic properties was collected from laboratory and field experiments. The successful results of the calibrated modelling exercise allow for the quantification of the soil volume interested by root water uptake (RWU). This volume is much smaller (with a surface area less than 2 m2, and about 40 cm thick) than expected and assumed in the design of classical drip irrigation schemes that prove to be losing at least half of the irrigated water which is not taken up by the plants

2016 Central Italy Earthquakes: comparison between GPS signals and low-cost distributed MEMS arrays

Abstract. Modern seismic ground-motion sensors have reached an excellent performance quality in terms of dynamic range and bandwidth resolution. The weakest point in the recording of seismic events remains spatial sampling and spatial resolution, due to the limited number of installed sensors. A significant improvement in spatial resolution can be achieved by the use of non-conventional motion sensors, such as low-cost distributed sensors arrays or positioning systems, capable of increasing the density of classical seismic recording networks. In this perspective, we adopted micro-electro mechanical system (MEMS) sensors to integrate the use of standard accelerometers for moderate-to-strong seismic events. In addition, we analyse high-rate distributed positioning system data that also record soil motion. In this paper, we present data from the 2016 Central Italy earthquakes as recorded by a spatially dense prototype MEMS array installed in the proximity of the epicentral area, and we compare the results to the signal of local 1s GPS stations. We discuss advantages and limitations of this joint approach, reaching the conclusion that such low-cost sensors and the use of high rate GPS signal could be an effective choice for integrate the spatial density of stations providing strong-motion parameters

Brief communication: Mountain permafrost acts as an aquitard during an infiltration experiment monitored with electrical resistivity tomography time-lapse measurements

Frozen layers within the subsurface of rock glaciers are generally assumed to act as aquicludes or aquitards. So far, this behavior has been mainly defined by analyzing the geochemical characteristics of spring waters. In this work, for the first time, we experimentally confirmed this assumption by executing an infiltration test in a rock glacier of the Southern Alps, Italy. Time-lapse electrical resistivity tomography (ERT) technique monitored the infiltration of 800 L of saltwater released on the surface of the rock glacier; 24 h ERT monitoring highlighted that the injected water was not able to infiltrate into the underlying frozen layer.</p

Use of small scale electrical resistivity tomography to identify soil-root interactions during deficit irrigation

Plant roots activity affect the exchanges of mass and energy between the soil and atmosphere. However, it is challenging to monitor the activity of the root-zone because roots are not visible from the soil surface, and root systems undergo spatial and temporal variations in response to internal and external conditions. Therefore, measurements of the activity of root systems are interesting to ecohydrologists in general, and are especially important for specific applications, such as irrigation water management. This study demonstrates the use of small scale three-dimensional (3-D) electrical resistivity tomography (ERT) to monitor the root-zone of orange trees irrigated by two different regimes: (i) full rate, in which 100% of the crop evapotranspiration (ETc) is provided; and (ii) partial root-zone drying (PRD), in which 50% of ETc is supplied to alternate sides of the tree. We performed time-lapse 3-D ERT measurements on these trees from 5 June to 24 September 2015, and compared the long-term and short-term changes before, during, and after irrigation events. Given the small changes in soil temperature and pore water electrical conductivity, we interpreted changes of soil electrical resistivity from 3-D ERT data as proxies for changes in soil water content. The ERT results are consistent with measurements of transpiration flux and soil temperature. The changes in electrical resistivity obtained from ERT measurements in this case study indicate that root water uptake (RWU) processes occur at the 0.1 m scale, and highlight the impact of different irrigation schemes. (C) 2017 Elsevier B.V. All rights reserved

Small-scale characterization of vine plant root water uptake via 3-D electrical resistivity tomography and mise-à-la-masse method

The investigation of plant roots is inherently difficult and often neglected. Being out of sight, roots are often out of mind. Nevertheless, roots play a key role in the exchange of mass and energy between soil and the atmosphere, in addition to the many practical applications in agriculture. In this paper, we propose a method for roots imaging based on the joint use of two electrical noninvasive methods: electrical resistivity tomography (ERT) and mise-à-la-masse (MALM). The approach is based on the key assumption that the plant root system acts as an electrically conductive body, so that injecting electrical current into the plant stem will ultimately result in the injection of current into the subsoil through the root system, and particularly through the root terminations via hair roots. Evidence from field data, showing that voltage distribution is very different whether current is injected into the tree stem or in the ground, strongly supports this hypothesis. The proposed procedure involves a stepwise inversion of both ERT and MALM data that ultimately leads to the identification of electrical resistivity (ER) distribution and of the current injection root distribution in the three-dimensional soil space. This, in turn, is a proxy to the active (hair) root density in the ground. We tested the proposed procedure on synthetic data and, more importantly, on field data collected in a vineyard, where the estimated depth of the root zone proved to be in agreement with literature on similar crops. The proposed noninvasive approach is a step forward towards a better quantification of root structure and functioning.</p

2D characterization of near-surface V P/V S: surface-wave dispersion inversion versus refraction tomography

International audienceThe joint study of pressure (P-) and shear (S-) wave velocities (Vp and Vs ), as well as their ratio (Vp /Vs), has been used for many years at large scales but remains marginal in near-surface applications. For these applications, and are generally retrieved with seismic refraction tomography combining P and SH (shear-horizontal) waves, thus requiring two separate acquisitions. Surface-wave prospecting methods are proposed here as an alternative to SH-wave tomography in order to retrieve pseudo-2D Vs sections from typical P-wave shot gathers and assess the applicability of combined P-wave refraction tomography and surface-wave dispersion analysis to estimate Vp/Vs ratio. We carried out a simultaneous P- and surface-wave survey on a well-characterized granite-micaschists contact at Ploemeur hydrological observatory (France), supplemented with an SH-wave acquisition along the same line in order to compare Vs results obtained from SH-wave refraction tomography and surface-wave profiling. Travel-time tomography was performed with P- and SH- wave first arrivals observed along the line to retrieve Vtomo p and Vtomo s models. Windowing and stacking techniques were then used to extract evenly spaced dispersion data from P-wave shot gathers along the line. Successive 1D Monte Carlo inversions of these dispersion data were performed using fixed Vp values extracted from Vtomo p the model and no lateral constraints between two adjacent 1D inversions. The resulting 1D Vsw s models were then assembled to create a pseudo-2D Vsw s section, which appears to be correctly matching the general features observed on the section. If the pseudo-section is characterized by strong velocity incertainties in the deepest layers, it provides a more detailed description of the lateral variations in the shallow layers. Theoretical dispersion curves were also computed along the line with both and models. While the dispersion curves computed from models provide results consistent with the coherent maxima observed on dispersion images, dispersion curves computed from models are generally not fitting the observed propagation modes at low frequency. Surface-wave analysis could therefore improve models both in terms of reliability and ability to describe lateral variations. Finally, we were able to compute / sections from both and models. The two sections present similar features, but the section obtained from shows a higher lateral resolution and is consistent with the features observed on electrical resistivity tomography, thus validating our approach for retrieving Vp/Vs ratio from combined P-wave tomography and surface-wave profiling

Shear wave profiles from surface wave inversion: the impact of uncertainty on seismic site response analysis

Inversion is a critical step in all geophysical techniques, and is generally fraught with ill-posedness. In the case of seismic surface wave studies, the inverse problem can lead to different equivalent subsoil models and consequently to different local seismic response analyses. This can have a large impact on an earthquake engineering design. In this paper, we discuss the consequences of non-uniqueness of surface wave inversion on seismic responses, with both numerical and experimental data. Our goal is to evaluate the consequences on common seismic response analysis in the case of different impedance contrast conditions. We verify the implications of inversion uncertainty, and consequently of data information content, on realistic local site responses. A stochastic process is used to generate a set of 1D shear wave velocity profiles from several specific subsurface models. All these profiles are characterized as being equivalent, i.e. their responses, in terms of a dispersion curve, are compatible with the uncertainty in the same surface wave data. The generated 1D shear velocity models are then subjected to a conventional one-dimensional seismic ground response analysis using a realistic input motion. While recent analyses claim that the consequences of surface wave inversion uncertainties are very limited, our test points out that a relationship exists between inversion confidence and seismic responses in different subsoils. In the case of regular and relatively smooth increase of shear wave velocities with depth, as is usual in sedimentary plains, our results show that the choice of a specific model among equivalent solutions strongly influences the seismic response. On the other hand, when the shallow subsoil is characterized by a strong impedance contrast (thus revealing a characteristic soil resonance period), as is common in the presence of a shallow bedrock, equivalent solutions provide practically the same seismic amplification, especially in the frequency range of engineering interest