Snow cover variability and snowmelt in a high‐altitude ungauged catchment

Snow variability is an integrated indicator of climate change, and it has important impacts on runoff regimes and water availability in high‐altitude catchments. Remote sensing techniques can make it possible to quantitatively detect the snow cover changes and associated hydrological effects in those poorly gauged regions. In this study, the spatial–temporal variations of snow cover and snow melting time in the Tuotuo River basin, which is the headwater of the Yangtze River, were evaluated based on satellite information from the Moderate Resolution Imaging Spectroradiometer snow cover product, and the snow melting equivalent and its contribution to the total runoff and baseflow were estimated by using degree–day model. The results showed that the snow cover percentage and the tendency of snow cover variability increased with rising altitude. From 2000 to 2012, warmer and wetter climate change resulted in an increase of the snow cover area. Since the 1960s, the start time for snow melt has become earlier by 0.9–3 days/10a and the end time of snow melt has become later by 0.6–2.3 days/10a. Under the control of snow cover and snow melting time, the equivalent of snow melting runoff in the Tuotuo River basin has been fluctuating. The average contributions of snowmelt to baseflow and total runoff were 19.6% and 6.8%, respectively. Findings from this study will serve as a reference for future research in areas where observational data are deficient and for planning of future water management strategies for the source region of the Yangtze River. Copyright © 2015 John Wiley & Sons, Ltd.     

Developing machine learning-based snow depletion curves and analysing their sensitivity over complex mountainous areas

A snow depletion curve (SDC), the relationship between snow mass (e.g., snow depth [SD]) and fractional snow cover area (SCF), is essential to parameterize the effect of snowpack within a physically based snow model. Existing SDCs are constructed using traditional statistic methods may not be applicable in complex mountainous areas. In this study, we developed an information fusion framework to define the relationship between SCF and SD as well as 12 auxiliary factors by using a traditional statistical method and four prevailing machine learning (ML) algorithms, which have comprehensively considered the variable conditions that cause spatiotemporal heterogeneity of snow cover. We also performed a single-dimensional sensitivity analysis to investigate the physical rationality of the newly developed SDCs. The Northern Xinjiang, Northwest China, is selected as the study area, and the data from 46 meteorological stations covering five snow seasons from 2010 to 2015 are used. The results illustrated that ML techniques can be used to establish high-accuracy and robust SDCs for complex mountainous areas. Compared with SDCs constructed by traditional statistical, the performance of the four ML-based SDCs is significantly improved, the RMSE values can be reduced by 50%, R² above 0.75, and an average relative variance close to 0. ML-based SDCs predicted SCF values showed a range of sensitivities to different input variables (e.g., Land surface temperature, aspect, longwave radiation and land cover type), in addition to SD, that were physically representative of effects that snow cover is sensitive to. Moreover, the complexity of SDCs can be reduced by removing insensitive input variables.     

Coupling the WRF model with a temperature index model based on remote sensing for snowmelt simulations in a river basin in the Altay Mountains,north‐west China

Snowmelt water is an important freshwater resource in the Altay Mountains in north‐west China; however, warming climate and rapid spring snowmelt can cause floods that endanger both public and personal property and safety. This study simulates snowmelt in the Kayiertesi River catchment using a temperature index model based on remote sensing coupled with high‐resolution meteorological data obtained from National Centers for Environmental Prediction (NCEP) reanalysis fields that were downscaled using the Weather Research Forecasting model and then bias corrected using a statistical downscaled model. Validation of the forcing data revealed that the high‐resolution meteorological fields derived from the downscaled NCEP reanalysis were reliable for driving the snowmelt model. Parameters of the temperature index model based on remote sensing were calibrated for spring 2014, and model performance was validated using Moderate Resolution Imaging Spectroradiometer snow cover and snow observations from spring 2012. The results show that the temperature index model based on remote sensing performed well, with a simulation mean relative error of 6.7% and a Nash–Sutcliffe efficiency of 0.98 in spring 2012 in the river of Altay Mountains. Based on the reliable distributed snow water equivalent simulation, daily snowmelt run‐off was calculated for spring 2012 in the basin. In the study catchment, spring snowmelt run‐off accounts for 72% of spring run‐off and 21% of annual run‐off. Snowmelt is the main source of run‐off for the catchment and should be managed and utilized effectively. The results provide a basis for snowmelt run‐off predictions, so as to prevent snowmelt‐induced floods, and also provide a generalizable approach that can be applied to other remote locations where high‐density, long‐term observational data are lacking. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluation of the MODIS snow cover fraction product

Eleven years of daily 500 m gridded Terra Moderate Resolution Imaging Spectroradiometer (MODIS) (MOD10A1) snow cover fraction (SCF) data are evaluated in terms of snow presence detection in Colorado and Washington states. The SCF detection validation study is performed using in‐situ measurements and expressed in terms of snow and land detection and misclassification frequencies. A major aspect addressed in this study involves the shifting of pixel values in time due to sensor viewing angles and gridding artifacts of MODIS sensor products. To account for this error, 500 m gridded pixels are grouped and aggregated to different‐sized areas to incorporate neighboring pixel information. With pixel aggregation, both the probability of detection (POD) and the false alarm ratios increase for almost all cases. Of the false negative (FN) and false positive values (referred to as the total error when combined), FN estimates dominate most of the total error and are greatly reduced with aggregation. The greatest POD increases and total error reductions occur with going from a single 500 m pixel to 3×3‐pixel averaged areas. Since the MODIS SCF algorithm was developed under ideal conditions, SCF detection is also evaluated for varying conditions of vegetation, elevation, cloud cover and air temperature. Finally, using a direct insertion data assimilation approach, pixel averaged MODIS SCF observations are shown to improve modeled snowpack conditions over the single pixel observations due to the smoothing of more error‐prone observations and more accurately snow‐classified pixels. Copyright © 2012 John Wiley & Sons, Ltd.     

Observations of late winter Canadian tundra snow cover properties

Tundra snow cover is important to monitor as it influences local, regional, and global‐scale surface water balance, energy fluxes, as well as ecosystem and permafrost dynamics. Observations are already showing a decrease in spring snow cover duration at high latitudes, but the impact of changing winter season temperature and precipitation on variables such as snow water equivalent (SWE) is less clear. A multi‐year project was initiated in 2004 with the objective to quantify tundra snow cover properties over multiple years at a scale appropriate for comparison with satellite passive microwave remote sensing data and regional climate and hydrological models. Data collected over seven late winter field campaigns (2004 to 2010) show the patterns of snow depth and SWE are strongly influenced by terrain characteristics. Despite the spatial heterogeneity of snow cover, several inter‐annual consistencies were identified. A regional average density of 0.293 g/cm³ was derived and shown to have little difference with individual site densities when deriving SWE from snow depth measurements. The inter‐annual patterns of SWE show that despite variability in meteorological forcing, there were many consistent ratios between the SWE on flat tundra and the SWE on lakes, plateaus, and slopes. A summary of representative inter‐annual snow stratigraphy from different terrain categories is also presented. © 2013 Her Majesty the Queen in Right of Canada. Hydrological Processes. © 2013 John Wiley & Sons, Ltd.     

Representation of canopy snow interception,unloading and melt in a parsimonious snowmelt model

Recent improvements in the Utah Energy Balance (UEB) snowmelt model are focused on snow–vegetation–atmosphere interactions to understand how different types of vegetation affect snow processes in the mountains of Western USA. This work presents field work carried out in the Rocky Mountains of Northern Utah to evaluate new UEB model algorithms that represent the processes of canopy snow interception, sublimation, mass unloading and melt. Four years' continuous field observations showed generally smaller accumulations of snow beneath the forest canopies in comparison with open (sage and grass) areas, a difference that is attributed to interception and subsequent sublimation and redistribution of intercepted snow by wind, much of it into surrounding open areas. Accumulations beneath the denser forest (conifer) canopies were found to be less than the accumulation beneath the less dense forest (deciduous) canopies. The model was able to represent the accumulation of snow water equivalent in the open and beneath the deciduous forest quite well but without accounting for redistribution tended to overestimate the snow water equivalent beneath the conifer forest. Evidence of redistribution of the intercepted snow from the dense forest (i.e. conifer forest) to the adjacent area was inferred from observations. Including a simple representation of redistribution in the model gave satisfactory prediction of snow water equivalent beneath the coniferous forest. The simulated values of interception, sublimation and unloading were also compared with previous studies and found in agreement. Copyright © 2013 John Wiley & Sons, Ltd.     

Evaluation of snow products over the Tibetan Plateau

Four satellite‐based snow products are evaluated over the Tibetan Plateau for the 2007–2010 snow seasons. The Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua snow cover daily L3 Global 500‐m grid products (MOD10A1 and MYD10A1), the National Oceanic and Atmospheric Administration Interactive Multisensor Snow and Ice Mapping System (IMS) daily Northern Hemisphere snow cover product and the Advanced Microwave Scanning Radiometer – Earth Observing System Daily Snow Water Equivalent were validated against Thematic Mapper (TM) snow cover maps of Landsat‐5 and meteorological station snow depth observations. The overall accuracy of MOD10A1, MYD10A1 and IMS is higher than 91% against stations observations and than 79% against Landsat TM images. In general, the daily MODIS snow cover products show better performance than the multisensor IMS product. However, the IMS snow cover product is suitable for larger scale (~4km) analysis and applications, with the advantage over MODIS to allow for mitigation for cloud cover. The accuracy of the three products decreases with decreasing snow depth. Overestimation errors are most common over forested regions; the IMS and Advanced Microwave Scanning Radiometer – Earth Observing System Snow Water Equivalent products also show poorer performance that the MODIS products over grassland. By identifying weaknesses in the satellite products, this study provides a focus for the improvement of snow products over the Tibetan plateau. The quantitative evaluation of the products proposed here can also be used to assess their relative weight in data assimilation, against other data sources, such as modelling and in situ measurement networks. Copyright © 2015 John Wiley & Sons, Ltd.     

Quantifying snow water equivalent using terrestrial ground penetrating radar and unmanned aerial vehicle photogrammetry

This study demonstrates the potential value of a combined unmanned aerial vehicle (UAV) Photogrammetry and ground penetrating radar (GPR) approach to map snow water equivalent (SWE) over large scales. SWE estimation requires two different physical parameters (snow depth and density), which are currently difficult to measure with the spatial and temporal resolution desired for basin-wide studies. UAV photogrammetry can provide very high-resolution spatially continuous snow depths (SD) at the basin scale, but does not measure snow densities. GPR allows nondestructive quantitative snow investigation if the radar velocity is known. Using photogrammetric snow depths and GPR two-way travel times (TWT) of reflections at the snow-ground interface, radar velocities in snowpack can be determined. Snow density (RSN) is then estimated from the radar propagation velocity (which is related to electrical permittivity of snow) via empirical formulas. A Phantom-4 Pro UAV and a MALA GX450 HDR model GPR mounted on a ski mobile were used to determine snow parameters. A snow-free digital surface model (DSM) was obtained from the photogrammetric survey conducted in September 2017. Then, another survey in synchronization with a GPR survey was conducted in February 2019 whilst the snowpack was approximately at its maximum thickness. Spatially continuous snow depths were calculated by subtracting the snow-free DSM from the snow-covered DSM. Radar velocities in the snowpack along GPR survey lines were computed by using UAV-based snow depths and GPR reflections to obtain snow densities and SWEs. The root mean square error of the obtained SWEs (384 mm average) is 63 mm, indicating good agreement with independent SWE observations and the error lies within acceptable uncertainty limits.     

Improvement of snowpack simulations in a regional climate model

To improve simulations of regional‐scale snow processes and related cold‐season hydroclimate, the Community Land Model version 3 (CLM3), developed by the National Center for Atmospheric Research (NCAR), was coupled with the Pennsylvania State University/NCAR fifth‐generation Mesoscale Model (MM5). CLM3 physically describes the mass and heat transfer within the snowpack using five snow layers that include liquid water and solid ice. The coupled MM5–CLM3 model performance was evaluated for the snowmelt season in the Columbia River Basin in the Pacific Northwestern United States using gridded temperature and precipitation observations, along with station observations. The results from MM5–CLM3 show a significant improvement in the SWE simulation, which has been underestimated in the original version of MM5 coupled with the Noah land‐surface model. One important cause for the underestimated SWE in Noah is its unrealistic land‐surface structure configuration where vegetation, snow and the topsoil layer are blended when snow is present. This study demonstrates the importance of the sheltering effects of the forest canopy on snow surface energy budgets, which is included in CLM3. Such effects are further seen in the simulations of surface air temperature and precipitation in regional weather and climate models such as MM5. In addition, the snow‐season surface albedo overestimated by MM5–Noah is now more accurately predicted by MM5–CLM3 using a more realistic albedo algorithm that intensifies the solar radiation absorption on the land surface, reducing the strong near‐surface cold bias in MM5–Noah. The cold bias is further alleviated due to a slower snowmelt rate in MM5–CLM3 during the early snowmelt stage, which is closer to observations than the comparable components of MM5–Noah. In addition, the over‐predicted precipitation in the Pacific Northwest as shown in MM5–Noah is significantly decreased in MM5–CLM3 due to the lower evaporation resulting from the longer snow duration. Copyright © 2011 John Wiley & Sons, Ltd.     

近50年云南省降雪的气候变化特征

利用云南省1961-2008年120个测站逐日降雪资料,分析了云南省降雪频次和范围的时空特征和气候变化.结果表明:近50年云南省的年降雪频次和范围总体呈减少趋势,平均每年频次约减少4.5频次.各月的降雪频次和范围都呈负趋势.12月降雪频次减少趋势最显著,4月降雪范围减少趋势最显著.降雪频次长期趋势变化有明显的空间变化.对于年降雪频次西北部比东北部和东部减少得多,滇西北降雪频次每年约减少0.44频次.进一步对云南省年降雪量和积雪深度的长期趋势变化进行分析.云南省近50年,降雪范围逐步减少,年降雪量和平均最大积雪深度呈增加的趋势.说明近年来在云南气候趋于暖湿背景下,年降雪频次和范围呈逐渐减少趋势,但强降雪的频次却增加了.