In ombrotrophic peatlands, the moisture content of the vadose zone (acrotelm) controls oxygen diffusion rates, redox state, and the turnover of organic matter. Whether peatlands act as sinks or sources of atmospheric carbon thus relies on variably saturated flow processes. The Richards equation is the standard model for water flow in soils, but it is not clear whether it can be applied to simulate water flow in live Sphagnum moss. Transient laboratory evaporation experiments were conducted to observe evaporative water fluxes in the acrotelm, containing living Sphagnum moss, and a deeper layer containing decomposed moss peat. The experimental data were evaluated by inverse modeling using the Richards equation as process model for variably‐saturated flow. It was tested whether water fluxes and time series of measured pressure heads during evaporation could be simulated. The results showed that the measurements could be matched very well providing the hydraulic properties are represented by a suitable model. For this, a trimodal parametrization of the underlying pore‐size distribution was necessary which reflects three distinct pore systems of the Sphagnum constituted by inter‐, intra‐, and inner‐plant water. While the traditional van Genuchten‐Mualem model led to great discrepancies, the physically more comprehensive Peters‐Durner‐Iden model which accounts for capillary and noncapillary flow, led to a more consistent description of the observations. We conclude that the Richards equation is a valid process description for variably saturated moisture fluxes over a wide pressure range in peatlands supporting the conceptualization of the live moss as part of the vadose zone.

A generalized modular framework for partitioning soil hydraulic property (SHP) functions into a capillary and a noncapillary part is developed. The full water retention curve (WRC) is modeled as a weighted sum of a parametric capillary saturation function and a new general model for the noncapillary saturation function. This model is directly computed from any selected capillary saturation function. With it, a physically complete, continuous, and flexible representation of the WRC is achieved, ensuring zero water content at oven dryness. In a modular and hierarchical framework, the expressions for the capillary and noncapillary saturation function are used to calculate the respective hydraulic conductivity curves (HCC). This is achieved by adopting Mualem's integral for the capillary part of the HCC only and calculating the noncapillary HCC directly from the new noncapillary saturation function. This leads to consistent descriptions of measured HCC data, including the often observed change in slope beyond −100 cm pressure head. Compared to the classical van Genuchten-Mualem approach, it requires only one additional model parameter. The SHP framework model describes both WRC and HCC adequately and coherently. We demonstrate the suitability of the SHP framework and versatility by describing measured WRC and HCC data across the full moisture range using soil samples from a wide range of textures and origins. The modular framework was implemented in the soil physics and soil hydrology (spsh) R-package, available from the Comprehensive R Archive Network. It contains several SHP models, model parameter estimation, and features options for goodness of fit statistics, and model selection.

No abstract

Abstract. In ombrotrophic peatlands, the moisture content of the acrotelm (vadoze zone) controls oxygen diffusion rates, redox state, and the turnover of organic matter. Thus, variably saturated flow processes determine whether peatlands act as sinks or sources of atmospheric carbon, and modelling these processes is crucial to assess effects of changed environmental conditions on the future development of these ecosystems. We show that the Richards equation can be used to accurately describe the moisture dynamics under evaporative conditions in variably saturated peat soil, encompassing the transition from the topmost living moss layer to the decomposed peat as part of the vadose zone. Soil hydraulic properties (SHP) were identified by inverse simulation of evaporation experiments on samples from the entire acrotelm. To obtain consistent descriptions of the observations, the traditional van GenuchtenMualem model was extended to account for non-capillary water storage and flow. We found that the SHP of the uppermost moss layer reflect a pore-size distribution (PSD) that combines three distinct pore systems of the Sphagnum moss. For deeper samples, acrotelm pedogenesis changes the shape of the SHP due to the collapse of inter-plant pores and an infill with smaller particles. This leads to gradually more homogeneous and bi-modal PSDs with increasing depth, which in turn can serve as a proxy for increasing state of pedogenesis in peatlands. From this, we derive a nomenclature and size classification for the pore spaces of Sphagnum mosses and define inter-, intra-, and inner-plant pore spaces, with effective pore diameters of > 300, 300-30, and 30-10 µm, respectively.

There is sparse information on reactive solute transport in peat; yet, with increasing development of peatland dominated landscapes, purposeful and accidental contaminant releases will occur, so it is important to assess their mobility. Previous experiments with peat have only evaluated single-component solutions, such that no information exists on solute transport of potentially competitively adsorbing ions to the peat matrix. Additionally, recent studies suggest chloride (Cl) might not be conservative in peat, as assumed by many past peat solute transport studies. Based on measured and modelled adsorption isotherms, this study illustrates concentration dependent adsorption of Cl to peat occurred in equilibrium adsorption batch (EAB) experiments, which could be described with a Sips isotherm. However, Cl adsorption was insignificant for low concentrations (<500 mg L) as used in breakthrough curve experiments (BTC). We found that competitive adsorption of Na, K, and NH transport could be observed in EAB and BTC, depending on the dissolved ion species present. Na followed a Langmuir isotherm, K a linear isotherm within the tested concentration range (~10 - 1500 mg L), while the results for NH are inconclusive due to potential microbial degradation. Only Na showed clear evidence of competitive behaviour, with an order of magnitude decrease in maximum adsorption capacity in the presence of NH (0.22 to 0.02 mol kg), which was confirmed by the BTC data where the Na retardation coefficient differed between the experiments with different cations. Thus, solute mobility in peatlands is affected by competitive adsorption.

Agroecosystem models need to reliably simulate all biophysical processes that control crop growth, particularly the soil water fluxes and nutrient dynamics. As a result of the erosion history, truncated and colluvial soil profiles coexist in arable fields. The erosion-affected field-scale soil spatial heterogeneity may limit agroecosystem model predictions. The objective was to identify the variation in the importance of soil properties and soil profile modifications in agroecosystem models for both agronomic and environmental performance. Four lysimeters with different soil types were used that cover the range of soil variability in an erosion-affected hummocky agricultural landscape. Twelve models were calibrated on crop phenological stages, and model performance was tested against observed grain yield, aboveground biomass, leaf area index, actual evapotranspiration, drainage, and soil water content. Despite considering identical input data, the predictive capability among models was highly diverse. Neither a single crop model nor the multi-model mean was able to capture the observed differences between the four soil profiles in agronomic and environmental

Abstract. The underlying processes governing solute transport in peat from an experimentally constructed fen peatland were analyzed by performing saturated and unsaturated solute breakthrough experiments using Na + and Cl − as reactive and non-reactive solutes, respectively. We tested the performance of three solute transport models, including the classical equilibrium convection-dispersion equation (CDE), a chemical non-equilibrium one-site adsorption model (OSA) and a model to account for physical non-equilibrium, the mobile-immobile (MIM) phases. The selection was motivated by the fact that the applicability of the MIM in peat soils finds a wide consensus. However, results from inverse modeling and a robust statistical evaluation of this peat provide evidence that the measured breakthrough of the conservative tracer, Cl − , could be simulated well using the CDE. Furthermore, the very high Damköhler number (which approaches infinity) suggests instantaneous equilibration between the mobile and immobile phases underscoring the redundancy of the MIM approach for this particular peat. Scanning electron microscope images of the peat show the typical multi-pore size distribution structures have been homogenized sufficiently by decomposition, such that physical non-equilibrium solute transport no longer governs the transport process. This result is corroborated by the fact the soil hydraulic properties were adequately described using a unimodal van Genuchten-Mualem model between saturation and a pressure head of ∼ −1000 cm of water. Hence, MIM was not the most suitable choice, and the long tailing of the Na + breakthrough curve was caused by chemical non-equilibrium. Successful description was possible using the OSA model. To test our results for the unsaturated case, we conducted an unsaturated steady-state evaporation experiment to drive Na + and Cl − transport. Using the parameterized transport models from the saturated experiments, we could numerically simulate the unsaturated transport using Hydrus-1-D. The simulation showed a good prediction of observed values, confirming the suitability of the parameters for use in a slightly unsaturated transport simulation. The findings improve the understanding of solute redistribution in the constructed fen and imply that MIM should not be automatically assumed for solute transport in peat but rather should be evidence based.

Modeling soil hydraulic properties requires an effective representation of capillary and noncapillary storage and conductivity. This is made possible by using physically comprehensive yet flexible soil hydraulic property models. Such a model (Brunswick [BW] model) was introduced by Weber et al. (2019, https://doi.org/10.1029/2018WR024584), and it overcomes some core deficiencies present in the widely used van Genuchten-Mualem (VGM) model. We first compared the performance of the BW model to that of the VGM model in its ability to describe water retention and hydraulic conductivity data on a set of measurements of 402 soil samples with textures covering the entire range of classes. Second, we developed a simple transfer function to predict BW parameters based on VGM parameters. Combined with our new function, any existing pedotransfer function for the prediction of the VGM parameters can be extended to predict BW model parameters. Based on information criteria, the smaller variance of the residuals, and a 40% reduction in mean absolute error in the hydraulic conductivity over all samples, the BW model clearly outperforms VGM. This is possible as the BW model explicitly accounts for hydraulic properties of dry soils. With the new pedotransfer function developed in this study, better descriptions of water retention and hydraulic conductivities are possible. We are convinced that this will strengthen the utility of the new model and enable improved field-scale simulations, climate change impact assessments on water, energy and nutrient fluxes, as well as crop productivity in agroecosystems by soil-crop and land-surface modeling. The models and the pedotransfer function are included in an R package spsh (https://cran.r-project.org/package¼spsh). Plain Language Summary Soil hydraulic property models are mathematical functions, which describe the relationship between the soil water pressure head and the state of soil water saturation, on the one hand, and the soil water pressure head and the unsaturated soil hydraulic conductivity, on the other. These types of mathematical functions are flexible by adjustable parameters. With one set of model equations, the hydraulic properties of soils which may have very different properties due to their vast natural variability can be described. The models treated in this work are (i) the van Genuchten-Mualem model, a model with well-known problems, but still frequently applied, and (ii) a relatively new physical comprehensive model, named the Brunswick model. First of all, in a data-based comparison of model performance, we demonstrate that the Brunswick model has systematic advantages. Second, knowledge about these above-mentioned parameters can be determined through other mathematical functions, so-called hydro-pedotransfer functions, which empirically relate these parameters to observed soil properties. The information about these soil properties can be measured in the laboratory and is also recorded in soil maps. We created a new pedotransfer function to facilitate the predictio...