Spatio‐temporal variability of Great Lakes basin snow cover ablation events

In the Great Lakes basin of North America, annual run‐off is dominated by snowmelt. This snowmelt‐induced run‐off plays an important role within the hydrologic cycle of the basin, influencing soil moisture availability and driving the seasonal cycle of spring and summer lake levels. Despite this, relatively little is understood about the patterns and trends of snow ablation event frequency and magnitude within the Great Lakes basin. This study uses a gridded dataset of Canadian and United States surface snow depth observations to develop a regional climatology of snow ablation events from 1960 to 2009. An ablation event is defined as an interdiurnal snow depth decrease within an individual grid cell. A clear seasonal cycle in ablation event frequency exists within the basin and peak ablation event probability is latitudinally dependent. Most of the basin experiences peak ablation frequency in March, while the northern and southern regions of the basin experience respective peaks in April and February. An investigation into the interannual frequency of ablation events reveals ablation events significantly decrease within the northeastern and northwestern Lake Superior drainage basins and significantly increase within the eastern Lake Huron and Georgian Bay drainage basins. In the eastern Lake Huron and Georgian Bay drainage basins, larger ablation events are occurring more frequently, and a larger impact to the hydrology can be expected. Trends in ablation events are attributed primarily to changes in snowfall and snow depth across the region.     

A spreadsheet‐based (Microsoft Excel) point surface energy balance model for glacier and snow melt studies

This paper describes a point surface energy balance model which runs within the Microsoft Excel spreadsheet package. The study incorporates a large amount of previous energy balance work and presents it in a useable form. The core model calculates the net shortwave and longwave radiation fluxes, the turbulent sensible and latent heat fluxes and the surface melt rate at a point on a melting ice or snow surface, from hourly inputs of incoming shortwave radiation, vapour pressure, air temperature and wind speed data. The latitude, longitude, slope angle, aspect, elevation, local temperature lapse rate, albedo and aerodynamic roughness of the study site, and the elevation of the meteorological station, can all be specified in the model. An output file containing the hourly and daily rates, and the totals of the energy fluxes is generated. The main advantages of the model are: first, that it requires only a PC or laptop computer running standard Microsoft Windows software, enabling it to be used at a desktop or in the field; and second, that it can be adapted quickly to different sites, meteorological data formats and other application requirements. Model calculations are compared with measured surface melt rates at five points on Haut Glacier d'Arolla, Switzerland, over a 115 day ablation period. Allowing for differences in shading between the meteorological station and the glacier, the root mean square error of the calculated melt rates is 2·0 mm day⁻¹ water equivalent melt (mean error +1·2 mm day⁻¹), for measured melt rates in the range 23 to 42 mm day⁻¹ water equivalent melt. Copyright © 2000 John Wiley & Sons, Ltd.     

Climate variability,snow, and physiographic controls on storm hydrographs in small forested basins,western Cascades,Oregon

Large floods are often attributed to the melting of snow during a rain event. This study tested how climate variability, snowpack presence, and basin physiography were related to storm hydrograph shape in three small (<1 km²) basins with old‐growth forest in western Oregon. Relationships between hydrograph characteristics and precipitation were tested for approximately 800 storms over a nearly 30‐year period. Analyses controlled for (1) snowpack presence/absence, (2) antecedent soil moisture, and (3) hillslope length and gradient. For small storms (<150 mm precipitation), controlling for precipitation, the presence of a snowpack on near‐saturated soil increased the threshold of precipitation before hydrograph rise, extended the start lag, centroid lag, and duration of storm hydrographs, and increased the peak discharge. The presence of a snowpack on near‐saturated soil sped up and steepened storm hydrographs in a basin with short steep slopes, but delayed storm hydrographs in basins with longer or more gentle slopes. Hydrographs of the largest events, which were extreme regional rain and rain‐on‐snow floods, were not sensitive to landform characteristics or snowpack presence/absence. Although the presence of a snowpack did not increase peak discharge in small, forested basins during large storms, it had contrasting effects on storm timing in small basins, potentially synchronizing small basin contributions to the larger basin hydrograph during large rain‐on‐snow events. By altering the relative timing of hydrographs, snowpack melting could produce extreme floods from precipitation events whose size is not extreme. Further work is needed to examine effects of canopy openings, snowpack, and climate warming on extreme rain‐on‐snow floods at the large basin scale. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling the spatial–temporal variability of spring snowmelt in an arctic catchment

Arctic spring landscapes are usually characterized by a mosaic of coexisting snow‐covered and bare ground patches. This phenomenon has major implications for hydrological processes, including meltwater production and runoff. Furthermore, as indicated by aircraft observations, it affects land‐surface–atmosphere exchanges, leading to a high degree of variability in surface energy terms during melt. The heterogeneity and related differences when certain parts of the landscape become snow free also affects the length of the growing season and the carbon cycle. Small‐scale variability in arctic snowmelt is addressed here by combining a spatially distributed end‐of‐winter snow cover with simulations of variable snowmelt energy balance factors for the small arctic catchment of Trail Valley Creek (63 km²). Throughout the winter, snow in arctic tundra basins is redistributed by frequent blowing snow events. Areas of above‐ or below‐average end‐of‐winter snow water equivalents were determined from land‐cover classifications, topography, land‐cover‐based snow surveys, and distributed surface wind‐field simulations. Topographic influences on major snowmelt energy balance factors (solar radiation and turbulent fluxes of sensible and latent heat) were modelled on a small‐scale (40 m) basis. A spatially variable complete snowmelt energy balance was subsequently computed and applied to the distributed snow cover, allowing the simulation of the progress of melt throughout the basin. The emerging patterns compared very well visually to snow cover observations from satellite images and aerial photographs. Results show the relative importance of variable end‐of‐winter snow cover, spatially distributed melt energy fluxes, and local advection processes for the development of a patchy snow cover. This illustrates that the consideration of these processes is crucial for an accurate determination of snow‐covered areas, as well as the location, timing, and amount of meltwater release from arctic catchments, and should, therefore, be included in hydrological models. Furthermore, the study shows the need for a subgrid parameterization of these factors in the land surface schemes of larger scale climate models. Copyright © 2005 John Wiley & Sons, Ltd.     

Estimating sublimation of intercepted and sub‐canopy snow using eddy covariance systems

Direct measurements of winter water loss due to sublimation were made in a sub‐alpine forest in the Rocky Mountains of Colorado. Above‐and below‐canopy eddy covariance systems indicated substantial losses of winter‐season snow accumulation in the form of snowpack (0·41 mm d^?1) and intercepted snow (0·71 mm d^?1) sublimation. The partitioning between these over and under story components of water loss was highly dependent on atmospheric conditions and near‐surface conditions at and below the snow/atmosphere interface. High above‐canopy sensible heat fluxes lead to strong temperature gradients between vegetation and the snow‐surface, driving substantial specific humidity gradients at the snow surface and high sublimation rates. Intercepted snowfall resulted in rapid response of above‐canopy latent heat fluxes, high within‐canopy sublimation rates (maximum = 3·7 mm d^?1), and diminished sub‐canopy snowpack sublimation. These results indicate that sublimation losses from the sub‐canopy snowpack are strongly dependent on the partitioning of sensible and latent heat fluxes in the canopy. This compels comprehensive studies of snow sublimation in forested regions that integrate sub‐canopy and over‐story processes. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact of an extreme melt event on the runoff and hydrology of a high Arctic glacier

On 28–30 July 2000, an extreme melt event was observed at John Evans Glacier (JEG), Ellesmere Island (79° 40′N, 74° 00′W). Hourly melt rates during this event fell in the upper 4% of the distribution of melt rates observed at the site during the period 1996–2000. Synoptic conditions during the event resulted in strong east‐to‐west flow over the northern sector of the Greenland Ice Sheet, with descending flow on the northwest side reaching Ellesmere Island. On JEG, wind speeds during the event averaged 8·1 m s^?1 at 1183 m a.s.l., with hourly mean wind speeds peaking at 11·6 m s^?1. Air temperatures reached 8°C, and rates of surface lowering measured by an ultrasonic depth gauge averaged 56 mm day^?1. Calculations with an energy balance model suggest that increased turbulent fluxes contributed to melt enhancement at all elevations on the glacier, while snow albedo feedback resulted in increased melting due to net radiation at higher elevations. The event was responsible for 30% of total summer melt at 1183 m a.s.l. and 15% at 850 m a.s.l. Conditions similar to those during the event occurred on only 0·1% of days in the period 1948–2000, but 61% of events occurred in the summer months and there was an apparent clustering of events in the 1950s and 1980s. Such events have the potential to impact significantly on runoff, mass balance and drainage system development at high Arctic glaciers, and changes in their incidence could play a role in determining how high Arctic glaciers respond to climate change and variability. Copyright © 2003 John Wiley & Sons, Ltd.     

Energy balance above a boreal coniferous forest: a difference in turbulent fluxes between snow‐covered and snow‐free canopies

To evaluate the interactive effects of snow and forest on turbulent fluxes between the forest surface and the atmosphere, the surface energy balance above a forest was measured by the eddy correlation method during the winter of 1995–1996. The forest was a young coniferous plantation comprised of spruce and fir. The study site, in Sapporo, northern Japan, had heavy and frequent snowfalls and the canopy was frequently covered with snow during the study period. A comparison of the observed energy balance above the forest for periods with and without a snow‐covered canopy and an analysis using a single‐source model gave the following results: during daytime when the canopy was covered with snow, the upward latent heat flux was large, about 80% of the net radiation, and the sensible heat flux was positive but small. On the other hand, during daytime when the canopy was dry and free from snow, the sensible heat flux was dominant and the latent heat flux was minor, about 10% of the net radiation. To explain this difference of energy partition between snow‐covered and snow‐free conditions, not only differences in temperature but also differences in the bulk transfer coefficients for latent heat flux were necessary in the model. Therefore, the high evaporation rate from the snow‐covered canopy can be attributed largely to the high moisture availability of the canopy surface. Evaporation from the forest during a 60‐day period in midwinter was estimated on a daily basis as net radiation minus sensible heat flux. The overall average evaporation during the 60‐day period was 0·6 mm day⁻¹, which is larger than that from open snow fields. Copyright © 1999 John Wiley & Sons, Ltd.     

An analysis of snow cover patterns in a small alpine catchment

Snow cover patterns in a 9.4 km² basin in the Austrian Alps are examined during spring and summer 1989. Digital mono-plotting from oblique aerophotographs is used for mapping. on the basis of a square grid with 25 m spacing, snow cover as mapped during nine surveys is analysed as a function of elevation and slope. During winter conditions the snow cover is found to be much better related to these terrain features than during the late ablation period.     

Assessing the significance of wet-canopy evaporation from forests during extreme rainfall events for flood mitigation in mountainous regions of the United Kingdom

There is increased interest in the potential of tree planting to help mitigate flooding using nature-based solutions or natural flood management. However, many publications based upon catchment studies conclude that, as flood magnitude increases, benefit from forest cover declines and is insignificant for extreme flood events. These conclusions conflict with estimates of evaporation loss from forest plot observations of gross rainfall, through fall and stem flow. This study explores data from existing studies to assess the magnitudes of evaporation and attempts to identify the meteorological conditions under which they would be supported. This is achieved using rainfall event data collated from publications and data archives from studies undertaken in temperate environments around the world. The meteorological conditions required to drive the observed evaporation losses are explored theoretically using the Penman–Monteith equation. The results of this theoretical analysis are compared with the prevailing meteorological conditions during large and extreme rainfall events in mountainous regions of the United Kingdom to assess the likely significance of wet canopy evaporation loss. The collated dataset showed that event Ewc losses between approximately 2 and 38% of gross rainfall (1.5 to 39.4 mm day⁻¹) have been observed during large rainfall events (up to 118 mm day⁻¹) and that there are few data for extreme events (>150 mm day⁻¹). Event data greater than 150 mm (reported separately) included similarly high percentage evaporation losses. Theoretical estimates of wet-canopy evaporation indicated that, to reproduce the losses towards the high end of these observations, relative humidity and the aerodynamic resistance for vapour transport needed to be lower than approximately 97.5% and 0.5 to 2 s m⁻¹ respectively. Surface meteorological data during large and extreme rainfall events in the United Kingdom suggest that conditions favourable for high wet-canopy evaporation are not uncommon and indicate that significant evaporation losses during large and extreme events are possible but not for all events and not at all locations. Thus the disparity with the results from catchment studies remains.     

Snow-distribution and melt modelling for glaciers in Zackenberg river drainage basin,north-eastern Greenland

A physically based snow-evolution modelling system (SnowModel) that includes four sub-models: MicroMet, EnBal, SnowPack, and SnowTran-3D, was used to simulate eight full-year evolutions of snow accumulation, distribution, sublimation, and surface melt from glaciers in the Zackenberg river drainage basin, in north-east Greenland. Meteorological observations from two meteorological stations were used as model inputs, and spatial snow depth observations, snow melt depletion curves from photographic time lapse, and a satellite image were used for model testing of snow and melt simulations, which differ from previous SnowModel tests methods used on Greenland glaciers. Modelled test-period-average end-of-winter snow water equivalent (SWE) depth for the depletion area differs by a maximum of 14 mm w.eq., or ∼6%, more than the observed, and modelled test-period-average snow cover extent differs by a maximum of 5%, or 0·8 km², less than the observed. Furthermore, comparison with a satellite image indicated a 7% discrepancy between observed and modelled snow cover extent for the entire drainage basin. About 18% (31 mm w.eq.) of the solid precipitation was returned to the atmosphere by sublimation. Modelled mean annual snow melt and glacier ice melt for the glaciers in the Zackenberg river drainage basin from 1997 through 2005 (September–August) averaged 207 mm w.eq. year⁻¹ and 1198 mm w.eq. year⁻¹, respectively, yielding a total averaging 1405 mm w.eq. year⁻¹. Total modelled mean annual surface melt varied from 960 mm w.eq. year⁻¹ to 1989 mm w.eq. year⁻¹. The surface-melt period started between mid-May and the beginning of June and lasted until mid-September. Annual calculated runoff averaged 1487 mm w.eq. year⁻¹ (∼150 × 10⁶ m³) (1997–2005) with variations from 1031 mm w.eq. year⁻¹ to 2051 mm w.eq. year⁻¹. The model simulated a total glacier recession averaging − 1347 mm w.eq. year⁻¹ (∼136 × 10⁶ m³) (1997–2005), which was almost equal to previous basin average hydrological water balance storage studies − 244 mm w.eq. year⁻¹ (∼125 × 10⁶ m³) (1997–2003). Copyright © 2007 John Wiley & Sons, Ltd.     

Distributed snow accumulation and ablation modeling in the American River basin

A spatially distributed snow model procedure for estimating snow melt, snow water equivalent and snow cover area is formulated and tested with data from the American River basin in California’s Sierra Nevada. An adaptation of the operational National Weather Service snow accumulation and ablation model is used for each model grid cell forced by spatially distributed precipitation and temperature data. The model was implemented with 6-hourly time steps on 1 km² grid cells for the snow season of 1999–2003. Temperature is spatially interpolated using the prevailing lapse rate and digital terrain elevation data. Precipitation is spatially interpolated using regional climatological analyses obtained from PRISM. Parameters that control snow melt are distributed using ground surface aspect. The model simulations are compared with data from 12 snow-sensors located in the basin and the daily 500-m snow cover extent product from the MODIS/Terra satellite mission. The results show that the distribution of snow pack over the area is generally captured. The snow pack quantity compared to snow gauges is well estimated in high elevations with increasing uncertainty in the snow pack at lower elevations. Sensitivity and uncertainty analyses indicate that the significant input uncertainty for precipitation and temperature is primarily responsible for model errors in lower elevations and near the snow line. The model is suitable for producing spatially resolved realistic snow pack simulations when forced with operationally available observed or predicted data.     

Effect of ENSO events on sediment production in a large coastal basin in northern Peru

Although the importance of ENSO on hydrological anomalies has been recognized, variations in sediment fluxes caused by these extreme events are poorly documented. The effect of ENSO is not limited to changes in sediment mobilization. Since ENSO events can affect terrestrial ecosystems, they may have important effects on sediment production and transport in river basins over time spans that are longer than the duration of the event itself. The Catamayo‐Chira basin is an interesting casestudy for investigating these geomorphic implications. The objectives were: (i) to study the effect of ENSO on stream flow and sediment yields in the basin, (ii) to investigate if ENSO events affect sediment yields in the post‐ENSO period and (iii) to understand which factors control the ENSO and post‐ENSO basin response. During strong negative ENSO periods, mean annual stream flow discharge at the inlet of the Poechos reservoir in the lower basin was 5.4 times higher than normal annual discharges, while average sediment fluxes exceeded those of normal years by a factor of about 11. In two heavily affected periods, 45.9% of the total sediment yield in the 29 years observation period was generated. Sediment fluxes in the post‐ENSO period are lower than expected, which proves post‐ENSO event dynamics are significantly different from pre‐event dynamics. Our analysis indicates the increase of vegetation growth in the lower basin is not the main reason explaining considerable sediment flux decrease in post‐ENSO periods. During strong ENSO events, sediment in alluvial stores in the lower part of the basin is removed due to enlarging and deepening of channels. In post‐ENSO periods, normal discharges and persisting sediment supplies from the middle/upper basin lead to river aggradation and storage of large amounts of sediment in alluvial plains. The decrease in sediment export will last for several years until the equilibrium is re‐established. Copyright © 2011 John Wiley & Sons, Ltd.     

Short‐term discharge and suspended sediment fluctuations in the proglacial Skeldal River,north‐east Greenland

The summer discharge pattern of the Skeldal River, which drains a 560 km² partly glacierized catchment in north‐east Greenland, is dominated by diurnal oscillations reflecting variations in the melt rate of snow and ice in the basin. Superimposed on this diurnal pattern are numerous short‐lived discharge fluctuations of irregular periodicity and magnitude. The larger fluctuations are described and attributed to both rainfall events and periodic collapse of the glacier margin damming flow from beneath the Skelbrae glacier. Other minor fluctuations are less readily explained but are associated with changes in the channelized and distributed reservoirs and possibly temporary blockage of subglacial conduits caused by ice melt with subsequent damming. Fluctuations in suspended sediment concentration (SSC) are normally associated with discharge fluctuations, although examples of ‘transient flushes’ were observed where marked increases in SSC occurred in the absence of corresponding discharge variations. A strong relationship between the event discharge increase and event SSC increase for rainfall‐induced events was established, but no such relationship existed for non‐rainfall‐induced events. There is some evidence for an exhaustion effect in the SSC patterns both at the event time‐scale and as the month proceeds. A mean suspended sediment load of 1765 ± 0·26 t day^?1 was estimated for the study period, which would be equivalent to a suspended sediment yield of 732 ± 4 t km^?2 year^?1. Copyright © 2001 John Wiley & Sons, Ltd.     

Lithological and hydrological controls on water composition: evaporite dissolution and glacial weathering in the south central Andes of Argentina (33°–34°S)

Lithological and hydrological influence on fluvial physical and chemical erosion was studied in a glacierized sedimentary basin with high evaporite presence. Suspended particulate matter (SPM), total dissolved solids (TDS) and major ion concentrations were analysed for 2 years of different hydrologic condition: (i) 2009–2010, Q = 100% average; and (ii) 2010–2011, Q = 60% average. Annual hydrograph was simple regime‐type with one peak in summer related to snow melting. The intra‐annual SPM and TDS variations were directly and inversely associated to Q, respectively. Snow chemistry showed continental influence (Na⁺/Ca²⁺ = 0.17), and atmospheric input of TDS was <1% of the total exported flux. River water was highly concentrated in Ca²⁺ and SO₄²⁻ (~4 mmol l⁻¹) and in Na⁺ and Cl⁻ (~3 mmol l⁻¹). Ca²⁺/SO₄²⁻ and Na⁺/Cl⁻ molar ratios were ~1 and related to Q, directly and inversely, respectively. Major ion relationships suggest that river chemistry is controlled by evaporite (gypsum and halite) dissolution having a summer input from sulfide oxidation and carbonate dissolution, and a winter input from subsurface flow loaded with silicate weathering products. This variation pattern resulted in nearly chemostatic behaviour for Ca⁺, Mg²⁺ and SO₄²⁻, whereas Na⁺, Cl⁻ and SiO₂ concentrations showed to be controlled by dilution/concentration processes. During the 2009–2010 hydrological year, the fluxes of water, SPM and TDS registered in the snow melting–high Q season were, respectively, 71%, 92% and 67% of the annual total, whereas for equal period in 2010–2011, 56% of water, 86% of SPM and 54% of TDS annual fluxes were registered. The SPM fluxes for 2009–2010 and 2010–2011 were 1.19 × 10⁶ and 0.79 × 10⁶ t year⁻¹, whereas TDS fluxes were 0.68 × 10⁶ and 0.55 × 10⁶ t year⁻¹, respectively. Export rates for 2009–2010 were 484 t km² year⁻¹ for SPM and 275 t km² year⁻¹ for TDS. These rates are higher than those observed in glacierized granite basins and in non‐glacierized evaporite basins, suggesting a synergistic effect of lithology and glaciers on physical and chemical erosion. Copyright © 2014 John Wiley & Sons, Ltd.     

Soil temperature and soil moisture dynamics in winter and spring under heavy snowfall conditions in North-Eastern Japan

Warm winters and high precipitation in north-eastern Japan generate snow covers of more than three meters depth and densities of up to 0.55 g cm⁻³. Under these conditions, rain/snow ratio and snowmelt have increased significantly in the last decade under increasing warm winters. This study aims at understanding the effect of rain-on-snow and snowmelt on soil moisture under thick snow covers in mid-winter, taking into account that snowmelt in spring is an important source of water for forests and agriculture. The study combines three components of the Hydrosphere (precipitation, snow cover and soil moisture) in order to trace water mobility in winter, since soil temperatures remained positive in winter at nearly 0.3°C. The results showed that soil moisture increased after snowmelt and especially after rain-on-snow events in mid-winter 2018/2019. Rain-on-snow events were firstly buffered by fresh snow, increasing the snow water equivalent (SWE), followed by water soil infiltration once the water storage capacity of the snowpack was reached. The largest increase of soil moisture was 2.35 vol%. Early snowmelt increased soil moisture with rates between 0.02 and 0.035 vol% hr⁻¹ while, rain-on-snow events infiltrated snow and soil faster than snowmelt and resulted in rates of up to 1.06 vol% hr⁻¹. These results showed the strong connection of rain, snow and soil in winter and introduce possible hydrological scenarios in the forest ecosystems of the heavy snowfall regions of north-eastern Japan. Effects of rain-on-snow events and snowmelt on soil moisture were estimated for the period 2012–2018. Rain/snow ratio showed that only 30% of the total precipitation in the winter season 2011/2012 was rain events while it was 50% for the winter 2018/2019. Increasing climate warming and weakening of the Siberian winter monsoons will probably increase rain/snow ratio and the number of rain-on-snow events in the near future.     

MELTWATER FLUXES AT AN ARCTIC FOREST-TUNDRA SITE

Models of surface energy balance and snow metamorphism are utilized to predict the energy and meltwater fluxes at an Arctic site in the forest–tundra transition zone of north-western Canada. The surface energy balance during the melt period is modelled using an hourly bulk aerodynamic approach. Once a snowcover becomes patchy, advection from the bare patches to the snow-covered areas results in a large spatial variation in basin snowmelt. In order to illustrate the importance of small-scale, horizontal advection, a simple parameterization scheme using sensible heat fluxes from snow free areas was tested. This scheme estimates the maximum horizontal advection of sensible heat from the bare patches to the snow-covered areas. Calculated melt was routed through the measured snowcover in each landscape type using a variable flow path, meltwater percolation model. This allowed the determination of the spatial variability in the timing and magnitude of meltwater release for runoff. Model results indicate that the initial release of meltwater first occurred on the shallow upland tundra sites, but meltwater release did not occur until nearly two weeks later on the deep drift snowcovers. During these early periods of melt, not all meltwater is available for runoff. Instead, there is a period when some snowpacks are only partially contributing to runoff, and the spatial variation of runoff contribution corresponds to landscape type. Comparisons of melt with and without advection suggests that advection is an important process controlling the timing of basin snowmelt.     

Quantifying the role of the snowpack in generating water available for run‐off during rain‐on‐snow events from snow pillow records

Rain‐on‐snow events have generated major floods around the world, particularly in coastal, mountainous regions. Most previous studies focused on a limited number of major rain‐on‐snow events or were based primarily on model results, largely due to a lack of long‐term records from lysimeters or other instrumentation for quantifying event water balances. In this analysis, we used records from five automated snow pillow sites in south coastal British Columbia, Canada, to reconstruct event water balances for 286 rain‐on‐snow events over a 10‐year period. For large rain‐on‐snow events (event rainfall >40 mm), snowmelt enhanced the production of water available for run‐off (WAR) by approximately 25% over rainfall alone. For smaller events, a range of antecedent and meteorological factors influenced WAR generation, particularly the antecedent liquid water content of the snowpack. Most large events were associated with atmospheric rivers. Rainfall dominated WAR generation during autumn and winter events, whereas snowmelt dominated during spring and summer events. In the majority of events, the sensible heat of rain contributed less than 10% of the total energy consumed by snowmelt. This analysis illustrated the importance of understanding the amount of rainfall occurring at high elevations during rain‐on‐snow events in mountainous regions.     

Assessing the impact of distributed snow water equivalent calibration and assimilation of Copernicus snow water equivalent on modelled snow and streamflow performance

Accuracy of the Copernicus snow water equivalent (SWE) product and the impact of SWE calibration and assimilation on modelled SWE and streamflow was evaluated. Daily snowpack measurements were made at 12 locations from 2016 to 2019 across a 4104 km² mixed-forest basin in the Great Lakes region of central Ontario, Canada. Sub-basin daily SWE calculated from these sites, observed discharge, and lake levels were used to calibrate a hydrologic model developed using the Raven modelling framework. Copernicus SWE was bias corrected during the melt period using mean bias subtraction and was compared to daily basin average SWE calculated from the measured data. Bias corrected Copernicus SWE was assimilated into the models using a range of parameters and the parameterizations from the model calibration. The bias corrected Copernicus product agreed well with measured data and provided a good estimate of mean basin SWE demonstrating that the product shows promise for hydrology applications within the study region. Calibration to spatially distributed SWE substantially improved the basin scale SWE estimate while only slightly degrading the flow simulation demonstrating the value of including SWE in a multi-objective calibration formulation. The particle filter experiments yielded the best SWE estimation but moderately degraded the flow simulation. The particle filter experiments constrained by the calibrated snow parameters produced similar results to the experiments using the upper and lower bounds indicating that, in this study, model calibration prior to assimilation was not valuable. The calibrated models exhibited varying levels of skill in estimating SWE but demonstrated similar streamflow performance. This indicates that basin outlet streamflow can be accurately estimated using a model with a poor representation of distributed SWE. This may be sufficient for applications where estimating flow is the primary water management objective. However, in applications where understanding the physical processes of snow accumulation, melt and streamflow generation are important, such as assessing the impact of climate change on water resources, accurate representations of SWE are required and can be improved via multi-objective calibration or data assimilation, as demonstrated in this study.     

Securing the future of ground water resources in the Great Lakes basin

Ground water is a vital, but underappreciated, natural resource in the Great Lakes basin. It meets many human needs and contributes significantly to the hydrology of the Great Lakes and the health of ecosystems. This paper provides an overview of ground water in the Great Lakes and the institutional and legal setting that governs the use, protection, diversion, and removal of water from the basin and proposes a citizen-centered vision for management of ground water in the 21st century.