Vadose Research Articles

Application of electrical and electromagnetic data to delineate shallow aquifer zones: a case study of Vea catchment, Ghana

Application of electrical and electromagnetic data to delineate shallow aquifer zones: a case study of Vea catchment, Ghana

Study on Soil Selenium-Induced Copper Deficiency in Yudong Black Goats.

Due to the degradation of pasture and strict restrictions on grazing ranges in recent years, copper (Cu) deficiency in Yudong black goats has been occurring, mainly manifested as emaciation, anemia, loss of appetite and lack of spirit. To explore the main causes of Cu deficiency in Yudong black goats, 40 black goats (1 year old, 25.11 ± 0.52 kg) were selected for this experiment; among them, 20 Yudong black goats with Cu deficiency from the experimental pasture were used as the experimental group, and 20 healthy Yudong black goats from the control pasture were used as the control group. In the pre-experiment, the mineral contents of the soil, forage, blood, and liver of black goats in both groups were determined, and in formal experiments, blood hematological, biochemical, antioxidant, and hemorheological parameters were analyzed. An experiment on the treatment of Cu deficiency in black goats was also conducted. This study showed that selenium (Se) levels in the soil, forage, blood, and liver from the experimental group were significantly lower than those from the control group (p < 0.01). The content of sulfur (S) in the forage was considerably higher than that of the control group (p < 0.01). The contents of Cu in the blood and liver from the experimental group were significantly lower than that from the control group (p < 0.01), and the content of S was considerably higher than that from the control group (p < 0.01). The blood hematology of the experimental group was affected, as evidenced by a decrease in hemoglobin, hematocrit value, mean corpuscular volume and mean corpuscular hemoglobin. The immunity and antioxidant capacity of black goats in the experiment group were impaired to varying degrees, with significant decreases in ceruloplasmin, immunoglobulin M, immunoglobulin G, glutathione peroxidase, and superoxide dismutase, and substantial increases in malondialdehyde. In addition, the experimental group showed a decrease in blood viscosity as evidenced by the rise in high shear viscosity, low shear viscosity, erythrocyte rigidity index, erythrocyte aggregation index, and erythrocyte deformation index, and a decrease in plasma viscosity. In the treatment experiment, oral administration of copper sulfate solution was carried out on 10 black goats with Cu deficiency. All the Cu deficiency goats were cured, and the Cu content in their bodies rebounded. In summary, low Se soil caused an increase in S content in the forage, and Yudong black goats feeding on high S forage resulted in a decrease in Cu absorption, which led to a secondary Cu deficiency.

Deep soil desiccation hinders nitrate leaching to groundwater in the global largest apple cultivation area

Massive use of N fertilizer in apple orchards aids China’s Loess Plateau (CLP) to become the largest region of apple production in the world, but it also leads to substantial nitrate accumulation in the vadose zone, posing a potential environmental risk. However, a holistic understanding of the spatial patterns and controls of soil nitrate in the vadose zone of apple orchards remains understudied, which is unfavorable for the formulation of a science-based N application strategy. To this end, we collected a total of 2852 paired data across apple orchards on the CLP from 52 published literature and 1000 paired data from our field sampling. We found that residual nitrate accumulation primarily occurred in the top 4–8.8 m soil layer, and it varied according to mean annual precipitation (MAP) and stand age; the accumulation peaked in sites with MAP of 500–600 mm and declined with a decrease in MAP; and among various stand ages, the highest nitrate content was observed in the 25-yr-old orchard. It was found that MAP and N fertilizer application amounts dominated nitrate accumulation in the soil. Furthermore, the desiccation index of the soil layer below the maximum cumulative depth was identified as the main factor influencing soil nitrate maximum accumulation depth . This means that deep soil desiccation hinders nitrate leaching towards deeper layers and reduces the risk of groundwater contamination. The findings here may provide insights into sustainable nitrogen and water management in apple orchards on the CLP and potentially in other similar regions.

Mechanism of vertical displacement beneath shallow compression zone in coastal area

Previous studies showed that there is a discrepancy between total subsidence (measured on the ground surface) and shallow compression (measured within several hundred meters beneath the ground surface) in Taiwan. This difference is referred to as deep displacement in this study. The variations of deep displacement are opposite to those of ground surface displacement and groundwater level within the depth of several hundred meters. The mechanism is unknown and requires further investigation. This study adopts two kinds of geodetic observation data and separates land subsidence into shallow compression and deep displacement to investigate the mechanism of deep displacement. A tectonically active coastal area in Taiwan is selected as the study area. The assessment results show that the associated variations are likely due to cyclic hydraulic loading and unloading. The variations of deep displacement are opposite to those of ground surface displacement. This study proposes that these variations are due to hydraulic expansion and contraction. The mechanism is demonstrated using a hydromechanical model. The results of a cross-correlation analysis show that hydraulic expansion and contraction occur at certain depths. The study results provide important information on the mechanism of deep displacement that can be used in tectonophysical and land subsidence investigations.Graphical

Rapid morphosedimentary development of a shallow subtropical embayment in response to extreme episodic flooding

Variability in supply of terrestrial sediments to the coastal zone impacts sediment transport processes and ecological dependencies. This study investigated surface sediment distribution and transportation within a 2020 km2 shallow embayment following a period of flood induced sediment supply. Extensive field sampling of 221 common sites between 2015 and the 2019 follow-up survey presented, identified significant reduction in fine sediment fractions, evident through a 367 km2 decrease in surface mud coverage. Wind-wave sediment transport processes controlled the spatial distribution of fine sediments, exporting 23 million tonnes of mud from the shallow sedimentation zones. In a broader context, sensitive benthic habitats such as seagrass meadows exhibited signs habitat recolonisation coinciding with areas of reduced mud content. Despite these promising short-term outlooks, the developed fine sediment budget suggests ongoing terrestrial sedimentation will exhaust key sediment sinks within Moreton Bay, effectively reducing the system's ability to buffer future sedimentation events.

Groundwater recharge estimation in Mediterranean mountain environments by isotope profiles – Partitioning of macropore and matrix flow

Developing and continuously improving effective and feasible methods to estimate groundwater recharge is essential to ensure sustainable resource assessment. Mountain environments represent key areas for the regional water balance. In this study, the isotope profile method was used to investigate water fluxes through the vadose zone and subsequent groundwater recharge in a Mediterranean mountain environment. Six high-resolution soil water isotope profiles were taken at two different locations in the Troodos Massif (Cyprus) in the eastern Mediterranean, focusing on estimating recharge in the soil of hillslopes. Analysis of the isotope profiles yielded recharge rates ranging from 30–70 mm/year, leading to a re-evaluation of previous studies on groundwater recharge. The comparison of results obtained by the peak-shift method and mathematical modeling was used to partition matrix and macropore flow, reaching up to 39 %. Potential methodological limits of the environmental isotope profile method in skeleton-rich soils were found.

Subsurface Transport of Plutonium in Organic and Aqueous Acidic Processing Wastes at the Hanford Site, USA.

Over 4 million liters of mixed acidic (∼pH 2.5), high ionic strength (∼5 M nitrate) plutonium (Pu) processing waste were released into the 216-Z-9 (Z-9) trench at the Hanford Site, USA, and trace Pu has migrated 37 m below the trench. In this study, we used flowthrough columns to investigate Pu transport in simplified processing waste through uncontaminated Hanford sediments to determine the conditions that led to Pu migration. In low pH aqueous fluids, some Pu breakthrough is observed at pH < 4, and increased Pu transport (14% total Pu breakthrough) is observed at pH < 2. However, Pu migrates in organic processing solvents through low pH sediments virtually uninhibited with approximately 94 and 86% total Pu breakthrough observed at pH 1 and pH 3, respectively. This study demonstrates that Pu migration can occur both with and without organic solvents at pH < 4, but significantly more Pu can be transported when partitioned into organic processing solvents. Our data suggest that under acidic conditions (pH < 4) in the vadose zone beneath the Z-9 trench, Pu present in organic processing solvents moved relatively unhindered and may explain the historical downward migration of Pu tens of meters below the Z-9 trench.

Inverse modeling of frequency domain-based one-dimensional soil water flow in layered soils

Soils are heterogeneous at multiple scales. Estimates of spatially varying saturated hydraulic conductivity (Ks) of soils dictate the resolution of modeling water flux based on the Richardson-Richards equation (RRE), but direct measurements of Ks everywhere in a field are practically impossible. As a result, this study develops an inverse approach to estimate spatially varying Ks utilizing soil moisture observations driven by periodic weather variability at different frequencies. Using a numerical solution to vertical, one-dimensional RRE based on Gardner-Kozeny (GK) model in a heterogeneous soil with periodic flux, we investigate the spatial cross-correlation between Ks and observations (the amplitude or phase shift) of soil moisture fluctuation at different frequencies without high computational cost in the conventional time-domain model. These correlation patterns vary with the spatial location and frequency, suggesting that soil moisture fluctuations at different frequencies at one monitoring location carry non-redundant information about the heterogeneous distribution of Ks. Guided by cross-correlation maps, several inverse experiments reveal the effectiveness of multi-frequency data of soil moisture fluctuation for mapping the heterogeneous distribution of Ks. Monte Carlo simulations suggest that multi-frequency data can provide non-redundant and useful information about the estimation of heterogeneous Ks. Synthetic data based on a commonly-applied nonlinear van Genuchten-Mualem (VG) model is then tested in our inversion algorithm. The results indicate that the spatial pattern of the inverted GK Ks is very similar to that of the VG Ks, although these two parameters are inherently not identical since they are used in two different models. In addition, the inverted GK Ks can successfully predict the soil moisture fluctuation of frequencies that were not used in the inversion. Despite having no field experiments to validate this model, we believe this is the first attempt to conduct frequency-domain inverse modeling in vadose zone hydrology using soil moisture fluctuation.

Amended Vegetation Filters as Nature-Based Solutions for the Treatment of Pharmaceuticals: Infiltration Experiments Coupled to Reactive Transport Modelling.

In small populations and scattered communities, wastewater treatment through vegetation filters (VFs), a nature-based solution, has proved to be feasible, especially for nutrient and organic matter removal. However, the presence of pharmaceuticals in wastewater and their potential to infiltrate through the vadose zone and reach groundwater is a drawback in the evaluation of VF performances. Soil amended with readily labile carbon sources, such as woodchips, enhances microbial activity and sorption processes, which could improve pharmaceutical attenuation in VFs. The present study aims to assess if woodchip amendments to a VF's soil are able to abate concentrations of selected pharmaceuticals in the infiltrating water by quantitatively describing the occurring processes through reactive transport modelling. Thus, a column experiment using soil collected from an operating VF and poplar woodchips was conducted, alongside a column containing only soil used as reference. The pharmaceuticals acetaminophen, naproxen, atenolol, caffeine, carbamazepine, ketoprofen and sulfamethoxazole were applied daily to the column inlet, mimicking a real irrigation pattern and periodically measured in the effluent. Ketoprofen was the only injected pharmaceutical that reached the column outlet of both systems within the experimental timeframe. The absence of acetaminophen, atenolol, caffeine, carbamazepine, naproxen and sulfamethoxazole in both column outlets indicates that they were attenuated even without woodchips. However, the presence of 10,11-epoxy carbamazepine and atenolol acid as transformation products (TPs) suggests that incomplete degradation also occurs and that the effect of the amendment on the infiltration of TPs is compound-specific. Modelling allowed us to generate breakthrough curves of ketoprofen in both columns and to obtain transport parameters during infiltration. Woodchip-amended columns exhibited Kd and μw values from one to two orders of magnitude higher compared to soil column. This augmentation of sorption and biodegradation processes significantly enhanced the removal of ketoprofen to over 96%.

Wastewater irrigation beneath the water table: analytical model of crop contamination risks

Wastewater irrigation alleviates freshwater scarcity. However, conventional (near)surface irrigation techniques directly expose crops to contaminants. Irrigating wastewater into shallow phreatic zones to raise the water table enhances groundwater evapotranspiration, while using the vadose zone as a bioreactor that attenuates contaminants through dilution, adsorption, and biodegradation. Nevertheless, contaminants may spread across the groundwater, soil, and vegetation. In this study, we focus on the crop contamination risks, and derive a simple analytical model to estimate crop solute uptake. Although crops are not directly exposed to the irrigated wastewater, contaminants (and nutrients) may spread to the root zone. Results show that crop contamination is primarily determined by the root zone water balance, and by solute dispersion and biogeochemical reaction parameters. The model contributes towards identifying hydrogeologically and climatically suitable locations for phreatic zone wastewater irrigation, determining acceptable levels of irrigation water quality, and evaluating crop contamination hazards against the fertigative value of wastewater.

The Impact of the Three-Dimensional Structure of a Subduction Zone on Time-dependent Crustal Deformation Measured by HR-GNSS

Accurately modeling time-dependent coseismic crustal deformation as observed on high-rate Global Navigation Satellite System (HR-GNSS) lends insight into earthquake source processes and improves local earthquake and tsunami early warning algorithms. Currently, time-dependent crustal deformation modeling relies most frequently on simplified 1D radially symmetric Earth models. However, for shallow subduction zone earthquakes, even low-frequency shaking is likely affected by the many strongly heterogeneous structures such as the subducting slab, mantle wedge, and the overlying crustal structure. We demonstrate that including 3D structure improves the estimation of key features of coseismic HR-GNSS time series, such as the peak ground displacement (PGD), the time to PGD (tPGD), static displacements (SD), and waveform cross-correlation values. We computed synthetic 1D and 3D, 0.25 Hz and 0.5 Hz waveforms at HR-GNSS stations for four M7.3+ earthquakes in Japan using MudPy and SW4, respectively. From these synthetics, we computed intensity-measure residuals between the synthetic and observed GNSS waveforms. Comparing 1D and 3D residuals, we observed that the 3D simulations show better fits to the PGD and SD in the observed waveforms than the 1D simulations for both 0.25 Hz and 0.5 Hz simulations. We find that the reduction in PGD residuals in the 3D simulations is a combined effect of both shallow and deep 3D structures; hence incorporating only the upper 30 km of 3D structure will still improve the fit to the observed PGD values. Our results demonstrate that 3D simulations significantly improve models of GNSS waveform characteristics and will not only help understand the underlying processes, but also improve local tsunami warning.

Evaluating the Relation of Cave Passage Formation to Stress-Field: Spatio-Temporal Correlation of Speleogenesis with Active Tectonics in Asprorema Cave (Mt. Pinovo, Greece)

Caves serve as time capsules, preserving significant markers of tectonic activity and offering insights into geological history. Fault geometries and past activations found in caves can be correlated with known deformational events in the broader area, temporally delimiting the speleogenesis. More specifically, cave passage formation is suggested to be affected by the regional stress-field. The Asprorema Cave in Northern Greece is a typical example of a fracture guided cave, with passage geometry influenced by relative sidewall movements, revealing these discontinuities as faults. This study constructs the timeframe and conceptual model of speleogenesis in relation to tectonic events, geomorphological evolution and hydrological zones, and verifies its relation to the stress-field. Active tectonics, mineralogy and cave geomorphology are investigated. Results suggest syntectonic speleogenesis under phreatic and epiphreatic conditions. The absence of corrosion on fault slip surfaces implies recent activations post cave’s shift to the vadose zone. Structural analysis identifies three main neotectonic phases: NNW-SSE striking faults (oldest group of structures), NE-SW striking faults with dextral strike-slip movement (post-middle Miocene), and NE-SW striking normal faults indicating extensional stress-regime (Quartenary). The last two phases affect cave passage shape causing wall displacement, highlighting passage formation along discontinuities perpendicular to the horizontal minimum stress axis.

Variations in Soil Organic Carbon after Farmland Conversion to Apple Orchard

The Chinese Loess Plateau has undergone extensive revegetation to restore degraded land and enhance carbon sequestration. However, soil organic carbon (SOC) sequestrated in the soil profiles of deep-rooted plants has not been fully studied. Here, we investigated the SOC within a 0–23 m profile in farmlands and apple orchards converted from farmlands with different ages (A5, &lt;5 years; A10, ~10 years; A15, ~15 years; A20, &gt;20 years) and the controlling factors on three loess tablelands (Changwu, Qingyang, and Luochuan). The results show that SOC stocks among farmlands and orchards showed no significant difference (p = 0.88); however, SOC stocks showed a trend with tree ages, i.e., a decrease for A5 and A10 but an increase for A15 and A20. For the vertical variability, the SOC stock was the highest within 0–1 m, regardless of the standing age; however, the SOC stock in this layer only accounted for 8.8% of the total SOC stock (97.93 ± 9.18 kg m−2). Climate accounted for 82% of the variations and controlled the changes in SOC in the 0–1 m range, while soil texture dominated the SOC in the soil below 1 m, accounting for 57% of the SOC variations. The variations in SOC in the thick, unsaturated zones provide implications for future land use management and the sustainability of apple orchards in arid regions.

An underground drip water monitoring network to characterize rainfall recharge of groundwater at different geologies, environments, and climates across Australia

Abstract. Understanding when and why groundwater recharge occurs is of fundamental importance for the sustainable use of this essential freshwater resource for humans and ecosystems. However, accurately capturing this component of the water balance is widely acknowledged to be a major challenge. Direct physical measurements identifying when groundwater recharge is occurring are possible by utilizing a sensor network of hydrological loggers deployed in underground spaces located in the vadose zone. Through measurements of water percolating into these spaces from above, we can record the potential groundwater recharge process in action. By using automated sensors, it is possible to precisely determine when recharge occurs (which event, month, or season and for which climate condition). Combined with daily rainfall data, it is possible to quantify the “rainfall recharge threshold”, the amount of rainfall needed to generate groundwater recharge, and its temporal and spatial variability. Australia's National Groundwater Recharge Observing System (NGROS) provides the first dedicated sensor network for observing groundwater recharge at an event scale across a wide range of geologies, environments, and climate types representing a wide range of Australian hydroclimates. Utilizing tunnels, mines, caves, and other subsurface spaces located in the vadose zone, the sensors effectively record “deep drainage”, water that can move beyond the shallow subsurface and root zone to generate groundwater recharge. The NGROS has the temporal resolution to capture individual recharge events, with multiple sensors deployed at each site to constrain the heterogeneity of recharge between different flow paths, and to quantify (including uncertainty bounds) rainfall recharge thresholds. Established in 2022, the network is described here together with examples of data being generated.

Estimation of virus transport parameters in both unsaturated-saturated zone using numerical simulation and flower pollination algorithm

To develop a virus transport model considering both the unsaturated and saturated zone, it is essential to consider separate parameters for both the zone. This is because the virus fate and transport process depends on the degree of saturation. As such, this paper introduces a parameter estimation model that simultaneously identifies virus transport parameters in an unconfined aquifer for both the unsaturated and saturated zones. The proposed model utilizes a linked simulation-optimization technique to minimize errors between observed and simulated virus concentrations. The simulation model used in this study is termed as “variable detachment model,” this is because the simulation component of this approach incorporates varying adsorption rates, which are influenced by the degree of saturation. This study also highlights the complexities involved in parameter estimation for virus transport models, such as non-uniqueness, instability, and non-identifiability. These challenges arise when multiple parameter combinations can produce similar model predictions, leading to difficulty in determining the actual parameter values. To address these challenges, a meta-heuristic optimization algorithm, Flower Pollination Algorithm (FPA), is used, effectively solving such non-convex problems. The model could estimate all the parameters accurately except the inactivation rates of viruses. Later when the global sensitivity analysis was carried out to understand the relative importance of each parameter, inactivation rates were found to be less influential to the model output. Thus, FPA based parameter estimation model can be used as an alternative for estimating the virus transport parameters in groundwater aquifers.

Ubiquitous but unique: Water depth and oceanographic attributes shape methane seep communities

AbstractIn the past decade, thousands of previously unknown methane seeps have been identified on continental margins around the world. As we have come to appreciate methane seep habitats to be abundant components of marine ecosystems, we have also realized they are highly dynamic in nature. With a focus on discrete depth ranges across the Cascadia Margin, we work to further unravel the drivers of seep‐associated microbial community structure. We found highly heterogenous environments, with depth as a deterministic factor in community structure. This was associated with multiple variables that covaried with depth, including surface production, prevailing oxygen minimum zones (OMZs), and geologic and hydrographic context. Development of megafaunal seep communities appeared limited in shallow depth zones (~ 200 m). However, this effect did not extend to the structure or function of microbial communities. Siboglinid tubeworms were restricted to water depths &gt; 1000 m, and we posit this deep distribution is driven by the prevailing OMZ limiting dispersal. Microbial community composition and distribution covaried most significantly with depth, but variables including oxygen concentration, habitat type, and organic matter, as well as iron and methane concentration, also explained the distribution of the microbial seep taxa. While members of the core seep microbiome were seen across sites, there was a high abundance of microbial taxa not previously considered within the seep microbiome as well. Our work highlights the multifaceted aspects that drive community composition beyond localized methane flux and depth, where environmental diversity adds to margin biodiversity in seep systems.

Alteration mineralogy, characteristics and shortwave infrared spectroscopy of white mica in the Zijinshan ore field, Fujian Province: Implications for porphyry Cu prospecting

The Zijinshan ore field is situated in the southeastern part of China. in which contains the large Zijinshan high-sulfidation Cu-Au deposit, the Luoboling porphyry Cu-Mo deposit, the Yueyang low-sulfidation Au-Ag deposit, the Longjiangting and the Wuziqilong tansitional Cu deposit, and the Xi’nan porphyry Cu-Mo prospect. Debates on genetic relationship among these deposits, especially the relationship between the Zijinshan Cu-Au deposit and the adjacent Luoboling deposit is a question of importance for further exploration in this ore field. In this contribution, the alteration mineralogy, alteration texture, zonation and mineral chemistry, in combination with shortwave infrared spectroscopy (SWIR) analysis of white mica were used to help focus exploration and assessment of hydrothermal deposits in the Zijinshan ore field.In this study, diaspore-pyrophyllite alteration zone, chlorite-sericite alteration zone and propylitic alteration zone were first identified in the deep drill core of the Zijinshan Cu-Au deposit, which gives way upward to quartz-alunite-kaolinite zone, quartz-alunite zone and the silicic alteration zone. In the Xi’nan Cu-Mo prospect, diaspore-pyrophyllite alteration zone was first recognized in the shallow level with the chlorite-sericite alteration zone and propylitic alteration zone in the deeper level. Andalusite- sericite ± diaspore alteration zone was first recognized in the upper part of the Luoboling Cu-Mo deposit and the deeper part of the Wuziqilong Cu deposit. It indicates that a pervasive occurrence of diaspore -pyrophyllite ± andalusite alteration zones in the Zijinshan ore field, which may potentially represent a transitional environment between the top of porphyry alteration zones and the base of advanced argillic alteration zones. Minerals of the white mica group are extensively distributed in various alteration zones of the Zijinshan ore field. The SWIR analysis of white mica indicated that the wavelength position of the Al-OH (∼2,200 nm; wAlOH) absorption feature increases gradually from the shallow kaolinite-dickite alteration zone to the deep potassic alteration zone. In the advanced argillic alteration zones, the range of wAlOH values of white mica is less than 2205 nm. Conversely, in the porphyry alteration zones, it always exceeds 2205 nm. It is evident that the andalusite-sericite alteration zone with higher wAlOH values (average of 2203 nm) is always proximal to the porphyry mineralization, which may act as indicators for deep porphyry prospecting. The results of electron probe microanalysis (EPMA) reveal that the early enrichment of Mg and depletion of Fe in white mica can be attributed to the preferential incorporation of Fe into iron oxides and iron sulfides in the potassic and propylitic alteration zones. Overall, there is a negative correlation between wAlOH and Al, while wAlOH exhibits a positive correlation with the sum of Fe + Mg + Mn. The distribution characteristics of alteration zones and the variation patterns of the wAlOH of white mica indicate that there are probably several paleo-thermal centers in the deep parts of the Zijinshan ore field, located respectively in the deep parts of the Xi’nan Cu-Mo prospect and the deep part of the Zijinshan Cu-Au deposit and the deep parts of the Wuziqilong Cu deposit, and the latest seems to be more fertile and deserves further exploration.

Vadose Zone Journal Special Section: Soil physics in agricultural production, water resources, and waste management

<i>Vadose Zone Journal</i> Special Section: Soil physics in agricultural production, water resources, and waste management

The dimensions of deep-layer soil desiccation and its impact on xylem hydraulic conductivity in dryland tree plantations

Abstract. In water-limited areas, planted trees can extract substantial amounts of soil water from deep layers (&gt; 200 cm) to meet their high water demand, resulting in deep soil desiccation, which influences not only regional water cycling but also the sustainability of trees per se in drylands. However, the specific dimensions of deep soil desiccation in relation to both the soil moisture limitation and the maximum root water uptake (RWU) depth are still not well determined. Whether the dimensions depend on tree species and how they will affect trees' xylem hydraulic conductivity are also unclear, restricting our ability to predict the fate of dryland tree plantations. Therefore, we studied the spatiotemporal distribution of deep soil moisture deficit (DSMD) for two typical planted trees, apple (Malus pumila Mill.) and black locust (Robinia pseudoacacia L.), based on published data and multiple field samplings on China's Loess Plateau. The results indicated that the lowest deep soil moisture (DSM; units of gravimetric percent, i.e., grav-%) occurred under the planted trees aged 24–28 years at all sites. The lowest DSMD varied around −0.6, which was close to the DSMD at the permanent wilting point (PWP, grav-%), regardless of tree species and site, although shallow (&lt; 200 cm) soil moisture was not reduced to the point of limitation. This suggests that PWP is a reliable indicator of the moisture limitation of deep-layer soil desiccation for the tree species examined. The corresponding depth of soil moisture use reached 18.0–22.0 m for these old planted trees at different sites, while it was more than 25 m for R. pseudoacacia in the drier site of Mizhi. Furthermore, the mean values of native percentage loss of hydraulic conductivity of planted trees' branches xylem reached 74.9 %–96.5 % in the plantations sampled in this study, indicating that tree mortality may occur. The findings help predict the sustainability of planted trees in semi-arid regions with a thick vadose zone.

Characterization of a Contaminated Site Using Hydro-Geophysical Methods: From Large-Scale ERT Surface Investigations to Detailed ERT and GPR Cross-Hole Monitoring

This work presents the results of an advanced geophysical characterization of a contaminated site, where a correct understanding of the dynamics in the unsaturated zone is fundamental to evaluate the effective management of the remediation strategies. Large-scale surface electrical resistivity tomography (ERT) was used to perform a preliminary assessment of the structure in a thick unsaturated zone and to detect the presence of a thin layer of clay supporting an overlying thin perched aquifer. Discontinuities in this clay layer have an enormous impact on the infiltration processes of both water and solutes, including contaminants. In the case here presented, the technical strategy is to interrupt the continuity of the clay layer upstream of the investigated site in order to prevent most of the subsurface water flow from reaching the contaminated area. Therefore, a deep trench was dug upstream of the site and, in order to evaluate the effectiveness of this approach in facilitating water infiltration into the underlying aquifer, a forced infiltration experiment was carried out and monitored using ERT and ground-penetrating radar (GPR) measurements in a cross-hole time-lapse configuration. The results of the forced infiltration experiment are presented here, with a particular emphasis on the contribution of hydro-geophysical methods to the general understanding of the subsurface water dynamics at this complex site.