Explore the vast range of titles available.

Assessing past distributions, variability and trends in the mountain snow cover and its first-order drivers, temperature and precipitation, is key for a wide range of studies and applications. In this study, we compare the results of various modeling systems (global and regional reanalyses ERA5, ERA5-Land, ERA5-Crocus, CERRA-Land, UERRA MESCAN-SURFEX and MTMSI and regional climate model simulations CNRM-ALADIN and CNRM-AROME driven by the global reanalysis ERA-Interim) against observational references (in situ, gridded observational datasets and satellite observations) across the European Alps from 1950 to 2020. The comparisons are performed in terms of monthly and seasonal snow cover variables (snow depth and snow cover duration) and their main atmospherical drivers (near-surface temperature and precipitation). We assess multi-annual averages of regional and subregional mean values, their interannual variations, and trends over various timescales, mainly for the winter period (from November through April). ERA5, ERA5-Crocus, MESCAN-SURFEX, CERRA-Land and MTMSI offer a satisfying description of the monthly snow evolution. However, a spatial comparison against satellite observation indicates that all datasets overestimate the snow cover duration, especially the melt-out date. CNRM-AROME and CNRM-ALADIN simulations and ERA5-Land exhibit an overestimation of the snow accumulation during winter, increasing with elevations. The analysis of the interannual variability and trends indicates that modeling snow cover dynamics remains complex across multiple scales and that none of the models evaluated here fully succeed to reproduce this compared to observational reference datasets. Indeed, while most of the evaluated model outputs perform well at representing the interannual to multi-decadal winter temperature and precipitation variability, they often fail to address the variability in the snow depth and snow cover duration. We discuss several artifacts potentially responsible for incorrect long-term climate trends in several reanalysis products (ERA5 and MESCAN-SURFEX), which we attribute primarily to the heterogeneities of the observation datasets assimilated. Nevertheless, many of the considered datasets in this study exhibit past trends in line with the current state of knowledge. Based on these datasets, over the last 50 years (1968–2017) at a regional scale, the European Alps have experienced a winter warming of 0.3 to 0.4 ∘C per decade, stronger at lower elevations, and a small reduction in winter precipitation, homogeneous with elevation. The decline in the winter snow depth and snow cover duration ranges from −7 % to −15 % per decade and from −5 to −7 d per decade, respectively, both showing a larger decrease at low and intermediate elevations. Overall, we show that no modeling strategy outperforms all others within our sample and that upstream choices (horizontal resolution, heterogeneity of the observations used for data assimilation in reanalyses, coupling between surface and atmosphere, level of complexity, configuration of the snow scheme, etc.) have great consequences on the quality of the datasets and their potential use. Despite their limitations, in many cases they can be used to characterize the main features of the mountain snow cover for a range of applications.

Accurate snow phenology detection, including snow cover days (SCD), snow start date (SSD), and snow end date (SED), is increasingly important for understanding mountain hydrology such as snow heterogeneity and snowmelt seasonality. Multiple cloud‐free daily snow cover products have recently been developed in China, employing diverse retrieval algorithms and cloud‐gap‐filling methods, resulting in varying accuracy levels. However, comprehensive analysis of differences among products and their impact on snow phenology detection is lacking. This study systematically evaluates eight state‐of‐the‐art snow cover products in China, focusing on the challenging Tibetan Plateau (TP). We introduce a novel metric, the consistency‐weighted correlation coefficient (CWR), customized for SSD and SED detection, and propose product‐combining schemes like “ensemble voting” and “sensor preference” to enhance reliability. Our findings highlight the prime influence of retrieval algorithms under clear‐sky conditions on accuracy, surpassing the importance of cloud‐gap‐filling methods. Specifically, a product optimizing normalized difference snow index thresholds for diverse landcover types consistently outperforms others in detecting all three snow phenology parameters, with correlation coefficients for SCD of 0.82 and 0.69, and CWR values for SSD of 0.54 and 0.40, and for SED of 0.53 and 0.37 in both China and the TP, respectively. Moreover, our proposed scheme combining three high‐accuracy products significantly enhances snow cover identification and SCD detection, especially when the best‐performing product alone faces substantial uncertainty. These findings provide immediate, crucial implications for optimizing the use of multiple cloud‐free products to enhance snow phenology detection, ultimately advancing the applicability of derived snow parameters in mountain hydrology research. Key Points Performances of eight state‐of‐the‐art cloud‐free products and their combinations in detecting snow phenology parameters are evaluated Impacts of retrieval algorithm, cloud‐gap‐filling method, and sensor differences are considered for analyzing products' accuracy differences Proposed scheme combining three high‐accuracy products notably enhances snow cover days detection but including less accurate ones degrades

We used the spatially distributed and physically based snow cover model SNOWGRID-CL to derive daily grids of natural snow conditions and snowmaking potential at a spatial resolution of 1 × 1 km for Austria for the period 1961–2020 validated against homogenized long-term snow observations. Meteorological driving data consists of recently created gridded observation-based datasets of air temperature, precipitation, and evapotranspiration at the same resolution that takes into account the high variability of these variables in complex terrain. Calculated changes reveal a decrease in the mean seasonal (November–April) snow depth (HS), snow cover duration (SCD), and potential snowmaking hours (SP) of 0.15 m, 42 days, and 85 h (26%), respectively, on average over Austria over the period 1961/62–2019/20. Results indicate a clear altitude dependence of the relative reductions (−75% to −5% (HS) and −55% to 0% (SCD)). Detected changes are induced by major shifts of HS in the 1970s and late 1980s. Due to heterogeneous snowmaking infrastructures, the results are not suitable for direct interpretation towards snow reliability of individual Austrian skiing resorts but highly relevant for all activities strongly dependent on natural snow as well as for projections of future snow conditions and climate impact research.

The European Alps stretch over a range of climate zones which affect the spatial distribution of snow. Previous analyses of station observations of snow were confined to regional analyses. Here, we present an Alpine-wide analysis of snow depth from six Alpine countries – Austria, France, Germany, Italy, Slovenia, and Switzerland – including altogether more than 2000 stations of which more than 800 were used for the trend assessment. Using a principal component analysis and k-means clustering, we identified five main modes of variability and five regions which match the climatic forcing zones: north and high Alpine, north-east, north-west, south-east, and south and high Alpine. Linear trends of monthly mean snow depth between 1971 and 2019 showed decreases in snow depth for most stations from November to May. The average trend among all stations for seasonal (November to May) mean snow depth was −8.4 % per decade, for seasonal maximum snow depth −5.6 % per decade, and for seasonal snow cover duration −5.6 % per decade. Stronger and more significant trends were observed for periods and elevations where the transition from snow to snow-free occurs, which is consistent with an enhanced albedo feedback. Additionally, regional trends differed substantially at the same elevation, which challenges the notion of generalizing results from one region to another or to the whole Alps. This study presents an analysis of station snow depth series with the most comprehensive spatial coverage in the European Alps to date.

Mean snow cover duration was derived for the entire continent of Europe based on the MODIS daily snow cover products MOD10A1 and MYD10A1 for the period from 2000 to 2011. Dates of snow cover start and snow cover melt were also estimated. Polar darkness north of ~62°N and extensive cloud coverage affected the daily snow cover, preventing a direct derivation of the desired parameters. Combining sensor data from both MODIS platforms and applying a temporal cloud filter, cloud coverage and polar darkness were removed from the input data and accuracy remained above 90% for 87% of the area. The typical snow cover characteristics of the whole continent are illustrated and constitute a unique dataset with respect to spatial and temporal resolution. Abnormal events, glacier inventories or studies on possible impacts of climate change on snow cover characteristics are only some examples for applications where the presented results may be utilized.

Analysis of in situ and satellite data shows evidence of different regional snow cover responses to the widespread warming and increasing winter precipitation that has characterized the Arctic climate for the past 40–50 years. The largest and most rapid decreases in snow water equivalent (SWE) and snow cover duration (SCD) are observed over maritime regions of the Arctic with the highest precipitation amounts. There is also evidence of marked differences in the response of snow cover between the North American and Eurasian sectors of the Arctic, with the North American sector exhibiting decreases in snow cover and snow depth over the entire period of available in situ observations from around 1950, while widespread decreases in snow cover are not apparent over Eurasia until after around 1980. However, snow depths are increasing in many regions of Eurasia. Warming and more frequent winter thaws are contributing to changes in snow pack structure with important implications for land use and provision of ecosystem services. Projected changes in snow cover from Global Climate Models for the 2050 period indicate increases in maximum SWE of up to 15% over much of the Arctic, with the largest increases (15–30%) over the Siberian sector. In contrast, SCD is projected to decrease by about 10–20% over much of the Arctic, with the smallest decreases over Siberia (<10%) and the largest decreases over Alaska and northern Scandinavia (30–40%) by 2050. These projected changes will have far-reaching consequences for the climate system, human activities, hydrology, and ecology.

High spatiotemporal resolution remote sensing data is crucial for monitoring heterogeneous mountainous snow cover. Although spatiotemporal fusion presents a promising approach for high‐resolution snow monitoring, cloud contamination and sparse observations remain a critical constraint on its large‐scale and long‐term implementation. To address this issue, we propose an adaptive time‐series fusion framework to generate 30‐m daily gapless snow cover data based on the Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM). We integrate multi‐source coarse‐resolution data and multi‐source fine‐resolution data to increase the number of valid pixels and enhance data density to capture the rapid spatiotemporal variations of snow cover. Additionally, we introduce time‐series image pairs to adapt the ESTARFM method, which overcomes the spatial completeness limitation of fine‐resolution data and dynamically selects the spatiotemporal information closest to the target time for each pixel. Comprehensive evaluations confirm the high accuracy of the fused results, as demonstrated by the consistency with reference data (R = 0.776–0.964). Furthermore, validation with ground‐based snow observations shows that the fused 30‐m daily snow cover data not only outperforms the widely used 500‐m data in capturing the temporal dynamics of snow cover, as evidenced by its strong alignment with ground‐based snow phenology metrics, but also provides new insights into the spatial distribution of mountainous snow cover. In areas with elevations below 3,500 m, slopes under 25°, or shaded slopes, the 30‐m data captures small‐scale, sparse, and fragmented snow cover, offering significant potential for hydrological research and practical applications that require accurate snow cover estimation.

Dust particles from high latitudes have a potentially large local, regional, and global significance to climate and the environment as short-lived climate forcers, air pollutants, and nutrient sources. Identifying the locations of local dust sources and their emission, transport, and deposition processes is important for understanding the multiple impacts of high-latitude dust (HLD) on the Earth’s systems. Here, we identify, describe, and quantify the source intensity (SI) values, which show the potential of soil surfaces for dust emission scaled to values 0 to 1 concerning globally best productive sources, using the Global Sand and Dust Storms Source Base Map (G-SDS-SBM). This includes 64 HLD sources in our collection for the northern (Alaska, Canada, Denmark, Greenland, Iceland, Svalbard, Sweden, and Russia) and southern (Antarctica and Patagonia) high latitudes. Activity from most of these HLD sources shows seasonal character. It is estimated that high-latitude land areas with higher (SI ≥ 0.5), very high (SI ≥ 0.7), and the highest potential (SI ≥ 0.9) for dust emission cover > 1 670 000 km2 , > 560 000 km2 , and > 240 000 km2 , respectively. In the Arctic HLD region (≥ 60◦ N), land area with SI ≥ 0.5 is 5.5 % (1 035 059 km2), area with SI ≥ 0.7 is 2.3 % (440 804 km2), and area with SI ≥ 0.9 is 1.1 % (208 701 km2). Minimum SI values in the northern HLD region are about 3 orders of magnitude smaller, indicating that the dust sources of this region greatly depend on weather conditions. Our spatial dust source distribution analysis modeling results showed evidence supporting a northern HLD belt, defined as the area north of 50◦ N, with a “transitional HLD-source area” extending at latitudes 50–58◦ N in Eurasia and 50–55◦ N in Canada and a “cold HLD-source area” including areas north of 60◦ N in Eurasia and north of 58◦ N in Canada, with currently “no dust source” area between the HLD and low-latitude dust (LLD) dust belt, except for British Columbia. Using the global atmospheric transport model SILAM, we estimated that 1.0 % of the global dust emission originated from the high-latitude regions. About 57 % of the dust deposition in snow- and ice-covered Arctic regions was from HLD sources. In the southern HLD region, soil surface conditions are favorable for dust emission during the whole year. Climate change can cause a decrease in the duration of snow cover, retreat of glaciers, and an increase in drought, heatwave intensity, and frequency, leading to the increasing frequency of topsoil conditions favorable for dust emission, which increases the probability of dust storms. Our study provides a step forward to improve the representation of HLD in models and to monitor, quantify, and assess the environmental and climate significance of HLD.

This study aims to examine the temporal impacts of atmospheric-oceanic patterns, on snow cover variations of the widest snow-prone zone in Iran (Central Alborz mountains ∼ 400,000 km 2 ). For this purpose, the snow cover area was derived from Moderate Resolution Imaging Spectroradiometer (MODIS) sensor images during 2000–2016 years. Thirty indices related to atmospheric-oceanic patterns were provided from the National Oceanic and Atmospheric Administration (NOAA) website. The strongest significant teleconnection effects occur with time lags. Atlantic Meridional Mode (AMM), North Pacific pattern (NP), Hurricane Activity (HA), Western Hemisphere Warm Pool (WHWP), Extreme Eastern ENSO (Niño [1 + 2]), and Eastern ENSO (Niño [3]), were the most effective signals on the snow cover variations in Central Alborz region. The AMM, HA, and WHWP indices, were the directly related ones with 2-, 4- and 4-month time lags, respectively. Indices NP, Niño (1 + 2), and Niño (3) were the inversely related signals and their most severe effects appeared after a 1-, 4-, and 4-month lags, respectively. A wavelet coherence analysis was also used for teleconnection investigation. It has illustrated that the mentioned time lags of these signals, have a significant 9–15-months period’ relation with snow cover variations, during the whole study years. In the following, the time lags of the indices were applied as inputs to forecast the snow cover area. The models including the least square support vector machine (LSSVM) and group method of data handling (GMDH) were utilized for this purpose. According to the values of normalized root mean square error (between 0.1 and 0.2) and Nash-Sutcliff (> 0.75), the GMDH had acceptable predictions, especially for long-term forecasting horizons (12- and 24-month-ahead). The applied input approach can have research value for forecasting snow cover area in other mountainous regions of the world, and also to reconstruct the snow data, in the years when satellite imagery was not available.

Snow cover mapping algorithms utilizing multispectral satellite data at various spatial resolutions are available, each treating subpixel variation differently. Past evaluations of snow mapping accuracy typically relied on satellite data collected at a higher spatial resolution than the data in question. However, these optical data cannot characterize snow cover mapping performance under forest canopies or at the meter scale. Here, we use 3 m spatial resolution snow depth maps collected on 116 d by an aerial laser scanner to validate band ratio and spectral-mixture snow cover mapping algorithms. Such a comprehensive evaluation of sub-canopy snow mapping performance has not been undertaken previously. The following standard (produced operationally by an agency) products are evaluated: NASA gap-filled Moderate Resolution Imaging Spectroradiometer (MODIS) MOD10A1F, NASA gap-filled Visible Infrared Imaging Radiometer Suite (VIIRS) VNP10A1F, and United States Geological Survey (USGS) Landsat 8 Level-3 Fractional Snow Covered Area. Two spectral-unmixing approaches are also evaluated: Snow-Covered Area and Grain Size (SCAG) and Snow Property Inversion from Remote Sensing (SPIReS), both of which are gap-filled MODIS products and are also run on Landsat 8. We assess subpixel snow mapping performance while considering the fractional snow-covered area (fSCA), canopy cover, sensor zenith angle, and other variables within six global seasonal snow classes. Metrics are calculated at the pixel and basin scales, including the root-mean-square error (RMSE), bias, and F statistic (a detection measure). The newer MOD10A1F Version 61 and VNP10A1F Version 1 product biases (− 7.1 %, −9.5 %) improve significantly when linear equations developed for older products are applied (2.8 %, −2.7 %) to convert band ratios to fSCA. The F statistics are unchanged (94.4 %, 93.1 %) and the VNP10A1F RMSE improves (18.6 % to 15.7 %), while the MOD10A1F RMSE worsens (12.7 % to 13.7 %). Consistent with previous studies, spectral-unmixing approaches (SCAG, SPIReS) show lower biases (−0.1 %, −0.1 %) and RMSE (12.1 %, 12.0 %), with higher F statistics (95.6 %, 96.1 %) relative to the band ratio approaches for MODIS. Landsat 8 products are all spectral-mixture methods with low biases (−0.4 % to 0.3 %), low RMSE (11.4 % to 15.8 %), and high F statistics (97.3 % to 99.1 %). Spectral-unmixing methods can improve snow cover mapping at the global scale.