Variability Of Snow Depth Research Articles

Estimating high-resolution snow depth over the North Hemisphere mountains utilizing active microwave backscatter and machine learning

While ground meteorological stations provide accurate snow depth data, their limited spatial coverage results in observational gaps. Satellites offer long-term, large-scale observations, addressing these gaps. Existing snow depth retrieval algorithms mainly use passive microwave remote sensing data with a 25 km resolution, insufficient for capturing snow depth variability in mountainous areas. This paper introduces active microwave backscatter data and machine learning techniques for high-resolution snow depth estimation. We conducted a preliminary exploration of the relationship between Sentinel-1 backscatter coefficient σ0 and snow depth. Due to factors such as vegetation coverage and underlying soil properties, the relationship between σ0 and snow depth is complex and nonlinear. Consequently, six machine learning models were trained to learn this relationship using σ0 and auxiliary data as input features, with in-situ snow depth serving as the target variable. After extensive validation, the Extreme Random Trees (ERT) model was selected for its high accuracy and stability. Using the ERT model, we generated 500 m-resolution snow depth data for northern hemisphere mountains, then analyzed temporal snow depth variations and altitudinal stratification.

Summer snow on Arctic sea ice modulated by the Arctic Oscillation

Since the 1970s, Arctic sea ice has undergone unprecedented change, becoming thinner, less extensive and less resilient to summer melt. Snow’s high albedo greatly reduces solar absorption in sea ice and the upper ocean, which mitigates sea–ice melt and ocean warming. However, the drivers of summertime snow depth variability are unknown. The Arctic Oscillation is a mode of natural climate variability, influencing Arctic snowfall and air temperatures. Thus, it may affect summertime snow conditions on Arctic sea ice. Here we examine the role of the Arctic Oscillation in summer snow depth variability on Arctic sea ice in 1980–2020 using atmospheric reanalysis, snow modelling and satellite data. The positive phase leads to greater snow accumulation, ranging up to ~4.5 cm near the North Pole, and higher surface albedo in summer. There are more intense, frequent Arctic cyclones, cooler temperatures aloft and greater snowfall relative to negative and neutral phases; these conditions facilitate a more persistent summer snow cover, which may lessen sea-ice melt and ocean warming. The Arctic Oscillation influence on summertime snow weakens after 2007, which suggests that future warming and Arctic sea-ice loss might modify the relationship between the Arctic Oscillation and snow on Arctic sea ice.

A machine learning approach for estimating snow depth across the European Alps from Sentinel-1 imagery

Seasonal snow plays a crucial role in society and understanding trends in snow depth and mass is essential for making informed decisions about water resources and adaptation to climate change. However, quantifying snow depth in remote, mountainous areas with complex topography remains a significant challenge. The increasing availability of high-resolution synthetic aperture radar (SAR) observations from active microwave satellites has prompted opportunistic use of the data to retrieve snow depth remotely across large spatial and frequent temporal scales and at a high spatial resolution. Nevertheless, these novel SAR-based snow depth retrieval methods face their own set of limitations, including challenges for shallow snowpacks, high vegetation cover, and wet snow conditions. In response, here we introduce a machine learning approach to enhance SAR-based snow depth estimation over the European Alps. By integrating Sentinel-1 SAR imagery, optical snow cover observations, and topographic, forest cover and snow class information, our machine learning retrieval method more accurately estimates snow depth at independent in-situ measurement sites than current methods. Further, our method provides estimates at 100 m horizontal resolution and is capable of better capturing local-scale topography-driven snow depth variability. Through detailed feature importance analysis, we identify optimal conditions for SAR data utilization, thereby providing insight into future use of C-band SAR for snow depth retrieval.

Spatio‐temporal analysis of snow depth and snow water equivalent in a mountainous catchment: Insights from in‐situ observations and statistical modelling

AbstractThis research, conducted in the mountainous catchment near Abant Lake in the Western Black Sea region of Türkiye, aimed to investigate the spatiotemporal variations of snow depth (SD) and snow water equivalent (SWE) throughout the snow season from December 2019 to March 2020, encompassing both accumulation and melting periods. In total, 14 snow surveys were conducted, covering 58 permanent snow measurement points (PSMP) marked with snow poles. The classification and regression tree (CART) method was employed to statistically analyse their relationships with eight variables: snow period, forest canopy, aspect, slope, elevation, slope position, plan and profile curvature. The root mean square error (RMSE) for SD and SWE was determined to be 0.15 m and 46 mm, respectively. The study findings revealed that mean SD and SWE values were higher in forest gaps compared with under‐forest and open areas. Although the snow cover disappeared earliest in under‐forest areas, the melting rate was observed to be 43% and 17% slower compared with forest gaps and open areas, respectively. Wind redistribution resulted in minimum snow accumulation on western aspects, upper slope positions and ridges, while maximum accumulation was observed on southern aspects, valleys and lower slope positions. Higher elevations (>1580 meters) experienced faster snow melting rates, leading to earlier disappearance of snow cover. PSMPs located on slopes with lower degrees (<15°) exhibited lesser accumulation and earlier snow disappearance. The CART model identified the snow period as the most significant factor in predicting SD and SWE, based on variations in snowfall and air temperature. Other significant variables included forest canopy, aspect and elevation. The study suggests that the CART method is well‐suited for modelling complex snow dynamics, providing valuable insights into spatiotemporal variations in SD and SWE in mountainous regions.

The Varied Role of Atmospheric Rivers in Arctic Snow Depth Variations

AbstractThe state and fate of snow on sea ice are crucial in the mass and energy balance of sea ice. The function of atmospheric rivers (ARs) on snow depth over sea ice has not been measured thus far, limiting the understanding of the mechanism of snow depth changes. Here, the effect of ARs on snow depth changes was explored. We found that increased AR frequency is responsible for winter‐autumn snow accumulation and spring‐summer snow melting. The 2 m air temperature (T2m), rainfall, snowfall, mean net longwave radiation (NLR), mean net shortwave radiation (NSR) and cloud radiative effect (CRE) during ARs explain the changes in snow depth triggered by AR occurrence. This work helps us understand how ARs affect snow depth changes through related physical processes, promotes an understanding of climate systems and provides a theoretical basis for snow treatment in sea ice models.

Implementing robust outlier detection to enhance estimation accuracy of GNSS-IR based seasonal snow depth retrievals

ABSTRACT Monitoring snow depth variations aligned with the seasonal cycle is crucial for studies on climate change and its impacts due to global warming. Especially in terms of snow hydrology, the determination and continuous monitoring of snow depth to determine snow water equivalent is a priority in climate, water science, drought, flood, and inundation prediction studies, particularly in the water cycle. Recently, the so-called Global Navigation Satellite Systems – Interferometric Reflectometry (GNSS-IR) method enables to extract environmental radiometric and geometric characteristics of surface where the signals transmitted from the satellites reflect. The reflected signal follows an extra path than the direct one and causes a major error source for accurate point positioning but a useful tool for sensing the environment. However, outliers remain a challenge in long-term GNSS-IR estimations. In this study, we proposed a median-based robust outlier detection (ROD) approach for identifying outliers in long-term GNSS-IR snow depth estimations. To validate our approach, we analysed 5-year GNSS L1 SNR and L2 SNR data from 1 January 2015 to 31 December 2019 provided by AB33 and AB39 GNSS stations in Alaska, U.S.A. using the GNSS-IR method. We validated GNSS-IR estimations using snow depth measurements from the Coldfoot and Fort Yukon stations in the SNOTEL network. Applying ROD to long-term snow depth estimates increased the highest correlation from 91.58% to 94.93%, and reduced the lowest RMSE from 8.6 cm to 6.7 cm. In addition, the improvement rates calculated to assess the contribution of ROD to the results showed improvements of up to 6.9% in correlation and 26.1% in RMSE. Overall, the results demonstrate that ROD can be effectively used to detect outliers in long-term GNSS-IR snow depth time series.

The MOSAiC Distributed Network: Observing the coupled Arctic system with multidisciplinary, coordinated platforms

Central Arctic properties and processes are important to the regional and global coupled climate system. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Distributed Network (DN) of autonomous ice-tethered systems aimed to bridge gaps in our understanding of temporal and spatial scales, in particular with respect to the resolution of Earth system models. By characterizing variability around local measurements made at a Central Observatory, the DN covers both the coupled system interactions involving the ocean-ice-atmosphere interfaces as well as three-dimensional processes in the ocean, sea ice, and atmosphere. The more than 200 autonomous instruments (“buoys”) were of varying complexity and set up at different sites mostly within 50 km of the Central Observatory. During an exemplary midwinter month, the DN observations captured the spatial variability of atmospheric processes on sub-monthly time scales, but less so for monthly means. They show significant variability in snow depth and ice thickness, and provide a temporally and spatially resolved characterization of ice motion and deformation, showing coherency at the DN scale but less at smaller spatial scales. Ocean data show the background gradient across the DN as well as spatially dependent time variability due to local mixed layer sub-mesoscale and mesoscale processes, influenced by a variable ice cover. The second case (May–June 2020) illustrates the utility of the DN during the absence of manually obtained data by providing continuity of physical and biological observations during this key transitional period. We show examples of synergies between the extensive MOSAiC remote sensing observations and numerical modeling, such as estimating the skill of ice drift forecasts and evaluating coupled system modeling. The MOSAiC DN has been proven to enable analysis of local to mesoscale processes in the coupled atmosphere-ice-ocean system and has the potential to improve model parameterizations of important, unresolved processes in the future.

Multi-Source Dataset Assessment and Variation Characteristics of Snow Depth in Eurasia from 1980 to 2018

Snow is an indicator of climate change. Its variation can affect surface energy, water balance, and atmospheric circulation, providing important feedback on climate change. There is a lack of assessment of the spatial characteristics of multi-source snow data in Eurasia, and these data exhibit high spatial variability and other differences. Therefore, using data obtained from the Global Historical Climatology Network Daily (GHCND) from 1980 to 2018, snow depth information from ERA5, MERRA2, and GlobSnow is assessed in this study. The spatiotemporal variation characteristics and the primary spatial modes of seasonal variations in snow depth are analyzed. The results show that the snow depth, according to GlobSnow data, is closer to that of the measured site data, while the ERA5_Land and MERRA2 data are overestimated. The annual variations in snow depth are consistent with seasonal variations in winter and spring, with an increasing trend in the mountains of Central Asia and Siberia and a decreasing trend in most of the rest of Eurasia. The dominant patterns of snow depth in late autumn, winter, and spring are all north–south dipole patterns, and there is overall consistency in summer.

Using Icepack to reproduce ice mass balance buoy observations in landfast ice: improvements from the mushy-layer thermodynamics

Abstract. Icepack (v1.1.0) – the column thermodynamics model of the Community Ice CodE (CICE) version 6 – is used to assess how changing the thermodynamics from the Bitz and Lipscomb (1999) physics (hereafter BL99) to the mushy-layer physics impacts the model performance in reproducing in situ landfast ice observations from two ice mass balance (IMB) buoys co-deployed in the landfast ice close to Nain (Labrador) in February 2017. To this end, a new automated surface retrieval algorithm is used to determine the in situ ice thickness, snow depth, basal ice congelation and snow-ice formation from the measured vertical temperature profiles. Icepack simulations are run to reproduce these observations using each thermodynamics scheme, with a particular interest in how the different physics influence the representation of snow-ice formation and ice congelation. Results show that the BL99 parameterization represents well the ice congelation but underrepresents the snow-ice contribution to the ice mass balance. In particular, defining snow-ice formation based on the hydrostatic balance alone does not reproduce the negative freeboards observed for several days in the IMB observations, resulting in an earlier snow-flooding onset, a positive ice thickness bias and reduced snow depth variations. We find that the mushy-layer thermodynamics with default parameters significantly degrades the model performance, overestimating both the congelation growth and snow-ice formation. The simulated thermodynamics response to flooding, however, better represents the observations, and the best results are obtained when allowing for negative freeboards in the mushy-layer physics. We find that the mushy-layer thermodynamics produces a larger variability in congelation rates at the ice bottom interface, alternating between periods of exceedingly fast growth and periods of unrealistic basal melt. This pattern is related to persistent brine dilution in the lowest ice layer by the congelation and brine drainage parameterizations. We also show that the mushy-layer congelation parameterization produces significant frazil formation, which is not expected in a landfast ice context. This behavior is attributed to the congelation parameterization not fully accounting for the conductive heat flux imbalance at the ice–ocean boundary. We propose a modification of the mushy-layer congelation scheme that largely reduces the frazil formation and allows for better tuning of the congelation rates to match the observations. Our results demonstrate that the mushy-layer physics and its parameters can be tuned to closely match the in situ observations, although more observations are needed to better constrain them.

A new interpolation method to resolve under-sampling of UAV-lidar snow depth observations in coniferous forests

Obtaining accurate snow depth estimates under dense canopies using airborne lidar (light detection and ranging) techniques is challenging due to the under-sampling of ground and snow surfaces. Existing interpolation techniques do not adequately address this problem and they often result in an overestimation of under-canopy snow depths. To address this issue, we introduce and evaluate a new interpolation method that incorporates intra-canopy snow depth variability to provide more accurate estimations at unsampled locations. Four interpolation methods were tested, considering systematic trends (landscape trend, canopy vs. gap trend, and intra-canopy trend) along with spatial interpolation of the residuals. Our results show that spatial interpolation methods without consideration of trends are sufficient to capture and reconstruct the small-scale variability of snow depths below a separation distance of 1 m between sampled and unsampled locations, (i.e., ground surface point density > 1 pt. m−2). However, beyond a separation distance of 2.5–3 m (point density < 0.33–0.40 pt. m−2), spatial interpolation based on proximity alone becomes unreliable because point separation becomes larger than the snow depth spatial correlation scale. Within these limiting distances, the method that incorporates trends along with spatial interpolation techniques can resolve the small-scale variability and thereby reduce the likely overestimation of snow depths under the canopy.

Long‐Term Variability and Trends in Snow Depth and Cover Days Throughout Iranian Mountain Ranges

AbstractIn Iran, the mountain snow cover generally feeds major rivers and thereby largely provides water resources required for improving human lives and protecting nature. Hence, understanding historical variability and trends in mountainous snowpack water resources in Iran in response to global warming and climate change can play a critical role in the sustainable development of this country. Accordingly, this study investigated long‐term (1982–2018) snowpack climatology at 13 hydrometeorological measurement stations scattered throughout the Iranian mountain ranges, with a focus on Elburz, Azerbaijan, Zagros, and Khorasan mountainous regions. The non‐parametric Mann‐Kendall test was used to detect statistically significant (p &lt; 0.05) trends, the Pettitt test to identify possible abrupt shift years, the Pearson's correlation coefficient to measure relationships among different time series, and the partial correlation to determine the most important climate factor influencing snowpack dynamics The annual snow depth (maximum snow depth) significantly declined throughout Iranian mountain ranges during 1982–2018, with an average rate of 1.0 (3.4) cm decade−1. The annual snow cover days (SCDs) also showed significant decreasing trends, ranging from 3 to 15 days decade−1 during 1982–2018, in 69% of the stations studied. Such considerable reductions in snow depth and cover days were mainly related to the compound effects of substantial increases in temperature, sunshine, and wind speed as well as decreases in precipitation and cloudiness during the SCDs across the Iranian mountain ranges. However, precipitation was the most influential climate factor controlling snow resources throughout both the Elburz and Zagros mountains in Iran.

Measuring the spatiotemporal variability in snow depth in subarctic environments using UASs – Part 2: Snow processes and snow–canopy interactions

Abstract. Detailed information on seasonal snow cover and depth is essential to the understanding of snow processes, to operational forecasting, and as input for hydrological models. Recent advances in uncrewed or unmanned aircraft systems (UASs) and structure from motion (SfM) techniques have enabled low-cost monitoring of spatial snow depth distribution in resolutions of up to a few centimeters. Here, we study the spatiotemporal variability in snow depth and interactions between snow and vegetation in different subarctic landscapes consisting of a mosaic of conifer forest, mixed forest, transitional woodland/shrub, and peatland areas. To determine the spatiotemporal variability in snow depth, we used high-resolution (50 cm) snow depth maps generated from repeated UAS–SfM surveys in the winter of 2018/2019 and a snow-free bare-ground survey after snowmelt. Due to poor subcanopy penetration with the UAS–SfM method, tree masks were utilized to remove canopy areas and the area (36 cm) immediately next to the canopy before analysis. Snow depth maps were compared to the in situ snow course and a single-point continuous ultrasonic snow depth measurement. Based on the results, the difference between the UAS–SfM survey median snow depth and single-point measurement increased for all land cover types during the snow season, from +5 cm at the beginning of the accumulation to −16 cm in coniferous forests and −32 cm in peatland during the melt period. This highlights the poor representation of point measurements in selected locations even on the subcatchment scale. The high-resolution snow depth maps agreed well with the snow course measurement, but the spatial extent and resolution of maps were substantially higher. The snow depth range (5th–95th percentiles) within different land cover types increased from 17 to 42 cm in peatlands and from 33 to 49 cm in the coniferous forest from the beginning of the snow accumulation to the melt period. Both the median snow depth and its range were found to increase with canopy density; this increase was greatest in the conifer forest area, followed by mixed forest, transitional woodland/shrub, and open peatlands. Using the high-spatial-resolution data, we found a systematic increase (2–20 cm) and then a decline in snow depth near the canopy with increasing distance (from 1 to 2.5 m) of the peak value through the snow season. This study highlights the applicability of the UAS–SfM in high-resolution monitoring of snow depth in multiple land cover types and snow–vegetation interactions in subarctic and remote areas where field data are not available or where the available data are collected using classic point measurements or snow courses.

Measuring the spatiotemporal variability in snow depth in subarctic environments using UASs – Part 1: Measurements, processing, and accuracy assessment

Abstract. Snow conditions in the Northern Hemisphere are rapidly changing, and information on snow depth is critical for decision-making and other societal needs. Uncrewed or unmanned aircraft systems (UASs) can offer data resolutions of a few centimeters at a catchment-scale and thus provide a low-cost solution to bridge the gap between sparse manual probing and low-resolution satellite data. In this study, we present a series of snow depth measurements using different UAS platforms throughout the winter in the Finnish subarctic site Pallas, which has a heterogeneous landscape. We discuss the different platforms, the methods utilized, difficulties working in the harsh northern environment, and the UAS snow depth results compared to in situ measurements. Generally, all UASs produced spatially representative estimates of snow depth in open areas after reliable georeferencing by using the structure from motion (SfM) photogrammetry technique. However, significant differences were observed in the accuracies produced by the different UASs compared to manual snow depth measurements, with overall root mean square errors (RMSEs) varying between 13.0 and 25.2 cm, depending on the UAS. Additionally, a reduction in accuracy was observed when moving from an open mire area to forest-covered areas. We demonstrate the potential of low-cost UASs to efficiently map snow surface conditions, and we give some recommendations on UAS platform selection and operation in a harsh subarctic environment with variable canopy cover.

Snow cover duration in northern Finland and the influence of key variables through a conceptual framework based on observed variations in snow depth

Seasonal snow cover duration is the net result from many processes acting on snow fallen on the Earth's surface. Several of these processes feed back into the atmosphere-cryosphere system causing non-linear interactions. The timing of snow retreat is of essential importance, but the duration of snow cover has large spatiotemporal variabilities. However, from a large data set of observed snow depth changes in northern Finland, systematic similar evolutions are identified that allow for a considerable simplification and reduction of the complexity in snow depth changes. Here, a novel conceptual framework is designed based on dividing the season into two main periods (dark and bright period, based on solar irradiance), for which snow depth decrease is parameterized based on three variables, average temperature, incoming shortwave radiation, and light-absorbing particles (LAP) in the snow. The processes are simplified into two linear relations, and a new formulation for concentration enhancement of LAP, which is dependent on snow depth decrease, is given. The results show that the seasonal snow cover duration is shifted by about one day for every 10 mm snow water equivalent of precipitation. This effect is comparable in scale to that of doubling of the amount of LAP concentration in snow. We also found that the combined shift in snow cover duration from interannual variability in ambient temperature and shortwave radiation (warm and bright vs. cold and dark season) is large enough to explain the variability of a couple of weeks for a given precipitation amount in Northern Finland.

Assessment of L-band InSAR snow estimation techniques over a shallow, heterogeneous prairie snowpack

Snow water equivalent (SWE) is a critical input for weather, climate, and water resource management models at local to global scales. Despite its importance, global SWE measurements that are accurate, consistent, and at sufficiently high spatiotemporal resolutions are not currently available. L-band interferometric synthetic aperture radar (InSAR) techniques have been used to measure SWE at local to regional scales, and two upcoming L-band SAR satellite missions have renewed interest in these techniques to provide regular SWE measurements at the global scale. However, previous research demonstrating the capabilities of L-band InSAR-SWE measurement has been limited to mountain or tundra snowpack regimes. Here we examine the feasibility of applying the same techniques over a prairie snowpack, which are typically characterized by shallow snow depths (mean snow depth of 0.22 m in this study), exposed agricultural vegetation, and high spatial variability over short distances. Our study area in central Montana, USA (47.060, -109.951) was a validation site for NASA SnowEx 2021, as part of the UAVSAR snow timeseries. Airborne L-band SAR imagery was acquired by the UAVSAR platform while concurrent snow measurements were collected using uncrewed aerial vehicle (UAV)-based LiDAR, UAV-based photogrammetry, and ground-based manual techniques. This validation dataset enables an investigation of the effects of sub-pixel snow cover heterogeneity and exposed agricultural vegetation stubble on SAR data and the resulting SWE estimations. Results based on repeated application of the Kolmogorov–Smirnov test show that UAVSAR VV phase change is sensitive to differences in snow cover but relatively unaffected by differences in agricultural stubble height. However, we did not find similarly definitive results when we used the same phase change data to estimate SWE. Although broad spatial patterns were similar in both LiDAR-derived and InSAR-derived SWE estimates, considerable differences in the two estimates were apparent in areas with large sub-pixel snow depth variability. Our results indicate that additional work is necessary to derive accurate SWE estimates in prairie environments. Regular measurements from L-band SAR satellites will provide an excellent opportunity to refine InSAR-based snow estimation techniques over shallow, heterogeneous snowpacks.

Consistent imprints of elevation, soil temperature and moisture on plant and arthropod communities across two subarctic landscapes

Abstract Factors shaping arthropod and plant community structure at fine spatial scales are poorly understood. This includes microclimate, which likely plays a large role in shaping local community patterns, especially in heterogeneous landscapes characterised by high microclimatic variability in space and in time. We explored differences in local microclimatic conditions and regional species pools in two subarctic regions: Kilpisjärvi in north‐west Finland and Varanger in north‐east Norway. We then investigated the relationship between fine‐scale climatic variation and local community characteristics (species richness and abundance) among plants and arthropods, differentiating the latter into two groups: flying and ground‐dwelling arthropods collected by Malaise and pitfall traps, respectively. Arthropod taxa were identified through DNA metabarcoding. Finally, we examined if plant richness can be used to predict patterns in arthropod communities. Variation in soil temperature, moisture and snow depth proved similar between regions, despite differences in absolute elevation. For each group of organisms, we found that about half of the species were shared between Kilpisjärvi and Varanger, with a quarter unique to each region. Plants and arthropods responded largely to the same drivers. The richness and abundance of both groups decreased as elevation increased and were positively correlated with higher soil moisture and temperature values. Plant species richness was a poor predictor of local arthropod richness, in particular for ground‐dwelling arthropods. Our results reveal how microclimatic variation within each region carves pronounced, yet consistent patterns in local community richness and abundance out of a joint species pool.

Characteristics of snow depth and period variations in snowy areas in Japan for 60 years from 1962 to 2021 winter season using annual average snow depth and annual maximum snow depth

Characteristics of snow depth and period variations in snowy areas in Japan for 60 years from 1962 to 2021 winter season using annual average snow depth and annual maximum snow depth

Quality of Snow Cover Characteristics Derived from ERA 5-Land Reanalysis for the Territory of Perm Krai

Received March 31, 2023; revised June 6, 2023; accepted June 27, 2023Agreement between values of the mean monthly snow depth provided by the ERA 5-Land reanalysis and similar characteristics of snow cover calculated on the basis of the extended hydrometeorological monitoring performed in the Perm Region for 1990–2020 is analyzed. It was found that ERA 5-Land in 73% of cases reproduces the presence/absence of snow during the onset period, and in 53% – during the period of snow loss. The conclusions made in the authors' previous studies based on more limited material were generally confirmed. It is shown that the reanalysis values of the snow depth are overestimated in relation to instrumental measurements for most of the hydrometeorological stations of the Perm Region. In space, the magnitude of the reanalysis error increases from the southwest to the northeast of the region, with the exception of its central part, where the influence of the Kama water reservoir is perceptible. But the interannual variability of the average snow depth in the Perm Region was reproduced by the ERA 5-Land reanalysis adequately. For 30 years, the magnitude of the reanalysis error decreased as it was compared with 61% observation points. The analysis of seasonal variability showed that in ERA 5-Land time of the maximum snow depth was shifted to earlier onset. The complete coincidence of the seasonal course was recorded only in 5% of hydrometeorological monitoring sites. The value of the average monthly discrepancies between the data of the reanalysis and the information of the posts as a whole exceeds the similar characteristic for the stations, which is especially evident during the period of active snowmelt.

Topographic and vegetation controls of the spatial distribution of snow depth in agro-forested environments by UAV lidar

Abstract. Accurate knowledge of snow depth distributions in forested regions is crucial for applications in hydrology and ecology. In such a context, understanding and assessing the effect of vegetation and topographic conditions on snow depth variability is required. In this study, the spatial distribution of snow depth in two agro-forested sites and one coniferous site in eastern Canada was analyzed for topographic and vegetation effects on snow accumulation. Spatially distributed snow depths were derived by unmanned aerial vehicle light detection and ranging (UAV lidar) surveys conducted in 2019 and 2020. Distinct patterns of snow accumulation and erosion in open areas (fields) versus adjacent forested areas were observed in lidar-derived snow depth maps at all sites. Omnidirectional semi-variogram analysis of snow depths showed the existence of a scale break distance of less than 10 m in the forested area at all three sites, whereas open areas showed comparatively larger scale break distances (i.e., 11–14 m). The effect of vegetation and topographic variables on the spatial variability in snow depths at each site was investigated with random forest models. Results show that the underlying topography and the wind redistribution of snow along forest edges govern the snow depth variability at agro-forested sites, while forest structure variability dominates snow depth variability in the coniferous environment. These results highlight the importance of including and better representing these processes in physically based models for accurate estimates of snowpack dynamics.

Study of snow cover/depth evolution characteristics in Tianshan region of China based on geographical partition

Based on the digital elevation data, snow depth and snow cover remote sensing data, this paper divides six snow evolution areas and geographical partitions, extracts the geographical partitions of each evolution area and obtains the geographical characteristics of the evolution area for analysis. The results show that: (1) From 2003 to 2017, the average snow area decreased at a rate of − 0.004, and the average snow depth increased at a rate of 0.03. (2) The snow in the middle altitude hill with shady gentle slope area is the most obvious in the seasonal evolution, and the percentage of this region in the seasonal snow evolution area is 5.46%, the snow depth in the middle altitude hill with sunny and gentle slopes area increased and decreased significantly in the past 15 years, and the percentage of this region in the SD significant changes evolution area was 6.32%. The snow in the low relief middle altitude mountain with shady and moderate slope area not only shows obvious seasonal evolution, but also increases and decreases significantly in snow depth. And the percentage of this region in the seasonal snow significant evolution area is 5.82%. (3) The geographical partitions with the largest area in all evolution areas is the middle altitude hill with sunny and gentle slopes area (4.75%). (4) The geographical partition with the largest variation of snow depth in Tianshan region is the low relief middle altitude mountain with shady and moderate slope area (12.02 cm). (5) The snow accumulation and melting are obvious in the range of 1000–3500 m above altitude, different geomorphology types lead to obvious differences in snow characteristics. The snow melting is most obvious in the gentle slope area of the low topographic relief geomorphology types, and the snow accumulation is most obvious in the steep slope area of the middle relief geomorphology types.