New insights into the parameterization of the dry surface layer and its hydrogeochemical mechanism: An experimental study

Knowledge of the parameterization of the dry surface layer (DSL) is essential for evaluating near-surface water flow and water balance in arid and semi-arid areas. Existing studies have parameterized DSL thickness and vapor flow as functions of the soil moisture content (SMC) in the surface layer to predict soil evaporation. However, hydrochemical processes related to DSL development have been ignored, including changes in hydrochemistry, the underlying hydrochemical mechanism, and the role of dissolved substances in the DSL development. Herein, we performed a series of soil evaporation experiments for 260 days and explored the factors influencing DSL development (e.g., soil texture, atmospheric temperature, SMC, solutes). Evaporation experiments were performed using silty loess, sandy loess, and fine sand with a 60-cm water table. Results showed that the cumulative evaporation of silty loess, sandy loess, and fine sand over the experimental period were 1,391.52, 460.10, and 185.53 mm, respectively, which determined by the maximum height of liquid flow continuity. The content of total dissolved solids (TDS) and major ions at the surface soil were significantly higher than the values at deep depths of 5‒55 cm, which largely depend on evaporative water loss. Evolutionary trends of chemical facies in sand media along the liquid water migration were from HCO₃-Ca type to SO₄·Cl-Na type. This was attributed to mineral dissolution at a depth of 5–55 cm and their transport with liquid water, resulting in the precipitation of salt crystals at the surface soil. Furthermore, a consolidated DSL with a thickness of 3.0–3.5 cm in the sandy loess and a loose DSL with a thickness of 1.5–2.0 cm in fine sand were observed at the end of the experiments. The accumulation of solutes at the surface leads to a reduction in effective porosity and the aggregation of soil particles during continuous drying, which facilitates the consolidation of DSL in sandy loess. This overestimated the DSL thickness, resulting in a difference between the experimental and predicted evaporation rates by Fick's law. Overall, these results highlight the limitations of considering DSL thickness as a function of SMC only, providing new insights into hydrochemical processes and dissolved solutes involving DSL parameterization during continuous soil drying.