Impact of climate change on groundwater recharge in shallow young glacial aquifers in northern Poland

We investigated the influence of climate change in the period 1951-2020 on shallow aquifers in the Brda and Wda outwash plains (Pomeranian Region, Northern Poland). There was a significant temperature rise (0.3 °C/10 years), which accelerated after 1980 (0.66 °C/10 years). Precipitation became increasingly irregular - extremely rainy years occurred right after or before extremely dry years, and intensive rainfall events became more frequent after 2000. The groundwater level decreased over the last 20 years, even though the average annual precipitation was higher than in the previous 50 years. We carried out numerical simulations of water flow in representative soil profiles for the years 1970-2020 using the HYDRUS-1D model, developed and calibrated during our earlier work at an experimental site in the Brda outwash plain (Gumuła-Kawęcka et al., 2022). We used a relationship between the water head and flux at the bottom of the soil profiles (the third-type boundary condition) to reproduce groundwater table fluctuations caused by recharge variability in time. The calculated daily recharge showed a decreasing linear trend for the last 20 years (0.05-0.06 mm d^-1/10 years), and dropping trends in water table level and soil water content in the entire profile of vadose zone. Field tracer experiments were performed to estimate impact of extremely rain events on water flux in vadose zone. The results suggest that tracer travel times are strongly determined by water content in the unsaturated zone which is determined by precipitation amount in span of weeks, rather than extremely high precipitation events.

Coupling time-lapse ground penetrating radar surveys and infiltration experiments to characterize two types of non-uniform flow

Understanding linkages between heterogeneous soil structures and non-uniform flow is fundamental for interpreting infiltration processes and improving hydrological simulations. Here, we utilized ground-penetrating radar (GPR) as a non-invasive technique to investigate those linkages and to complement current traditional methods that are labor-intensive, invasive, and non-repeatable. We combined time-lapse GPR surveys with different types of infiltration experiments to create three-dimensional (3D) diagrams of the wetting dynamics. We carried out the GPR surveys and validated them with in situ observations, independent measurements and field excavations at two experimental sites. Those sites were selected to represent different mechanisms that generate non-uniform flow: (1) preferential water infiltration initiated by tree trunk and root systems; and (2) lateral subsurface flow due to soil layering. Results revealed links between different types of soil heterogeneity and non-uniform flow. The first experimental site provided evidence of root-induced preferential flow paths along coarse roots, emphasizing the important role of coarse roots in facilitating preferential water movement through the subsurface. The second experimental site showed that water infiltrated through the restrictive layer mainly following the plant root system. The presented approach offers a non-invasive, repeatable and accurate way to detect non-uniform flow.

Field measurement of effects of individual and combined application of biochar and polyacrylamide on erosion variables in loess and marl soils

Controlling soil erosion, especially in its initial stages, is greatly important in natural resources management. Consequently, the present research aimed to control splash and interrill erosion in two soil types (marl at Marzan-Abad and loess at Maraveh-Tapeh sites in northern Iran) using biochar (BC) and polyacrylamide (PAM). We established 0.5 × 0.5-m plots and applied BC (800 g·m^-2), PAM (2 g·m^-2), and BC + PAM (800 g·m^-2 + 2 g·m^-2) with control plots and three replications on a slope of ~25%. We used a rainfall simulator to achieve rainfall intensity of 50 mm·h^-1 with 30-min duration in the experiments. Analysis of the results obtained from the variables of splash and interrill erosion during the rainfall-runoff process showed that the PAM significantly (p ≤ 0.05) increased all study variables of splash erosion. For interrill erosion, it reduced the variables of soil loss and sediment concentration. However, the difference was not significant (p > 0.05) compared to the control plot and runoff from the two treatment sites increased relative to that from the control plots. The plot treated with BC showed decreased runoff volume, runoff coefficient, and soil loss compared to the control plot at the Marzan-Abad site, but the differences were not statistically significant (p > 0.05). However, the plot in which loess soil was treated with BC at the Maraveh-Tapeh site exhibited considerably (p ≤ 0.05) increased runoff and soil loss compared to the control plot. The entire results verified a wide range for benefit reduction of study treatments from +25.09 to -37.49% for runoff and from +38.59 to -231% for soil loss with more effectiveness for Maraveh-Tapeh Loess soil as well as combined application of BC and PAM. These findings contribute to improved understanding of proper application of soil amendments to control runoff and soil loss in loam and loess soils.

Detecting infiltrated water and preferential flow pathways through time-lapse ground-penetrating radar surveys

The objective of this paper was to identify the incidence and extent of preferential flow at two experimental areas located in Lyon, France. We used time-lapse ground-penetrating radar (GPR) surveys in conjunction with automatized single-ring infiltration experiments to create three-dimensional (3D) representations of infiltrated water. In total we established three 100 cm × 100 cm GPR grids and used differenced radargrams from pre- and post-infiltration surveys to detect wetting patterns. The analyzed time-lapse GPR surveys revealed the linkage between nonuniform flow and heterogeneous soil structures and plant roots. At the first experimental area, subsurface coarse gravels acted as capillary barriers that concentrated flow into narrow pathways via funneled flow. At the second experimental area, the interpolated 3D patterns closely matched direct observation of dyed patterns, thereby validating the applied protocol. They also highlighted the important role of plant roots in facilitating preferential water movement through the subsurface. The protocol presented in this study represents a valuable tool for improving the hydraulic characterization of highly heterogeneous soils, while also alleviating some of the excessive experimental efforts currently needed to detect preferential flow pathways in the field.

Experimental and modeling of the unsaturated transports of S-metolachlor and its metabolites in glaciofluvial vadose zone solids

The transport of pesticides to groundwater is assumed to be impacted by flow processes and geochemical interactions occurring in the vadose zone. In this study, the transport of S-metolachlor (SMOC) and its two metabolites ESA-metolachlor (MESA) and OXA-metolachlor (MOXA) in vadose zone materials of a glaciofluvial aquifer is studied at laboratory scale. Column experiments are used to study the leaching of a conservative tracer (bromide) and SMOC, MESA and MOXA under unsaturated conditions in two lithofacies, a bimodal gravel (Gcm,b) and a sand (S-x). Tracer experiments showed water fractionation into mobile and immobile compartments more pronounced in bimodal gravel columns. In both lithofacies columns, SMOC outflow is delayed (retardation factor>2) and mass balance reveals depletion (mass balance of 0.59 and 0.77 in bimodal gravel and sand, respectively). However, complete mass elution associated with retardation factors close to unity shows that there is no adsorption of MESA and MOXA in either lithofacies. SMOC transport is characterized by non-equilibrium sorption and sink term in both bimodal gravel and sand columns. Batch experiments carried out using agitation times consistent with column water residence times confirmed a time-dependence of SMOC sorption and high adsorption rates (>80%) of applied concentrations. Desorption experiments confirm the irreversibility of a major part of the SMOC adsorption onto particles, corresponding to the sink term in columns. In the bimodal gravel column, SMOC adsorption occurs mainly on reactive particles in contact with mobile water because of flow regionalization whereas in the sand column, there is pesticide diffusion to the immobile water. Such results clearly show that sorption mechanisms in the vadose zone solids below the soil are both solute and contact-time-dependent and are impacted by hydrodynamic conditions. The more rapid transport of MESA and MOXA to the aquifer would be controlled mainly by water flow through the unsaturated zone whereas SMOC transport is retarded by sorption processes within the vadose zone.

Glyphosate and AMPA adsorption in soils: laboratory experiments and pedotransfer rules

Adsorption of the herbicide glyphosate and its main metabolite AMPA (aminomethylphosphonic acid) was investigated on 17 different agricultural soils. Batch equilibration adsorption data are shown by Freundlich adsorption isotherms. Glyphosate adsorption is clearly affected by equilibration concentrations, but the nonlinear AMPA adsorption isotherms indicate saturation of the adsorption sites with increasing equilibrium concentrations. pHCaCl2 (i.e. experimental pH) is the major parameter governing glyphosate and AMPA adsorption in soils. However, considering pHCaCl2 values, available phosphate amount, and amorphous iron and aluminium oxide contents by using a nonlinear multiple regression equation, obtains the most accurate and powerful pedotransfer rule for predicting the adsorption constants for these two molecules. As amorphous iron and aluminium oxide contents in soil are not systematically determined, we also propose a pedotransfer rule with two variables-pHCaCl2 values and available phosphate amount-that remains acceptable for both molecules. Moreover, the use of the commonly measured pHwater or pHKCl values gives less accurate results compared to pHCaCl2 measurements. To our knowledge, this study is the first AMPA adsorption characterization for a significant number of temperate climate soils.

Combined effect of capillary barrier and layered slope on water, solute and nanoparticle transfer in an unsaturated soil at lysimeter scale

It is well recognized that colloidal nanoparticles are highly mobile in soils and can facilitate the transport of contaminants through the vadose zone. This work presents the combined effect of the capillary barrier and soil layer slope on the transport of water, bromide and nanoparticles through an unsaturated soil. Experiments were performed in a lysimeter (1×1×1.6m(3)) called LUGH (Lysimeter for Urban Groundwater Hydrology). The LUGH has 15 outputs that identify the temporal and spatial evolution of water flow, solute flux and nanoparticles in relation to the soil surface conditions and the 3D system configuration. Two different soil structures were set up in the lysimeter. The first structure comprises a layer of sand (0-0.2cm, in diameter) 35cm thick placed horizontally above a layer of bimodal mixture also 35cm thick to create a capillary barrier at the interface between the sand and bimodal material. The bimodal material is composed of a mixture 50% by weight of sand and gravel (0.4-1.1cm, in diameter). The second structure, using the same amount of sand and bimodal mixture as the first structure represents an interface with a 25% slope. A 3D numerical model based on Richards equation for flow and the convection dispersion equations coupled with a mechanical module for nanoparticle trapping was developed. The results showed that under the effect of the capillary barrier, water accumulated at the interface of the two materials. The sloped structure deflects flow in contrast to the structure with zero slope. Approximately 80% of nanoparticles are retained in the lysimeter, with a greater retention at the interface of two materials. Finally, the model makes a good reproduction of physical mechanisms observed and appears to be a useful tool for identifying key processes leading to a better understanding of the effect of capillary barrier on nanoparticle transfer in an unsaturated heterogeneous soil.