A land surface model incorporated with soil freeze/thaw and its application in GAME/Tibet

Land surface process is of great importance in global climate change, moisture and heat exchange in the interface of the earth and atmosphere, human impacts on the environment and eco- system, etc. Soil freeze/thaw plays an important role in cold land surface processes. In this work the diurnal freeze/thaw effects on energy partition in the context of GAME/Tibet are studied. A sophisti- cated land surface model is developed, the particular aspect of which is its physical consideration of soil freeze/thaw and vapor flux. The simultaneous water and heat transfer soil sub-model not only reflects the water flow from unfrozen zone to frozen fringe in freezing/thawing soil, but also demon- strates the change of moisture and temperature field induced by vapor flux from high temperature zone to low temperature zone, which makes the model applicable for various circumstances. The modified Picard numerical method is employed to help with the water balance and convergence of the numerical scheme. Finally, the model is applied to analyze the diurnal energy and water cycle char- acteristics over the Tibetan Plateau using the Game/Tibet datasets observed in May and July of 1998. Heat and energy transfer simulation shows that: (i) There exists a negative feedback mechanism between soil freeze/thaw and soil temperature/ground heat flux; (ii) during freezing period all three heat fluxes do not vary apparently, in spite of the fact that the negative soil temperature is higher than that not considering soil freeze; (iii) during thawing period, ground heat flux increases, and sensible heat flux decreases, but latent heat flux does not change much; and (iv) during freezing period, soil temperature decreases, though ground heat flux increases.     

Evaporation from seasonally frozen bare and vegetated ground at various groundwater table depths in the Ordos Basin,Northwest China

In cold climates, the process of freezing–thawing significantly affects the ground surface heat balance and water balance. To better understand the mechanism of evaporation from seasonally frozen soils, we performed field experiments at different water table depths on vegetated and bare ground in a semiarid region in China. Soil moisture and temperature, air temperature, precipitation, and water table depths were measured over a 5‐month period (November 1, 2016, to March 14, 2017). The evaporation, which was calculated by a mass balance method, was high in the periods of thawing and low in the periods of freezing. Increased water table depth in the freezing period led to high soil moisture in the upper soil layer, whereas lower initial groundwater levels during freezing–thawing decreased the cumulative evaporation. The extent of evaporation from the bare ground was the same in summer as in winter. These results indicate that a noteworthy amount of evaporation from the bare ground is present during freezing–thawing. Finally, the roots of Salix psammophila could increase the soil temperature. This study presents an insight into the joint effects of soil moisture, temperature, ground vegetation, and water table depths on the evaporation from seasonally frozen soils. Furthermore, it also has important implications for water management in seasonally frozen areas.     

A contemporary winter ground thermal profile in the Lesotho highlands and implications for active and relict soil frost phenomena

A 2 m deep ground thermal pro?le is constructed from temperature data collected over the winter and spring of 2000 at 3220 m a.s.l. near the Thabana Ntlenyana summit (3482 m) in Lesotho, southern Africa. The zero isotherm is found to have penetrated to 0·16 m soil depth. Ground remained frozen at 0·05 m for a total of 79 days and for shorter periods at 0·02 m and 0·10 m. Diurnal freezing and thawing is restricted to the upper 0·10 m and conforms to the observed depth of active micro‐patterned ground found in the region. Holocene temperature depressions projected along the thermal pro?le can account for freezing down to 0·45 to 0·65 m. Deeper sorting to 1 m, evident from relict patterned ground near the logger site, corresponds to at least a 2·5 °C temperature depression and such landforms are evidently pre‐Holocene. Projections indicate a seasonal freezing depth exceeding 2 m during the Pleistocene Last Glacial Maximum although the existence of permafrost appears unlikely. Copyright © 2003 John Wiley & Sons, Ltd.     

Hysteretic freezing characteristics of riparian peatlands in the Western Boreal Forest of Canada

Freezing characteristics were investigated for a sedge covered floating fen and spruce covered swamp located beside a shallow lake in the Western Boreal Forest of Canada. Thermal properties were measured in situ for one freeze‐thaw cycle, and for two freeze‐thaw cycles in laboratory columns. Thermal conductivity and liquid water content were related to a range of subsurface temperatures above and below the freezing thresholds, and clearly illustrate hysteresis between the freezing and thawing process. Thermal hysteresis occurs because of the large change in thermal conductivity between water and ice, high water content of the peat, and wide variation in pore sizes that govern ice formation. Field and laboratory results were combined to develop linear freezing functions, which were tested in a heat transfer model. For surface temperature boundary conditions, subsurface temperatures were simulated for the over‐winter period and compared with field measurements. Replication of the transient subsurface thermal regime required that freezing functions transition gradually from thawed to frozen state (spanning the ?0·25 to ?2 °C range) as opposed to a more abrupt step function. Subsurface temperatures indicate that the floating fen underwent complete phase change (from water to ice) and froze to approximately the same depth as lake ice thickness. Therefore, the floating fen peatland froze as a ‘shelf’ adjacent to the lake, whereas the spruce covered swamp had a higher capacity for thermal buffering, and subsurface freezing was both more gradual and limited in depth. These thermal properties, and the timing and duration of frozen state, are expected to control the interaction of water and nutrients between surface water and groundwater, which will be affected by changes in air temperature associated with global climate change. Copyright © 2009 John Wiley & Sons, Ltd.     

The influence of freeze–thaw cycles of active soil layer on surface runoff in a permafrost watershed

As a result of global warming, the discharges from rivers in permafrost regions have varied significantly. However, its mechanism remains unclear. One of possible factors is active soil freeze–thaw cycle, which may influence surface runoff in the variation of permafrost water cycle processes. In this study, a typical permafrost watershed in the Qinghai-Tibet plateau was selected, its hydrological processes were monitored from 2004 to 2007, and the effects of the freezing and thawing depth of the soil active layer on runoff processes were assessed. The runoff modulus, runoff coefficient, direct runoff ratio, recession gradient and their seasonal variations were estimated and analyzed. The active soil dynamics and water budget were analyzed to prove the features of the surface runoff and the influences of active soil freeze–thaw processes. The primary factors influencing surface runoff processes during different seasons were analyzed by Principal Component Analysis (PCA) and statistical regression methods. The results showed that the high runoff coefficient and low direct runoff ratio were the main characteristics during the spring flood period (May–June) and during the autumn recession period (September). The runoff modulus and its year-to-year variability were the greatest in the summer flood period. The direct runoff ratio decreased from 0.43 in May to 0.29 in September, with the exception of the highest ratio, which occurred during the summer recession period (July). The active soil thawing in the upper layer of depth of 60 cm had contributed to increase in discharge, but the increase in thawing depth deeper than 60 cm led to a decrease in surface runoff and slowness in the recession process. Precipitation played a small role in the spring flood runoff and the autumn runoff. The soil active layer freeze–thaw variation, which affected seasonal soil water dynamic and water budget and reformed seasonal runoff characteristics, along with vegetation cover changes, is considered the potential major factor in control of the hydrological processes in the permafrost region.     

Distribution and changes of active layer thickness (ALT) and soil temperature (TTOP) in the source area of the Yellow River using the GIPL model

Active layer thickness (ALT) is critical to the understanding of the surface energy balance, hydrological cycles, plant growth, and cold region engineering projects in permafrost regions. The temperature at the bottom of the active layer, a boundary layer between the equilibrium thermal state (in permafrost below) and transient thermal state (in the atmosphere and surface canopies above), is an important parameter to reflect the existence and thermal stability of permafrost. In this study, the Geophysical Institute Permafrost Model (GIPL) was used to model the spatial distribution of and changes in ALT and soil temperature in the Source Area of the Yellow River (SAYR), where continuous, discontinuous, and sporadic permafrost coexists with seasonally frozen ground. Monthly air temperatures downscaled from the CRU TS3.0 datasets, monthly snow depth derived from the passive microwave remote-sensing data SMMR and SSM/I, and vegetation patterns and soil properties at scale of 1:1000000 were used as input data after modified with GIS techniques. The model validation was carried out carefully with ALT in the SAYR has significantly increased from 1.8 m in 1980 to 2.4 m in 2006 at an average rate of 2.2 cm yr^?1. The mean annual temperature at the bottom of the active layer, or temperature at the top of permafrost (TTOP) rose substantially from ?1.1°C in 1980 to ?0.6°C in 2006 at an average rate of 0.018°C yr^?1. The increasing rate of the ALT and TTOP has accelerated since 2000. Regional warming and degradation of permafrost has also occurred, and the changes in the areal extent of regions with a sub-zero TTOP shrank from 2.4×10⁴ to 2.2×10⁴ km² at an average rate of 74 km² yr^?1. Changes of ALT and temperature have adversely affected the environmental stability in the SAYR.     

Characteristics of land surface heat and water exchange under different soil freeze/thaw conditions over the central Tibetan Plateau

Freezing and thawing processes at the soil surface play an important role in determining the nature of Tibetan land and atmosphere interactions. In this study, land surface water and heat exchanges under different freezing and thawing conditions over the central Tibetan Plateau were investigated using observations from the Coordinated Enhanced Observing Period/Asia‐Australia Monsoon Project on the Tibetan Plateau, and the Simultaneous Heat and Water Model. During the freezing and thawing stages, significant diurnal variation of soil temperature resulted in a diurnal cycle of unfrozen water content at the surface. Radiation and energy components and evapotranspiration averaged over four freeze/thaw stages also changed diurnally. On average, the surface albedo (0·68) during the completely frozen stage was sharply higher than those during the freezing, thawing, and completely thawed stages due to the snow cover. The Bowen ratios were 3·1 and 2·5 in the freezing and thawing stages, respectively, but the ratio was only 0·5 in the completely thawed stage. Latent heat flux displayed distinctly better correlation with unfrozen soil water content during the freezing and thawing stages than during the completely frozen and thawed stages. This implies that the diurnal cycle of unfrozen soil water, resulting from diurnal freeze/thaw cycles at the surface, has a significant impact on latent heat flux. A surface energy imbalance problem was encountered, and the possible sources of error were analysed. Copyright © 2011 John Wiley & Sons, Ltd.     

Interaction study of permafrost and highway along Qinghai-Xizang Highway

Qinghai-Xizang Highway crosses permafrost areas, sporadic permafrost areas and seasonal frozen soil areas. It is very important to study the engineering geology of frozen soil in construc-tion of cold regions engineering, and the choices of engineering st…     

A new soil freeze/thaw discriminant algorithm using AMSR‐E passive microwave imagery

The soil freeze–thaw controls the hydrological and carbon cycling and thus affects water and energy exchanges at land surface. This article reported a newly developed algorithm for distinguishing the freeze/thaw status of surface soil. The algorithm was based on information from Advanced Microwave Scanning Radiometer Enhanced (AMSR‐E) which records brightness temperature (Tb) in the afternoon and after midnight. The criteria and discriminant functions were obtained from both radiometer observations and model simulations. First of all, the microwave radiation from freeze–thaw soil was examined by carrying out experimental measurements at 18·7 and 36·5 GHz using a Truck‐mounted Multi‐frequency Microwave Radiometer (TMMR) in the Heihe River of China. The experimental results showed that the soil moisture is a key component that differentiates the microwave radiation behaviours during the freeze–thaw process, and the differences in soil temperature and emissivity between frozen and thawed soils were found to be the most important criteria. Secondly, a combined model was developed to consider the impacts of complex ground surface conditions on the discrimination. The model simulations quite followed the trend of in situ observations with an overall relation coefficient (R) of approximately 0·88. Finally, the ratio of Tb18·7H (horizontally polarized Tb at 18·7 GHz) to Tb36·5V was considered primarily as the quasi‐emissivity, which is more reasonable and explicit in measuring the microwave radiation changes in soil freezing and thawing than the spectral gradient. By combining Tb36·5V to indicate the soil temperature variety, a Fisher linear discrimination analysis was used to establish the discriminant functions. After being corrected by TMMR measurements, the new discriminant algorithm had an overall accuracy of 86% when validated by 4‐cm soil temperature. The multi‐year discriminant results also provided a good agreement with the classification map of frozen ground in China. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrological impacts of seismic lines in the wetland‐dominated zone of thawing,discontinuous permafrost,Northwest Territories,Canada

Intensive seismic exploration in the Northwest Territories began in the late 1960s. Since that time, the legacy of seismic surveys – i.e. straight lines cutting through boreal forest and tundra – has remained visible throughout northern Canada and Alaska. The removal of trees and compaction of the ground surface alter the thermophysical properties of the active (i.e. seasonally thawed) layer to such an extent that the underlying permafrost seriously degrades or even disappears completely. Such a transformation along linear corridors that cut indiscriminately across different terrain types with contrasting hydrological functions has potentially serious implications to the redistribution of water and energy within and among landscape units with feedbacks to permafrost thaw, land cover change and run‐off generation. This paper characterizes the flow and storage of water and energy along a seismic cut line in the high boreal zone of discontinuous permafrost in order to improve the understanding of these processes, their interactions and hydrological implications. As such, this paper lays a conceptual foundation for the development of numerical models needed to predict the hydrological and thermal impact of seismic lines in this sensitive region. We used ground‐penetrating radar and multi‐year ground temperatures and water levels along a seismic line to estimate the degree of permafrost degradation below it. The seismic line studied extends from a permafrost‐free wetland (flat bog), over a permafrost body (peat plateau) and into another permafrost‐free wetland (channel fen). It was found that once thaw had lowered the permafrost table below the ground surface elevation of the flat bog and channel fen, the seismic line forms a hydrological connection between them. It was also shown that during the permafrost thaw process, seismic lines develop a perennially thawed layer (talik) between the overlying active layer and underlying permafrost and that the talik conveys water as a conduit throughout the year. The implications of such drainage through seismic lines and networks on basin drainage in peatland‐dominated regions with discontinuous permafrost are also discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

Thermal regimes and degradation modes of permafrost along the Qinghai-Tibet Highway

Permafrost on the Qinghai-Tibet Plateau (QTP) is widespread, thin, and thermally unstable. Under a warming climate during the past few decades, it has been degrading extensively with generally rising ground temperatures, the deepening of the maximum summer thaw, and with lessening of the winter frost penetration. The permafrost has degraded downward, upward and laterally. Permafrost has thinned or, in some areas, has totally disappeared. The modes of permafrost degradation have great significance in geocryology, in cold regions engineering and in cold regions environmental management. Permafrost in the interior of the QTP is well represented along the Qing-hai-Tibet Highway (QTH), which crosses the Plateau through north to south and traverses 560 km of permafrost-impacted ground. Horizontally, the degradation of permafrost occurs more visibly in the sporadic permafrost zone in the vicinity of the lower limit of permafrost (LLP), along the margins of taliks, and around permafrost islands. Downward degradation develops when the maximum depth of seasonal thaw exceeds the maximum depth of seasonal frost, and it generally results in the formation of a layered talik disconnecting the permafrost from the seasonal frost layer. The downward degrada- tion is divided into four stages: 1) initial degradation, 2) accelerated degradation, 3) layered talik and 4) finally the conversion of permafrost to seasonally frozen ground (SFG). The upward degradation occurs when the geothermal gradient in permafrost drops to less than the geothermal gradients in the underlying thawed soil layers. Three types of permafrost temperature curves (stable, degrading, and phase-changing transitory permafrost) illustrate these modes. Although strong differentiations in local conditions and permafrost types exist, the various combinations of the three degradation modes will ultimately transform permafrost into SFG. Along the QTH, the downward degradation has been proceeding at annual rates of 6 to 25 cm, upward degradation at 12 to 30 cm, and lateral degradation in the sporadic permafrost zone at 62 to 94 cm during the last quarter century. These rates exceed the 4 cm per year for the past 20 years reported for the discontinuous permafrost zone in subarctic Alaska, the 3 to 7 cm per year reported in Mongolia, and that of the thaw-stable permafrost in subarctic Yakutia and Arctic Alaska.     

The Response of Shallow Groundwater Levels to Soil Freeze–Thaw Process on the Qinghai-Tibet Plateau

The Qinghai-Tibet plateau has the world's largest area of seasonally frozen ground. Here, shallow groundwater displays behavior that is distinct from that elsewhere in the world. In the present study, we explore the seasonal and interannual variation of the shallow groundwater levels from 2012 to 2016, and attempt to quantitatively evaluate the relative influences of individual driving factors on the shallow groundwater levels based on boosted regression trees. The results show that: (1) on a seasonal scale, the groundwater levels were characterized by a double peak and double valley relationship, while on an interannual scale the groundwater levels showed a slightly downwards trend from 2012 to 2016; and (2) during the frozen period, the seasonal variation of groundwater levels was determined by mean air temperature through its effect on the soil thaw–freeze process, accounting for 53.15% of total variation. Meanwhile, ET₀ and rainfall exerted little impact on the seasonal variation of groundwater levels, which might be attributed to the aquitard of frozen soil that impedes the exchange between surface water and groundwater. Moreover, there was a lag between groundwater levels and soil freezing–thawing. During the non-frozen period, the mean air temperature was again the most important factor impacting the variation of groundwater levels, through its effect on ET₀, and accounted for 40.75% of total variation, while rainfall had little effect on groundwater levels when rainfall intensity was less than 12 mm/day. These results will benefit predictions of future trends in groundwater levels within the context of global warming.     

Impact of permafrost change on the Qinghai-Tibet Railroad engineering

Permafrost along the Qinghai-Tibet Railroad produces the great change under the influence of climate change, such as the decreasing of permafrost table, the rising of permafrost temperatures, etc. Climate effect on permafrost is the long-term process. Engineering action makes rapidly permafrost the large extent change. On the basis of analyzing the permafrost change under the climate change and engineering action, the thermal regime and spatial distribution of permafrost are predicted for air temperature rising 1℃ and 2℃ after 50 years in this paper. The results show that climate change results in the larger change for the thermal regime and spatial distribution of permafrost. Permafrost change will produce the great effect on the Qinghai-Tibet Railroad engineering, not only resulting in the decreasing of permafrost table beneath the roadbed, but also resulting in thawing settlement due to the thawing of ground ice near permafrost table. The idea of cooling roadbed and active protecting permafrost for the Qinghai-Tibet Railroad engineering could adjust and control the permafrost thermal state, some better methods are provided to ensure the engineering stability in the areas of warm permafrost and high ice content.     

Synergistic effect of vegetation and air temperature changes on soil water content in alpine frost meadow soil in the permafrost region of Qinghai‐Tibet

Seasonal changes over 2 years (2004–2006) in soil moisture content (θ_v) of frozen alpine frost meadow soils of the Qinghai‐Tibet plateau permafrost region under three different levels of vegetation cover were investigated. Vegetation cover and air temperature changes had significant effects (synergistic effect) on θ_v and its distribution in the soil profile. During periods of soil freezing or thawing, the less the vegetation cover, the quicker the temperature drop or rise of soil water, and the shorter the duration of the soil water freeze–thaw response in the active soil layer. Under 30% and 65% vegetation cover the amplitude of variation in θ_v during the freezing period was 20–26% greater than that under 93% cover, while during the thawing period, it was 1·5‐ to 40·5‐fold greater. The freezing temperature of the surface soil layer, ^fT_s, was 1·6 °C lower under 30% vegetation cover than under 93% vegetation cover. Changes in vegetation cover of the alpine frost meadow affected θ_v and its distribution, as well as the relationship between θ_v and soil temperature (T_s). As vegetation cover decreased, soil water circulation in the active layer increased, and the response to temperature of the water distribution across the soil profile was heightened. The quantity of transitional soil phase water at different depths significantly increased as vegetation cover decreased. The influence of vegetation cover and soil temperature distribution led to a relatively dry soil layer in the middle of the profile (0·70–0·80 m) under high vegetation cover. Alpine meadow θ_v and its pattern of distribution in the permafrost region were the result of the synergistic effect of air temperature and vegetation cover. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact of heat advection on the thermal regime of roads built on permafrost

In northern regions, transportation infrastructure can experience severe structural damages due to permafrost degradation. Water infiltration and subsurface water flow under an embankment affect the energy balance of roadways and underlying permafrost. However, the quantification of the processes controlling these changes and a detailed investigation of their thermal impacts remain largely unknown due to a lack of available long-term embankment temperature data in permafrost regions. Here, we report observations of heat advection linked to surface water infiltration and subsurface flow based on a 9-year (from 2009 to 2017) thermal monitoring at an experimental road test site built on ice-rich permafrost conditions in southwestern Yukon, Canada. Our results show that snowmelt water infiltration in the spring rapidly increases temperature in the upper portion of the embankment. The earlier disappearance of snow deposited at the embankment slope increases the thawing period and the temperature gradient in the embankment compared with the natural ground. Infiltrated summer rainfall water lowered the near-surface temperatures and subsequently warmed embankment fill materials down to 3.6-m depth. Heat advection caused by the flow of subsurface water produced warming rates at depth in the embankment subgrade up to two orders of magnitude faster than by atmospheric warming (heat conduction). Subsurface water flow promoted permafrost thawing under the road embankment and led to an increase in active layer thickness. We conclude that the thermal stability of roadways along the Alaska Highway corridor is not maintainable in situations where water is flowing under the infrastructure unless mitigation techniques are used. Severe structural damages to the highway embankment are expected to occur in the next decade.     

Mercury and methylmercury biogeochemistry in a thawing permafrost wetland complex,Northwest Territories,Canada

In arctic and sub‐arctic environments, mercury (Hg), more specifically toxic methylmercury (MeHg), is of growing concern to local communities because of its accumulation in fish. In these regions, there is particular interest in the potential mobilization of atmospherically deposited Hg sequestered in permafrost that is thawing at unprecedented rates. Permafrost thaw and the resulting ground surface subsidence transforms forested peat plateaus into treeless and permafrost‐free thermokarst wetlands where inorganic Hg released from the thawed permafrost and draining from the surrounding peat plateaus may be transformed to MeHg. This study begins to characterize the spatial distribution of MeHg in a peat plateau–thermokarst wetland complex, a feature that prevails throughout the wetland‐dominated southern margin of thawing discontinuous permafrost in Canada's Northwest Territories. We measured pore water total Hg, MeHg, dissolved organic matter characteristics and general water chemistry parameters to evaluate the role of permafrost thaw on the pattern of water chemistry. A gradient in vegetation composition, water chemistry and dissolved organic matter characteristics followed a toposequence from the ombrotrophic bogs near the crest of the complex to poor fens at its downslope margins. We found that pore waters in poor fens contained elevated levels of MeHg, and the water draining from these features had dissolved MeHg concentrations 4.5 to 14.5 times higher than the water draining from the bogs. It was determined through analysis of historical aerial images that the poor fens in the toposequence had formed relatively recently (early 1970s) as a result of permafrost thaw. Differences between the fens and bogs are likely to be a result of their differences in groundwater function, and this suggests that permafrost thaw in this landscape can result in hotspots for Hg methylation that are hydrologically connected to downstream ecosystems. Copyright © 2016 John Wiley & Sons, Ltd.     

青藏高原多年冻土区天然气水合物形成条件模拟研究

基于野外气体地球化学调查研究,以及前人有关冻土表层温度、冻土层内地温梯度、冻土层下地温梯度等的资料,对青藏高原多年冻土区天然气水合物的形成条件开展了模拟研究. 结果显示：研究区冻土条件能够满足天然气水合物形成的基本要求;气体组成、冻土特征（如冻土厚度或冻土表层温度、冻土层内地温梯度、冻土层下地温梯度等）是影响研究区天然气水合物稳定带厚度的最重要因素,其在不同点位上的差异性可能导致天然气水合物分布的不均匀性的主要原因;研究区最可能的天然气水合物为甲烷与重烃（乙烷和丙烷）的混合气体型天然气水合物;在天然气水合物分布的区域,其产出的上临界点深度在几十至一百多米间,下临界点深度在几百至近一千米间,厚度可达到几百米. 与Canadian Mallik三角洲多年冻土区相比,青藏高原多年冻土区除了冻土厚度小些外,其他条件,如冻土层内地温梯度、冻土层下地温梯度、气体组成等条件较为相近,具有一定的可比性,预示着良好的天然气水合物潜力.     

Effects of sand-mulch thickness on soil evaporation during the freeze–thaw period

Control of evaporation from seasonally frozen soil is an important method for alleviating water shortages in arid and semi-arid areas. To investigate the inhibition of soil evaporation by sand and the major factors that influence soil evaporation, a series of field experiments with five sand-mulch thicknesses (0 cm, bare soil [BS], 1 cm [T1], 2 cm [T2], 3 cm [T3] and 4 cm [T4], with an average diameter of 1 mm) were conducted during the freeze–thaw period in Northern China. Soil evaporation characteristics in the three freeze–thaw stages were revealed and the major factors influencing soil evaporation were analysed using grey correlation analysis. The results showed that the cumulative soil evaporation decreased with increasing sand-mulch thickness during the freeze–thaw period, and only small differences in soil evaporation were observed between the T3 and T4 treatments. The reduction in soil evaporation under different sand-mulch thicknesses was 19.2–62.6% in the unstable freezing stage (P1), 2.0–28.3% in the stable freezing stage (P2) and 4.8–20.4% in the thawing stage (P3). In P1, solar radiation was a major factor influencing soil evaporation in all treatments and vapour pressure was a major factor in the sand-mulch treatments, and the influence of relative humidity on soil evaporation decreased in the T4 treatment. During the coldest P2, solar radiation was lowest so that relative humidity and wind speed became the more dominant influence factors on soil evaporation in all treatments, and surface soil water content was a major factor in the sand-mulch treatments. In P3, average air temperature and solar radiation were major factor influencing soil evaporation in all treatments and vapour pressure was a major factor in the BS and T1 treatments, whereas water surface evaporation was the major factor in the T2, T3 and T4 treatments. The results suggest that the addition of sand mulch in agricultural fields may be a beneficial practice to reduce water stress in arid and semi-arid areas.