Saturated Conductivity Research Articles

Study of the effect of the compaction level on the hydrodynamic properties of loamy sand soil in an agricultural context

Agricultural soil compaction adversely affects crop water use and yield performance and should be avoided or remediated through appropriate soil management strategies. The investigation of the impact of different levels of soil compaction on its hydrodynamic properties remains a crucial step in improving water use and crop yields. We examined five compaction levels of silty sand soil sampled from a potato field in the agricultural regions of northern Quebec (Canada). Soil hydraulic characteristics (saturated and unsaturated hydraulic conductivity, soil water retention capacity) were measured using the constant head method, the HYPROP device, and a WP4C dew point potentiometer. The sixteen hydraulic models integrated into the HYPROP software were fitted to the soil water retention curve (SWRC) data for the studied compaction levels. Statistical parameters such as the mean bias error, mean absolute error, correlation coefficient, and root mean square error were used to measure the performance of the models. The results show that saturated and unsaturated conductivity decreases with increasing soil compaction. The lowest saturated hydraulic conductivity (Ks) value is observed for the highest level of soil compaction, reflecting a solid medium with less pore space and connectivity. Among the hydraulic models, the Peters-Durner-Iden (PDI) variant of van Genuchten's unconstrained bimodal model (VGm-b-PDI) outperformed all other models for SWRC simulation of different soil compaction levels and was, accordingly, selected as the optimal model. This model was implemented in HYDRUS-1D to estimate the amount of irrigation for different compaction levels. We simulated irrigation scenarios with the dual-porosity model. The results indicated that soil compaction can strongly influence soil hydraulic properties and water differently. However, the amount of irrigation for the potato crop was optimal at a moderate level of soil compaction. Overall, combined HYPROP and HYDRUS 1D can provide helpful information on the soil hydraulics properties dynamics and a rigorous simulation for irrigation planning and management in potato fields.

Modeling on effective thermal conductivity of hydrate-bearing sediments considering the shape of sediment particle

The estimation of effective thermal conductivity (ETC) of hydrate-bearing sediment is of great importance for the exploitation of natural gas hydrate. At present, existing ETC models have been lack of consideration of the influences of sediment particle shape and stacking patterns. Thus, in this work, new models considering the shape (octahedral, plate-shaped, spindle-shaped and elongated) of sediment particles and stacking patterns were established and verified by a series of existed laboratory data. The verification results indicate that the octahedral model, plate-shape model and the average values of the four models have high accuracy in predicting ETC of hydrate-bearing sediments. The ETC of hydrate-bearing sediment under water-saturated state were calculated by the models in this work. The average ETC predicted by four models decreased from 2.87 W⋅m-1⋅K-1 to 1.47 W⋅m-1⋅K-1 with porosity increasing from 0.22 to 0.47 and decreased from 1.61 W⋅m-1⋅K-1 to 1.58 W⋅m-1⋅K-1 with hydrate saturation increasing from 0 to 1. Finally, the effects of porosity, hydrate saturation and intrinsic thermal conductivity of sediment particles on the ETC were discussed and analyzed. This modeling work could contribute to the understanding and prediction of the ETC of hydrate-bearing sediment with different shapes of sediment particles in the complex porous system.

Excluding quartz content from the estimation of saturated soil thermal conductivity: Combined use of differential effective medium theory and geometric mean method

Saturated soil thermal conductivity (λsat) is the maximum soil thermal conductivity value of a given soil. Although it can be determined accurately with a heat pulse sensor, there are challenges to prepare fully saturated soil samples. Numerous models have been developed to estimate λsat, and among these, the geometric mean method (GMM) generally performs well. The GMM requires soil mineral composition or quartz content information, which is unavailable for most soils. Earlier studies commonly used assumed that quartz content (fquartz) was equal to sand content (fsand) or to 0.5 × fsand, which led to significant λsat estimation errors especially on coarse-textured soils. We derived a novel method to estimate λsat from soil porosity (ϕ) based on a combination of the GMM and differential effective medium theory (DEM). The new DEM-GMM approach has a single parameter, cementation exponent (m). Using a calibration dataset of 43 soils, we determined best fit m values for soils in three groups: 1.66 for Group I (fsand < 0.4), 1.62 for Group II (0.4 <= fsand < 1) and m = -1.34ϕ+1.70 for Group III (fsand = 1). Using best fit m values for different groups, the new model can estimate λsat values from ϕ. Independent validation results on another 46 soils showed that the new model outperformed the GMM method with the assumption that fquartz = fsand or fquartz = 0.5 × fsand. The mean RMSE, Bias and R2 values of the DEM-GMM approach were 0.202 W m−1 K−1, 0.013 W m−1 K−1 and 0.89, respectively, and corresponding values of the GMM with the two assumptions were 0.295 and 0.476 W m−1 K−1, 0.056 and -0.28 W m−1 K−1, 0.80 and 0.82, respectively. The robust performance of the DEM-GMM approach suggests that it can be incorporated into thermal conductivity models to accurately estimate the thermal conductivity of unsaturated soils.

Soil water retention and hydraulic conductivity measured in a wide saturation range

Abstract. Soil hydraulic properties (SHPs), particularly soil water retention capacity and hydraulic conductivity of unsaturated soils, are among the key properties that determine the hydrological functioning of terrestrial systems. Some large collections of SHPs, such as the UNSODA and HYPRES databases, have already existed for more than 2 decades. They have provided an essential basis for many studies related to the critical zone. Today, sample-based SHPs can be determined in a wider saturation range and with higher resolution by combining some recently developed laboratory methods. We provide 572 high-quality SHP data sets from undisturbed, mostly central European samples covering a wide range of soil texture, bulk density and organic carbon content. A consistent and rigorous quality filtering ensures that only trustworthy data sets are included. The data collection contains (i) SHP data, which consist of soil water retention and hydraulic conductivity data, determined by the evaporation method and supplemented by retention data obtained by the dewpoint method and saturated conductivity measurements; (ii) basic soil data, which consist of particle size distribution determined by sedimentation analysis and wet sieving, bulk density and organic carbon content; and (iii) metadata, which include the coordinates of the sampling locations. In addition, for each data set, we provide soil hydraulic parameters for the widely used van Genuchten–Mualem model and for the more advanced Peters–Durner–Iden model. The data were originally collected to develop and test SHP models and associated pedotransfer functions. However, we expect that they will be very valuable for various other purposes such as simulation studies or correlation analyses of different soil properties to study their causal relationships. The data are available at https://doi.org/10.5880/fidgeo.2023.012 (Hohenbrink et al., 2023).

Inspecting the role of vapour loss and other model strategies in the mod-elling of a bentonite thermo-hydraulic cell

Inspecting the role of vapour loss and other model strategies in the mod-elling of a bentonite thermo-hydraulic cell

Experimental research on the thermal conductivity of unsaturated rocks in geothermal engineering

The thermal conductivity of rock is an essential thermophysical parameter in geothermal engineering. Experiments have shown that when rock changes from an unsaturated to a saturated state, thermal conductivity can be increased by 50% at most. In this paper, the variation of thermal conductivity with saturation of 101 samples of six different rock types from Songliao Basin, Gonghe Basin, and Ordos Basin in China was measured in the laboratory. Correlations between saturation, porosity, internal structural characteristics of rocks and thermal conductivity were investigated. The results show that the effect of saturation on the thermal conductivity of rock increases exponentially with porosity. Based on three classical models and the relationship between porosity, saturation, and thermal conductivity, improved models and a fitted model of thermal conductivity are proposed. By comparing the improved models to the fitted model, the applicable scope of each model is obtained. The errors of each model are analyzed in detail. The average MAE value of the improved models is 0.405 and the average RMSE value is 1.207 over the applicable model scopes. The R2 of the fitted model is 0.94. Compared with experimental data from literature, MAE values of the improved models and the fitted model are 0.040 and 0.003 respectively. Unsaturated rock with a single fracture heat transfer model is established by simulation. A thermal-Hydrological coupling model is used to simulate seepage heat transfer through fracture. The effect of variable thermal conductivity simulation is the best with the experimental value, and the MAPE value was 8.3%. This study fills a gap in available models of unsaturated thermal conductivity of rock, and provide theoretical support for the accurate calculation of heat transfer efficiencies in related geothermal engineering applications.

Characterization of DNAPL source zones in clay-sand media via joint inversion of DC resistivity, induced polarization and borehole data

Toxic organic contaminants in groundwater are pervasive at many industrial sites worldwide. These contaminants, such as chlorinated solvents, often appear as dense non-aqueous phase liquids (DNAPLs). To design efficient remediation strategies, detailed characterization of DNAPL Source Zone Architecture (SZA) is required. Since invasive borehole-based investigations suffer from limited spatial coverage, a non-intrusive geophysical method, direct current (DC) resistivity, has been applied to image the DNAPL distribution; however, in clay-sand environments, the ability of DC resistivity for DNAPLs imaging is limited since it cannot separate between DNAPLs and surrounding clay-sand soils. Moreover, the simplified parameterization of conventional inversion approaches cannot preserve physically realistic patterns of SZAs, and tends to smooth out any sharp spatial variations. In this paper, the induced polarization (IP) technique is combined with DC resistivity (DCIP) to provide plausible DNAPL characterization in clay-sand environments. Using petrophysical models, the DCIP data is utilized to provide tomograms of the DNAPL saturation (SN) and hydraulic conductivity (K). The DCIP-estimated K/SN tomograms are then integrated with borehole measurements in a deep learning-based joint inversion framework to accurately parameterize the highly irregular SZA and provide a refined DNAPL image. To evaluate the performance of the proposed approach, we conducted numerical experiments in a heterogeneous clay-sand aquifer with a complex SZA. Results demonstrate the standalone DC resistivity method fails to infer the DNAPL in complex clay-sand environments. In contrast, the combined DCIP technique provides the necessary information to reconstruct the large-scale features of K/SN fields, while integrating DCIP data with sparse but accurate borehole data results in a high resolution characterization of the SZA.

Numerical simulation of the effective thermal conductivity of cement-based composites: Microscopic insights into diatomite/silica aerogel

Silica aerogel used in cement composites to prepare novel thermal superinsulation material has been gradually applied in buildings and can effectively reduce building energy consumption. In this paper, the geometric voxel construction method is used to heterogeneity model the thermal conductivity of a cement-based diatomite/silica aerogel composite (DSAC), and thermal conductivity simulations are carried out using a heat conduction model framework analogous to an electrical resistance model framework. The simulation results show that the simulated thermal conductivities of cementitious DSAC materials with different DSAC composite contents agree well with the corresponding experimental values. The saturated thermal conductivity of the cementitious DSAC material is 2.5 times the dry thermal conductivity at 28 d. This shows that the thermal conductivity model proposed in this paper can effectively predict the effective thermal conductivity of cementitious DSAC materials and can provide strong support for the material design and performance prediction of cementitious DSAC materials.

A new model to predict soil thermal conductivity

Thermal conductivity is a basic parameter of soil heat transferring, playing an important role in many fields including groundwater withdrawal, ground source heat pump, and heat storage in soils. However, it usually requires a lot of time and efforts to obtain soil thermal conductivity. To conveniently obtain accurate soil thermal conductivity, a new model describes the relationship between soil thermal conductivity (λ) and degree of saturation (Sr) was proposed in this study. Dry soil thermal conductivity (λdry) and saturated soil thermal conductivity (λsat) were described using a linear expression and a geometric mean model, respectively. A quadratic function with one constant was added to calculate λ beyond the lower λdry and upper λsat limit conditions. The proposed model is compared with five other frequently used models and measured data for 51 soil samples ranging from sand to silty clay loam. Results show that the proposed model match the measured data well. The proposed model can be used to determine soil thermal conductivity of a variety of soil textures over a wide range of water content.

Dynamic monitoring of dielectric properties during two phase immiscible displacements in sand packs using frequency domain electromagnetic sweeps

Electromagnetic heating is a promising soil remediation method especially in thin formations. The lack of a wide-spread adoption of this method stems from insufficient knowledge of how the complex dielectric properties, that govern propagation of the electromagnetic waves through porous media, change with changing frequency, water saturation, displacement types and flow regimes. To breach these gaps several sets of spontaneous deionized (DI) water imbibition experiments, followed by the primary drainage floods, that were followed by the secondary DI water imbibition floods in confined uniform sand packs were performed. The frequency domain relative dielectric constant and conductivities were extracted from the two-port complex S-parameter measurements taken with the vector network analyzer during these immiscible displacements at various water saturation levels at ambient conditions. A novel coaxial transmission line core holder was designed and commissioned, and a modified version of a plane-invariant dielectric extraction algorithm was developed for this purpose. Series, parallel and semi-disperse mixing models were applied to fit the water saturation dependent relative dielectric constant and conductivity values sampled at 500 MHz from the extracted frequency domain spectra. The Maxwell-Garnett parallel model was proved to be the most flexible because it could capture the sampled conductivity values in all secondary imbibition floods before and after the breakthroughs, where the inflection points were observed. These inflection points were attributed to silica production and a potential shear-stripping flow. This observation was further confirmed by conducting a single-phase Darcy's law analysis of two DI water imbibition floods.

Global transonic solution of hot accretion flow with thermal conduction

ABSTRACT We examine the effect of thermal conduction on the low-angular momentum hot accretion flow (HAF) around non-rotating black holes accreting mass at very low rate. While doing so, we adopt the conductive heat flux in the saturated form, and solve the set of dynamical equations corresponding to a steady, axisymmetric, viscous, advective accretion flow using numerical methods. We study the dynamical and thermodynamical properties of accreting matter in terms of the input parameters, namely energy (ε0), angular momentum (ℓ0), viscosity parameter (α), and saturation constant (Φs) regulating the effect of thermal conduction. We find that Φs plays a pivotal role in deciding the transonic properties of the global accretion solutions. In general, when Φs is increased, the critical point (rc) is receded away from the black hole, and flow variables are altered particularly in the outer part of the disc. To quantify the physically acceptable range of Φs, we compare the global transonic solutions with the self-similar solutions, and observe that the maximum saturation constant ($\Phi ^{\rm max}_{\rm s}$) estimated from the global solutions exceeds the saturated thermal conduction limit (Φsc) derived from the self-similar formalism. Moreover, we calculate the correlation between α and $\Phi ^{\rm max}_{\rm s}$ and find ample disagreement between global solutions and self-similar solutions. Further, using the global flow variables, we compute the Bernoulli parameter (Be) which remains positive all throughout the disc, although flow becomes loosely unbound for higher Φs. Finally, we indicate the relevance of this work in the astrophysical context in explaining the possibility of mass-loss/outflows from the unbound disc.

Analysis of Water Infiltration Characteristics and Hydraulic Parameters of Sierozem Soil under Humic Acid Addition

The farmland in Yinchuan is composed of sierozem soil, which is characterized by high sand content and low organic matter content, resulting in poor water-holding capacity and weak soil structure. Humic acid is a natural organic polymer soil amendment. It is critical to study how humic acid affects soil water infiltration in sierozem soil at the microlevel. A one-dimensional vertical infiltration experiment was conducted to explore how adding different amounts of humic acid (0, 1%, 2%, 3% and 4%) affected the infiltration characteristics and hydraulic parameters of the sierozem soil. The results revealed that the wetting front and cumulative infiltration decreased with the increase in humic acid addition. When the infiltration time was 90 min, the wetting front of the 1%, 2%, 3% and 4% treatments was 6.50%, 10.00%, 15.00% and 21.00% lower than CK (0 for CK), and the cumulative infiltration volume was 4.50%, 11.14%, 18.42% and 23.60% lower than CK, respectively. Among the three infiltration models created by Philip, Horton and Kostiakov, the Kostiakov model (R2 &gt; 0.95) could more accurately describe the soil water infiltration process in the study area. After infiltration, the moisture content of each soil layer increased with the increase in humic acid, which improved the water-holding capacity of the sierozem soil. Using Hydrus-1D to calculate soil hydraulic parameters, we found that the humic acid addition affected the hydraulic parameters. With the increase in the amount of humic acid addition, the retention water content θr and saturated water θs were positively correlated with the humic acid addition amount and negatively correlated with the saturated water conductivity Ks and the reciprocal of air-entry α. The results showed that humic acid could increase the water-holding capacity of soil and improve the rapid water loss and poor water-holding capacity of sierozem soil.

Prediction of the absolute hydraulic conductivity function from soil water retention data

Abstract. For modeling flow and transport processes in the soil–plant–atmosphere system, knowledge of the unsaturated hydraulic properties in functional form is mandatory. While much data are available for the water retention function, the hydraulic conductivity function often needs to be predicted. The classical approach is to predict the relative conductivity from the retention function and scale it with the measured saturated conductivity, Ks. In this paper we highlight the shortcomings of this approach, namely, that measured Ks values are often highly uncertain and biased, resulting in poor predictions of the unsaturated conductivity function. We propose to reformulate the unsaturated hydraulic conductivity function by replacing the soil-specific Ks as a scaling factor with a generally applicable effective saturated tortuosity parameter τs and predicting total conductivity using only the water retention curve. Using four different unimodal expressions for the water retention curve, a soil-independent general value for τs was derived by fitting the new formulation to 12 data sets containing the relevant information. τs was found to be approximately 0.1. Testing of the new prediction scheme with independent data showed a mean error between the fully predicted conductivity functions and measured data of less than half an order of magnitude. The new scheme can be used when insufficient or no conductivity data are available. The model also helps to predict the saturated conductivity of the soil matrix alone and thus to distinguish between the macropore conductivity and the soil matrix conductivity.

Few large and many small, how to control the thermal behavior of composite materials

Here we show that the distribution of the pore sizes within a solid matrix can contribute to improving the thermal insulation of composite materials. Building on the analogy between species diffusion through saturated porous materials and heat conduction through composite materials made of a solid matrix and air spheres, we develop a methodology to allocate the air space in such a way that the thermal conductivity decreases. The methodology relies on the few-large-and-many-small approach in which the largest air spheres are step-by-step replaced by smaller ones while maintaining the overall porosity constant. Finally, a numerical example highlights the potential of such composite material to act as a thermal barrier in high temperature applications.

Coupled pressure-driven flow and spontaneous imbibition in shale oil reservoirs

Coupled pressure-driven (viscous) flow and spontaneous imbibition are the main regimes during shale oil production. Revealing the unclear mechanisms of this coupled flow is a major concern for scholars and field engineers. In this work, the oil–water flow mechanisms within shale pore structures are investigated by pore-scale modeling methods in focused ion beam scanning electron microscopy digital rocks enhanced by applying super-resolution reconstruction (SRR). More small pores are identified with SRR, and the connectivity is improved. The enhanced pore size distribution is consistent with the nitrogen adsorption measurement; hence, more representative capillary pressure and relative permeability curves are obtained with essential experimental measurements. Then, an analytical solution of coupled pressure-driven (viscous) flow and spontaneous imbibition is derived, and a corresponding algorithm is proposed. Based on the pore-scale calculated relative permeability and capillary pressure curves, the analytical solution is applied to investigate the variations in water saturation profiles and conductance of the oil phase during the shale reservoir development. The results demonstrate that most of the shale oil is recovered by pressure dropdown-induced viscous flow and that imbibition is a minor factor. The overall oil-relative permeability decreases due to imbibition invasion. When the fracture spacing increases, the impairment of the overall oil-relative permeability decreases.

Effect of Different Soil Organic Carbon Content in Different Soils on Water Holding Capacity and Soil Health

An experiment was carried out to study the effect of soil organic carbon (SOC) and soil texture on the distance of the wetting front, cumulative water infiltration (I), infiltration rate (IR), saturated water conductivity (Ks), and water holding capacity (WHC). Three levels ( 0, 10, 20, and 30 g OC kg-1 ) from organic carbon (OC) were mixed with different soil materials sandy, loam, and clay texture soils. Field capacity (FC) and permanent wilting point (PWP) were estimated. Soil materials were placed in transparent plastic columns(12 cm soil column ), and water infiltration(I) was measured as a function of time, the distance of the wetting front and Ks. Results showed that advance wetting front as a function of time for soil column was 6 minutes and with no differences between OC levels for sandy soils, while it ranged between 90 minutes (0% OC) - 130 minutes (3% OC) for loam soils, and between 470 minutes (0 %OC) and 590 minutes (1%OC) for clay soils, at the same time cumulative water infiltration(I) increases at the beginning of infiltration and decreases with time and levels of OC. The highest infiltration values were in sandy soils, giving data of 0.05 and 0.12 cm min-1, with no significant differences with OC rates. IR values decreased when OC increased in loam soils, and IR increased exponentially in clay soils with increasing OC levels. The values of Ks decrease with increasing OC for sandy and loam soils, and increase when OC increases above 3% for clay soils. FC and WP values were increased for sandy, loam and clay soils when OC was increased. The AW values decreased for both sandy and clay soils compared to loam soils. It can be concluded that AW can be estimated from FC values regardless of texture and OC by the linear function: AW=0.51(FC)+0.005.

Impacts of liquid phase distribution on the effective thermal conductivity of closed-cell thermal insulation

This paper proposes a multiscale effective thermal conductivity model combined with liquid phase distribution to explore in detail the influence of the internal structure of closed-cell thermal insulation on the thermal behavior in the presence of moisture. By simulating the liquid phase distribution in the cracks between the insulation cell skeleton and combining it with the representative elementary volume (REV) scale and macroscopic scale thermal conductivity models, the correlations between the structural parameters, saturation, and effective thermal conductivity of phenolic thermal insulation are analyzed. A critical saturation point is proposed for the first time, and the related impact factors are discussed. For the heat and moisture transfer in porous media, the larger the critical saturation, the better the insulation performance. Based on the newly developed improved effective thermal conductivity model, nine groups of phenolic samples with different structural parameters were investigated to compare the impact of porosity, closed-open pore volume ratio, relative deviation ratio, contact angle, and closed-pore diameter on the effective thermal conductivity. Simulation results indicate that the influence of the liquid phase distribution is mainly reflected at low saturation at the REV scale; at the macroscale, the influence of the liquid phase distribution is significantly reduced.

Salinity-specific stomatal conductance model parameters are reduced by stomatal saturation conductance and area via leaf nitrogen

Modeling stomatal behavior is necessary for accurate stomatal simulation and predicting the terrestrial water‑carbon cycle. Although the Ball-Berry and Medlyn stomatal conductance (gs) models have been widely used, variations and the drivers of their key slope parameters (m and g1) remain poorly understood under salinity stress. We measured leaf gas exchange, physiological and biochemical traits, soil water content and electrical conductivity of saturation extract (ECe), and fitted slope parameters of two genotypes of maize growing in two water and two salinity levels. We found m was different between the genotypes, but no difference in g1. Salinity stress reduced m and g1, saturated stomatal conductance (gsat), the fraction of leaf epidermis area allocation to stomata (fs), and leaf nitrogen (N) content, and increased ECe, but no marked decrease in slope parameters under drought. Both m and g1 were positively correlated with gsat, fs, and leaf N content, and negatively correlated with ECe in the same fashion among the two genotypes. Salinity stress altered m and g1 by modulating gsat and fs via leaf N content. The prediction accuracy of gs was improved using salinity-specific slope parameters, with root mean square error (RMSE) being decreased from 0.056 to 0.046 and 0.066 to 0.025 mol m−2 s−1 for the Ball-Berry and Medlyn models, respectively. This study provides a modeling approach to improving the simulation of stomatal conductance under salinity.

Conversion of Forest Hillslopes into Tea Fields Increases Soil Nutrient Losses through Surface Runoff

The conversion of forest hillslopes into tea fields changes the slope structure, soil properties, and vegetation, which may increase the great variations in the nitrogen (N) and phosphorus (P) losses through the surface runoff and induce agricultural non-point source pollution, particularly during the first years of conversion. Therefore, a natural forest system and a newly reclaimed tea system were compared to determine the land use conversion effects on the surface runoff and resulting N and P losses during 2013–2016. The results showed that the soil bulk density (BD) at all the depths and the saturated soil conductivities (Ks) in the surface soil were higher in the tea system than in the forest system, indicating that great land use conversion impacts soil’s physical and hydraulic properties. The annual surface runoff volumes and ratios were significantly higher in the tea system than in the forest system (p &lt; 0.05), and this increased surface runoff was largely attributed to the slope structure and soil property changes. The total P (TP), dissolved P (DP), and particle P (PP) concentrations in the surface runoff volumes were similar in the two systems, while the total N (TN) concentrations were higher in the tea system than in the forest system due to the high N fertilizer application rates in the tea system. The percentages of dissolved N (DN) in the TN and particle N (PN) in the TN were, respectively, lower and higher in the tea system than in the forest system, which was associated with increased soil erodibility due to the land use conversion. The mean annual N and P losses through the surface runoff were higher in the tea system (11.07 ± 4.03 kg N ha−1 yr−1 and 1.73 ± 0.19 kg P ha−1 yr−1) than in the forest system (4.51 ± 0.66 kg N ha−1 yr−1 and 0.78 ± 0.18 kg P ha−1 yr−1). The N and P losses exponentially correlated with the surface runoff in both systems; however, the coefficients of determination (R2) were higher in the forest system than in the tea system, indicating that the exponential relationship between the surface runoff and the N and P losses should be more stable in the forest system than in the newly reclaimed tea system. Therefore, more effort should be dedicated to reducing the N and P losses though the surface runoff when forest hillslopes are converted into tea fields, particularly during the first three years of conversion.

Synergistic effects of biochar and carboxymethyl cellulose sodium (CMC) applications on improving water retention and aggregate stability in desert soils

Making improvements to the water-holding characteristics and water-erosion resistance of desert soils, particularly in inland extremely arid areas, is vital for achieving both sustainable water resource utilisation and food security. The aim of this study is to evaluate the effects of the co-application of biochar and carboxymethyl cellulose sodium (CMC) on the physical properties of sandy desert soil, including infiltration rate, saturated water conductivity, field water-holding capacity and aggregate stability. Sandy desert soil samples were collected from jujube plantations on the southern edge of the Taklimakan Desert in the Hotan Prefecture, Xinjiang, China. Five CMC application ratios (C0:0, C1:0.01 g/kg, C2:0.02 g/kg, C3:0.04 g/kg and C4:0.08 g/kg) and five biochar application ratios (B0:0, B1:1.0 g/kg, B2:2.0 g/kg, B3:4.0 g/kg and B4:8.0 g/kg) were designed and a total of 11 experimental treatments were performed, which were labelled as CK (control group), B2C0, B2C1, B2C2, B2C3, B2C4, B4C4, B0C2, B1C2, B3C2 and B4C2. Compared with CK, the combined application of biochar and CMC reduced the soil bulk density (BD) by 1.29–9.41% and the saturated hydraulic conductivity (Ks) by 29.64–94.98%, and increased the soil saturated water content (SSWC) by 8.81–30.74% and the water holding capacity (WHC) by 13.91–36.87%. Similarly, the water-stable aggregates that were co-applied with biochar and CMC increased by 29.10–256.86%. This resulted in significant improvement in the stability of sandy desert soil against water erosion. The principal component analysis (PCA) results found B4C4 to have the best comprehensive improvement effect. Therefore, 0.08 g/kg of CMC and 8.0 g/kg of biochar were used as recommended for improving the hydraulic properties of desert soils. Generally, CMC and biochar have a mutually complementary effect on improving sandy desert soil, providing new ideas and approaches for the improvement of soil and the sustainable development of agriculture in desert areas.