Diagnosing snow accumulation errors in a rain‐snow transitional environment with snow board observations

Diagnosing the source of errors in snow models requires intensive observations, a flexible model framework to test competing hypotheses, and a methodology to systematically test the dominant snow processes. We present a novel process‐based approach to diagnose model errors through an example that focuses on snow accumulation processes (precipitation partitioning, new snow density, and snow compaction). Twelve years of meteorological and snow board measurements were used to identify the main source of model error on each snow accumulation day. Results show that modeled values of new snow density were outside observational uncertainties in 52% of days available for evaluation, while precipitation partitioning and compaction were in error 45% and 16% of the time, respectively. Precipitation partitioning errors mattered more for total winter accumulation during the anomalously warm winter of 2014–2015, when a higher fraction of precipitation fell within the temperature range where partition methods had the largest error. These results demonstrate how isolating individual model processes can identify the primary source(s) of model error, which helps prioritize future research.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Exploration of remotely sensed forest structure and ultrasonic range sensor metrics to improve empirical snow models

Current methods to estimate snow accumulation and ablation at the plot and watershed levels can be improved as new technologies offer alternative approaches to more accurately monitor snow dynamics and their drivers. Here we conduct a meta‐analysis of snow and vegetation data collected in British Columbia to explore the relationships between a wide range of forest structure variables – obtained from Light Detection and Ranging (LiDAR), hemispherical photography (HP) and Landsat Thematic Mapper – and several indicators of snow accumulation and ablation estimated from manual snow surveys and ultrasonic range sensors. By merging and standardizing all the ground plot information available in the study area, we demonstrate how LiDAR‐derived forest cover above 0.5 m was the variable explaining the highest percentage of absolute peak snow water equivalent (SWE) (33%), while HP‐derived leaf area index and gap fraction (45° angle of view) were the best potential predictors of snow ablation rate (explaining 57% of variance). This study reveals how continuous SWE data from ultrasonic sensors are fundamental to obtain statistically significant relationships between snow indicators and structural metrics by increasing mean r² by 20% when compared to manual surveys. The relationships between vegetation and spectral indices from Landsat and snow indicators, not explored before, were almost as high as those shown by LiDAR or HP and thus point towards a new line of research with important practical implications. While the use of different data sources from two snow seasons prevented us from developing models with predictive capacity, a large sample size helped to identify outliers that weakened the relationships and suggest improvements for future research. A concise overview of the limitations of this and previous studies is provided along with propositions to consistently improve experimental designs to take advantage of remote sensing technologies, and better represent spatial and temporal variations of snow. Copyright © 2013 John Wiley & Sons, Ltd.     

Rain on snow induced urban floods in cold maritime climate: Risk,indicators and trends

Urban floods pose a societal and economical risk. This study evaluated the risk and hydro-meteorological conditions that cause pluvial flooding in coastal cities in a cold climate. Twenty years of insurance claims data and up to 97 years of meteorological data were analysed for Reykjavík, Iceland (64.15°N; <100 m above sea level). One third of the city's wastewater collection system is combined, and pipe grades vary from 0.5% to 10%. Results highlight semi-intensive rain (<7 mm/h; ≤3 year return period) in conjunction with snow and frozen ground as the main cause for urban flood risk in a climate which undergoes frequent snow and frost cycles (avg. 13 and 19 per season, respectively). Floods in winter were more common, more severe and affected a greater number of neighbourhoods than during summer. High runoff volumes together with debris remobilized with high winds challenged the capacity of wastewater systems regardless of their age or type (combined vs. separate). The two key determinants for the number of insurance claims were antecedent frost depth and total precipitation volume per event. Two pluvial regimes were particularly problematic: long duration (13–25 h), late peaking rain on snow (RoS), where snowmelt enhanced the runoff intensity, elongated and connected independent rainfall into a singular, more voluminous (20–76 mm) event; shorter duration (7–9 h), more intensive precipitation that evolved from snow to rain. Closely timed RoS and cooling were believed to trigger frost formation. A positive trend was detected in the average seasonal snow depth and volume of rain and snowmelt during RoS events. More emphasis, therefore, needs to be placed on designing and operating urban drainage infrastructure with regard to RoS co-acting with frozen ground. Furthermore, more detailed, routine monitoring of snow and soil conditions is important to predict RoS flood events.     

Inference and uncertainty of snow depth spatial distribution at the kilometre scale in the Colorado Rocky Mountains: the effects of sample size,random sampling,predictor quality,and validation procedures

Historically, observing snow depth over large areas has been difficult. When snow depth observations are sparse, regression models can be used to infer the snow depth over a given area. Data sparsity has also left many important questions about such inference unexamined. Improved inference, or estimation, of snow depth and its spatial distribution from a given set of observations can benefit a wide range of applications from water resource management, to ecological studies, to validation of satellite estimates of snow pack. The development of Light Detection and Ranging (LiDAR) technology has provided non‐sparse snow depth measurements, which we use in this study, to address fundamental questions about snow depth inference using both sparse and non‐sparse observations. For example, when are more data needed and when are data redundant? Results apply to both traditional and manual snow depth measurements and to LiDAR observations. Through sampling experiments on high‐resolution LiDAR snow depth observations at six separate 1.17‐km² sites in the Colorado Rocky Mountains, we provide novel perspectives on a variety of issues affecting the regression estimation of snow depth from sparse observations. We measure the effects of observation count, random selection of observations, quality of predictor variables, and cross‐validation procedures using three skill metrics: percent error in total snow volume, root mean squared error (RMSE), and R². Extremes of predictor quality are used to understand the range of its effect; how do predictors downloaded from internet perform against more accurate predictors measured by LiDAR? Whereas cross validation remains the only option for validating inference from sparse observations, in our experiments, the full set of LiDAR‐measured snow depths can be considered the ‘true’ spatial distribution and used to understand cross‐validation bias at the spatial scale of inference. We model at the 30‐m resolution of readily available predictors, which is a popular spatial resolution in the literature. Three regression models are also compared, and we briefly examine how sampling design affects model skill. Results quantify the primary dependence of each skill metric on observation count that ranges over three orders of magnitude, doubling at each step from 25 up to 3200. Whereas uncertainty (resulting from random selection of observations) in percent error of true total snow volume is typically well constrained by 100–200 observations, there is considerable uncertainty in the inferred spatial distribution (R²) even at medium observation counts (200–800). We show that percent error in total snow volume is not sensitive to predictor quality, although RMSE and R² (measures of spatial distribution) often depend critically on it. Inaccuracies of downloaded predictors (most often the vegetation predictors) can easily require a quadrupling of observation count to match RMSE and R² scores obtained by LiDAR‐measured predictors. Under cross validation, the RMSE and R² skill measures are consistently biased towards poorer results than their true validations. This is primarily a result of greater variance at the spatial scales of point observations used for cross validation than at the 30‐m resolution of the model. The magnitude of this bias depends on individual site characteristics, observation count (for our experimental design), and sampling design. Sampling designs that maximize independent information maximize cross‐validation bias but also maximize true R². The bagging tree model is found to generally outperform the other regression models in the study on several criteria. Finally, we discuss and recommend use of LiDAR in conjunction with regression modelling to advance understanding of snow depth spatial distribution at spatial scales of thousands of square kilometres. Copyright © 2012 John Wiley & Sons, Ltd.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

A rain on snow climatology and temporal analysis for the eastern United States

ABSTRACT

Rain-on-snow (ROS) has the potential to produce devastating floods by enhancing runoff from snowmelt. Although a common phenomenon across the eastern United States, little research has focused on ROS in this region. This study used a gridded observational snow dataset from 1960–2009 to establish a comprehensive seasonal climatology of ROS for this region. Additionally, different rain and snow thresholds were compared while considering temporal trends in ROS occurrence at four grid cells representing individual locations. Results show most ROS events occur in MAM (March-April-May). ROS events identified with rainfall >1 cm are more frequent near the east coast and events identified with >1 cm snow loss are more common in higher latitudes and/or elevations. Decreasing trends in DJF (December-January-February) ROS events were identified near the coastal areas, with increasing trends in the northern portion of the domain. Significant decreasing trends in MAM ROS are likewise present on a regional scale. Factors playing a role in snowpack depth and rainfall, such as movement of storm tracks in this region, should be considered with future work to discern mechanisms causing the changes in ROS frequency.     

Development of a rain‐on‐snow detection algorithm using passive microwave radiometry

Currently observed climate warming in the Arctic has numerous consequences. Of particular relevance, the precipitation regime is modified where mixed and liquid precipitation can occur during the winter season leading to rain‐on‐snow (ROS) events. This phenomenon is responsible for ice crust formation, which has a significant impact on ecosystems (such as biological, hydrological, ecological and physical processes). The spatially and temporally sporadic nature of ROS events makes the phenomenon difficult to monitor using meteorological observations. This paper focuses on the detection of ROS events using passive microwave (PMW) data from a modified brightness temperature (T_B) gradient approach at 19 and 37 GHz. The approach presented here was developed empirically for observed ROS events with coincident ground‐based PMW measurements in Sherbrooke, Quebec, Canada. It was then tested in Nunavik, Quebec, with the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR‐E). We obtained a detection accuracy of 57, 71 and 89% for ROS detection for three AMSR‐E grid cells with a maximum error of 7% when considering all omissions and commissions with regard to the total number of AMSR‐E passes throughout the winter period. Copyright © 2016 John Wiley & Sons, Ltd.     

Impact of climate indicators on continental‐scale potential groundwater recharge in Africa

In the last decades, human activity has been contributing to climate change that is closely associated with an increase in temperatures, increase in evaporation, intensification of extreme dry and wet rainfall events, and widespread melting of snow and ice. Understanding the intricate linkage between climate warming and the hydrological cycle is crucial for sustainable management of groundwater resources, especially in a vulnerable continent like Africa. This study investigates the relationship between climate‐change drivers and potential groundwater recharge (PGR) patterns across Africa for a long‐term record (1960–2010). Water‐balance components were simulated by using the PCR‐GLOBWB model and were reproduced in both gridded maps and latitudinal trends that vary in space with minima on the Tropics and maxima around the Equator. Statistical correlations between temperature, storm occurrences, drought, and PGR were examined in six climatic regions of Africa. Surprisingly, different effects of climate‐change controls on PGR were detected as a function of latitude in the last three decades (1980–2010). Temporal trends observed in the Northern Hemisphere of Africa reveal that the increase in temperature is significantly correlated to the decline of PGR, especially in the Northern Equatorial Africa. The climate indicators considered in this study were unable to explain the alarming negative trend of PGR observed in the Sahelian region, even though the Standardized Precipitation‐Evapotranspiration Index (SPEI) values report a 15% drought stress. On the other hand, increases in temperature have not been detected in the Southern Hemisphere of Africa, where increasing frequency of storm occurrences determine a rise of PGR, particularly in southern Africa. Time analysis highlights a strong seasonality effect, while PGR is in‐phase with rainfall patterns in the summer (Northern Hemisphere) and winter (Southern Hemisphere) and out‐of‐phase during the fall season. This study helps to elucidate the mechanism of the processes influencing groundwater resources in six climatic zones of Africa, even though modelling results need to be validated more extensively with direct measurements in future studies. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of snowpack dynamics during the spring melt season for a sub‐alpine site using point measurements and numerical modeling

Snowpack dynamics through October 2014–June 2017 were described for a forested, sub‐alpine field site in southeastern Wyoming. Point measurements of wetness and density were combined with numerical modeling and continuous time series of snow depth, snow temperature, and snowpack outflow to identify 5 major classes of distinct snowpack conditions. Class (i) is characterized by no snowpack outflow and variable average snowpack temperature and density. Class (ii) is characterized by short durations of liquid water in the upper snowpack, snowpack outflow values of 0.0008–0.005 cm hr^?1, an increase in snowpack temperature, and average snow density between 0.25–0.35 g cm^?3. Class (iii) is characterized by a partially saturated wetness profile, snowpack outflow values of 0.005–0.25 cm hr^?1, snowpack temperature near 0 °C, and average snow density between 0.25–0.40 g cm^?3. Class (iv) is characterized by strong diurnal snowpack outflow pattern with values as high as 0.75 cm hr^?1, stable snowpack temperature near 0 °C, and stable average snow density between 0.35–0.45 g cm^?3. Class (v) occurs intermittently between Classes (ii)–(iv) and displays low snowpack outflow values between 0.0008–0.04 cm hr^?1, a slight decrease in temperature relative to the preceding class, and similar densities to the preceding class. Numerical modeling of snowpack properties with SNOWPACK using both the Storage Threshold scheme and Richards' equation was used to quantify the effect of snowpack capillarity on predictions of snowpack outflow and other snowpack properties. Results indicate that both simulations are able to predict snow depth, snow temperature, and snow density reasonably well with little difference between the 2 water transport schemes. Richards' equation more accurately simulates the timing of snowpack outflow over the Storage Threshold scheme, especially early in the melt season and at diurnal timescales.