Fractional Snow Cover Research Articles

Retrieval of Fractional Snow Cover over High Mountain Asia Using 1 km and 5 km AVHRR/2 with Simulated Mid-Infrared Reflective Band

Accurate long-term snow-covered-area mapping is essential for climate change studies and water resource management. The NOAA AVHRR/2 provides a unique data source for long-term, large-spatial-scale monitoring of snow-covered areas at a daily scale. However, the value of AVHRR/2 in mapping snow-covered areas is limited, due to its lack of a shortwave infrared band for snow/cloud discrimination. We simulated the reflectance in the 3.75 µm mid-infrared band with a radiative transfer model and then developed three fractional-snow-cover retrieval algorithms for AVHRR/2 imagery at 1 km and 5 km resolutions. These algorithms are based on the multiple endmember spectral mixture analysis algorithm (MESMA), snow index (SI) algorithm, and non-snow/snow two endmember model (TEM) algorithm. Evaluation and comparison of these algorithms were performed using 313 scenarios that referenced snow-cover maps from Landsat-5/TM imagery at 30 m resolution. For all the evaluation data, the MESMA algorithm outperformed the other two algorithms, with an overall accuracy of 0.84 (0.85) and an RMSE of 0.23 (0.21) at the 1 km (5 km) scale. Regarding the effect of land cover type, we found that the three AVHRR/2 fractional-snow-cover retrieval algorithms have good accuracy in bare land, grassland, and Himalayan areas; however, the accuracy decreases in forest areas due to the shading of snow by the canopy. Regarding the topographic effect, the accuracy evaluation indices showed a decreasing and then increasing trend as the elevation increased. The accuracy was worst in the 4000–5000 m range, which was due to the severe snow fragmentation in the High Mountain Asia region; the early AVHRR/2 sensors could not effectively monitor the snow cover in this region. In this study, by increasing the number of bands of AVHRR/2 1 km data for fractional-snow-cover retrieval, a good foundation for subsequent long time series kilometre- resolution snow-cover monitoring has been laid.

Evaluation of Snow and Streamflows Using Noah-MP and WRF-Hydro Models in Aroostook River Basin, Maine

Snow influences land–atmosphere interactions in snow-dominated areas, and snow melt contributes to basin streamflows. However, estimating snowpack properties such as the snow depth (SD) and snow water equivalent (SWE) from land surface model simulations remains a challenge. This is, in part, due to uncertainties in the atmospheric forcing variables, which propagate into hydrological model predictions. This study implements the Weather Research and Forecasting (WRF)-Hydro framework with the Noah-Multiparameterization (Noah-MP) land surface model in the NOAA’s National Water Model (NWM) version 2.0 configuration to estimate snow in a single column and subsequently the streamflow across the Aroostook River’s sub-basins in Maine for water years (WY) 2014–2016. This study evaluates how differences between two atmospheric forcing datasets, the North American Land Data Assimilation version 2 (NLDAS-2) and in situ (Station), translate into differences in the simulation of snow. NLDAS-2 was used as the meteorological forcing in the retrospective NWM 2.0 simulations. The results from the single-column study showed that differences in the simulated SWE and SD were linked to differences in the 2 m air temperature (T2m), which influenced the precipitation partitioning of rain and snow, as parameterized in Noah-MP. The negative mean bias of −0.7 K (during the accumulation period) in T2m for NLDAS-2, compared to the Station forcing, was a major factor that contributed to the positive mean bias of +52 mm on average in the peak SWE in the NLDAS-2-forced Noah-MP simulation during the study period. The higher T2m values at the Station led to higher sensible heat fluxes towards the snowpack, which led to a higher amount of net energy at the snow’s surface and melt events during the accumulation season in Station-forced Noah-MP simulations. The results from the retrospective NWM version 2.0′s simulation in the basin showed that the streamflow estimates were closer to the United States Geological Survey gage observations at the two larger sub-basins (NSE = 0.9), which were mostly forested, compared to the two smaller sub-basins (NSE ≥ 0.4), which had more agricultural land-use. This study also showed that the spring snowmelt timing was captured quite well by the timing of the decline in the simulated SWE and SD, providing an early indication of melt in most sub-basins. The simulated fractional snow cover area (fSCA) however provided less information about the changes in snow or onset of snowmelt as it was mostly binary (full snow cover in winter), which differed from the more realistic fSCA values shown by the Moderate Resolution Imaging Spectroradiometer.

Snow Cover in the Three Stable Snow Cover Areas of China and Spatio-Temporal Patterns of the Future

In the context of global warming, relevant studies have shown that China will experience the largest temperature rise in the Qinghai–Tibet Plateau and northwestern regions in the future. Based on MOD10A2 and MYD10A2 snow products and snow depth data, this study analyzes the temporal and spatial evolution characteristics of the snow cover fraction, snow depth, and snow cover days in the three stable snow cover areas in China, and combines 15 modes in CMIP6 snow cover data in four different scenarios with three kinds of variables, predicting the spatiotemporal evolution pattern of snow cover in China’s three stable snow cover areas in the future. The results show that (1) the mean snow cover fraction, snow depth, and snow cover days in the snow cover area of Northern Xinjiang are all the highest. Seasonal changes in the snow cover areas of the Qinghai–Tibet Plateau are the most stable. The snow cover fraction, snow depth, and snow cover days of the three stable snow cover areas are consistent in spatial distribution. The high values are mainly distributed in the southeast and west of the Qinghai–Tibet Plateau, the south and northeast of Northern Xinjiang, and the north of the snow cover area of Northeast China. (2) The future snow changes in the three stable snow cover areas will continue to decline with the increase in development imbalance. Snow cover fraction and snow depth decrease most significantly in the Qinghai–Tibet Plateau and the snow cover days in Northern Xinjiang decrease most significantly under the SSPs585 scenario. In the future, the southeast of the Qinghai–Tibet Plateau, the northwest of Northern Xinjiang, and the north of Northeast China will be the center of snow cover reduction. (3) Under the four different scenarios, the snow cover changes in the Qinghai–Tibet Plateau and Northern Xinjiang are the most significant. Under the SSPs126 and SSPs245 scenarios, the Qinghai–Tibet Plateau snow cover has the most significant change in response. Under the SSPs370 and SSPs585 scenarios, the snow cover in Northern Xinjiang has the most significant change.

Bias adjustment and downscaling of snow cover fraction projections from regional climate models using remote sensing for the European Alps

Abstract. Mountain seasonal snow cover is undergoing major changes due to global climate change. Assessments of future snow cover usually rely on physically based models, and often include post-processed meteorology. Alternatively, we here propose a direct statistical adjustment of snow cover fraction from regional climate models by using long-term remote-sensing observations. We compared different bias-adjustment routines (delta change, quantile mapping, and quantile delta mapping) and explored a downscaling based on historical observations for the Greater Alpine Region in Europe. All bias-adjustment methods account for systematic biases, for example due to topographic smoothing, and reduce model spread in future projections. The trend-preserving methods delta change and quantile delta mapping were found to be more suitable for snow cover fraction than quantile mapping. Averaged over the study region and whole year, snow cover fraction decreases from 12.5 % in 2001–2020 to 10.4 % (8.9 %, 11.5 %; model spread) in 2071–2100 under RCP2.6 (representative concentration pathway), and to 6.4 % (4.1 %, 7.8 %) under RCP8.5 (bias-adjusted estimates from quantile delta mapping). In addition, changes strongly depended on season and elevation. The comparison of the statistical downscaling to a high-resolution physically based model yields similar results for the elevation range covered by the climate models, but different elevation gradients of change above and below. Downscaling showed overall potential but requires further research. Since climate model and remote-sensing observations are available globally, the proposed methods are potentially widely applicable but are limited to snow cover fraction.

Distribution and Attribution of Earlier Start of the Growing Season over the Northern Hemisphere from 2001–2018

The start of the growing season (SOS) is a vital ecological indicator for climate change and the terrestrial ecosystem. Previous studies have reported that the SOS over the Northern Hemisphere (NH) has experienced remarkable changes in the past few decades. However, because of the different spatial and temporal coverages of existing SOS studies, a coherent and robust account for SOS changes in the NH has been lacking. Using satellite-retrieved vegetation-phenology datasets, ground observations, and several auxiliary datasets, this study evaluated the performance of the latest MODIS vegetation-dynamics dataset (MCD12Q2-C6) and explored the distribution and attribution of the SOS to climate change over the NH for the period 2001–2018. The validation results using the Chinese Ecosystem Research Network (CERN) and Lilac-leafing observations (Lilac) displayed that the MCD12Q2-C6 has a good performance in SOS monitoring over the NH mid-latitudes. Meanwhile, evidence from MCD12Q2-C6 pointed out that the SOS was advanced by 2.08 days on average over the NH during 2001–2018, especially for Europe, China, and Alaska, United States. In addition, detailed-sensitivity analysis showed that the increased surface air temperature (Ts) (−1.21 ± 0.34 days °C−1) and reduced snow-cover fraction (Sc) (0.62 ± 0.29 days%−1) were the key driving factors of the observed SOS changes over the NH during 2001–2018. Compared with Ts and Sc, the role of total precipitation (Pt) was minor in dominating the spring vegetation-phenology changes at the same period. The findings of this study contribute to our understanding of the responses of SOS to the competing changes of Ts, Pt, and Sc over the NH.

Deciphering the Role of Aerosol‐Induced Snow Albedo Feedback on Dust Emission Over the Tibetan Plateau

AbstractThe aerosol‐induced snow albedo effect (SAE) and its feedback on the dust emission over the Tibetan Plateau are investigated using the regional climate model (RegCM4.6) coupled with the SNow, ICe, and Aerosol Radiative model. The deposition of absorbing aerosols (black carbon and dust) decreases the albedo of snow over the Himalayan‐Tibetan region, which results in significant surface warming of 1–4 K. This aerosol‐induced surface warming and accelerated snow melting decreased fractional snow cover by 8% and number of snow cover days by 20 days over the region. Compared to the control simulation, this early and more exposure of nonvegetated land resulted in the increase of the dust emission flux by 60% over central Tibet and therefore the columnar dust loading by 10%–40%. This effect alone results in a surface cooling of −4 Wm−2 and also increases the occurrence of dust events by 10%–20%. Therefore, the snow darkening effect due to the deposition of black carbon and dust increases the local emission of dust aerosols, which further increases the aerosol‐induced snow melting and thus perturbs the air quality and radiation balance over the high mountains of Asia.

Subpixel Snow Mapping Using Daily AVHRR/2 Data over Qinghai–Tibet Plateau

Based on a linear spectral mixture model and multiple endmember spectral mixture analysis, using daily advanced very-high-resolution radiometer (AVHRR/2) data of the Qinghai–Tibet Plateau, a subpixel snow mapping algorithm was proposed in this paper, for prolonging the historical time series of the fractional snow cover data to 40 years. In particular, the normalized difference vegetation index (NDVI), and channels 1 and 2 of AVHRR/2 data were used to automatically select the end-members directly, from a certain AVHRR/2 image. A look-up table of sample spectra of mixed pixels and their respective snow cover percentages was introduced for one AVHRR/2 image. According to the established look-up tables, the fractional snow cover of each mixed pixel can then be extracted from the AVHRR/2 images. Before the subpixel snow mapping, the cloud pollution of the AVHRR/2 images was mitigated, with both the thick and thin clouds almost removed from the AVHRR/2 images. It turns out that the processing speed of the subpixel snow mapping is three times faster than the process not using the look-up table. The mapping algorithm was validated against the snow-covered area from Thematic Mapper (TM) data, with the root-mean-square errors (RMSEs) well below 0.12. Results show that the proposed algorithm for subpixel snow mapping is both effective and efficient, especially in such a mountainous region as the Qinghai–Tibet Plateau.

Snow Albedo Feedbacks Enhance Snow Impurity‐Induced Radiative Forcing in the Sierra Nevada

This study employs a fully coupled meteorology‐chemistry‐snow model to investigate the impacts of light‐absorbing particles (LAPs) on snow darkening in the Sierra Nevada. After comprehensive evaluation with spatially and temporally complete satellite retrievals, the model shows that LAPs in snow reduce snow albedo by 0.013 (0–0.045) in the Sierra Nevada during the ablation season (April‐July), producing a midday mean radiative forcing of 4.5 W m−2 which increases to 15–22 W m−2 in July. LAPs in snow accelerate snow aging processes and reduce snow cover fraction, which doubles the albedo change and radiative forcing caused by LAPs. The impurity‐induced snow darkening effects decrease snow water equivalent and snow depth by 20 and 70 mm in June in the Sierra Nevada bighorn sheep habitat. The earlier snowmelt reduces root‐zone soil water content by 20%, deteriorating the forage productivity and playing a negative role in the survival of bighorn sheep.

A conceptual model for simulating streamflow in a changing snow-covered catchment: Application to the data-sparse upper Brahmaputra River basin

Understanding the hydrological processes in the Tibetan Plateau (TP) is increasingly important as the snow resources feed the demand for freshwater for a vast downstream population. Given the limited information available, the task of formulating a hydrological model that characterises future streamflow at downstream locations is challenging. A flexible conceptual hydrological model is proposed relating snow cover fraction (SCF) and snow water equivalent (SWE) relationship based on remote sensing data with the aim of characterising flow in the region into the future where temperature profiles will be considerably warmer than present. The model’s novelty lies in its ability to simulate dynamically varying SCF, SWE, and streamflow using limited available data. Results show that NSE of modelled SCF, SWE and streamflow are 0.82, 0.66 and 0.75 respectively in the calibration period (1981 to 1992) whereas 0.78, 0.56 and 0.70 respectively in the validation period (1993 to 2000).

Trend and variability in the long‐term relationship between Eurasian snow cover and Indian summer monsoon rainfall

AbstractThe predictability of the Indian summer monsoon rainfall (ISMR) at any given time period depends on the strength of its relationship with predictable drivers like the El Nino–Southern Oscillation (ENSO) that are known to undergo significant epochal variations. While a similar relationship between the Eurasian snow cover fraction (SCF) and ISMR has also shown epochal variation in recent decades, its stationarity on centennial or longer timescales remains unknown. Here, on the basis of the observed relationship between the dominant modes of ISMR and SCF over a period of 115 years (1901–2015), we unravel two indices of SCF that encapsulate its space–time variability. We find that the SCF indices and ISMR relationship have a statistically significant increasing trend with a weak multidecadal variability superimposed on it, making significant positive correlations between the two highly probable in the coming decades. With SCF driving the North Atlantic Oscillation (NAO) that in turn driving the ISMR, we find that the NAO plays a pivotal role in modulating the SCF–ISMR teleconnection on a multidecadal time scale. Further, the SCF also regulates the ISMR low‐frequency variability through modulation of the tropospheric temperature gradient of the Indian region.

A Novel Machine Learning–Based Gap-Filling of Fine-Resolution Remotely Sensed Snow Cover Fraction Data by Combining Downscaling and Regression

Abstract Satellite-based remotely sensed observations of snow cover fraction (SCF) can have data gaps in spatially distributed coverage from sensor and orbital limitations. We mitigate these limitations in the example fine-resolution Moderate Resolution Imaging Spectroradiometer (MODIS) data by gap-filling using auxiliary 1-km datasets that either aid in downscaling from coarser-resolution (5 km) MODIS SCF wherever not fully covered by clouds, or else by themselves via regression wherever fully cloud covered. This study’s prototype predicts a 1-km version of the 500-m MOD10A1 SCF target. Due to noncollocatedness of spatial gaps even across input and auxiliary datasets, we consider a recent gap-agnostic advancement of partial convolution in computer vision for both training and predictive gap-filling. Partial convolution accommodates spatially consistent gaps across the input images, effectively implementing a two-dimensional masking. To overcome reduced usable data from noncollocated spatial gaps across inputs, we innovate a fully generalized three-dimensional masking in this partial convolution. This enables a valid output value at a pixel even if only a single valid input variable and its value exist in the neighborhood covered by the convolutional filter zone centered around that pixel. Thus, our gap-agnostic technique can use significantly more examples for training (∼67%) and prediction (∼100%), instead of only less than 10% for the previous partial convolution. We train an example simple three-layer legacy super-resolution convolutional neural network (SRCNN) to obtain downscaling and regression component performances that are better than baseline values of either climatology or MOD10C1 SCF as relevant. Our generalized partial convolution can enable multiple Earth science applications like downscaling, regression, classification, and segmentation that were hindered by data gaps.

How do snow cover fraction change and respond to climate in Altai Mountains of China?

AbstractInvestigating the spatial and temporal changes in snow cover over mountain areas is significant for understanding the impact of regional climate variability. In this study, using cloud‐removed snow cover data, which are generated based on Moderate Resolution Imaging Spectroradiometer (MODIS) daily snow cover products, the spatiotemporal changes of snow cover fraction (SCF) and its relationship with temperature and precipitation changes from 2002 to 2020 were examined over Altai Mountains, China. The results demonstrate that the distribution and changes of SCF are highly spatiotemporally heterogeneous. Within a year, the maximum SCF occurs in January at 98.6%, and the minimum appears in July at 8.7%. The annual‐mean SCF shows an increasing trend at 0.09%·annum−1, owing to the significantly increasing SCF in the snow accumulation period at 0.5%·annum−1 and the decreasing SCF in the snow melting period at −0.2%·annum−1. The SCF distribution, as well as its interannual change, is greatly influenced by elevation. During the snow cover period, a positive linear correlation between SCF and elevation is found at 0.02%·m−1 (p < .01). The annual‐mean SCF decreases in the area below 1,200 m, whereas it increases in the area above 1,200 m. Accordingly, the elevation‐dependent SCF results in various SCF distributions on different slopes and watersheds. The SCF shows an apparent pattern in different aspects, with similar SCFs between the north and east aspects and between the west and south aspects but a difference between the northeast aspects and the westsouth aspects. The SCF is negatively correlated with air temperature (r = −0.74, p < 0.01) and positively correlated with precipitation (r = 0.74, p < 0.01). In addition, temperature shows a significant and larger correlation with SCF in both the snow accumulation and melting periods, indicating the major factor of temperature for the changes in SCF.

Brown Carbon Fuel and Emission Source Attributions to Global Snow Darkening Effect

AbstractSnow and ice albedo reduction due to deposition of absorbing particles (snow darkening effect [SDE]) warms the Earth system and is largely attributed to black carbon (BC) and dust. Absorbing organic aerosol (BrC) also contributes to SDE but has received less attention due to uncertainty and challenges in model representation. This work incorporates the SDE of absorbing organic aerosol (BrC) from biomass burning and biofuel sources into the Snow Ice and Aerosol Radiative (SNICAR) model within a variant of the Community Earth System Model. Additionally, 12 different emission regions of BrC and BC from biomass burning and biofuel sources are tagged to quantify the relative contribution to global and regional SDE. BrC global SDE (0.021–0.056 Wm−2 over land area and 0.0061–0.016 Wm−2 over global area) is larger than other model estimates, corresponding to 37%–98% of the SDE from BC. When compared to observations, BrC simulations have a range in median bias (−2.5% to +21%), with better agreement in the simulations that include BrC photochemical bleaching. The largest relative contributions to global BrC SDE are traced to Northern Asia (23%–31%), Southeast Asia (16%–21%), and South Africa (13%–17%). Transport from Southeast Asia contributes nearly half of the regional BrC SDE in Antarctica (0.084–0.3 Wm−2), which is the largest regional input to global BrC SDE. Lower latitude BrC SDE is correlated with snowmelt, in‐snow BrC concentrations, and snow cover fraction, while polar BrC SDE is correlated with surface insolation and snowmelt. This indicates the importance of in‐snow processes and snow feedbacks on modeled BrC SDE.

Benchmarking algorithm changes to the Snow CCI+ snow water equivalent product

The European Space Agency (ESA) Snow Climate Change Initiative (CCI+) provides long-term, global time series of daily snow cover fraction and snow water equivalent (SWE). The Snow CCI+ SWE Version 1 (CCIv1) product is built on the GlobSnow algorithm, which combines passive microwave (PMW) data with in situ snow depth (SD) measurements to estimate SWE. While CCIv1 remains algorithmically similar to the most recent GlobSnow product (GlobSnow Version 3), Snow CCI+ SWE Version 2 (CCIv2) incorporates two notable differences. CCIv2 uses updated PMW data from the NASA MEaSUREs Calibrated Passive Microwave Daily EASE-Grid 2.0 Earth Science Data Record and is generated in EASE-Grid 2.0 with 12.5 km grid spacing. It also adjusts SWE retrievals in post-processing by incorporating spatially and temporally varying snow density information. Due to the phased product development framework CCI+ employs, proposed changes between CCIv1 and CCIv2 were implemented in a series of step-wise developmental datasets. Using these developmental datasets, we analyze how changes to input PMW and SD data and the snow density parameterization affect the resulting SWE product. Using in situ snow courses as reference data, we demonstrate that the correlation and RMSE of the CCIv2 developmental product improved 18% (0.10) and 12% (5 mm), respectively, relative to CCIv1. The timing of peak snow mass is shifted two weeks later and a temporal discontinuity in the monthly northern hemisphere snow mass time series associated with the shift from the Special Sensor Microwave/Imager (SSM/I) to the Special Sensor Microwave Imager/Sounder (SSMIS) in 2009 is also removed.

Sentinel-1 time series for mapping snow cover depletion and timing of snowmelt in Arctic periglacial environments: case study from Zackenberg and Kobbefjord, Greenland

Abstract. Snow cover (SC) and timing of snowmelt are key regulators of a wide range of Arctic ecosystem functions. Both are strongly influenced by the amplified Arctic warming and essential variables to understand environmental changes and their dynamics. This study evaluates the potential of Sentinel-1 (S-1) synthetic aperture radar (SAR) time series for monitoring SC depletion and snowmelt with high spatiotemporal resolution to capture their understudied small-scale heterogeneity. We use 97 dual-polarized S-1 SAR images acquired over northeastern Greenland and 94 over southwestern Greenland in the interferometric wide swath mode from the years 2017 and 2018. Comparison of S-1 intensity against SC fraction maps derived from orthorectified terrestrial time-lapse imagery indicates that SAR backscatter can increase before a decrease in SC fraction is observed. Hence, the increase in backscatter is related to changing snowpack properties during the runoff phase as well as decreasing SC fraction. We here present a novel empirical approach based on the temporal evolution of the SAR signal to identify start of runoff (SOR), end of snow cover (EOS) and SC extent for each S-1 observation date during melt using backscatter thresholds as well as the derivative. Comparison of SC with orthorectified time-lapse imagery indicates that HV polarization outperforms HH when using a global threshold. The derivative avoids manual selection of thresholds and adapts to different environmental settings and seasonal conditions. With a global configuration (threshold: 4 dB; polarization: HV) as well as with the derivative, the overall accuracy of SC maps was in all cases above 75 % and in more than half of cases above 90 %. Based on the physical principle of SAR backscatter during snowmelt, our approach is expected to work well in other low-vegetation areas and, hence, could support large-scale SC monitoring at high spatiotemporal resolution (20 m, 6 d) with high accuracy.

Improving Fractional Snow Cover Retrieval From Passive Microwave Data Using a Radiative Transfer Model and Machine Learning Method

Optical sensors are subject to cloud obscuration and sunlight dependence, resulting in large proportions of missing snow cover information. Microwave sensors are a good alternative to snow cover monitoring in all weather conditions. Thus far, few studies in the literature have directly derived the fractional snow cover (FSC) from passive microwave data, and none have considered the relationship between FSC and brightness temperature (TB). This study first explores the FSC–TB relationship with a radiation transfer model, exhibiting that no generic function can properly describe the nonlinear and complex FSC–TB relationship. Therefore, a new algorithm based on machine learning method was designed to improve FSC retrieval from TB data, considering other auxiliary information, including soil property, land surface, and geography information. Benchmarked against the Moderate Resolution Imaging Spectroradiometer (MODIS) reference FSC, our FSC retrieval model performed well with an average correlation coefficient of 0.70, the mean absolute error ranging from 0.15 to 0.17, and the root-mean-square error ranging from 0.19 to 0.21. The generated FSC maps reasonably characterized the seasonal dynamics and spatial distribution patterns of snow cover; time series analysis with three AmeriFlux stations observation indicated effective capture of snowpack evolution process by the generated FSC. In addition, the verification of snow mapping capability using snow depth measurements from 13 521 stations indicates that it was relatively stable with overall accuracy greater than 0.88. For precise monitoring of snow cover extent in all weather conditions, particularly for subpixel snow cover areas, the development of the FSC estimation scheme with TB data should be extensively encouraged and implemented.

Characterization of NDSI Variation: Implications for Snow Cover Mapping

The normalized difference snow index (NDSI) plays an important role in mapping snow cover with spaceborne visible and shortwave-infrared imagery. The NDSI variation depends on illuminating-viewing geometry and snow physical properties, including equivalent grain size (EGS), snow depth (SD) and impurity concentration, as well as fractional snow cover (FSC) within a mixed pixel; however, it is still not fully understood. To quantifiably characterize the pattern of snow NDSI variation, we use a light scattering model of snow to calculate bidirectional reflectance and consequent NDSI values for a wide range of illuminating-viewing geometries, SD, and EGS values. In the model designated bicontinuous snow model using Geometric Optics theory and Radiative Transfer Equation (bicontinuous-GO/RTE), snowpack is represented by a bicontinuous microstructure, and bidirectional reflectance is simulated based on geometric optics and vector radiative transfer equation. The discrete ordinates radiative transfer (DISORT) algorithm is used to simulate the soot concentration effect on snow NDSI. A soil spectral reflectance model (SOILSPECT) is utilized with the assumption of the linear spectral mixture of snow and soil to quantify the effect of FSC on NDSI. As for discontinuous forests, an analytical hybrid geometric-optical and radiative transfer (GORT) model in conjunction with the bicontinuous-GO/RTE model and the PROpriétésSPECTrales (PROSPECT) model is used to examine the effect of canopy cover, which is related to the maximal FSC viewable to satellites. Modeling results indicate that: 1) snow NDSI is comparably low at off-nadir viewing angles, and this effect is exacerbated by the decline in solar elevation but limited by large EGS; 2) snow NDSI increases with EGS yet becomes saturated at EGS of 500 <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\mu \text{m}$ </tex-math></inline-formula> ; 3) the effect of SD that is as low as 1 cm on NDSI is rarely noticeable; 4) the concentration of internally mixed soot up to 1 ppm has little reducing effect on snow NDSI (<0.05); 5) the angular effects on the NDSI thresholding method to identify binary snow cover for soil background are nonsignificantly pronounced in the cross plane (CP); 6) NDSI increases linearly with FSC for soil background at nadir as well as in the CP; 7) the SNOw-MAPping algorithm (SNOMAP) method could omit snow cover at large zenith angles for solar zenith angles <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$&gt; 20^{\circ }$ </tex-math></inline-formula> in forests; and 8) forests complicate the nonlinear relationship between NDSI and canopy cover with fully snow-covered ground beneath canopies. These findings imply important uncertainty sources of binary and FSC mapping with NDSI.

Attribution of the spatial heterogeneity of Arctic surface albedo feedback to the dynamics of vegetation, snow and soil properties and their interactions

The Arctic warming rate is triple the global average, which is partially caused by surface albedo feedback (SAF). Understanding the varying pattern of SAF and the mechanisms is therefore critical for predicting future Arctic climate under anthropogenic warming. To date, however, how the spatial pattern of seasonal SAF is influenced by various land surface factors remains unclear. Here, we aim to quantify the strengths of seasonal SAF across the Arctic and to attribute its spatial heterogeneity to the dynamics of vegetation, snow and soil as well as their interactions. The results show a large positive SAF above −5% K−1 across Baffin Island in January and eastern Yakutia in June, while a large negative SAF beyond 5% K−1 is observed in Canada, Chukotka and low latitudes of Greenland in January and Nunavut, Baffin Island and Krasnoyarsk Krai in July. Overall, a great spatial heterogeneity of Arctic land warming induced by positive SAF is found with a coefficient of variation (CV) larger than 61.5%, and the largest spatial difference is detected in wintertime with a CV > 643.9%. Based on the optimal parameter-based geographic detector model, the impacts of snow cover fraction (SCF), land cover type (LC), normalized difference vegetation index (NDVI), soil water content (SW), soil substrate chemistry (SC) and soil type (ST) on the spatial pattern of positive SAF are quantified. The rank of determinant power is SCF > LC > NDVI > SW > SC > ST, which indicates that the spatial patterns of snow cover, land cover and vegetation coverage dominate the spatial heterogeneity of positive SAF in the Arctic. The interactions between SCF, LC and SW exert further influences on the spatial pattern of positive SAF in March, June and July. This work could provide a deeper understanding of how various land factors contribute to the spatial heterogeneity of Arctic land warming at the annual cycle.

Ground- and Satellite-Based Evaluation of WRF Snowfall Prediction

This study performed 4-day numerical integration in 1-hour intervals using the Weather Research and Forecasting (WRF) model for four major cases of heavy snowfall that occurred from 2020 to 2021. The model-predicted snow depth data were compared with the ground-observed snow depth and the satellite-observed snow cover data and then were statistically verified. The scalar verification results for ground data from the four cases showed a root–mean–square error of 2.55-16.67 cm and a correlation coefficient of 0.48-0.80, whereas the verification results with satellite data showed the correlation coefficients of 0.38-0.60. For categorical verification, using a threshold value of a snow depth exceeding 5 cm, the proportion correct was 90% or higher for ground observations of each case. In addition, in the satellite categorical verification, when the threshold value of the Snow Cover Fraction (SCF) exceeds 0.5, the proportion correct was 50% or more. These results are meaningful because the model snow depth verification methods were devised strategically for the first time using both the snow depth data of the mesoscale ground observation networks and ultra-high-resolution Sentinel-2 satellite data currently available in Korea. The findings of this study will contribute to the development of a high-resolution numerical prediction model and its verification methodology for snowfalls in the Korean Peninsula, eventually leading to increased prediction accuracy and reduced snow damage.

Binary and Fractional MODIS Snow Cover Mapping Boosted by Machine Learning and Big Landsat Data

It is promising to improve MODIS snow cover mapping by training effective machine learning models. However, considering the strong spatiotemporal heterogeneity of snow cover, the representativeness of training samples becomes a critical issue to ensure the model effectiveness over large areas and long terms. To deal with this issue, we propose the <i>m</i>-day dynamic training strategy. The core of this strategy is to divide a long-term snow cover mapping task into multiple short-term tasks with consecutive <i>m</i> days (<i>m</i> is greater or equal to the number of days that higher-resolution remote sensing images used as equivalent ground references cover the study area once). It can ensure the spatial representativeness of training samples by utilizing all available equivalent ground references data, and reduce the issue caused by snow cover changing over time. We apply this strategy to random forest for both binary (BSC) and fractional (FSC) snow cover mapping using Landsat-8 data as equivalent ground reference and validate its effect on the Tibetan Plateau. Results show that this strategy achieves higher accuracies than other training strategies, with <i>F1</i> of 0.9596 for BSC mapping, <i>R</i> and <i>MSE</i> of 0.9037 and 0.0260 for FSC mapping, respectively. In addition, spatiotemporal analysis of results further demonstrates that this strategy holds advantages for snow cover mapping in areas of complex terrain and in periods when snow cover changes rapidly. In conclusion, this strategy can effectively ensure the representativeness of training samples generated from Landsat-8 data and thus improve MODIS snow cover mapping through machine learning.