Pedotransfer functions (PTFs) are widely used to estimate soil hydraulic properties (SHPs) from easily measurable characteristics. However, most existing PTFs rely on unimodal hydraulic models, which fail to accurately represent the bimodal SHPs caused by soil structure common in field conditions. In this study, we developed new PTFs using two bimodal soil hydraulic models and introduced soil physics-informed neural networks (SPINN) to embed the models into PTF training. The results showed that the new PTFs effectively captured bimodality in hydraulic conductivity curves, achieving an RMSE of 0.578 in the test set, compared to 0.709 for unimodal models. The PTFs also improved soil water retention curve (SWRC) predictions but struggled with bimodal SWRCs for some samples, likely due to the limited number of bimodal SWRCs in the dataset. An independent dataset evaluation revealed that the RMSE for hydraulic conductivity predicted by the new PTFs was approximately one-third of that of classic PTFs. This underscores the significant role of soil structure in SHPs, which classic PTFs fail to capture. Additionally, PTFs developed using the SPINN method outperformed those optimized fitted hydraulic parameters via machine learning, a common approach in the literature. We also found that separate versus simultaneous optimization of water retention and hydraulic conductivity greatly affects PTF performance. Finally, we provided global 1 km-resolution maps of soil hydraulic parameters for the bimodal model.

Here, we review in depth how soils can remember moisture anomalies across spatial and temporal scales, embedded in the concept of soil moisture memory (SMM), and we explain the mechanisms and factors that initiate and control SMM. Specifically, we explore external and internal drivers that affect SMM, including extremes, atmospheric variables, anthropogenic activities, soil and vegetation properties, soil hydrologic processes, and groundwater dynamics. We analyze how SMM considerations should affect sampling frequency and data source collection. We discuss the impact of SMM on weather variability, land surface energy balance, extreme events (drought, wildfire, and flood), water use efficiency, and biogeochemical cycles. We also discuss the effects of SMM on various land surface processes, focusing on the coupling between soil moisture, water and energy balance, vegetation dynamics, and feedback on the atmosphere. We address the spatiotemporal variability of SMM and how it is affected by seasonal variation, location, and soil depth. Regarding the representation and integration of SMM in land surface models, we provide insights on how to improve predictions and parameterizations in LSMs and address model complexity issues. The possible use of satellite observations for identifying and quantifying SMM is also explored, emphasizing the need for greater temporal frequency, spatial resolution, and coverage of measurements. We provide guidance for further research and practical applications by providing a comprehensive definition of SMM, considering its multifaceted perspective.

A document by Shuying Chen. Click on the document to view its contents.

Soil hydraulic properties (SHPs) are impacted by various mechanisms such as soil structure, capillarity, and adsorption forces, often showing a bimodal shape. Developing soil hydraulic models that describe SHPs over the entire saturation range often involves balancing the representation of multiple processes while minimizing the number of free-fitted parameters. Existing soil hydraulic models rarely capture bimodal soil hydraulic properties across the full moisture range or introduce too many free-fitted parameters. In this study, we propose a novel framework to describe SHPs over the entire moisture range, accounting for the effects of soil structure, capillarity, adsorption forces, and vapor diffusion. In its four free-fitted parameters form, the proposed models can capture unimodal soil water retention curves (SWRC) and bimodal hydraulic conductivity curves (HCC). With one additional free-fitted parameter, the proposed models can capture both bimodal SWRC and HCC. Testing with 107 and 52 soil samples from two public datasets demonstrated that the proposed models performed exceptionally well in describing SHPs across the entire moisture range. The reported lowest root-mean-square error values were 0.005 and 0.009 cm3 cm−3 for fitting SWRCs, and 0.465 and 0.666 for predicting HCCs, respectively. Due to the minimal introduction of free-fitted parameters, the proposed framework showed significant application potential.

The need to develop and provide integrated observation systems to better understand and manage global and regional environmental change is one of the major challenges facing Earth system science today. In 2008, the German Helmholtz Association took up this challenge and launched the German research infrastructure TERrestrial ENvironmental Observatories (TERENO). The aim of TERENO is the establishment and maintenance of a network of observatories as a basis for an interdisciplinary and long-term research programme to investigate the effects of global environmental change on terrestrial ecosystems and their socio-economic consequences. State-of-the-art methods from the field of environmental monitoring, geophysics, remote sensing, and modelling are used to record and analyze states and fluxes in different environmental disciplines from groundwater through the vadose zone, surface water, and biosphere, up to the lower atmosphere. Over the past 15 years we have collectively gained experience in operating a long-term observing network, thereby overcoming unexpected operational and institutional challenges, exceeding expectations, and facilitating new research. Today, the TERENO network is a key pillar for environmental modelling and forecasting in Germany, an information hub for practitioners and policy stakeholders in agriculture, forestry, and water management at regional to national levels, a nucleus for international collaboration, academic training and scientific outreach, an important anchor for large-scale experiments, and a trigger for methodological innovation and technological progress. This article describes TERENO’s key services and functions, presents the main lessons learned from this 15-year effort, and emphasises the need to continue long-term integrated environmental monitoring programmes in the future.

The assessment of soil health has evolved from focusing primary on agricultural productivity to an integrated evaluation of soil biota and biotic processes that impact soil properties. Consequently, soil health assessment has shifted from a predominantly physico-chemical approach to incorporating ecological, biological and molecular microbiology methods. These methods enable a comprehensive exploration of soil microbial community properties and their responses to environmental changes arising from climate change and anthropogenic disturbances. Despite the increasing availability of soil health indicators (physical, chemical, and biological), a holistic mechanistic linkage between indicators and soil functions across multiple spatiotemporal scales has not yet been fully established. This article reviews the state-of-the-art of soil health monitoring, focusing on understanding how soil-microbiome-plant processes contribute to feedback mechanisms and causes of changes in soil properties, as well as the impact these changes have on soil functions. Furthermore, we survey the opportunities afforded by the soil-plant digital twin approach, an integrative framework that amalgamates process-based models, Earth Observation data, data assimilation, and physics-informed machine learning, to achieve a nuanced comprehension of soil health. This review delineates the prospective trajectory for monitoring soil health by embracing a digital twin approach to systematically observe and model the soil-plant system. We further identify gaps and opportunities, and provide perspectives for future research for an enhanced understanding of the intricate interplay between soil properties, soil hydrological processes, soil-plant hydraulics, soil microbiomes, and landscape genomics.

Cosmic ray neutron sensors (CRNS) allow to determine field-scale soil moisture content non-invasively due to the dependence of aboveground measured epithermal neutrons on the amount of hydrogen. Because other pools besides soil contain hydrogen (e.g. biomass), it is necessary to consider these for accurate soil moisture content measurements, especially when they are changing dynamically (e.g., arable crops, de- and reforestation). In this study, we compare four approaches for the correction of biomass effects on soil moisture content measurements with CRNS using experiments with three crops (sugar beet, winter wheat and maize) on similar soils: I) site-specific functions based on in-situ measured biomass, II) a generic approach, III) the thermal-to-epithermal neutron ratio (Nr) and IV) the thermal neutron intensity. Calibration of the CRNS during bare soil conditions resulted in root mean square errors (RMSE) of 0.097, 0.041 and 0.019 m3/m3 between estimated and reference soil moisture content of the cropped soils, respectively. Considering in-situ measured biomass for correction reduced the RMSE to 0.015, 0.018 and 0.009 m3/m3. When thermal neutron intensity was considered for correction, similarly accurate results were obtained. Corrections based on Nr and the generic approach were less accurate. We also explored the use of CRNS for biomass estimation. The use of Nr only provided accurate biomass estimates for sugar beet. However, significant site-specific relationships between biomass and thermal neutron intensity were obtained for all three crops. It was concluded that thermal neutron intensity can be used to correct soil moisture content estimates from CRNS and to estimate biomass.

Modelling of the land surface water-, energy-, and carbon balance provides insight into the behaviour of the Earth System, under current and future conditions. Currently, there exists a substantial variability between model outputs, for a range of model types, whereby differences between model input parameters could be an important reason. For large-scale land surface, hydrological, and crop models, soil hydraulic properties (SHP) are required as inputs, which are estimated from pedotransfer functions (PTFs). To analyse the functional sensitivity of widely used PTFs, the water fluxes for different scenarios using HYDRUS-1D was simulated and predictions compared. The results showed that using different PTFs causes substantial variability in predicted fluxes. In addition, an in-depth analysis of the soil SHPs and derived soil characteristics was performed to analyse why the SHPs estimated from the different PTFs cause the model to behave differently. The results obtained provide guidelines for the selection of PTFs in large scale models. The model performance in terms of numerical stability, time-integrated behaviour of cumulative fluxes, as well as instantaneous fluxes was evaluated, in order to compare the suitability of the PTFs. Based on this, the Rosetta, Wösten, and Tóth PTF seem to be the most robust PTFs for the Mualem van Genuchten SHPs and the PTF of Cosby et al. (1984) for the Brooks Corey functions. Based on our findings, we strongly recommend to harmonize the PTFs used in model inter-comparison studies to avoid artefacts originating from the choice of PTF rather from different model structures.

The movement and spreading of contaminated groundwater plumes and their mixing with non-contaminated water is strongly influenced by the heterogeneity of the aquifer properties, which may vary strongly over small spatial scales. Thus, imaging these small-scale features and monitoring transport of tracer plumes at a fine resolution is of interest to characterize transport processes in aquifers. Full-waveform inversion (FWI) of crosshole ground penetrating radar (GPR) measurements can provide an aquifer characterization at decimeter-scale resolution. The method produces images of both relative dielectric permittivity (εr) and bulk electrical conductivity (σb), which related to hydraulic aquifer properties and tracer distributions. To test the potential of time-lapse GPR FWI for imaging tracer plumes, we conducted a numerical experiment of tracer transport in a heterogeneous aquifer. Concentration was converted to saline and desalinated tracers, which changed σb, and to ethanol, which changed both εr and σb. The simulated εr and σb distributions in a crosshole plane were considered to simulate GPR data. These data were subsequently used to reconstruct εr and σb distributions using the crosshole 2D GPR FWI. Tracer concentrations were retrieved from the inverted εr and σb models using information about petrophysical parameters. GPR FWI εr images could recover preferential paths of ~0.2 m width, while the σb images resolved structures up to ~ 0.2-0.3 m. The results highlight that changes in εr, e.g., ethanol and hot water, can be used to image transport processes with high resolution by time-lapse GPR FWI, while the accuracy of the recovery of σb is limited.

In his seminal paper on solution of the infiltration equation, Philip (1957) proposed a gravity time, tgrav, to estimate practical convergence time of his infinite time series expansion, TSE. The parameter tgrav refers to a point in time where infiltration is dominated equally by capillarity and gravity derived from the first two (dominant) terms of the TSE expansion. Evidence that higher order TSE terms describe the infiltration process better for longer times. Since the conceptual definition of tgrav is valid regardless of the infiltration model used, we opted to reformulate tgrav using the analytic approximation proposed by Parlange et al. (1982) valid for all times. In addition to the roles of soil sorptivity (S) and saturated (Ks) and initial (Ki) hydraulic conductivities, we explored effects of a soil specific shape parameter β on the behavior of tgrav. We show that the reformulated tgrav (notably tgrav= F(β) S^2/(Ks - Ki)^2 where F(β) is a β-dependent function) is about 3 times larger than the classical tgrav given by tgrav, Philip= S^2/(Ks - Ki)^2. The differences between original tgrav, Philip and the revised tgrav increase for fine textured soils. Results show that the proposed tgrav is a better indicator for convergence time than tgrav, Philip. For attainment of the steady-state infiltration, both time parameters are suitable for coarse-textured soils, but not for fine-textured soils for which tgrav is too conservative and tgrav, Philip too short. Using tgrav will improve predictions of the soil hydraulic parameters (particularly Ks) from infiltration data as compared to tgrav, Philip.

The advance of the cosmic ray neutron (CRN) sensing method for estimating field scale soil moisture relied largely on simulations of the footprint properties of epithermal neutrons (~0.5 eV - 100 keV). Commercially available CRN probes are usually additionally equipped with a thermal neutron (< 0.5 eV) detector. The potential of these measurements is rarely explored because relevant features of thermal neutrons, such as the footprint and the sensitivity to soil moisture are unknown. Here, we used neutron transport modeling and a river crossing experiment to assess the thermal neutron footprint. We found that the horizontal thermal neutron footprint ranges between 43 and 48 m distance from the probe and that the vertical footprint extends to soil depths between 10 and 65 cm depending on soil moisture. Furthermore, we derived weighting functions that quantify the footprint characteristics of thermal neutrons. These results will enable new applications of thermal neutrons.