Abstract:Accurately quantifying the mass (water, nutrients, and carbon) and energy exchange processes between the Earth's atmosphere, biosphere, and lithosphere requires the accurate parameterization of soil hydraulic properties (SHPs) and their spatial heterogeneity. Because direct measurements of SHPs are difficult, time-consuming, and impossible at larger spatial scales, various pedotransfer functions (PTFs) have been developed in the last few decades, providing divergent estimates of SHPs from readily measurable variables. However, existing PTFs are mostly developed for specific regions and may not be suitable for other pedoclimatic conditions. Here, PTFs were developed using multiple machine-learning algorithms to estimate SHPs and examined the weaknesses and strengths of each method in estimating the average and variability of SHPs across China's arid region. The optimal PTFs were applied to a 1 × 1 km 2 regional map of texture and bulk density, thus producing maps of the saturated hydraulic conductivity ( K s ), the parameters of the van Genuchten (VG) formulation, field capacity ( θ fc ), wilting point ( θ wp ), plant available water ( θ pa ), and soil macroporosity ( φ m ) in the 0–2 m soil profile throughout the region. The results indicate that the ensemble model with the averaging method (EMA) is the most robust for estimating all SHPs. The EMA-PTFs for K s yielded the best performance for sand textures, followed by sandy loam, loam, silty loam, silty clay loam, loamy sand, and clay loam textures. There were no significant differences in estimating soil water retention curve parameters among the different soil texture classes. Using the same observed data set, we demonstrated that the new EMA-PTFs outperformed those of existing PTFs, such as Rosetta and HYPRES, with RMSE values decreasing by 25–83 % depending on the SHPs. Furthermore, our SHP datasets exhibited significantly deeper soil profiles and higher accuracy than other available regional and global SHP products. The VG retention parameter α shows the greatest variation vertically throughout the 0–2 m soil profile, followed by ln( K s ), θ r , θ pa , θ fc , θ wp , φ m , θ s and n . All SHPs exhibit much lower variability in the desert than in other areas, likely due to the homogeneous particle size distribution of desert areas. This study provides more accurate SHP estimates in China's arid region than the existing SHP products by developing advanced PTFs and highlights the inconsistency in SHPs among different global or regional products; the uncertainties induced by PTFs should thus be considered in future terrestrial biosphere modeling.

Development of pedotransfer functions for predicting the bulk density in the critical zone on the Loess Plateau, China

Ensemble Modelling-Based Pedotransfer Functions for Predicting Soil Bulk Density in China

Machine Learning Based Pedotransfer Function Improves Soil Bulk Density Prediction but Not for Soil Organic Carbon Stock

Spatial non‐stationary effects of explanatory variables on soil bulk density in the critical zone of the Chinese Loess Plateau

Development of pedotransfer functions for predicting hydraulic parameters of van Genuchten model by incorporating environmental variables on the Qinghai-Tibet Plateau

European topsoil bulk density and organic carbon stock database (0–20 cm) using machine-learning-based pedotransfer functions

Pedogenic knowledge-aided modelling of soil inorganic carbon stocks in an alpine environment

Analytical modelling of soil porosity and bulk density across the soil organic matter and land-use continuum

A method for estimating the bulk density and particle density of granite residual soil based on the construction of pedotransfer functions

Using Organic Matter and pH to Estimate the Bulk Density of Afforested/Reforested Soils in Northwest and Northeast China

Developing novel ensemble models for predicting soil hydraulic properties in China's arid region

Distributions of glycerol dialkyl glycerol tetraethers in surface soils of Qinghai Tibetan Plateau: implications of GDGT based proxies in cold and dry regions.

Quantification of the Soil Stiffness Constants Using Physical Properties of Paddy Soils in Yangtze Delta Plain, China

Evaluation of pedotransfer functions for estimating saturated hydraulic conductivity in coastal salt-affected mud farmland

Estimating soil stress distribution by using depth‐dependent soil bulk‐density data

Vertical distribution of soil carbon, nitrogen, and phosphorus in typical Chinese terrestrial ecosystems

Variations in soil physico-chemical properties, soil stocks, and soil stoichiometry under different soil layers, the major forest region Liupan Mountains of Northwest China

Bulk Density of Soils in Relation to Soil Physical and Chemical Properties

Vertical distribution of soil available phosphorus and soil available potassium in the critical zone on the Loess Plateau, China

C:N:P stoichiometric characteristics and drivers of surface soil in the drylands of China

Topography Dominates the Spatial and Temporal Variability of Soil Bulk Density in Typical Arid Zones