Search

Water uptake by plant roots is an important component of the soil water balance. Predicting to what extent potential transpiration from the canopy, that is, transpiration demand, can be met by supply of water from the soil through the root system is crucial to simulate the actual transpiration and assess vegetation water stress. In models that simulate the dynamics of vertical soil water content profiles as a function of water fluxes and soil water potential gradients, the root water uptake (RWU) distribution is represented by macroscopic sink terms. We present RWU functions that calculate sink terms based on a mechanistic model of water flow in the soil–root system. Based on soil–root hydraulics, we define α‐supply functions representing the maximal uptake by the root system from a certain soil depth when the root collar water potential equals the wilting point, ω‐supply factors representing the maximal supply from the entire root system, and a critical ωc factor representing the potential transpiration demand. These functions and factors are subsequently used to calculate RWU distributions directly from potential transpiration or demand and the soil water potentials. Unlike currently used approaches, which define α‐stress functions and ω factors representing ratios of actual uptake to uptake demand, the supply‐based formulations can be derived directly from soil and root hydraulic properties and can represent processes like root hydraulic redistribution and hydraulic lift. Core Ideas: Functions to calculate root water uptake (RWU) from different soil depths were derived.RWU is calculated directly using soil water potentials and the potential transpiration rate.The functions are derived using soil and root hydraulic properties.Phenomena like RWU compensation, hydraulic redistribution, and hydraulic lift are reproduced.Similarities and differences with currently used uptake functions are discussed. [ABSTRACT FROM AUTHOR]

The two‐term infiltration equation I(t)=St+At$I(t) = S\sqrt {t} + A t$ is commonly used to determine the sorptivity, S$S$(LT−1/2)$({\text{LT}}^{-1/2})$, and product, A=cKs$A = c K_{\text{s}}$(LT−1)$({\text{LT}}^{-1})$, of the dimensionless multiple c$c$ and saturated soil hydraulic conductivity Ks$K_{\text{s}}$(LT−1)$({\text{LT}}^{-1})$ from cumulative vertical infiltration measurements I∼1,...,I∼n$\tilde{I}_{1},\ldots,\tilde{I}_{n}$ (L) at times t1,...,tn$t_{1},\ldots,t_{n}$ (T). This reduced form of the quasi‐analytical power series solution of Richardson's equation of Philip enjoys a solid physical underpinning but at the expense of a limited time validity. Using simulated infiltration data, Jaiswal et al. have shown this time validity to equal about 2.5 cm of cumulative infiltration. The goals of this work are twofold. First, we investigate the extent to which cumulative infiltration measurements larger than 2.5 cm bias the estimates of S$S$ and Ks$K_{\text{s}}$. Second, we investigate the impact of epistemic errors on the inferred time validities and parameters. Partial infiltration curves up to 2.5 cm of cumulative vertical infiltration improve substantially the agreement between actual and least squares estimates of S$S$ and Ks$K_{\text{s}}$. But this only holds if the data generating infiltration process follows Richardson's equation and experimental conditions satisfy assumptions of soil homogeneity and a uniform initial water content. Otherwise, autocorrelated cumulative infiltration residuals will bias the least squares estimates of S$S$ and Ks$K_{\text{s}}$. Our findings reiterate and reinvigorate earlier conclusions of Haverkamp et al. and show that epistemic errors deteriorate the physical significance of the coefficients of infiltration functions. As a result, the parameters of infiltration functions cannot simply be used in storm water and vadose zone flow models to forecast runoff and recharge at field and landscape scales unless these predictions are accompanied by realistic uncertainty bounds. We conclude that the time validity of Philip's two‐term equation is an elusive theoretical quantity with arbitrary physical meaning. Core Ideas: Structural model errors of infiltration functions bias the values of S$S$ and Ks$K_s$.Least squares estimates of Ks$K_s$ compare poorly to measured values.Epistemic uncertainty turns coefficients of infiltration functions into fitting parameters.The time validity of Philip's two‐term equation is an elusive theoretical quantity.Treatment of epistemic errors enhances usefulness of plot‐scale infiltration experiments for large‐scale prediction. [ABSTRACT FROM AUTHOR]

Electromagnetic (EM) coupling effects including both inductive and capacitive coupling have long been an essential problem in broad-band spectral electrical impedance tomography (sEIT) measurements at the field scale. Efforts have been made to remove EM coupling numerically or to suppress the effects by modified data acquisition strategies. For near-surface applications with relatively small survey layouts, inductive coupling can be well removed in the mHz to kHz frequency range. With the use of shielded coaxial cables and so-called active electrodes where the amplifiers are mounted at the electrodes, capacitive coupling in sEIT measurements can also be reduced. However, it remains challenging to cope with capacitive coupling between the cable shield and the ground, especially in resistive field conditions. The aim of this study is to deal with this type of capacitive coupling effect by identifying and filtering out sEIT measurements that are strongly affected by capacitive coupling. Based on a correction method for capacitive coupling proposed in a previous study, an approach to estimate measurement errors due to capacitive coupling is presented first. In the second step, a workflow was proposed to calculate the capacitive coupling strength (CCS) for each electrode configuration, which is defined as the ratio of the imaginary part of the impedance induced by capacitive coupling and the imaginary part of the impedance due to the subsurface electrical conductivity. In the final step, measurements with low CCS were selected for inversion and the results were compared with inversion results obtained using the previously developed correction approach. It was found that the filtering method based on CCS is more capable in tackling capacitive coupling compared to using model-based corrections. Spectrally consistent sEIT results up to kHz were obtained using the newly developed filtering method, which were not achieved in previous work using model-based correction. [ABSTRACT FROM AUTHOR]

Soil moisture is an essential climate variable of the Earth system. Understanding its spatiotemporal dynamics is essential for predicting weather patterns and climate variability, monitoring and mitigating the effects and occurrence of droughts and floods, improving irrigation in agricultural areas, and sustainably managing water resources. Here we review in depth how soils can remember information on soil moisture anomalies over time, as embedded in the concept of soil moisture memory (SMM). We explain the mechanisms underlying SMM and explore its external and internal drivers; we also discuss the impacts of SMM on different land surface processes, focusing on soil‐plant‐atmosphere coupling. We explore the spatiotemporal variability, seasonality, locality, and depth‐dependence of SMM and provide insights into both improving its characterization in land surface models and using satellite observations to quantify it. Finally, we offer guidance for further research on SMM. Plain Language Summary: Our review paper takes an in‐depth look at soil moisture memory, which is how soil records its moisture history over time and space. Analogous to human psychology, which seeks to understand how a person's/society's memory influences his/her present and future behavior, understanding soil moisture memory encourages consideration of how such memory determines present state and might determine future behavior of soils exposed to environmental disturbances. Soil moisture memory can be affected by a variety of factors, both external (e.g., weather extremes) and internal (soil's unique properties). It affects everything from the air to the way our landscapes respond to disasters like droughts, wildfires, and floods. We also studied how this phenomenon affects the balance of water and energy in our environment, the health of our plants, and even how it communicates with the atmosphere. We show how it can change depending on where you are on the planet, the time of year, and how deep you dig into the soil. We offer scientists insights into how weather and land surface models can become more accurate by accounting for soil moisture memory. Its understanding not only helps us predict and manage our environment, but also provides opportunities for exciting scientific discoveries. Key Points: Atmospheric forcings, land use and management, and soil processes and mechanisms explain how and why soil moisture memory emerges in ecosystemsNonlocality of moisture memory, its spread across different regions, and its interaction with large‐scale climate phenomena are underexploredFurther advances in land surface models and closer integration of model simulations and observations are needed to better characterize moisture memory [ABSTRACT FROM AUTHOR]

Plant water uptake and plant and soil water status are important for the soil water balance and plant growth. They depend on atmospheric water demand and the accessibility of soil water to plant roots, which is in turn related to the hydraulic properties of the root system and the soil around root segments. We present a simulation model that describes water flow in the soil–plant system mechanistically considering both root and soil hydraulic properties. We developed an approach to upscale three‐dimensional (3D) flow in the soil toward root segments of a 3D root architecture to a model that considers one‐dimensional flow between horizontal soil layers and radial flow to root segments in that layer. The upscaled model couples upscaled linear flow equations in the root system with an analytical solution of the nonlinear radial flow equation between the soil and roots. The upscaled model avoids simplifying assumptions about root hydraulic properties and water potential drops near roots made in, respectively, soil‐ and root‐centered models. Xylem water potentials and soil–root interface potentials are explicitly simulated and show, respectively, large variations with depth and large deviations from bulk soil water potentials under dry soil conditions. Accounting for hydraulic gradients in the soil around root segments led to an earlier but slower reduction of transpiration during a drought period and a better plant water status with higher nighttime plant water potentials. Core Ideas: Root water uptake depends on root and soil hydraulic properties.Water uptake at root element scale was upscaled to the root system scale.The upscaled model can be implemented in one‐dimensional soil water flow models.Low conductance of dry soil prevents low nighttime plant water potentials. [ABSTRACT FROM AUTHOR]

Renewable energy is recognized in Africa as a means for climate change mitigation, but also to provide access to electricity in sub‐Saharan Africa, where three‐quarters of the global population without electricity resides. Reliable and highly resolved renewable energy potential (REP) information is indispensable to support power plants expansion. Existing atmospheric data sets over Africa that are used for REP estimates are often characterized by data gaps, or coarse resolution. With the aim to overcome these challenges, the ICOsahedral Nonhydrostatic (ICON) Numerical Weather Prediction (ICON‐NWP) model in its Limited Area Mode (ICON‐LAM) is implemented and run over southern Africa in a hindcast dynamical downscaling setup at a convection‐permitting 3.3 km horizontal resolution. The simulation time span covers contrasting solar and wind weather years from 2017 to 2019. To assess the suitability of the novel simulations for REP estimates, the simulated hourly 10 m wind speed (sfcWind) and hourly surface solar irradiance (rsds) are extensively evaluated against a large compilation of in situ observations, satellite, and composite data products. ICON‐LAM reproduces the spatial patterns, temporal evolution, the variability, and absolute values of sfcWind sufficiently well, albeit with a slight overestimation and a mean bias (mean error (ME)) of 1.12 m s−1 over land. Likewise the simulated rsds with an ME of 50 W m−2 well resembles the observations. This new ICON simulation data product will be the basis for ensuing REP estimates that will be compared with existing lower resolution data sets. Plain Language Summary: There are still approximately 580 million people in Africa without reliable electricity supply. Renewable energy is broadly accepted as an important solution for Africa to fill the power supply gap and to mitigate climate change. Renewable energy potential (REP) information is thereby imperative for the expansion of renewable energy power planning. With conventional REP estimates, challenges are often linked to the meteorological input data, due to either the relatively coarse spatial resolution, data gaps in space and time, or data quality in general. In this study, we implemented and ran the atmospheric model ICOsahedral Nonhydrostatic (ICON) from the German Weather Service and partners at 3 km high‐resolution. The renewable energy variables wind speed and solar irradiance from these simulations are evaluated against an extensive in situ observations data set, as well as satellite, and other composite data products. In a comparison with more than 200 stations from three different in situ observation networks, it can be shown that ICON can reproduce REP‐related variables with a level of sophistication that the data is likely to offer added value over conventional inputs to REP assessments. The study is an example on how numerical models can fill in gaps in data‐scarce regions to produce useable information. Key Points: A new convection‐permitting regional ICOsahedral Nonhydrostatic (ICON) model setup over southern Africa is presented and evaluatedThe spatially and temporally highly resolved ICON wind and solar fields are inputs for improved renewable energy potential estimatesICON outputs agree well with hourly in situ observations and satellite data; wind speeds are slightly overestimated [ABSTRACT FROM AUTHOR]

This comment challenges Gao et al. (2023)'s perspective rejecting the role of soil processes in hydrology. We argue that the authors present a false dichotomy between soil-centric and ecosystem-centric views. These two views of hydrology are complementary and reflect on the inherent multiscale complexity of hydrology where soil processes dominate at certain scales but other processes may become important at catchment scale. We recognize the need for a new scale aware framework that reconciles the interplay between soil processes at small scales with emergent behaviors driven by vegetation, topography and climate at large scales. [ABSTRACT FROM AUTHOR]

This comment challenges Gao et al. (2023)'s perspective rejecting the role of soil processes in hydrology. We argue that the authors present a false dichotomy between soil-centric and ecosystem-centric views. These two views of hydrology are complementary and reflect on the inherent multiscale complexity of hydrology where soil processes dominate at certain scales but other processes may become important at catchment scale. We recognize the need for a new scale aware framework that reconciles the interplay between soil processes at small scales with emergent behaviors driven by vegetation, topography and climate at large scales. [ABSTRACT FROM AUTHOR]

Data assimilation plays a dual role in advancing the "scientific" understanding and serving as an "engineering tool" for the Earth system sciences. Land data assimilation (LDA) has evolved into a distinct discipline within geophysics, facilitating the harmonization of theory and data and allowing land models and observations to complement and constrain each other. Over recent decades, substantial progress has been made in the theory, methodology, and application of LDA, necessitating a holistic and in‐depth exploration of its full spectrum. Here, we present a thorough review elucidating the theoretical and methodological developments in LDA and its distinctive features. This encompasses breakthroughs in addressing strong nonlinearities in land surface processes, exploring the potential of machine learning approaches in data assimilation, quantifying uncertainties arising from multiscale spatial correlation, and simultaneously estimating model states and parameters. LDA has proven successful in enhancing the understanding and prediction of various land surface processes (including soil moisture, snow, evapotranspiration, streamflow, groundwater, irrigation and land surface temperature), particularly within the realms of water and energy cycles. This review outlines the development of global, regional, and catchment‐scale LDA systems and software platforms, proposing grand challenges of generating land reanalysis and advancing coupled land‒atmosphere DA. We lastly highlight the opportunities to expand the applications of LDA from pure geophysical systems to coupled natural and human systems by ingesting a deluge of Earth observation and social sensing data. The paper synthesizes current LDA knowledge and provides a steppingstone for its future development, particularly in promoting dual driven theory‐data land processes studies. Plain Language Summary: Land Data Assimilation (LDA) integrates numerical models with observation data to enhance predictions of key variables related to land surface processes, including soil moisture, snow, evapotranspiration, and groundwater. Consequently, LDA effectively tackles two of the fundamental scientific challenges in the Earth system: how to enhance the usability of exponentially increasing land observations and how to improve the predictive capability and accuracy of land models. Significant advancements have been made in the theory, methods, and applications of LDA since the 21st century. Breakthroughs have been attained in characterizing uncertainties, estimating multiple variables and parameters simultaneously, and incorporating big data analytics. LDA systems at global, regional, and catchment‐scales, along with widely used software platforms, have been developed. However, future challenges facing LDA comprise improving meteorological forcing data, creating long‐term land reanalysis, devising operational applications, and broadening its range to cover the critical zone. In the era of big data, LDA will evolve further by assimilating big Earth data and incorporating machine learning to develop digital twins from pure geophysical systems to coupled natural and human systems. Key Points: Land data assimilation (LDA) advances scientific understanding and serves as an engineering tool for Earth system sciencesLDA reflects the trend of harmonizing theory and data in the era of big data and artificial intelligenceFuture LDA research should expand the applications from pure geophysical systems to coupled natural and human systems [ABSTRACT FROM AUTHOR]

Land surface models (LSMs) are an important tool for advancing our knowledge of the Earth system. LSMs are constantly improved to represent the various terrestrial processes in more detail. High-quality data, freely available from various observation networks, are being used to improve the prediction of terrestrial states and fluxes of water and energy. To optimize LSMs with observations, data assimilation methods and tools have been developed in the past decades. We apply the coupled Community Land Model version 5 (CLM5) and Parallel Data Assimilation Framework (PDAF) system (CLM5-PDAF) for 13 forest field sites throughout Europe covering different climate zones. The goal of this study is to assimilate in situ soil moisture measurements into CLM5 to improve the modeled evapotranspiration fluxes. The modeled fluxes will be evaluated using the predicted evapotranspiration fluxes with eddy covariance (EC) systems. Most of the sites use point-scale measurements from sensors placed in the ground; however, for three of the forest sites we use soil water content data from cosmic-ray neutron sensors, which have a measurement scale closer to the typical land surface model grid scale and EC footprint. Our results show that while data assimilation reduced the root-mean-square error for soil water content on average by 56 % to 64 %, the root-mean-square error for the evapotranspiration estimation is increased by 4 %. This finding indicates that only improving the soil water content (SWC) estimation of state-of-the-art LSMs such as CLM5 is not sufficient to improve evapotranspiration estimates for forest sites. To improve evapotranspiration estimates, it is also necessary to consider the representation of leaf area index (LAI) in magnitude and timing, as well as uncertainties in water uptake by roots and vegetation parameters. [ABSTRACT FROM AUTHOR]

The continuous monitoring of soil water content is commonly carried out using low-frequency capacitance sensors that require a site-specific calibration to relate sensor readings to apparent dielectric bulk permittivity (Kb) and soil water content (θ). In fine-textured soils, the conversion of Kb to θ is still challenging due to temperature effects on the bound water fraction associated with clay mineral surfaces, which is disregarded in factory calibrations. Here, a multi-point calibration approach accounts for temperature effects on two soils with medium to high clay content. A calibration strategy was developed using repacked soil samples in which the Kb-θ relationship was determined for temperature (T) steps from 10 to 40 °C. This approach was tested using the GS3 and TEROS-12 sensors (METER Group, Inc. Pullman, WA, USA; formerly Decagon Devices). Kb is influenced by T in both soils with contrasting T-Kb relationships. The measured data were fitted using a linear function θ = a K b + b with temperature-dependent coefficients a and b. The slope, a(T), and intercept, b(T), of the loam soil were different from the ones of the clay soil. The consideration of a temperature correction resulted in low RMSE values, ranging from 0.007 to 0.033 cm3 cm−3, which were lower than the RMSE values obtained from factory calibration (0.046 to 0.11 cm3 cm−3). However, each experiment was replicated only twice using two different sensors. Sensor-to-sensor variability effects were thus ignored in this study and will be systematically investigated in a future study. Finally, the applicability of the proposed calibration method was tested at two experimental sites. The spatial-average θ from a network of GS3 sensors based on the new calibration fairly agreed with the independent area-wide θ from the Cosmic Ray Neutron Sensor (CRNS). This study provided a temperature-corrected calibration to increase the accuracy of commercial sensors, especially under dry conditions, at two experimental sites. [ABSTRACT FROM AUTHOR]

Land surface modelling runs with CLM5 over Africa at 3 km resolution were carried out and we assessed the impact of different sources of soil information and different upscaling strategies of the soil information, also in combination with different atmospheric forcings and different temporal resolutions of those atmospheric forcings. FAO and SoilGrids250m were used as soil information. SoilGrids information at 250 m resolution was upscaled to the 3 km grid scale by three different methods: (i) random selection of one of the small SoilGrids250m grid cells contained in the model grid cell; (ii) arithmetic averaging of SoilGrids soil texture values and (iii) selection of the dominant soil texture. These different soil model inputs were combined with different atmospheric forcings model inputs, which provide inputs at different temporal resolutions: CRUNCEPv7 (6-hourly input resolution), GSWPv3 (3-hourly) and WFDE5 (hourly). We found that varying the source of soil texture information (FAO or SoilGrids250m) influences model water balance outputs more than the upscaling methodology of the soil texture maps. However, for high temporal resolution of atmospheric forcings (WFDE5) the different soil texture upscaling methods result in large differences in simulated evapotranspiration, surface runoff and subsurface runoff at the local and regional scales related to the higher resolution representation of rainfall intensity in the model. The upscaling methodology of fine scale soil texture information influences land surface model simulation results, but only clearly in combination with high temporal resolution atmospheric forcings. [ABSTRACT FROM AUTHOR]

The dynamic interaction between subsurface flow and soil mechanics is often simplified in the stability assessment of variably saturated landslide-prone hillslopes. The aim of this study is to analyze the impact of conventional simplifications in coupling and modeling strategies on stability assessment of such hillslopes in response to precipitation using the local factor of safety (LFS) concept. More specifically, it investigates (1) the impact of neglecting poroelasticity, (2) transitioning from full coupling between hydrological and mechanical models to sequential coupling, and (3) reducing the two-phase flow system to a one-phase flow system (Richards' equation). Two rainfall scenarios, with the same total amount of rainfall but two different relatively high (4 mm h−1) and low (1 mm h−1) intensities are considered. The simulation results of the simplified approaches are compared to a comprehensive, fully coupled poroelastic hydromechanical model with a two-phase flow system. It was found that the most significant difference from the comprehensive model occurs in areas experiencing the most transient changes due to rainfall infiltration in all three simplified models. Among these simplifications, the transformation of the two-phase flow system to a one-phase flow system showed the most pronounced impact on the simulated local factor of safety (LFS), with a maximum increase of +21.5% observed at the end of the high-intensity rainfall event. Conversely, using a rigid soil without poroelasticity or employing a sequential coupling approach with no iteration between hydromechanical parameters has a relatively minor effect on the simulated LFS, resulting in maximum increases of +2.0% and +1.9%, respectively. In summary, all three simplified models yield LFS results that are reasonably consistent with the comprehensive poroelastic fully coupled model with two-phase flow, but simulations are more computationally efficient when utilizing a rigid porous media and one-phase flow based on Richards' equation. [ABSTRACT FROM AUTHOR]

Non‐invasive imaging of processes within the soil–plant continuum, particularly root and soil water distributions, can help optimize agricultural practices such as irrigation and fertilization. In this study, in‐situ time‐lapse horizontal crosshole ground penetrating radar (GPR) measurements and root images were collected over three maize crop growing seasons at two minirhizotron facilities (Selhausen, Germany). Root development and GPR permittivity were monitored at six depths (0.1–1.2 m) for different treatments within two soil types. We processed these data in a new way that gave us the information of the "trend‐corrected spatial permittivity deviation of vegetated field," allowing us to investigate whether the presence of roots increases the variability of GPR permittivity in the soil. This removed the main non‐root‐related influencing factors: static influences, such as soil heterogeneities and rhizotube deviations, and dynamic effects, such as seasonal moisture changes. This trend‐corrected spatial permittivity deviation showed a clear increase during the growing season, which could be linked with a similar increase in root volume fraction. Additionally, the corresponding probability density functions of the permittivity variability were derived and cross‐correlated with the root volume fraction, resulting in a coefficient of determination (R2) above 0.5 for 23 out of 46 correlation pairs. Although both facilities had different soil types and compaction levels, they had similar numbers of good correlations. A possible explanation for the observed correlation is that the presence of roots causes a redistribution of soil water, and therefore an increase in soil water variability. Core Ideas: Analying spatial and temporal belowground GPR signal variability and root images for maize crop roots.Root images and GPR data show differences for different treatments.Linking of root volume fraction and GPR permittivity variability. [ABSTRACT FROM AUTHOR]

An accurate assessment of the physical states of the Earth system is an essential component of many scientific, societal, and economical considerations. These assessments are becoming an increasingly challenging computational task since we aim to resolve models with high resolutions in space and time, to consider complex coupled partial differential equations, and to estimate uncertainties, which often requires many realizations. Machine learning methods are becoming a very popular method for the construction of surrogate models to address these computational issues. However, they also face major challenges in producing explainable, scalable, interpretable, and robust models. In this paper, we evaluate the perspectives of geoscience applications of physics-based machine learning, which combines physics-based and data-driven methods to overcome the limitations of each approach taken alone. Through three designated examples (from the fields of geothermal energy, geodynamics, and hydrology), we show that the non-intrusive reduced-basis method as a physics-based machine learning approach is able to produce highly precise surrogate models that are explainable, scalable, interpretable, and robust. [ABSTRACT FROM AUTHOR]

Electrical impedance tomography (EIT) is a promising method to image the frequency-dependent complex electrical conductivity distribution of the subsurface in the mHz to kHz frequency range. In contrast to the well-developed electrical resistivity tomography (ERT) method, the inversion approach for EIT data is less established. Different inversion strategies have been proposed, but the implications of the differences between these methods have not been investigated yet. In this study, we aim to compare four different inversion strategies for EIT measurements. The first strategy (CVI) formulates the inverse problem in the complex number domain and is mathematically the most elegant method. The second strategy (RVI) is the established real-valued inversion method, which decouples the inversion of the real and imaginary parts and completely ignores the complex nature. The third strategy (ALT) is very similar to the RVI strategy in case of small phase angles, but it considers the complex coupling in the forward operator and alternately updates the real and imaginary parts of the model in the case of large phase angles. The fourth and final strategy (CVI+) was newly formulated in this study. It fully considers the complex nature of EIT measurements but separates the treatment of the real and imaginary part in terms of the data weighting and regularization. The different inversion strategies were tested with two synthetic models. The first model has a small phase contrast and the second model has a large phase contrast. In the case of a small phase contrast, the CVI strategy was able to resolve the distribution of electrical conductivity amplitude, but the inversion result for the phase angle was less reliable. The other three strategies presented similar results and the models were well resolved within the expected data misfit. In the case of a model with large phase contrast, only the newly formulated CVI + strategy was able to produce reliable results. It was found that the extremely large phase angle can have a significant influence on the modelled amplitude of data. The cross-sensitivity (i.e. the imaginary part of the sensitivity) that describes the influence on the real part of data due to a change in the imaginary part of model, or that on the imaginary part of data due to a change in the real part of model, provided unique information during the inversion. It was concluded that the CVI + strategy is theoretically the most comprehensive and correct approach for EIT inversion, but that in the case of small phase angles the RVI strategy has the practical advantage that no complex calculations are required, which substantially reduces the required computational effort. [ABSTRACT FROM AUTHOR]

Integrated long-term, in-situ observations are needed to document ongoing environmental change, to "ground-truth" remote sensing and model outputs and to predict future Earth system behaviour. The scientific and societal value of in-situ observations increases with site representativeness, temporal duration, number of parameters measured and comparability within and across sites. Research Infrastructures (RIs) can support harmonised, cross-site data collection, curation and publication. Integrating RI networks through site co-location and standardised observation methods can help answers three questions about the terrestrial carbon sink: (i) What are present and future carbon sequestration rates in northern European forests? (ii) How are these rates controlled? (iii) Why do the observed patterns exist? Here, we present a conceptual model for RI co-location and highlight potential insights into the terrestrial carbon sink achievable when long-term in-situ Earth observation sites participate in multiple RI networks (e.g., ICOS and eLTER). Finally, we offer recommendations to promote RI co-location. [ABSTRACT FROM AUTHOR]

Copyright of Vadose Zone Journal is the property of Wiley-Blackwell and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

In most of the world, conditions conducive to wildfires are becoming more prevalent. Net carbon emissions from wildfires contribute to a positive climate feedback that needs to be monitored, quantified, and predicted. Here we use a causal inference approach to evaluate the influence of top-down weather and bottom-up fuel precursors on wildfires. The top-down dominance on wildfires is more widespread than bottom-up dominance, accounting for 73.3% and 26.7% of regions, respectively. The top-down precursors dominate in the tropical rainforests, mid-latitudes, and eastern Siberian boreal forests. The bottom-up precursors dominate in North American and European boreal forests, and African and Australian savannahs. Our study identifies areas where wildfires are governed by fuel conditions and hence where fuel management practices may be more effective. Moreover, our study also highlights that top-down and bottom-up precursors show complementary wildfire predictability across timescales. Seasonal or interannual predictions are feasible in regions where bottom-up precursors dominate. This paper shows that weather and fuel precursors show complementary predictability of wildfires extending across different timescales, which may be leveraged for seasonal or interannual wildfire prediction. [ABSTRACT FROM AUTHOR]

The production of crops secure the human food supply, but climate change is bringing new challenges. Dynamic plant growth and corresponding environmental data are required to uncover phenotypic crop responses to the changing environment. There are many datasets on above-ground organs of crops, but roots and the surrounding soil are rarely the subject of longer term studies. Here, we present what we believe to be the first comprehensive collection of root and soil data, obtained at two minirhizotron facilities located close together that have the same local climate but differ in soil type. Both facilities have 7m-long horizontal tubes at several depths that were used for crosshole ground-penetrating radar and minirhizotron camera systems. Soil sensors provide observations at a high temporal and spatial resolution. The ongoing measurements cover five years of maize and wheat trials, including drought stress treatments and crop mixtures. We make the processed data available for use in investigating the processes within the soil–plant continuum and the root images to develop and compare image analysis methods. [ABSTRACT FROM AUTHOR]

Long-range weather forecasts provide predictions of atmospheric, ocean and land surface conditions that can potentially be used in land surface and hydrological models to predict the water and energy status of the land surface or in crop growth models to predict yield for water resources or agricultural planning. However, the coarse spatial and temporal resolutions of available forecast products have hindered their widespread use in such modelling applications, which usually require high-resolution input data. In this study, we applied sub-seasonal (up to 4 months) and seasonal (7 months) weather forecasts from the latest European Centre for Medium-Range Weather Forecasts (ECMWF) seasonal forecasting system (SEAS5) in a land surface modelling approach using the Community Land Model version 5.0 (CLM5). Simulations were conducted for 2017–2020 forced with sub-seasonal and seasonal weather forecasts over two different domains with contrasting climate and cropping conditions: the German state of North Rhine-Westphalia (DE-NRW) and the Australian state of Victoria (AUS-VIC). We found that, after pre-processing of the forecast products (i.e. temporal downscaling of precipitation and incoming short-wave radiation), the simulations forced with seasonal and sub-seasonal forecasts were able to provide a model output that was very close to the reference simulation results forced by reanalysis data (the mean annual crop yield showed maximum differences of 0.28 and 0.36 t ha -1 for AUS-VIC and DE-NRW respectively). Differences between seasonal and sub-seasonal experiments were insignificant. The forecast experiments were able to satisfactorily capture recorded inter-annual variations of crop yield. In addition, they also reproduced the generally higher inter-annual differences in crop yield across the AUS-VIC domain (approximately 50 % inter-annual differences in recorded yields and up to 17 % inter-annual differences in simulated yields) compared to the DE-NRW domain (approximately 15 % inter-annual differences in recorded yields and up to 5 % in simulated yields). The high- and low-yield seasons (2020 and 2018) among the 4 simulated years were clearly reproduced in the forecast simulation results. Furthermore, sub-seasonal and seasonal simulations reflected the early harvest in the drought year of 2018 in the DE-NRW domain. However, simulated inter-annual yield variability was lower in all simulations compared to the official statistics. While general soil moisture trends, such as the European drought in 2018, were captured by the seasonal experiments, we found systematic overestimations and underestimations in both the forecast and reference simulations compared to the Soil Moisture Active Passive Level-3 soil moisture product (SMAP L3) and the Soil Moisture Climate Change Initiative Combined dataset from the European Space Agency (ESA CCI). These observed biases of soil moisture and the low inter-annual differences in simulated crop yield indicate the need to improve the representation of these variables in CLM5 to increase the model sensitivity to drought stress and other crop stressors. [ABSTRACT FROM AUTHOR]

The retrieval of soil moisture information with spatially and temporally high resolution from Synthetic Aperture Radar (SAR) observations is still a challenge. By using multi-orbit Sentinel-1 C-band time series, we present a novel approach for estimating volumetric soil moisture content for agricultural areas with a temporal resolution of one to two days, based on a short-term change detection method. By applying an incidence angle normalization and a Fourier Series transformation, the effect of varying incidence angles on the backscattering signal could be reduced. As the C-band co-polarized backscattering signal is prone to vegetational changes, it is used in this study for the vegetational correction of its related backscatter ratios. The retrieving algorithm was implemented in a cloud-processing environment, enabling a potential global and scalable application. Validated against eight in-situ cosmic ray neutron probe stations across the Rur catchment (Germany) as well as six capacitance stations at the Apulian Tavoliere (Italy) site for the years 2018 to 2020, the method achieves a correlation coefficient of R of 0.63 with an unbiased Root Mean Square Error of 0.063 m3/m3. [ABSTRACT FROM AUTHOR]

An accurate assessment of the physical states of the Earth system is an essential component of many scientific, societal and economical considerations. These assessments are becoming an increasingly challenging computational task since we aim to resolve models with high resolutions in space and time, to consider complex coupled partial differential equations, and to estimate uncertainties, which often requires many realizations. Machine learning methods are becoming a very popular method for the construction of surrogate models to address these computational issues. However, they also face major challenges in producing explainable, scalable, interpretable and robust models. In this manuscript, we evaluate the perspectives of geoscience applications of physics-based machine learning, which combines physics-based and data-driven methods to overcome the limitations of each approach taken alone. Through three designated examples (from the fields of geothermal energy, geodynamics, and hydrology), we show that the non-intrusive reduced basis method as a physics-based machine learning approach is able to produce highly precise surrogate models that are explainable, scalable, interpretable, and robust. [ABSTRACT FROM AUTHOR]

Land surface models (LSM) are an important tool for advancing our knowledge of the Earth system. LSM are constantly improved to represent the various terrestrial processes in more detail. High quality data, freely available from various observation networks, are providing being used to improve the prediction of terrestrial states and fluxes of water and energy. To optimize LSM with observations, data assimilation methods and tools have been developed in the past decades. We apply the coupled Community Land Model version 5 (CLM5) and Parallel Data Assimilation Framework (PDAF) system (CLM5-PDAF) for thirteen forest field sites throughout Europe covering different climate zones. The goal of this study is to assimilate in-situ soil moisture measurements into CLM5 to improve the modeled evapotranspiration fluxes. The modeled fluxes will be evaluated using the predicted evapotranspiration fluxes with eddy covariance (EC) systems. Most of the sites use point scale measurements from, however for three of the forest sites we use soil water content data from cosmic-ray neutron sensors, which have a measurement scale closer to the typical land surface model grid scale and EC footprint. Our results show that while data assimilation reduced the root-mean-square error for soil water content on average by 56 to 64 %, the root-mean-square error for the evapotranspiration estimation is increased by 4 %. This finding indicates that state-of-the-art LSM such as CLM5 still suffer from uncertainties in the representation of soil hydrological processes in forests, e.g. deep root water uptake, or highly uncertain vegetation parameters. [ABSTRACT FROM AUTHOR]

Long-range weather forecasts provide predictions of atmospheric, ocean and land surface conditions that can potentially be used in land surface and hydrological models to predict the water and energy status of the land surface or in crop growth models to predict yield for water resources or agricultural planning. However, the coarse spatial and temporal resolutions of available forecast products have hindered their widespread use in such modelling applications that usually require high resolution input data. In this study, we applied sub-seasonal (up to 4 months) and seasonal (7 months) weather forecasts from the latest European Centre for Medium-Range Weather Forecasts (ECMWF) seasonal forecasting system (SEAS5) in a land surface modelling approach using the Community Land Model version 5.0 (CLM5). Simulations were conducted for 2017-2020 forced with subseasonal and seasonal weather forecasts over two different domains with contrasting climate and cropping conditions: the German state of North Rhine-Westphalia and the Australian state of Victoria. We found that, after pre-processing of the forecast products (temporal downscaling of precipitation and incoming shortwave radiation), the simulations forced with seasonal and sub-seasonal forecasts were able to generate a model system response very close to reference simulation results forced by reanalysis data. Differences between seasonal and sub-seasonal experiments were insignificant. The forecast experiments were able to satisfactorily capture recorded inter-annual variations of crop yield. In addition, they also reproduced the generally higher inter-annual variability in crop yield across the Australian domain (approximately 50 % inter-annual variability in recorded yields and up to 17 % in simulated yields) compared to the German domain (approximately 15 % inter-annual variability in recorded yields and up to 5 % in simulated yields). The high and low yield seasons (2020 and 2018) among the four simulated years were clearly reproduced in forecast simulation results. Furthermore, sub-seasonal and seasonal simulations reflected the early harvest in the drought year of 2018 in the German domain. However, the simulated inter-annual yield variability was lower in all simulations compared to the official statistics. While general soil moisture trends, such as the European drought in 2018, were captured by the seasonal experiments, we found systematic over- and underestimations in both the forecast and the reference simulations compared to the Soil Moisture Active Passive Level-3 soil moisture product (SMAP L3) and the Soil Moisture Climate Change Initiative Combined dataset from the European Space Agency's (ESA CCI). These observed biases of soil moisture as well as the low inter-annual variability of simulated crop yield indicate the need to improve the representation of these variables in CLM5 to increase the model sensitivity to drought stress and other crop stressors. [ABSTRACT FROM AUTHOR]

Aims: We aimed at assessing the influence of above- and below-ground environmental conditions over the performance of Centaurea jacea L., a drought-resistant grassland forb species. Methods: Transpiration rate, CO2 assimilation rate, leaf water potential, instantaneous and intrinsic water use efficiency, temperature, relative humidity, vapor pressure deficit and soil water content in one plant and root length density in four plants, all grown in custom-made columns, were monitored daily for 87 days in the lab. The soil water isotopic composition in eleven depths was recorded daily in a non-destructive manner. The isotopic composition of plant transpiration was inferred from gas chamber measurements. Vertical isotopic gradients in the soil column were created by adding labeled water. Daily root water uptake (RWU) profiles were computed using the multi-source mixing model Stable Isotope Analysis in R (Parnell et al. PLoS ONE 5(3):1–5, 2010). Results: RWU occurred mainly in soil layer 0–15 cm, ranging from 79 to 44%, even when water was more easily available in deeper layers. In wet soil, the transpiration rate was driven mainly by vapor pressure deficit and light intensity. Once soil water content was less than 0.12 cm3 cm− 3, the computed canopy conductance declined, which restricted leaf gas exchange. Leaf water potential dropped steeply to around − 3 MPa after soil water content was below 0.10 cm3 cm− 3. Conclusion: Our comprehensive data set contributes to a better understanding of the effects of drought on a grassland species and the limits of its acclimation in dry conditions. [ABSTRACT FROM AUTHOR]

Soil water status, which refers to the wetness or dryness of soils, is crucial for the productivity of agroecosystems, as it determines nutrient cycling and uptake physically via transport, biologically via the moisture‐dependent activity of soil flora, fauna, and plants, and chemically via specific hydrolyses and redox reactions. Here, we focus on the dynamics of nitrogen (N), phosphorus (P), and sulfur (S) and review how soil water is coupled to the cycling of these elements and related stoichiometric controls across different scales within agroecosystems. These scales span processes at the molecular level, where nutrients and water are consumed, to processes in the soil pore system, within a soil profile and across the landscape. We highlight that with increasing mobility of the nutrients in water, water‐based nutrient flux may alleviate or even exacerbate imbalances in nutrient supply within soils, for example, by transport of mobile nutrients towards previously depleted microsites (alleviating imbalances), or by selective loss of mobile nutrients from microsites (increasing imbalances). These imbalances can be modulated by biological activity, especially by fungal hyphae and roots, which contribute to nutrient redistribution within soils, and which are themselves dependent on specific, optimal water availability. At larger scales, such small‐scale effects converge with nutrient inputs from atmospheric (wet deposition) or nonlocal sources and with nutrient losses from the soil system towards aquifers. Hence, water acts as a major control in nutrient cycling across scales in agroecosystems and may either exacerbate or remove spatial disparities in the availability of the individual nutrients (N, P, S) required for biological activity. [ABSTRACT FROM AUTHOR]

Due to surface water scarcity and variability in rainwater events in the Thies region, Senegal, farmers mostly rely on groundwater resources to improve yield production, although water quality in the hydrosystem and its suitability may be an issue. The objective of this study is to evaluate suitability of the Maastrichtian and Paleocene aquifers located the central western part of Senegal in the Thies region for irrigation purposes. For this purpose, chemical analyses were carried out on the major ions on the Maastrichtian and Paleocene aquifers through a network of 62 sample locations (39 from the Maastrichtian and 23 from the Paleocene) sampled in June 2019. Different water suitability assessment indices were used to evaluate the potential for irrigation, including sodium absorption ration (SAR), percentage of sodium (Na%), permeability index, Kelly ratio (RK), and salinity potential (SP) together with Wilcox and USSL diagrams as methods. The results indicate, that the Ca-Mg-HCO3 water type is dominant in the Maastrichtian aquifer, while in the Paleocene aquifer, HCO3-Ca-Mg and Cl-Ca-Mg are the main water types. The combination of these computed index values indicates that the percentage of suitable water for irrigation purposes accounts for 89.7% and 86.9% of the samples for the Maastrichtian and the Paleocene aquifers, respectively. On the other hand, water samples unsuitable for irrigation purposes range between only 10.3% and 13.1% for the two aquifers. [ABSTRACT FROM AUTHOR]

Current databases of soil hydraulic properties (SHPs) have typically been used to develop pedotransfer functions (PTFs) to estimate water retention [θ(h)] assuming a unimodal pore‐size distribution. However, natural soils often show the presence of bimodal to multimodal pore‐size distributions. Here, we used three widely spread databases for PTF development: UNsaturated SOil hydraulic DAtabase (UNSODA) 2.0, Vereecken, and European hydropedological data inventory (EU‐HYDI), to analyze the presence of structural effects in both θ(h) and hydraulic conductivity [K(h)]. Only undisturbed samples were included in the analysis that contained enough datapoints for both θ(h) and K(h) properties, especially in the wet range. One‐hundred ninety‐two samples were suitable for our analysis, which is only 1% of the total samples in the three databases. Results showed that 65% of the samples exhibited a bimodal pore‐size distribution, and bimodality was not limited to fine‐textured but also coarser‐textured soils. The Mualem–van Genuchten (MvG) expression for both unimodal and bimodal soils was not able to predict the observed unsaturated K. Only a joint fitting of measured θ(h) and K(h) functions provided parameter estimates that were able to describe unsaturated K for uni‐ and bimodal soils. In addition, we observed a negative relationship between α and n in the case of low sand content (<52%) for both unimodal and bimodal matrix domain properties, contradicting the classical notion. The ratio of α for the macropore and matrix domain was positively correlated with the fraction of macropores and sand content. We anticipate that the results will contribute to deriving PTF for structured soils and avoid unrealistic combinations of MvG parameters. Core Ideas: Only 1% of the samples in the three databases was suitable for analyzing the presence of bimodality.Bimodality was seen in 65% of the samples, both in fine‐textured and coarser‐textured soils.Joint fitting of θ(h) and K(h) data outperforms the fitting of only θ(h) and using MvG function to estimate K(h).Contradicting the classical notion, a negative relationship between the shape parameters α and n was observed for both uni‐ and bimodal properties.α for macropore and matrix domains was positively correlated with the fraction of macropores and sand content. [ABSTRACT FROM AUTHOR]

Projections of global climate models suggest that ongoing human-induced climate change will lead to an increase in the frequency of severe droughts in many important agricultural regions of the world. Eco-hydrological models that integrate current understanding of the interacting processes governing soil water balance and plant growth may be useful tools to predict the impacts of climate change on crop production. However, the validation status of these models for making predictions under climate change is still unclear, since few suitable datasets are available for model testing. One promising approach is to test models using data obtained in "space-for-time" substitution experiments, in which samples are transferred among locations with contrasting current climates in order to mimic future climatic conditions. An important advantage of this approach is that the soil type is the same, so that differences in soil properties are not confounded with the influence of climate on water balance and crop growth. In this study, we evaluate the capability of a relatively simple eco-hydrological model to reproduce 6 years (2013–2018) of measurements of soil water contents, water balance components and grass production made in weighing lysimeters located at two sites within the TERENO-SoilCan network in Germany. Three lysimeters are located at an upland site at Rollesbroich with a cool, wet climate, while three others had been moved from Rollesbroich to a warmer and drier climate on the lower Rhine valley floodplain at Selhausen. Four of the most sensitive parameters in the model were treated as uncertain within the framework of the GLUE (generalized likelihood uncertainty estimation) methodology, while the remaining parameters in the model were set according to site measurements or data in the literature. The model satisfactorily reproduced the measurements at both sites, and some significant differences in the posterior ranges of the four uncertain parameters were found. In particular, the results indicated greater stomatal conductance as well an increase in dry-matter allocation below ground and a significantly larger maximum root depth for the three lysimeters that had been moved to Selhausen. As a consequence, the apparent water use efficiency (above-ground harvest divided by evapotranspiration) was significantly smaller at Selhausen than Rollesbroich. Data on species abundance on the lysimeters provide one possible explanation for the differences in the plant traits at the two sites derived from model calibration. These observations showed that the plant community at Selhausen had changed significantly in response to the drier climate, with a significant decrease in the abundance of herbs and an increase in the proportion of grass species. The differences in root depth and leaf conductance may also be a consequence of plasticity or acclimation at the species level. Regardless of the reason, we may conclude that such adaptations introduce significant additional uncertainties into model predictions of water balance and plant growth in response to climate change. [ABSTRACT FROM AUTHOR]

Crop model intercomparison studies have mostly focused on the assessment of predictive capabilities for crop development using weather and basic soil data from the same location. Still challenging is the model performance when considering complex interrelations between soil and crop dynamics under a changing climate. The objective of this study was to test the agronomic crop and environmental flux‐related performance of a set of crop models. The aim was to predict weighing lysimeter‐based crop (i.e., agronomic) and water‐related flux or state data (i.e., environmental) obtained for the same soil monoliths that were taken from their original environment and translocated to regions with different climatic conditions, after model calibration at the original site. Eleven models were deployed in the study. The lysimeter data (2014–2018) were from the Dedelow (Dd), Bad Lauchstädt (BL), and Selhausen (Se) sites of the TERENO (TERrestrial ENvironmental Observatories) SOILCan network. Soil monoliths from Dd were transferred to the drier and warmer BL site and the wetter and warmer Se site, which allowed a comparison of similar soil and crop under varying climatic conditions. The model parameters were calibrated using an identical set of crop‐ and soil‐related data from Dd. Environmental fluxes and crop growth of Dd soil were predicted for conditions at BL and Se sites using the calibrated models. The comparison of predicted and measured data of Dd lysimeters at BL and Se revealed differences among models. At site BL, the crop models predicted agronomic and environmental components similarly well. Model performance values indicate that the environmental components at site Se were better predicted than agronomic ones. The multi‐model mean was for most observations the better predictor compared with those of individual models. For Se site conditions, crop models failed to predict site‐specific crop development indicating that climatic conditions (i.e., heat stress) were outside the range of variation in the data sets considered for model calibration. For improving predictive ability of crop models (i.e., productivity and fluxes), more attention should be paid to soil‐related data (i.e., water fluxes and system states) when simulating soil–crop–climate interrelations in changing climatic conditions. Core Ideas: We demonstrate the use of high precision weighable lysimeter for full model calibration and validation.Lysimeter data from translocated soils represent effects of changing climatic conditions.We compare calibration with blind forward simulations (fixed soil and calibrated crop parameter).We compare individual crop model predictions with multi‐model mean.We test the predictive ability of crop models and multi‐model mean. [ABSTRACT FROM AUTHOR]

Integrated terrestrial system models predict the coupled water, energy and biogeochemical cycles. Simulations with these models are affected by uncertainties of model parameters, initial and boundary conditions, atmospheric forcings and the biophysical processes. Data assimilation (DA) can quantify and reduce the uncertainty. This has been tested intensively for single compartment models, but far less for integrated models with multiple compartments. We constructed a virtual reality (VR) with a coupled land surface‐subsurface model under the Terrestrial Systems Modeling Platform, which mimics the Neckar catchment in southern Germany. Soil moisture and groundwater level (GWL) data extracted from the simulated VR are used as measurements to be assimilated with state‐only/state‐hydraulic parameter estimation. Soil moisture DA improves soil moisture characterization in the vertical profile and the neighboring grid cells, with a 40 ∼ 60% reduction of root mean square error (RMSE) over the observation points. In spite of a small ensemble size of 64 members, assimilating soil moisture data improved saturated hydraulic conductivity estimation around the measurement locations. The characterization of evapotranspiration and river discharge only show limited improvements (1% at observation points and less than 0.1% in RMSE at 3 selected gauge locations respectively). GWL DA not only improves the GWL characterization (76 ∼ 88% RMSE reduction at observation locations) but also soil moisture for some cases. In addition, a clear improvement in GWL characterization is observed up to 8 km from the observations, and updating the model states of the saturated zone only instead of the complete domain gives better performance. Key Points: Groundwater level (GWL) and soil moisture data assimilation (DA) improve characterization of subsurface statesGWL DA improves GWL characterization up to 8 km from observation locationsEvapotranspiration and discharge hardly better estimated, despite better estimates of soil moisture and GWL [ABSTRACT FROM AUTHOR]

Climate change amplifies the intensity and occurrence of dry periods leading to drought stress in vegetation. For monitoring vegetation stresses, sun-induced chlorophyll fluorescence (SIF) observations are a potential game-changer, as the SIF emission is mechanistically coupled to photosynthetic activity. Yet, the benefit of SIF for drought stress monitoring is not yet understood. This paper analyses the impact of drought stress on canopy-scale SIF emission and surface reflectance over a lettuce and mustard stand with continuous field spectrometer measurements. Here, the SIF measurements are linked to the plant's photosynthetic efficiency, whereas the surface reflectance can be used to monitor the canopy structure. The mustard canopy showed a reduction in the biochemical component of its SIF emission (the fluorescence emission efficiency at 760 nm—ϵ760) as a reaction to drought stress, whereas its structural component (the Fluorescence Correction Vegetation Index—FCVI) barely showed a reaction. The lettuce canopy showed both an increase in the variability of its surface reflectance at a sub-daily scale and a decrease in ϵ760 during a drought stress event. These reactions occurred simultaneously, suggesting that sun-induced chlorophyll fluorescence and reflectance-based indices sensitive to the canopy structure provide complementary information. The intensity of these reactions depend on both the soil water availability and the atmospheric water demand. This paper highlights the potential for SIF from the upcoming FLuorescence EXplorer (FLEX) satellite to provide a unique insight on the plant's water status. At the same time, data on the canopy reflectance with a sub-daily temporal resolution are a promising additional stress indicator for certain species. [ABSTRACT FROM AUTHOR]

The integrated distributed hydrological model ParFlow-CLM was used to predict water and energy transport between subsurface, land surface, and atmosphere for the Stettbach headwater catchment, Germany. Based on this model, global sensitivity analysis was performed using the Latin-Hypercube (LH) sampling strategy followed by the One-factor-Ata-Time (OAT) method to identify the most influential and interactive parameters affecting the main hydrologic processes. total 12 parameters were evaluated including soil hydraulic properties, storage, Manning coefficient, leaf area index, stem area index, and aerodynamic resistance that characterize water and energy fluxes in soil and vegetation. In addition, the sensitivity analysis was also carried out for different slopes and meteorological conditions to test the transferability of the results regions with other topographies and climates. Our results show that the simulated energy fluxes, i.e. latent heat flux and sensible heat flux are sensitive to the parameters such as wilting point, leaf area index, and stem area index, especially for steep slope and subarctic climate conditions. The simulated soil evaporation, plant transpiration, infiltration, and runoff, are most sensitive to soil porosity, the van Genuchten parameter n representing the soil pore size distribution, soil wilting point, and leaf area index. The subsurface soil water storage and groundwater storage are most sensitive to soil porosity, while the surface water storage was most sensitive to the soil roughness parameter. For the different slope and climate conditions, the rank order of input parameter sensitivity is consistent, but the magnitude of parameter sensitivity is very different. The strongest deviation in parameter sensitivity occurred for sensible heat flux under the different slope conditions as well as for transpiration under different climate conditions. Overall, this study provides an insight into the most important input parameters that control hydrological fluxes and how the simulated variables vary with the change in parameter values, which can improve our understanding of the key processes in the model and help us to reduce the computational demands of completing multiple simulations of expensive domains. [ABSTRACT FROM AUTHOR]

In his seminal paper on the solution of the infiltration equation, Philip (1969), https://doi.org/10.1016/b978-1-4831-9936-8.50010-6 proposed a gravity time, tgrav, to estimate practical convergence time and the time domain validity of his infinite time series expansion, TSE, for describing the transient state. The parameter tgrav refers to a point in time where infiltration is dominated equally by capillarity and gravity as derived from the first two (dominant) terms of the TSE. Evidence suggests that applicability of the truncated two‐term equation of Philip has a time limit requiring higher‐order TSE terms to better describe the infiltration process for times exceeding that limit. Since the conceptual definition of tgrav is valid regardless of the infiltration model used, we opted to reformulate tgrav using the analytic implicit model proposed by Parlange et al. (1982), https://doi.org/10.1097/00010694-198206000-00001 valid for all times and related TSE. Our derived gravity times ensure a given accuracy of the approximations describing transient states, while also providing insight about the times needed to reach steady state. In addition to the roles of soil sorptivity (S) and the saturated (Ks) and initial (Ki) hydraulic conductivities, we explored the effects of a soil specific shape parameter β, involved in Parlange's model and related to the type of soil, on the behavior of tgrav. We show that the reformulated tgrav (notably tgrav=F(β)S2/Ks−Ki2, ${t}_{\text{grav}}=\,F(\beta){S}^{2}/{\left({K}_{s}-{K}_{i}\right)}^{2},$ where F(β) is a β‐dependent function) is about three times larger than the classical tgrav given by tgrav,Philip=S2/Ks−Ki2 $\,{t}_{\text{grav},\text{Philip}}={S}^{2}/{\left({K}_{s}-{K}_{i}\right)}^{2}$. The differences between the classical tgrav,Philip and the reformulated tgrav increase for fine‐textured soils, attributed to the time needed to attain steady‐state infiltration and thus i + nfiltration for inferring soil hydraulic properties. Results show that the proposed tgrav is a better indicator of time domain validity than tgrav,Philip. For the attainment of steady‐state infiltration, the reformulated tgrav is suitable for coarse‐textured soils. Still neither the reformulated tgrav nor the classical tgrav,Philip are suitable 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 compared to tgrav,Philip. Key Points: A new formulation for infiltration characteristic time, tgrav, is providedThe reformulated tgrav seems to be a better criterion for convergence time of Philip's truncated infiltration equationsThe usage of reformulated tgrav improves predictions of soil hydraulic parameters [ABSTRACT FROM AUTHOR]

Soil bulk density (ρb) is an important indicator of soil quality, productivity, compaction and porosity. Despite its importance, ρb is often omitted from global datasets due to the costs of making many direct ρb measurements and the difficulty of direct measurement on rocky, sandy, very dry, or very wet soils. Pedotransfer functions (PTFs) are deployed to address these limitations. Using readily available soil properties, PTFs employ estimator equations to fit existing datasets to estimate properties like ρb. However, PTF performance often declines when applied to soils outside those in the training dataset. Potentially, recalibrating existing PTFs using new observations would leverage the power of large datasets used in the original PTF derivation, while updating information based on new soil observations. Here, we evaluate such a recalibration approach for ρb estimation, benchmarking its performance against two alternatives: the original, uncalibrated PTFs, and novel, local PTFs derived solely from new soil observations. Using a ρb dataset of N = 360 total observations obtained in West Azerbaijan, Iran, we varied the local dataset size (with N = 15, 30, 60, and 360) and recalibrated four existing PTFs with these data. Local PTFs were generated based on stepwise multiple linear regression for the same datasets. The same PTFs (original, recalibrated, and local) were also applied to the study area, and the resulting ρb estimates were compared with the global SoilGrids dataset. Recalibration of PTFs reduced errors relative to the original uncalibrated PTFs; for instance, the NSE increased from −22.07 to 0.30 (uncalibrated) to 0.20–0.41 (recalibrated), and RMSE decreased from 0.12 to 0.60 Mg m−3 (uncalibrated) to 0.10–0.13 Mg m−3 (recalibrated). The recalibrated PTFs performance was comparable to or better than local PTFs applied to the same data. Recalibration of existing PTFs with local/regional uses provides a viable alternative to the use of global datasets or the development of local PTFs in data‐scarce regions. Highlights: Existing global PTFs were calibrated and tested using a small dataset for local utilisation.Several new local PTFs were also developed using the same datasets.Recalibration of existing global PTFs is comparable to or more accurate than developing new PTFs. [ABSTRACT FROM AUTHOR]

Soil hydraulic conductivity (k soil ) drops significantly in dry soils, resulting in steep soil water potential gradients (ψ s ) near plant roots during water uptake. Coarse soil grid resolutions in root system scale (RSS) models of root water uptake (RWU) generally do not spatially resolve this gradient in drying soils which can lead to a large overestimation of RWU. To quantify this, we consider a benchmark scenario of RWU from drying soil for which a numerical reference solution is available. We analyze this problem using a finite volume scheme and investigate the impact of grid size on the RSS model results. At dry conditions, the cumulative RWU was overestimated by up to 300% for the coarsest soil grid of 4.0 cm and by 30% for the finest soil grid of 0.2 cm, while the computational demand increased from 19 s to 21 h. As an accurate and computationally efficient alternative to the RSS model, we implemented a continuum multi-scale model where we keep a coarse grid resolution for the bulk soil, but in addition, we solve a 1-dimensional radially symmetric soil model at rhizosphere scale around individual root segments. The models at the two scales are coupled in a mass-conservative way. The multi-scale model compares best to the reference solution (−20%) at much lower computational costs of 4 min. Our results demonstrate the need to shift to improved RWU models when simulating dry soil conditions and highlight that results for dry conditions obtained with RSS models of RWU should be interpreted with caution. [ABSTRACT FROM AUTHOR]

Cloud computing facilities can provide crucial computing support for processing the time series of satellite data and exploiting their spatio-temporal information content. However, dedicated efforts are still required to develop workflows, executable on cloud-based platforms, for ingesting the satellite data, performing the targeted processes, and generating the desired products. In this study, an operational workflow is proposed, based on monthly Evaporative Stress Index (ESI) anomaly, and implemented in cloud-based online Virtual Earth Laboratory (VLab) platform, as a demonstration, to monitor European agricultural water stress. To this end, daily time-series of actual and reference evapotranspiration (ETa and ET0), from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) sensor, were used to execute the proposed workflow successfully on VLab. The execution of the workflow resulted in obtaining one decade (2011–2020) of European monthly agricultural water stress maps at 0.04˚ spatial resolution and corresponding stress reports for each country. To support open science, all the workflow outputs are stored in GeoServer, documented in GeoNetwork, and made available through MapStore. This enables creating a dashboard for better visualization of the results for end-users. The results from this study demonstrate the capability of VLab platform for water stress detection from time series of SEVIRI-ET data. [ABSTRACT FROM AUTHOR]

Fast and accurate large‐scale localization and quantification of harmfully compacted soils in recultivated post‐mining landscapes are of particular importance for mining companies and the following farmers. The use of heavy machinery during recultivation imposes soil stress and can cause irreversible subsoil compaction limiting crop growth in the long term. To overcome or guide classical point‐scale methods to determine compaction, fast methods covering large areas are required. In our study, a recultivated field of the Garzweiler mine in North Rhine‐Westphalia, Germany, with known variability in crop performance was intensively studied using non‐invasive electromagnetic induction (EMI) and electrode‐based electrical resistivity tomography (ERT). Additionally, soil bulk density, volumetric soil water content and soil textures were analysed along two transects covering different compaction levels. The results showed that the measured EMI apparent electrical conductivity (ECa) along the transects was highly correlated (R2 >.7 for different dates and depths below 0.3 m) to subsoil bulk density. Finally, the correlations established along the transects were used to predict harmful subsoil compaction within the field, whereby a spatial probabilistic map of zones of harmful compaction was developed. In general, the results revealed the feasibility of using the EMI derived ECa to predict harmful compaction. They can be the basis for quick monitoring of the recultivation process and implementation of necessary melioration to return a well‐structured soil with good water and nutrient accessibility, and rooting depths for increased crop yields to the farmers. [ABSTRACT FROM AUTHOR]

Many different equations have been proposed to describe quantitatively one‐dimensional soil water infiltration. The unknown coefficients of these equations characterize soil hydraulic properties and may be estimated from a n record, {t∼i,I∼i}i=1n$\{ {\tilde t_i},{\tilde I_i}\} _{i = 1}^n$, of cumulative infiltration measurements using curve fitting techniques. The two‐term infiltration equation, I(t)=St+cKst$I(t) = S\sqrt t + c{K_{\rm{s}}}t$, of Philip has been widely used to describe measured infiltration data. This function enjoys a solid mathematical–physical underpinning and admits a closed‐form solution for the soil sorptivity, S [L T−1/2], and multiple, c [−], of the saturated hydraulic conductivity, Ks [L T−1]. However, Philip's two‐term equation has a limited time validity, tvalid [T], and thus cumulative infiltration data, I∼(t∼)$\tilde I(\tilde t)$, beyond t=tvalid$t = {t_{{\rm{valid}}}}$ will corrupt the estimates of S and Ks. This paper introduces a novel method for estimating S, c, Ks, and tvalid of Philip's two‐term infiltration equation. This method, coined parasite inversion, use as vehicle Parlange's three‐parameter infiltration equation. As prerequisite to our method, we present as secondary contribution an exact, robust and efficient numerical solution of Parlange's infiltration equation. This solution admits Bayesian parameter estimation with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm and yields as byproduct the marginal distribution of Parlange's β parameter. We evaluate our method for 12 USDA soil types using synthetic infiltration data simulated with HYDRUS‐1D. An excellent match is observed between the inferred values of S and Ks and their "true" values known beforehand. Furthermore, our estimates of c and tvalid correlate well with soil texture, corroborate linearity of the c(β)$c({{\beta}})$ relationship for 0≤t≤tvalid$0 \le t \le {t_{{\rm{valid}}}}$, and fall within reported ranges. A cumulative vertical infiltration of about 2.5 cm may serve as guideline for the time‐validity of Philip's two‐term infiltration equation. Core Ideas: We describe a method that can determine the time validity of infiltration equations.This method infers as well the hydraulic properties of unsaturated soils.We present a robust and efficient numerical solution of Parlange's three‐parameter infiltration equation.Posterior distribution determines how well parameters are defined by calibration to data.Two and a half centimeters of cumulative infiltration is a good proxy for time validity Philip's two‐term equation. [ABSTRACT FROM AUTHOR]

Effects of climate change on the ecosystem productivity and water fluxes have been studied in various types of experiments. However, it is still largely unknown whether and how the experimental approach itself affects the results of such studies. We employed two contrasting experimental approaches, using high-precision weighable monolithic lysimeters, over a period of 4 years to identify and compare the responses of water fluxes and aboveground biomass to climate change in permanent grassland. The first, manipulative, approach is based on controlled increases of atmospheric CO 2 concentration and surface temperature. The second, observational, approach uses data from a space-for-time substitution along a gradient of climatic conditions. The Budyko framework was used to identify if the soil ecosystem is energy limited or water limited. Elevated temperature reduced the amount of non-rainfall water, particularly during the growing season in both approaches. In energy-limited grassland ecosystems, elevated temperature increased the actual evapotranspiration and decreased aboveground biomass. As a consequence, elevated temperature led to decreasing seepage rates in energy-limited systems. Under water-limited conditions in dry periods, elevated temperature aggravated water stress and, thus, resulted in reduced actual evapotranspiration. The already small seepage rates of the drier soils remained almost unaffected under these conditions compared to soils under wetter conditions. Elevated atmospheric CO 2 reduced both actual evapotranspiration and aboveground biomass in the manipulative experiment and, therefore, led to a clear increase and change in seasonality of seepage. As expected, the aboveground biomass productivity and ecosystem efficiency indicators of the water-limited ecosystems were negatively correlated with an increase in aridity, while the trend was unclear for the energy-limited ecosystems. In both experimental approaches, the responses of soil water fluxes and biomass production mainly depend on the ecosystems' status with respect to energy or water limitation. To thoroughly understand the ecosystem response to climate change and be able to identify tipping points, experiments need to embrace sufficiently extreme boundary conditions and explore responses to individual and multiple drivers, such as temperature, CO 2 concentration, and precipitation, including non-rainfall water. In this regard, manipulative and observational climate change experiments complement one another and, thus, should be combined in the investigation of climate change effects on grassland. [ABSTRACT FROM AUTHOR]

Droughts can have important impacts on environment and economy like in the year 2018 in parts of Europe. Droughts can be analyzed in terms of meteorological drought, agricultural drought, hydrological drought and social-economic drought. In this paper, we focus on meteorological and agricultural drought and analyzed drought trends for the period 1965-2019 and assessed how extreme the drought year 2018 was in Germany and the Netherlands. The analysis was made on the basis of the following drought indices: standardized precipitation index (SPI), standardized soil moisture index (SSI), potential precipitation deficit (PPD) and ET deficit. SPI and SSI were computed at two time scales, the period April-September and a 12-months period. In order to analyze drought trends and the ranking of the year 2018, HYDRUS 1-D simulations were carried out for 31 sites with long-term meteorological observations and soil moisture, potential evapotranspiration (ET) and actual ET were determined for five soil types (clay, silt, loam, sandy loam and loamy sand). The results show that the year 2018 was severely dry, which was especially related to the highest potential ET in the time series 1965-2019, for most of the sites. For around half of the 31 sites the year 2018 had the lowest SSI, and largest PPD and ET-deficit in the 1965-2019 time series, followed by 1976 and 2003. The trend analysis reveals that meteorological drought (SPI) hardly shows significant trends over 1965-2019 over the studied domain, but agricultural droughts (SSI) are increasing, at several sites significantly, and at even more sites PPD and ET deficit show significant trends. The increasing droughts over Germany and Netherlands are mainly driven by increasing potential ET and increasing vegetation water demand. [ABSTRACT FROM AUTHOR]

Background and aims: Characterizing root system architectures of field-grown crops is challenging as root systems are hidden in the soil. We investigate the possibility of estimating root architecture model parameters from soil core data in a Bayesian framework. Methods: In a synthetic experiment, we simulated wheat root systems in a virtual field plot with the stochastic CRootBox model. We virtually sampled soil cores from this plot to create synthetic measurement data. We used the Markov chain Monte Carlo (MCMC) DREAM(ZS) sampler to estimate the most sensitive root system architecture parameters. To deal with the CRootBox model stochasticity and limited computational resources, we essentially added a stochastic component to the likelihood function, thereby turning the MCMC sampling into a form of approximate Bayesian computation (ABC). Results: A few zero-order root parameters: maximum length, elongation rate, insertion angles, and numbers of zero-order roots, with narrow posterior distributions centered around true parameter values were identifiable from soil core data. Yet other zero-order and higher-order root parameters were not identifiable showing a sizeable posterior uncertainty. Conclusions: Bayesian inference of root architecture parameters from root density profiles is an effective method to extract information about sensitive parameters hidden in these profiles. Equally important, this method also identifies which information about root architecture is lost when root architecture is aggregated in root density profiles. [ABSTRACT FROM AUTHOR]

Incorporation of crop residues into the soil has been widely recommended as an effective method to sustain soil fertility and improve soil carbon sequestration in arable lands. However, it may lead to an increase in the emission of nitrous oxide (N2O) and leaching of nitrate (NO3−) to groundwater due to higher nitrogen (N) availability after crop residue incorporation. Here, we conducted a meta-analysis based on 345 observations from 90 peer-reviewed studies to evaluate the effects of crop residue return on soil N2O emissions and NO3− leaching for different locations, climatic and soil conditions, and agricultural management strategies. On average, crop residue incorporation significantly stimulated N2O emissions by 29.7%, but decreased NO3− leaching by 14.4%. The increase in N2O emissions was negatively and significantly correlated with mean annual temperature and mean annual precipitation, and with the most significant changes occurring in the temperate climate zone. Crop residues stimulated N2O emission mainly in soils with pH ranging between 5.5 and 6.5, or above 7.5 in soils with low clay content. In addition, crop residue application decreased NO3− leaching significantly in soils with sandy loam, silty clay loam, and silt loam textures. Our analysis reveals that an appropriate crop residue management adapted to the site-specific soil and environmental conditions is critical for increasing soil organic carbon stocks and decreasing nitrogen losses. The most important novel finding is that residue return, despite stimulation of N2O emissions, is particularly effective in reducing NO3− leaching in soils with loamy texture, which are generally among the most productive arable soils. [ABSTRACT FROM AUTHOR]

The original article can be found online at https://doi.org/10.1007/s13280-023-01930-4. B Publisher Correction to: Ambio b https://doi.org/10.1007/s13280-023-01930-4 In the original publication of article, a typo was introduced in the co-author name during the creation of the article proof. [Extracted from the article]

The Committee on Earth Observation Satellites (CEOS) Land Product Validation (LPV) subgroup has been established to coordinate the development of standardized validation across the satellite-derived products from different platforms, sensors, and algorithms with reference measurements from the in situ networks. Soil moisture exhibits a high variability in space that challenges the in situ validation. One of the main drivers for this variability is the characteristic heterogeneity in the soil texture. By the machine learning methods using the soil profile measurements and the remotely sensed predictors, spatially continuous maps of basic soil properties such as soil texture and bulk density are available. Those can be used to estimate soil moisture variability within a satellite product grid cell, here exemplarily shown for the Soil Moisture Active Passive (SMAP) 36-km product. The soil moisture standard deviation is described as a function of the mean soil moisture, whereby the approach needs the mean and standard deviation of the hydraulic parameters as input. The resulting global data set helps identifying the number of in situ stations necessary to validate the coarse soil moisture products. For most SMAP grid cells, three to four stations are adequate to estimate the mean soil moisture for validation; however, also regions were identified where 80 stations are necessary. [ABSTRACT FROM AUTHOR]

Root water uptake is an important process in the terrestrial water cycle. How this process depends on soil water content, root distributions, and root properties is a soil–root hydraulic problem. We compare different approaches to implement root hydraulics in macroscopic soil water flow and land surface models. By upscaling a three-dimensional hydraulic root architecture model, we derived an exact macroscopic root hydraulic model. The macroscopic model uses the following three characteristics: the root system conductance, Krs , the standard uptake fraction, SUF , which represents the uptake from a soil profile with a uniform hydraulic head, and a compensatory matrix that describes the redistribution of water uptake in a non-uniform hydraulic head profile. The two characteristics, Krs and SUF , are sufficient to describe the total uptake as a function of the collar and soil water potential, and water uptake redistribution does not depend on the total uptake or collar water potential. We compared the exact model with two hydraulic root models that make a priori simplifications of the hydraulic root architecture, i.e., the parallel and big root model. The parallel root model uses only two characteristics, Krs and SUF , which can be calculated directly following a bottom-up approach from the 3D hydraulic root architecture. The big root model uses more parameters than the parallel root model, but these parameters cannot be obtained straightforwardly with a bottom-up approach. The big root model was parameterized using a top-down approach, i.e., directly from root segment hydraulic properties, assuming a priori a single big root architecture. This simplification of the hydraulic root architecture led to less accurate descriptions of root water uptake than by the parallel root model. To compute root water uptake in macroscopic soil water flow and land surface models, we recommend the use of the parallel root model with Krs and SUF computed in a bottom-up approach from a known 3D root hydraulic architecture. [ABSTRACT FROM AUTHOR]

Coupled numerical models, which simulate water and energy fluxes in the subsurface–land-surface–atmosphere system in a physically consistent way, are a prerequisite for the analysis and a better understanding of heat and matter exchange fluxes at compartmental boundaries and interdependencies of states across these boundaries. Complete state evolutions generated by such models may be regarded as a proxy of the real world, provided they are run at sufficiently high resolution and incorporate the most important processes. Such a simulated reality can be used to test hypotheses on the functioning of the coupled terrestrial system. Coupled simulation systems, however, face severe problems caused by the vastly different scales of the processes acting in and between the compartments of the terrestrial system, which also hinders comprehensive tests of their realism. We used the Terrestrial Systems Modeling Platform (TerrSysMP), which couples the meteorological Consortium for Small-scale Modeling (COSMO) model, the land-surface Community Land Model (CLM), and the subsurface ParFlow model, to generate a simulated catchment for a regional terrestrial system mimicking the Neckar catchment in southwest Germany, the virtual Neckar catchment. Simulations for this catchment are made for the period 2007–2015 and at a spatial resolution of 400 m for the land surface and subsurface and 1.1 km for the atmosphere. Among a discussion of modeling challenges, the model performance is evaluated based on observations covering several variables of the water cycle. We find that the simulated catchment behaves in many aspects quite close to observations of the real Neckar catchment, e.g., concerning atmospheric boundary-layer height, precipitation, and runoff. But also discrepancies become apparent, both in the ability of the model to correctly simulate some processes which still need improvement, such as overland flow, and in the realism of some observation operators like the satellite-based soil moisture sensors. The whole raw dataset is available for interested users. The dataset described here is available via the CERA database (Schalge et al., 2020): 10.26050/WDCC/Neckar_VCS_v1. [ABSTRACT FROM AUTHOR]