Vadose Research Articles

Hotspots of carbon release from rock weathering in the unsaturated zone

Hotspots of carbon release from rock weathering in the unsaturated zone

Assessing groundwater contamination risk in industrial zone of Ranipet district, Southern India: A modified DRASTIC and Fuzzy-AHP approach

Groundwater plays a substantial function in sustaining the water requirements of Ranipet district, Tamil Nadu, India, across a variety of domains, including drinking water, agriculture, and industry. However, rising demand and industrial activities have caused concern about groundwater pollution and susceptibility in the region. So, this study employed a modified DRASTIC model to find the regions that are most vulnerable to groundwater contamination. The model includes a novel parameter, "proximity to industry,” in addition to the seven traditional factors (depth to water level, net recharge, aquifer media, soil media, topography, impact of vadose zone and hydraulic conductivity) to account for the industrial nature of the study region. The fuzzy-analytical hierarchy process (FAHP) was used to compute parameter weights and ratings. All the parameters are overlayed in ArcGIS software to create the groundwater vulnerability map (GWVM). The map revealed that the central area along the Palar River is highly vulnerable to contamination. The southern and northern parts have medium to extremely low groundwater contamination risks. The depth of the water table was the most significant component impacting the vulnerability, as per the sensitivity analysis. The model is further validated by chromium concentrations, which were shown to be high in highly sensitive zones and low in low-risk zones. This study also suggests the need for extensive and periodic groundwater quality assessments in highly vulnerable zones for sustainable water resource management.

The added value of CoreVESS score and penetration resistance in predicting the soil–water retention curve

Soil hydraulic properties, such as the soil–water retention curve (SWRC), play a crucial role in simulating water transport within the vadose zone. They can be directly measured in the laboratory or field, or indirectly predicted by using Pedotransfer Functions (PTFs). It has been advocated that accounting for soil structure could marginally improve the prediction accuracy of PTFs in its wet range. The aim of this study was to test whether the use of an easy-to-determine soil structure-related variable could effectively increase the prediction accuracy of water retention curve PTFs. Additionally, we explored whether including soil strength, another soil property that is easy and quickly to determine in the field in many replicates, could improve PTF performance. Our investigation involved extensive sampling of 252 soil horizons across 42 cropped fields, each exhibiting within-field variations in soil structural degradation. Soil structure was represented by a soil structural quality score (Sq) taken from a semi-quantitative visual soil assessment (CoreVESS), which was included in the PTFs as predictor variable or as discriminator (Sq) for data grouping in addition to or instead of ‘classical’ predictors such as soil organic carbon content, clay and sand content, bulk density and soil layer. Penetration resistance (PR) was taken as a variable for soil strength. Various regression methods from classical regression to machine learning were evaluated on Train and Test subdatasets. To evaluate the effect of data grouping, the dataset was split in two based on structure and on texture. Both point and parametric PTFs were developed, with the latter predicting the parameters of the Van Genuchten (VG) water retention equation. Overall, the k-nearest neighbors PTFs (KNN-PTFs) achieved the best performance showing the highest R2 and lowest RMSE values. Moreover, the point PTFs showed much better and more stable predictive ability than parametric PTFs. The inclusion of PR as an additional predictor variable slightly improved the prediction performance in the wet range of the SWRC. Interestingly, Sq proved beneficial when substituted for bulk density but not when used in addition to it. Grouping data based on texture and structure improved the prediction accuracy of the PTFs, particularly evident in textural grouping and MLR, but not in KNN.

Deep desert aquifers as an archive for Mid- to Late Pleistocene hydroclimate: An example from the southeastern Mediterranean

Many efforts have been made to illuminate the nature of past hydroclimates in semi-arid and arid regions, where current and future shifts in water availability have enormous consequences on human subsistence. Deep desert aquifers, where groundwater is stored for prolonged periods, might serve as a direct record of major paleo-recharge events. To date, groundwater-based paleoclimate reconstructions have mainly focused on a relatively narrow timescale (up to ∼40 kyr), limited by the relatively short half-life of the widely used radiocarbon (5.73 kyr). Here we demonstrate the usage of deep regional aquifers in the arid southeastern Mediterranean as a hydroclimate archive for earlier Mid-to-Late Pleistocene epochs. State-of-the-art dating tools, primarily the 81Kr radioisotope (t1/2 = 229 kyr), were combined with other atmosphere-derived tracers to illuminate the impact of four distinguishable wetter episodes over the past 400 kyr, with differences in climatic conditions and paleo-recharge locations. Variations in stable water isotope composition suggest moisture transport from more proximal (Mediterranean) and distal (Atlantic) sources to different parts of the region at distinct times. Large variability in the computed noble gas-based recharge temperature (NGT), ranging ~15–30 °C, cannot be explained by climate variations solely, and points to different recharge pathways, including geothermal heating in the deep unsaturated zone and recharge from high-elevation (colder) regions. The obtained groundwater record complements and enhances the interpretation of other terrestrial archives in the arid region, including a contribution of valuable information regarding the moisture source origin as reflected in the deuterium-excess values, which is unattainable from the common practice analysis of calcitic cave deposits. We conclude that similar applications in other deep (hundred-m-order) regional groundwater systems (e.g., the Sahara desert aquifers) can significantly advance our understanding of long-term (up to 1 Myr) paleo-hydroclimate in arid regions, including places where no terrestrial remnants, such as cave, lake, and spring sediments, are available.

A self-sustaining effect induced by iron sulfide generation and reuse in pyrite-woodchip mixotrophic bioretention systems: An experimental and modeling study

Dual electron donor bioretention systems have emerged as a popular strategy to enhance dissolved nitrogen removal from stormwater runoff. Pyrite-woodchip mixotrophic bioretention systems showed a promoted and stabilized removal of dissolved nutrients under complex rainfall conditions, but the sulfate reduction process that can induce iron sulfide generation and reuse was overlooked. In this study, experiments and models were applied to investigate the effects of filler configuration and dissolved organic carbon (DOC) dissolution rate on treatment performance and iron sulfide generation in pyrite-woodchip bioretention systems. Key parameters govern that DOC dissolution and microbe-mediated processes were obtained by experiments. The water quality models that integrate one-dimensional constant flow, sorption and microbial transformation kinetics were used to predict the performance of bioretention systems. Results showed that the mixotrophic bioretention system with woodchip mixed in the vadose zone and pyrite in the saturated zone achieves a better performance in both nitrogen removal efficiency and by-product control. Comparably, woodchip and pyrite mixed in the saturated zone could encounter a high secondary pollution risk. The sensitivity coefficients of oxic/anoxic DOC dissolution rates to total nitrogen removal are 0.36 and -2.43 respectively. Iron sulfide generation was affected by DOC distribution and the competition between heterotrophic denitrifiers, autotrophic denitrifiers, and sulfate-reducing bacteria (SRB). DOC accumulation has an antagonistic effect on iron production and sulfate reduction. Extra DOC accumulation favors sulfate reduction while high DOC concentration inhibits pyrite-based denitrification and reduces Fe(III) production. The recycling of iron sulfide can improve the robustness and sustainability of bioretention systems.

Migration of depleted uranium from a corroded penetrator in soil vadose zone in Bosnia and Herzegovina

Depleted uranium (DU) from corroded armor penetrators can migrate through the soil vadose zone and cause environmental problems, yet studies on such migration at former theatres of war are scarce. Here, we investigated vertical DU migration in a soil profile due to a penetrator (3–8 cm beneath the soil surface) corroded over 7 years in Bosnia and Herzegovina. The highest concentration of DU was ∼45,300 mg/kg at 6–10 cm, with the concentration decreasing markedly with increasing depth. The majority of the DU accumulated within the top 20 cm and the DU front reached ∼42 cm beneath the penetrator. In addition, particles with varying U concentrations (3–65 wt%) were observed at 0–15 cm, with U primarily co-located with O, Si, Al, maghemite, and hematite. Particularly, metaschoepite was identified at 6–10 cm. Finally, X-ray absorption spectroscopy analysis found U was hexavalent in the soil profile. These findings suggest that the downward migration of DU was likely present as a soluble form adsorbed on clay minerals and Fe oxides. Overall, we show that the rate of DU migration within the vadose zone is comparatively slow, although if the penetrator is left in the soil for decades, it could pose a serious long-term risk. Environmental ImplicationsOver 90 % of the depleted uranium (DU) penetrators fired in previous conflicts missed their armored targets and were left in the soil to corrode. The corroded penetrators can not only contaminate soil but also pose a risk to groundwater. The present study examined the migration of DU in a soil profile that included a DU penetrator that had been corroding for over 7 years. Studying the dynamics of DU migration is essential to develop effective remediation strategies to mitigate long-term environmental risks and safeguard ecosystems and human health from DU contamination.

Fault rock properties and conditions produce variance in slip during earthquake rupture propagation at the Nankai Trough

Although drilled samples of fault rocks have yielded information on frictional features of shallow subduction zones, the relationship of rupture propagation to the levels of friction and pore-fluid pressure remains uncertain. To investigate this topic, we performed dynamic rupture simulations along the megasplay fault that slipped during the 1944 Mw 8.0 Tonankai earthquake in the Nankai Trough. We used actual data from friction experiments on rocks from the fault segment and pre-existing pore pressures deduced from geophysical surveys for the shallow portion of 0–10 km depth along the fault. Simulations of low friction (friction coefficient ca. 0.04) produced large slip (about 30 m), whereas simulations using higher friction (friction coefficient ca. 0.2) suppressed the rupture. In simulations with low friction in which the pore-fluid pressure was nearly equal to the lithostatic stress, the slip decreased to about 25 m. However, when the simulations included slip-strengthening at shallow depth and higher friction, the slip still reached roughly 20 m. Such variability in slip during rupture propagation is caused by differences in the friction features and fluid pressure conditions of fault rocks, in which the friction features might be related to the mineral composition. Spatiotemporal heterogeneity in fault-rock type and their physical and hydraulic properties may fundamentally produce the complexity and variability of earthquake rupture propagation along the Nankai plate-subduction boundary.

Unravelling the bridge of polyphenol composition and mineralomics in Se biofortified edible Allium species

Selenium is an essential element for the well-being of humans and animals. Micronutrient deficiencies affect over half of the global population. Biofortification of crops with selenium represents an economically sustainable solution to combat this problem. In this study, seven Allium species underwent Se biofortification, via Se soil enrichment with selenate in pots, at three distinct concentration levels and were used to explore the relationship between Se plant internalization, polyphenol composition, and mineral content. Allium species were chosen because they are recognized as secondary accumulators of Se. Plant morphology, as well as the presence of fifteen phenolic acids and flavonoids, were analysed using a UHPLC-DAD system. Mineral content was determined using both ICP-OES and FAES. Se levels were quantified using a validated method based on HG-HR-CS-QT-AAS. Even at the lowest biofortification level, Allium species demonstrated potential as sources of dietary Se, offering viable alternatives to Se supplements, without negatively impacting the main mineral content. Higher biofortification levels led to reduced levels of flavonoids and phenolic acids, accompanied by potential decreases in biomass for certain species. The study highlighted how Se concentration affects key polyphenols like quercetin, gallic acid, and luteolin, although the response varied among species.

Assessment of Soil Horizons and Their Matric Potential from Ground-Penetrating Radar Signal Attributes

Soil plays significant roles in different phases and in the continuous existence of human life. Its comprehensive knowledge, particularly as related to its physical characteristics, enhances its utilization, conservation, and management. The traditional methods of soil study are characterized with some pitfalls such as much time needed to perform such assessments. There are also issues of invasiveness that affect the soil structures and discrete sampling that may not reflect true spatial attributes in the outcome of such techniques. These problems are largely due to the concealing nature of soil layers that made its thorough evaluation difficult. In this study, an alternative geophysical approach has been adopted. The technique is the ground-penetrating method (GPR) that utilizes electromagnetic pulse energy via its equipment’s sensors, which can allow the investigation of soil properties, even in its concealing state. This study aimed at qualitatively evaluating the soil horizons and the matric potentials using the GPR signal attributes within the unsaturated zone with a view of having insight into the test field’s characterization. Field data measurements were obtained using MALA ProEX GPR equipment with its accessories manufactured by MALA Geosciences, Stockholm, Sweden. Evaluation of the processed field data results and computed attributes show soil characteristics variations with depth that was interpreted as the layers. This can be seen from the GPR data presentation as an image representing the subsurface of the zones of propagation of the pulse energy. Spectral analysis of the GPR signals allows for the delineation of two zones of contrasting features, which were tagged as high and low matric potentials. Although the conventional direct measurement of the matric potential was not made at the time of the study to complement and confirm the veracity of the approach, the results indicate the possibility of the approach towards a quick and in situ technique of soil investigations. Such evaluation may be valuable input in precision agriculture where accurate data are sought for implementation.

A stable localized weak strong form radial basis function method for modelling variably saturated groundwater flow induced by pumping and injection

The unsaturated zone profoundly affects groundwater (GW) flow induced by pumping and injection due to the capillary forces. However, current radial basis function (RBF) numerical models for GW pumping and injection mostly ignore the unsaturated zone. To bridge this gap, we developed a new three-dimensional weak strong form RBF model in this study, called CCHE3D-GW-RBF. Flow processes were modelled by the mixed-form Richards equation which was iteratively solved by the modified Picard iteration. Soil-water characteristic curves were represented by the widely applicable formulas, the van Genuchten (1980) model. Differential operators were approximated by the localized Gaussian RBF, and the weighted singular value decomposition method was used to construct stable bases. The injection/pumping wells and the flux boundaries were handled by the weak formulation using Meshless Local Petrov Galerkin method, and the strong-form equation using the collocation RBF method was enforced on the other points. Good agreement was found between the simulation results from our numerical model and the well-accepted solutions in all three verification cases, demonstrating the accuracy and applicability of this model. In addition, a smaller RBF shape parameter was found to produce a more accurate modelling resulting, indicating the necessity of implementing stable bases for RBF models.

Depth dimension in seepage detection: Insights for exploration and geological gas storage surveillance

The search for surface hydrocarbon seeps has historically served as a fundamental tool in the early stages of oil and gas exploration. However, as easily detectable surface hydrocarbon seeps become increasingly scarce, advancement in analytical capabilities have shifted the focus towards detecting ultra-trace levels of hydrocarbons, commonly referred to as microseepages. While it is critical to understand the distribution of hydrocarbon microseepage for assessing petroleum systems and optimizing exploration efforts, the role of sampling depth in the detection of microseepage remains poorly understood. In this study, we investigated the influence of sampling depth on microseepage detection, using passive soil-gas sampling techniques at depths ranging from 0.5 to 10 m at sites with hydrocarbon-bearing reservoirs and dry wells. Organic geochemical analysis, including thermal desorption-gas chromatography-mass spectrometry (GC-MS), was employed to identify and quantify hydrocarbon and nonhydrocarbon compounds.The study result demonstrates that sampling at the Minimal Interference Zone (MIZ) depths of 5 m or greater significantly mitigates the risk of interference signals, thereby enabling the identification and quantification of hydrocarbon microseepages at trace and ultra-trace concentrations. At these depths, the attenuation of non-hydrocarbon compounds was observed, making the hydrocarbon signal markedly more distinguishable. Geochemical characterization of the samples revealed spatial and vertical heterogeneity in the distribution of hydrocarbon compounds, indicative of variations in the microseepage mechanisms potentially influenced by factors such as sedimentary basin settings and secondary alterations associated with the vadose zone. By utilizing passive soil-gas sampling techniques coupled with GC-MS, this study emphasizes the importance of sampling depth to improve the detection of reliable microseepage signal while reducing interference and background noise. The proposed technique not only refines existing methodologies for petroleum system assessments but also paves the way for broader geological applications. Specifically, the approach offers critical insights for the development of high-resolution geochemical surveys adaptable for surface background profiling in CO2 sequestration projects to monitor for storage efficiency and seal integrity.

Climatic Effects on Vapor Flow and Behavior in the Vadose Zone

AbstractThe concentrations and transport of volatile organic compounds (VOCs) and other vapors in the vadose zone may exhibit some degree of temporal variability due to the effect of various climatic factors, including (1) Air temperature; (2) Barometric pressure; (3) Surface winds; and (4) Soil moisture, including the effects of any water infiltration and/or changes in groundwater level. These variables may directly affect the rates of gas transport through the vadose zone or may indirectly affect transport by changing the soil‐gas concentrations at a given location and depth. To understand the potential effect of these factors due to climate change, it is first necessary to understand their effect over typical time periods of one to several days, seasonally, and annually. In this paper, the effects of the above variables over various time periods are presented and the long‐term effects due to climate change are discussed. Standard approaches for soil‐gas measurement attempt to account for these variables, either to negate their potential influence or to capture data under reasonably worst‐case conditions. The appropriateness and adequacy of typical soil vapor measurement approaches are discussed.

Distribution of Thanetian-Ilerdian biozones in Lorestan Zone, West of Iran

The Taleh Zang Formation is a carbonate sequence of the late Paleocene–early Eocene age in the Lorestan zone, west of Iran. This Formation overlays siliciclastic rocks of the Amiran Formation and is overlain by red sandstone of the Kashkan Formation. A biostratigraphic scheme, based on the larger benthic foraminifera, is proposed for the Thanetian–Ilerdian deposits in Lorestan Zone. We compare our results with the available data from other parts of Tethys and the zonation of shallow-marine Paleocene–Eocene of Mediterranean Tethys is made. In total, 367 samples were collected from the six sections. From a biochronostratigraphic point of view, three biozones and one indeterminate zone are identified. The zones are shallow benthic zone 3 (early Thanetian), shallow benthic zone 4 (late Thanetian), and shallow benthic zone 6 (early Ilerdian). In the northern part of the Lorestan Zone (Kabutar Bula and Psht-e Jangal sections), the Taleh Zang Formation contains poor fauna. Due to the stratigraphic position, an early Thanetian age can be considered. Our data showed that the late Paleocene and Eocene biozones in Lorestan Zone developed parallel to the axis of the folding, in a northwest-southeast trend. They get younger from the northeast to the southwest. Keywords: biostratigraphy, Eocene, Foraminifera, Lorestan Zone, Paleocene.

Multiscale discrete fracture network modelling of shallow-water carbonates: East Agri Valley Basin, Southern Italy

We investigate the multiscale geometrical and dimensional properties of fracture and fault networks crosscutting Lower Jurassic to Cretaceous shallow-water carbonate rocks. The carbonates are exposed along the flanks of the Viggiano Mt., on the eastern side of the High Agri Valley Basin of southern Italy. Furthermore, we employ the results of field and digital structural analyses to derive the input parameters for subsequent Discrete Fracture Network modelling (DFN) of geocellular volumes whose dimension and architecture resemble those of the investigated sites. DFN modelling is carried out for geocellular volumes representing the studied bed packages (5m-side volumes), outcrops (50m-side volumes), and reservoir-scale carbonate cliffs (500m-side volumes). Notably, modelling is obtained considering the minimum mechanical aperture value obtained in the field for porosity computations and theoretical values of fracture hydraulic aperture for permeability computations. This is because the mechanical aperture values and roughness profiles collected in the field along single fractures are unreliable due to weathering and localized karst-related dissolution. Our findings reveal that a scale-dependent geometry characterizes the Cretaceous carbonates over three orders of magnitude. These results support previously published data, and contrast with those obtained for the Lower Jurassic carbonates, which suffer of some bias due to the quality of the exposures. After DFN modelling, we show that the amount of computed fracture porosity is mainly due to the SB and NSB fractures, and that equivalent permeability is greater within faulted rock volumes. In terms of horizontal permeability, which is the most reliable result after DFN modelling, we document near isotropic horizontal permeability ellipses at all scales of observations. At individual outcrops, we note that the permeability ellipses are elongated parallel to the dominant fault sets that denotes how localized strain affects the directionality of fluid flow. Finally, we summarize the original data in a conceptual model illustrating the modalities of meteoric fluid infiltration through the vadose zone, and then subsequent horizontal flow in the phreatic zone present within the fracture carbonates at shallow depths.

Intra-annual variations in isotopic composition of atmospheric water vapour in a semi-arid monsoon region: observations and inferences

The isotopic composition of atmospheric water vapour dynamically changes due to the complex interplay between the moisture-contributing sources (precipitation, evaporating soil moisture, transpiration, evaporating water bodies, etc.), operating ecohydrological processes, and the climatic conditions. In the present study, the isotopic composition of different sources of water was sampled; air moisture (n = 47) twice a month for the period December, 2021 to January, 2023; transpiration samples (n = 76); evaporation and rain samples (n = 12 each). The samples were analyzed using the mass spectrometry to investigate the dynamic nature of the isotopic composition of the atmospheric vapour and the factors influencing it. The research findings showed that the isotopic composition of atmospheric moisture rapidly depleted after the onset of the monsoon, indicating the replacement of pre-monsoon local vapour by monsoon moisture of depleted isotopic composition. The process of evaporation from the soil surface during the post-monsoon period was found to progress parallel to other soil hydrological processes, such as capillary action and diffusion-driven downward flux of enriched water molecules. This resulted in the broadening of an isotopically enriched layer of soil moisture in the vadose zone with time, and accordingly, there was transpiration of this enriched soil moisture from the deeper layers of soils with time by plants/trees. Evaporation of soil water and the transpiration of plants and trees of enriched water vapour gradually altered the isotopic composition of the atmospheric water vapour composition from the depleted isotopic composition of the monsoon period to an enriched isotopic composition. The transpiration of enriched water vapour from deep-rooted trees indicated more prolonged arid climate conditions. The findings of the study are useful in understanding the regional hydrological and hydroclimatic processes which can be implemented in agro-climatic studies and water resource management.

Microbial community structure and function in a sandy soil column under dynamic hydrologic regimes

The capillary fringe (CF) exhibits steep chemical and redox gradients, representing a critical reactive interface in the subsurface. The spatial extent of the CF is temporally dynamic because of water table and soil moisture fluctuations. To simulate the spatio-temporal responses of the CF’s microbial community structure and activity, a 91 cm long sandy soil column experiment was conducted by imposing three different sequential hydrologic regimes: a stable water table separating the unsaturated and saturated zones (static regime), periodic pulse infiltration events with a constant water table position (infiltration regime), and water table fluctuations under three drainage-imbibition cycles of 6, 12, and 18 days (fluctuating regime). The investigation focused on the microbially-mediated turnover of nitrogen (N) and carbon (C). Geochemical profiles in the soil column were monitored continuously, while microbial DNA was extracted for 16S rDNA sequencing at the start, middle and end of the experiment. The observed microbial community compositions corresponded to five categories: (1) initial soil, (2) near-surface and (3) CF and saturated soil under the infiltration regime, (4) unsaturated and (5) saturated soil in the fluctuating regime. Based on co-occurrence network analysis, microbial physiologies matched the predominant geochemical conditions, and relatively stable community structures prevailed throughout the experiment. NO3- and NH4+ were more available to the microbes in the fluctuating regime’s alternating oxic/anoxic zone, leading to a more flexible network. The prevalent microbial metabolic pathways predicted using FAPROTAX were mainly associated with N and C metabolisms and significantly correlated with the taxonomic compositions. However, the dominant microbial groups were shared among all the soil samples. The limited spatial differences in microbial community structure under the fluctuating regime were likely due to a more pronounced role of microbial activity than microbial diversity, and similar drying-wetting conditions between the soil column and the filed site. Internally consistent geochemical and microbial depth stratifications were identified by (1) the redox potential and nitrogen species distributions, (2) spatial moment analysis of NO3- and NH4+ concentrations, (3) microbial community composition and (4) function. The five imposed hydrologic regimes generated distinct and coupled geochemical and microbial signatures. These signatures resulted from variations in the soil microbial community’s activity, rather than diversity or abundance.

Effect of earth-air on water transport in the vadose zone of the loess plateau

Effect of earth-air on water transport in the vadose zone of the loess plateau

Using a comparative of DRASTIC and Bayesian weights of evidence approach to assess transboundary aquifer vulnerability in a data scarcity region: Tuli-Karoo aquifer

Study regionThe Tuli-Karoo Transboundary Aquifer system is shared among Botswana, South Africa and Zimbabwe. It hosts a complex aquifer predominately of sand, heavily tapped in South Africa and sandstone substantially used in Botswana threatening it’s quality. Study focusTwo vulnerability models, the classic DRASTIC (Depth to water table [D], Recharge to the aquifer [R], Aquifer media [A], Soil media [S], Topography [T], Impact of vadose zone [I], Conductivity [C]) and the bayesian based, Weights of Evidence (WofE), were used to delineate zones vulnerable to groundwater contamination under limited nitrate observations. New hydrological insights for the regionThe study revealed that the generic DRASTIC vulnerability results were comparable to the WofE, for the highly vulnerable zones, as both had an area under the curve (AUC) below 0.6. The high vulnerability was attributed to the shallow depth of water level, fractured geological condition (based on transmissivity) and high net recharge rate in and around the Tuli-Karoo Aquifer. However, there are differences in spatial variability in the low/medium vulnerability zones. WofE was superior in capturing spatial variability of groundwater vulnerability. Considerations must therefore be provided, under data scarcity, to incorporate proxies that include landuse-landcover and population density as evidential layers. Overall, WofE performs better due to better prediction power under data-scarcity conditions. This study recommends that under data scarcity, high predictive models like WofE should be used for environmental planning and groundwater protection strategies.

Quartz luminescence sensitivity enhanced by residence time in the critical zone

The emerging use of quartz luminescence properties to characterize Earth-surface processes shows promise, with optically stimulated luminescence (OSL) sensitivity proposed as a valuable tool for provenance or sediment history tracing. However, the geologic processes that lead to quartz sensitization remain unclear. Here we study the impact of source rock and surface processes on the luminescence properties of quartz sand from bedrock and modern and Late Pleistocene alluvium generated from a mountainous catchment in northern Utah, USA. Continuous wave and linear modulated OSL are used to characterize the luminescence sensitivity and intensity of the fast-decay component. We compare the OSL sensitivity with sand-grain provenance and with proxies for surface processes such as topographic metrics, cosmogenic 10Be-derived erosion rates, chemical weathering indices, and magnetic susceptibility. Late Pleistocene sediment has low OSL sensitivity and a weak fast-decay component, similar to bedrock samples from the source area. In contrast, modern alluvium is dominated by the fast-decay component and has higher and more variable OSL sensitivity, with no clear relationship to upstream bedrock source. There is, however, an inverse relationship between OSL sensitivity and catchment-averaged erosion rates and a positive relationship with chemical weathering indices and magnetic susceptibility. These metrics suggest that the modern alluvium has experienced increased residence time in the shallow critical zone compared to the Late Pleistocene sediments. We suggest that changes in hillslope processes between the effectively wetter, cooler Pleistocene and the dryer, warmer conditions of the Holocene enhanced the luminescence properties. The results suggest that climatic controls on rates and processes of chemical and mechanical weathering and sediment transport and residence within the critical zone are encoded in the luminescence properties of quartz sand.

Preliminary Validation of Persulfate Delivery Through Porous Pipe in the Vadose Zone: A Two-Dimensional Sandbox Experiment

ABSTRACT When remediating contamination in the vadose zone, which is a potential source of continuous groundwater and surface water pollution, technical challenges are often encountered. In this study, the feasibility of delivering sodium persulfate (SPS) oxidant to unsaturated soils was explored by integrating porous pipes, commonly used in agriculture, with in-situ chemical oxidation (ISCO). An evaluation of the operating characteristics of the pipe determined a recommended influent rate of < 0.1 L min−1 and an initial influent water pressure < 20 kPa to prevent percolation pore enlargement. Two-dimensional sandbox experiments revealed that pulse infiltration at an influent rate of 0.05 L min−1 was effective in creating matrix flow and lateral movement, minimizing preferential flow induced by continuous infiltration. Application of SPS (10% by wt.) under these conditions resulted in approximately 90% wetting of the sandbox area SPS flow at the experimental time that infiltrated from the pipe. After a 14-d post infiltration resting state, residual SPS concentrations ranged from 6,000 -12,000 mg/kg in soil, with soil pH levels of 3–4, ORP around 600 mV, and 40–60% decreases in soil organic carbon. This developed method helps retain the oxidant and promote chemical reactions, making ISCO a viable option for remediating unsaturated soil contamination.