Freezing precipitation has major consequences for ground and air transportation, the health of citizens, and power networks. Previous studies using coarse resolution climate models have shown a northward migration of freezing rain in the future. Increased model resolution can better define local topography leading to improved representation of conditions that are favorable for freezing rain. The goal of this study is to examine the climatology and characteristics of future freezing rain events using very-high resolution climate simulations. Historical and pseudo-global warming simulations with a 4-km horizontal grid length were used and compared with available observations. Simulations revealed a northerly shift of freezing rain occurrence, and an increase in the winter. Freezing rain was still shown to occur in the Saint-Lawrence River Valley in a warmer climate, primarily due to stronger wind channeling. Up to 50% of the future freezing rain events also occurred in present day climate within 12 h of each other. In northern Maine, they are typically shorter than 6 h in current climate and longer than 6 h in warmer conditions due to the onset of precipitation during low-pressure systems occurrences. The occurrence of freezing rain also locally increases slightly north of Québec City in a warmer climate because of freezing rain that is produced by warm rain processes. Overall, the study shows that high-resolution regional climate simulations are needed to study freezing rain events in warmer climate conditions, because high horizontal resolutions better define small-scale topographic features and local physical mechanisms that have an influence on these events.

Freezing precipitation has major consequences for ground and air transportation, the health of citizens, and power networks. Previous studies using coarse resolution climate models have shown a northward migration of freezing rain in the future. Increased model resolution can better define local topography leading to improved representation of conditions that are favorable for freezing rain. The goal of this study is to examine the climatology and characteristics of future freezing rain events using very-high resolution climate simulations. Historical and pseudo-global warming simulations with a 4-km horizontal grid length were used and compared with available observations. Simulations revealed a northerly shift of freezing rain occurrence, and an increase in the winter. Freezing rain was still shown to occur in the Saint-Lawrence River Valley in a warmer climate, primarily due to stronger wind channeling. Up to 50% of the future freezing rain events also occurred in present day climate within 12 h of each other. In northern Maine, they are typically shorter than 6 h in current climate and longer than 6 h in warmer conditions due to the onset of precipitation during low-pressure systems occurrences. The occurrence of freezing rain also locally increases slightly north of Québec City in a warmer climate because of freezing rain that is produced by warm rain processes. Overall, the study shows that high-resolution regional climate simulations are needed to study freezing rain events in warmer climate conditions, because high horizontal resolutions better define small-scale topographic features and local physical mechanisms that have an influence on these events.

Abstract Temperatures near 0°C represent a critical threshold for many environmental processes and socio‐economic activities. This study examines surface air temperatures ( T ) near 0°C (−2°C ≤ T ≤ 2°C) across much of southern Canada over a 13 year period (October 2000–September 2013). It utilized hourly data from 39 weather stations and from 4‐km resolution Weather Research and Forecasting model simulations that were both a retrospective simulation as well as a pseudo‐global warming simulation applicable near the end of the 21st century. Average annual occurrences of near‐0°C conditions increase by a relatively small amount of 5.1% from 985 hr in the current climate to 1,035 hr within the future one. Near‐0°C occurrences with precipitation vary from <5% to approximately 50% of these values. Near‐0°C occurrences are sometimes higher than values of neighboring temperatures. These near‐0°C peaks in temperature distributions can occur in both the current and future climate, in only one, or in neither. Only 4.3% of southern Canada is not associated with a near‐0°C peak and 65.8% is associated with a near‐0°C peak in both climates. It is inferred that latent heat exchanges from the melting and freezing of, for example, precipitation and the snowpack contribute significantly to some of these findings.

Solid precipitation falling near 0 °C, mainly snow, can adhere to surface features and produce major impacts. This study is concerned with characterizing this precipitation over the Canadian Prairie provinces of Manitoba and Saskatchewan in the current (2000–2013) and pseudo-global warming future climate, with an average 5.9 °C temperature increase, through the use of high resolution (4 km) model simulations. On average, simulations in the current climate suggest that this precipitation occurs within 11 events per year, lasting 33.6 h in total and producing 27.5 mm melted equivalent, but there are wide spatial variations that are partly due to enhancements arising from its relatively low terrain. Within the warmer climate, average values generally increase, and spatial patterns shift somewhat. This precipitation consists of four categories covering its occurrence just below and just above a wet-bulb temperature of 0 °C, and with or without liquid precipitation. It generally peaks in March or April, as well as in October, and these peaks move towards mid-winter by approximately one month within the warmer climate. Storms producing this precipitation generally produce winds with a northerly component during or shortly after the precipitation; these winds contribute to further damage. Overall, this study has determined the features of and expected changes to adhering precipitation across this region.

In mountains, the precipitation phase greatly varies in space and time and affects the evolution of the snow cover. Snowpack models usually rely on precipitation-phase partitioning methods (PPMs) that use near-surface variables. These PPMs ignore conditions above the surface thus limiting their ability to predict the precipitation phase at the surface. In this study, the impact on snowpack simulations of atmospheric-based PPMs, incorporating upper atmospheric information, is tested using the snowpack scheme Crocus. Crocus is run at 2.5-km grid spacing over the mountains of southwestern Canada and northwestern United States and is driven by meteorological fields from an atmospheric model at the same resolution. Two atmospheric-based PPMs were considered from the atmospheric model: the output from a detailed microphysics scheme and a post-processing algorithm determining the snow level and the associated precipitation phase. Two ground-based PPMs were also included as lower and upper benchmarks: a single air temperature threshold at 0°C and a PPM using wet-bulb temperature. Compared to the upper benchmark, the snow-level based PPM improved the estimation of snowfall occurrence by 5% and the simulation of snow water equivalent (SWE) by 9% during the snow melting season. In contrast, due to missing processes, the microphysics scheme decreased performances in phase estimate and SWE simulations compared to the upper benchmark. These results highlight the need for detailed evaluation of the precipitation phase from atmospheric models and the benefit for mountain snow hydrology of the post-processed snow level. The limitations to drive snowpack models at slope scale are also discussed.

In mountains, the precipitation phase greatly varies in space and time and affects the evolution of the snow cover. Snowpack models usually rely on precipitation-phase partitioning methods (PPMs) that use near-surface variables. These PPMs ignore conditions above the surface thus limiting their ability to predict the precipitation phase at the surface. In this study, the impact on snowpack simulations of atmospheric-based PPMs, incorporating upper atmospheric information, is tested using the snowpack scheme Crocus. Crocus is run at 2.5-km grid spacing over the mountains of southwestern Canada and northwestern United States and is driven by meteorological fields from an atmospheric model at the same resolution. Two atmospheric-based PPMs were considered from the atmospheric model: the output from a detailed microphysics scheme and a post-processing algorithm determining the snow level and the associated precipitation phase. Two ground-based PPMs were also included as lower and upper benchmarks: a single air temperature threshold at 0°C and a PPM using wet-bulb temperature. Compared to the upper benchmark, the snow-level based PPM improved the estimation of snowfall occurrence by 5% and the simulation of snow water equivalent (SWE) by 9% during the snow melting season. In contrast, due to missing processes, the microphysics scheme decreased performances in phase estimate and SWE simulations compared to the upper benchmark. These results highlight the need for detailed evaluation of the precipitation phase from atmospheric models and the benefit for mountain snow hydrology of the post-processed snow level. The limitations to drive snowpack models at slope scale are also discussed.

Abstract. The amount and phase of cold season precipitation accumulating in the upper Saint John River basin are critical factors in determining spring runoff, ice-jams, and flooding in downstream communities. To study the impact of winter and spring storms on the snowpack in the upper Saint John River (SJR) basin, the Saint John River Experiment on Cold Season Storms (SAJESS) utilized meteorological instrumentation, upper air soundings, human observations, and hydrometeor macrophotography during winter/spring 2020–21. Here, we provide an overview of the SAJESS study area, field campaign, and existing data networks surrounding the upper SJR basin. Initially, meteorological instrumentation was co-located with an Environment and Climate Change Canada station near Edmundston, New Brunswick, in early December 2020. This was followed by an intensive observation period that involved manual observations, upper-air soundings, a multi-angle snowflake camera, macrophotography of solid hydrometeors, and advanced automated instrumentation throughout March and April 2021. The resulting datasets include optical disdrometer size and velocity distributions of hydrometeors, micro rain radar output, near-surface meteorological observations, and wind speed, temperature, pressure and precipitation amounts from a K63 Hotplate precipitation gauge, the first one operating in Canada. These data are publicly available from the Federated Research Data Repository at https://doi.org/10.20383/103.0591 (Thompson et al., 2022). We also include a synopsis of the data management plan and data processing, and a brief assessment of the rewards and challenges of utilizing community volunteers for hydro-meteorological citizen science.

The southern Canadian Rockies are prone to extreme precipitation that often leads to high streamflow, deep snowpacks, and avalanche risks. Many of these precipitation events are associated with rain–snow transitions, which are highly variable in time and space due to the complex topography. A warming climate will certainly affect these extremes and the associated rain–snow transitions. The goal of this study is to investigate the characteristics and variability of rain–snow transitions aloft and how they will change in the future. Weather Research and Forecasting (WRF) simulations were conducted from 2000 to 2013 and these were repeated in a warmer pseudo-global warming (PGW) future. Rain–snow transitions occurred aloft throughout the year over the southern Canadian Rockies, but their elevations and depths were highly variable, especially across the continental divide. In PGW conditions, with future air temperatures up to 4–5°C higher on average over the Canadian Rockies, rain–snow transitions are projected to occur more often throughout the year, except during summer. The near-0°C conditions associated with rain–snow transitions are expected to increase in elevation by more than 500 m, resulting in more rain reaching the surface. Overall, this study illustrates the variability of rain–snow transitions, which often impact the location of the snowline. This study also demonstrates the non-uniform changes under PGW conditions, due in part to differences in the types of weather patterns that generate rain–snow transitions across the region.

ABSTRACT Winter storms in eastern Canada can bring heavy precipitation, including large amounts of freezing rain. The resulting ice accumulation on structures such as trees and power lines can lead to widespread power outages and damage to infrastructure. The objective of this study is to provide a better understanding of the processes that led to extreme freezing rain events over New Brunswick (NB), Canada, during past events and how they may change in the future. To accomplish this, freezing rain events that affected the power network over NB were identified and analysed using high-resolution convection-permitting simulations. These simulations were produced from 2000 to 2013 climate data and using the pseudo global warming (WRF-PGW) approach, assuming warmer climate conditions. Our results show that through the process of cold air damming, the Appalachians enhance the development of strong temperature inversions, leading to an increase in the amount of freezing rain in central and southern NB. The occurrence of freezing rain events generally decreases by 40% in southern and eastern NB, while the occurrence of long-duration events (>6 h) increases slightly in northwestern NB in the WRF-PGW simulation. Overall, key local orographic effects that influence atmospheric conditions favorable for freezing precipitation were identified. This knowledge will enable us to better anticipate the impact of climate change on similar storms.

ABSTRACT A devastating storm struck southern Manitoba, Canada on 10–13 October 2019, producing a large region of mainly sticky and wet snow. Accumulations reached 75 cm, wind gusts exceeded 100 km h−1, and surface temperature (T) remained near 0°C (−1°C ≤ T ≤ 1°C) for up to 88 h. It produced the largest October snowfall and was the earliest to produce at least 20 cm since 1872 in Winnipeg. These factors led to unparalleled damage and power restoration challenges for Manitoba Hydro and, with leaves still largely on vegetation, the most damaging storm to Winnipeg’s trees ever recorded. The storm’s track was uncommon, and produced elevated convection related to buoyancy-driven instability and conditional symmetric instability (CSI), with a moist absolutely unstable layer (MAUL) near 500 hPa. Instabilities were released via lift through lower-tropospheric warm advection and frontogenesis, differential cyclonic vorticity advection, and jet streak dynamics. Precipitation bands, elevated convection, and lake effect snow bands enhanced local snowfall. Snow adhering to structures was not always wet but, when present, it sometimes occurred because of incomplete freezing of particles partially melted aloft in a near-surface (<100 m deep) inversion. Although other storms over the historical record have produced a similar combination of severe precipitation, temperature and wind conditions, none have done this for such a long period.

Abstract Freezing precipitation has major consequences for ground and air transportation, the health of citizens, and power networks. Previous studies using coarse resolution climate models have shown a northward migration of freezing rain in the future. Increased model resolution can better define local topography leading to improved representation of conditions that are favorable for freezing rain. The goal of this study is to examine the climatology and characteristics of future freezing rain events using very-high resolution climate simulations. Historical and pseudo-global warming simulations with a 4-km horizontal grid length were used and compared with available observations. Simulations revealed a northerly shift of freezing rain occurrence, and an increase in the winter. Freezing rain was still shown to occur in the Saint-Lawrence River Valley in a warmer climate, primarily due to stronger wind channeling. Up to 50% of the future freezing rain events also occurred in present day climate within 12 h of each other. In northern Maine, they are typically shorter than 6 h in current climate and longer than 6 h in warmer conditions due to the onset of precipitation during low-pressure systems occurrences. The occurrence of freezing rain also locally increases slightly north of Québec City in a warmer climate because of freezing rain that is produced by warm rain processes. Overall, the study shows that high-resolution regional climate simulations are needed to study freezing rain events in warmer climate conditions, because high horizontal resolutions better define small-scale topographic features and local physical mechanisms that have an influence on these events.

Abstract Freezing precipitation has major consequences for ground and air transportation, the health of citizens, and power networks. Previous studies using coarse resolution climate models have shown a northward migration of freezing rain in the future. Increased model resolution can better define local topography leading to improved representation of conditions that are favorable for freezing rain. The goal of this study is to examine the climatology and characteristics of future freezing rain events using very-high resolution climate simulations. Historical and pseudo-global warming simulations with a 4-km horizontal grid length were used and compared with available observations. Simulations revealed a northerly shift of freezing rain occurrence, and an increase in the winter. Freezing rain was still shown to occur in the Saint-Lawrence River Valley in a warmer climate, primarily due to stronger wind channeling. Up to 50% of the future freezing rain events also occurred in present day climate within 12 h of each other. In northern Maine, they are typically shorter than 6 h in current climate and longer than 6 h in warmer conditions due to the onset of precipitation during low-pressure systems occurrences. The occurrence of freezing rain also locally increases slightly north of Québec City in a warmer climate because of freezing rain that is produced by warm rain processes. Overall, the study shows that high-resolution regional climate simulations are needed to study freezing rain events in warmer climate conditions, because high horizontal resolutions better define small-scale topographic features and local physical mechanisms that have an influence on these events.

Abstract Winter precipitation is the source of many inconveniences in many regions of North America, for both infrastructure and the economy. The ice storm that hit the Canadian Maritime Provinces on 24–26 January 2017 remains one of the most expensive in history for the province of New Brunswick. Up to 50 mm of freezing rain caused power outages across the province, depriving up to one-third of New Brunswick residences of electricity, with some outages lasting 2 weeks. This study aims to use high-resolution atmospheric modeling to investigate the meteorological conditions during this severe storm and their contribution to major power outages. The persistence of a deep warm layer aloft, coupled with the slow movement of the associated low pressure system, contributed to widespread ice accumulation. When combined with the strong winds observed, extensive damage to electricity networks was inevitable. A 2-m temperature cold bias was identified between the simulation and the observations, in particular during periods of freezing rain. In the northern part of New Brunswick, cold-air advection helped keep temperatures below 0°C, while in southern regions, the 2-m temperature increased rapidly to slightly above 0°C because of radiational heating. The knowledge gained in this study on the processes associated with either maintaining or stopping freezing rain will enhance the ability to forecast and, in turn, to mitigate the hazards associated with those extreme events. Significance Statement A slow-moving low pressure system produced up to 50 mm of freezing rain for 31 h along the east coast of New Brunswick, Canada, on 24–26 January 2017, causing unprecedented power outages. Warm-air advection aloft, along with a combination of higher wind speeds and large amounts of ice accumulation, created ideal conditions for severe freezing rain. The storm began with freezing rain along the entire north–south cross section of eastern New Brunswick and changed to rain only in the south, when local temperatures increased to >0°C. Near-surface cold-air advection kept temperatures below 0°C in the north. Warming from the latent heat produced by freezing contributed to persistent near-0°C conditions during freezing rain.

Abstract Given their potentially severe impacts, understanding how freezing rain events may change as the climate changes is of great importance to stakeholders including electrical utility companies and local governments. Identification of freezing rain in climate models requires the use of precipitation-type algorithms, and differences between algorithms may lead to differences in the types of precipitation identified for a given thermodynamic profile. We explore the uncertainty associated with algorithm selection by applying four algorithms (Cantin and Bachand, Baldwin, Ramer, and Bourgouin) offline to an ensemble of simulations of the fifth-generation Canadian Regional Climate Model (CRCM5) at 0.22° grid spacing. First, we examine results for the CRCM5 driven by ERA-Interim reanalysis to analyze how well the algorithms reproduce the recent climatology of freezing rain and how results vary depending on algorithm parameters and the characteristics of available model output. We find that while the Ramer and Baldwin algorithms tend to be better correlated with observations than Cantin and Bachand or Bourgouin, their results are highly sensitive to algorithm parameters and to the number of pressure levels used. We also apply the algorithms to four CRCM5 simulations driven by different global climate models (GCMs) and find that the uncertainty associated with algorithm selection is generally similar to or greater than that associated with choice of driving GCM for the recent past climate. Our results provide guidance for future studies on freezing rain in climate simulations and demonstrate the importance of accounting for uncertainty between algorithms when identifying precipitation type from climate model output. Significance Statement Freezing rain events and ice storms can have major consequences, including power outages and dangerous road conditions. It is therefore important to understand how climate change might affect the frequency and severity of these events. One source of uncertainty in climate studies of these events is related to the choice of algorithm used to detect freezing rain in model output. We compare the frequency of freezing rain identified using four different algorithms and find sometimes large differences depending on the algorithm chosen over some regions. Our findings highlight the importance of taking this source of uncertainty into account and will provide researchers with guidance as to which algorithms are best suited for climate studies of freezing rain.

Abstract The Canadian Rockies are a triple-continental divide, whose high mountains are drained by major snow-fed and rain-fed rivers flowing to the Pacific, Atlantic, and Arctic Oceans. The objective of the April–June 2019 Storms and Precipitation Across the continental Divide Experiment (SPADE) was to determine the atmospheric processes producing precipitation on the eastern and western sides of the Canadian Rockies during springtime, a period when upslope events of variable phase dominate precipitation on the eastern slopes. To do so, three observing sites across the divide were instrumented with advanced meteorological sensors. During the 13 observed events, the western side recorded only 25% of the eastern side’s precipitation accumulation, rainfall occurred rather than snowfall, and skies were mainly clear. Moisture sources and amounts varied markedly between events. An atmospheric river landfall in California led to moisture flowing persistently northward and producing the longest duration of precipitation on both sides of the divide. Moisture from the continental interior always produced precipitation on the eastern side but only in specific conditions on the western side. Mainly slow-falling ice crystals, sometimes rimed, formed at higher elevations on the eastern side (>3 km MSL), were lifted, and subsequently drifted westward over the divide during nonconvective storms to produce rain at the surface on the western side. Overall, precipitation generally crossed the divide in the Canadian Rockies during specific spring-storm atmospheric conditions although amounts at the surface varied with elevation, condensate type, and local and large-scale flow fields.

In mountains, the precipitation phase greatly varies in space and time and affects the evolution of the snow cover. Snowpack models usually rely on precipitation-phase partitioning methods (PPMs) that use near-surface variables. These PPMs ignore conditions above the surface thus limiting their ability to predict the precipitation phase at the surface. In this study, the impact on snowpack simulations of atmospheric-based PPMs, incorporating upper atmospheric information, is tested using the snowpack scheme Crocus. Crocus is run at 2.5-km grid spacing over the mountains of southwestern Canada and northwestern United States and is driven by meteorological fields from an atmospheric model at the same resolution. Two atmospheric-based PPMs were considered from the atmospheric model: the output from a detailed microphysics scheme and a post-processing algorithm determining the snow level and the associated precipitation phase. Two ground-based PPMs were also included as lower and upper benchmarks: a single air temperature threshold at 0°C and a PPM using wet-bulb temperature. Compared to the upper benchmark, the snow-level based PPM improved the estimation of snowfall occurrence by 5% and the simulation of snow water equivalent (SWE) by 9% during the snow melting season. In contrast, due to missing processes, the microphysics scheme decreased performances in phase estimate and SWE simulations compared to the upper benchmark. These results highlight the need for detailed evaluation of the precipitation phase from atmospheric models and the benefit for mountain snow hydrology of the post-processed snow level. The limitations to drive snowpack models at slope scale are also discussed.

In mountains, the precipitation phase greatly varies in space and time and affects the evolution of the snow cover. Snowpack models usually rely on precipitation-phase partitioning methods (PPMs) that use near-surface variables. These PPMs ignore conditions above the surface thus limiting their ability to predict the precipitation phase at the surface. In this study, the impact on snowpack simulations of atmospheric-based PPMs, incorporating upper atmospheric information, is tested using the snowpack scheme Crocus. Crocus is run at 2.5-km grid spacing over the mountains of southwestern Canada and northwestern United States and is driven by meteorological fields from an atmospheric model at the same resolution. Two atmospheric-based PPMs were considered from the atmospheric model: the output from a detailed microphysics scheme and a post-processing algorithm determining the snow level and the associated precipitation phase. Two ground-based PPMs were also included as lower and upper benchmarks: a single air temperature threshold at 0°C and a PPM using wet-bulb temperature. Compared to the upper benchmark, the snow-level based PPM improved the estimation of snowfall occurrence by 5% and the simulation of snow water equivalent (SWE) by 9% during the snow melting season. In contrast, due to missing processes, the microphysics scheme decreased performances in phase estimate and SWE simulations compared to the upper benchmark. These results highlight the need for detailed evaluation of the precipitation phase from atmospheric models and the benefit for mountain snow hydrology of the post-processed snow level. The limitations to drive snowpack models at slope scale are also discussed.

Abstract Freezing rain events have caused severe socioeconomic and ecosystem impacts. An understanding of how these events may evolve as the Earth warms is necessary to adequately adapt infrastructure to these changes. We present an analysis of projected changes to freezing rain events over North America relative to the 1980–2009 recent past climate for the periods during which +2, +3, and +4°C of global warming is attained. We diagnose freezing rain using four precipitation‐type algorithms (Cantin and Bachand, Bourgouin, Ramer, and Baldwin) applied to four simulations of the fifth‐generation Canadian Regional Climate Model (CRCM5) driven by four global climate models (GCMs). We find that the choice of driving GCM strongly influences the spatial pattern of projected change. The choice of algorithm has a comparatively smaller impact, and primarily affects the magnitude but not the sign of projected change. We identify several regions where all simulations and algorithms agree on the sign of change, with increases projected over portions of western Canada and decreases over the central, eastern, and southern United States. However, we also find large regions of disagreement on the sign of change depending on driving GCM and even ensemble member of the same GCM, highlighting the importance of examining freezing rain events in a multi‐member ensemble of simulations driven by multiple GCMs to sufficiently account for uncertainty in projections of these hazardous events.

Abstract The collection efficiency of a typical precipitation gauge-shield configuration decreases with increasing wind speed, with a high scatter for a given wind speed. The high scatter in the collection efficiency for a given wind speed arises in part from the variability in the characteristics of falling snow and atmospheric turbulence. This study uses weighing gauge data collected at the Marshall Field Site near Boulder, Colorado, during the WMO Solid Precipitation Intercomparison Experiment (SPICE). Particle diameter and fall speed data from a laser disdrometer were used to show that the scatter in the collection efficiency can be reduced by considering the fall speed of solid precipitation particles. The collection efficiency was divided into two classes depending on the measured mean-event particle fall speed during precipitation events. Slower-falling particles were associated with a lower collection efficiency. A new transfer function (i.e., the relationship between collection efficiency and other meteorological variables, such as wind speed or air temperature) that includes the fall speed of the hydrometeors was developed. The root-mean-square error of the adjusted precipitation with the new transfer function with respect to a weighing gauge placed in a double fence intercomparison reference was lower than using previously developed transfer functions that only consider wind speed and air temperature. This shows that the measured fall speed of solid precipitation with a laser disdrometer accounts for a large amount of the observed scatter in weighing gauge collection efficiency.

Abstract Accurate snowfall measurement is challenging because it depends on the precipitation gauge used, meteorological conditions, and the precipitation microphysics. Upstream of weighing gauges, the flow field is disturbed by the gauge and any shielding used usually creates an updraft, which deflects solid precipitation from falling in the gauge, resulting in significant undercatch. Wind shields are often used with weighing gauges to reduce this updraft, and transfer functions are required to adjust the snowfall measurements to consider gauge undercatch. Using these functions reduces the bias in precipitation measurement but not the root-mean-square error (RMSE). In this study, the accuracy of the Hotplate precipitation gauge was compared to standard unshielded and shielded weighing gauges collected during the WMO Solid Precipitation Intercomparison Experiment program. The analysis performed in this study shows that the Hotplate precipitation gauge bias after wind correction is near zero and similar to wind corrected weighing gauges. The RMSE of the Hotplate precipitation gauge measurements is lower than weighing gauges (with or without an Alter shield) for wind speeds up to 5 m s −1 , the wind speed limit at which sufficient data were available. This study shows that the Hotplate precipitation gauge measurement has a low bias and RMSE due to its aerodynamic shape, making its performance mostly independent of the type of solid precipitation.

Abstract. The continental divide along the spine of the Canadian Rockies in southwestern Canada is a critical headwater region for hydrological drainages to the Pacific, Arctic, and Atlantic oceans. Major flooding events are typically attributed to heavy precipitation on its eastern side due to upslope (easterly) flows. Precipitation can also occur on the western side of the divide when moisture originating from the Pacific Ocean encounters the west-facing slopes of the Canadian Rockies. Often, storms propagating across the divide result in significant precipitation on both sides. Meteorological data over this critical region are sparse, with few stations located at high elevations. Given the importance of all these types of events, the Storms and Precipitation Across the continental Divide Experiment (SPADE) was initiated to enhance our knowledge of the atmospheric processes leading to storms and precipitation on either side of the continental divide. This was accomplished by installing specialized meteorological instrumentation on both sides of the continental divide and carrying out manual observations during an intensive field campaign from 24 April–26 June 2019. On the eastern side, there were two field sites: (i) at Fortress Mountain Powerline (2076 m a.s.l.) and (ii) at Fortress Junction Service, located in a high-elevation valley (1580 m a.s.l.). On the western side, Nipika Mountain Resort, also located in a valley (1087 m a.s.l.), was chosen as a field site. Various meteorological instruments were deployed including two Doppler light detection and ranging instruments (lidars), three vertically pointing micro rain radars, and three optical disdrometers. The three main sites were nearly identically instrumented, and observers were on site at Fortress Mountain Powerline and Nipika Mountain Resort during precipitation events to take manual observations of precipitation type and microphotographs of solid particles. The objective of the field campaign was to gather high-temporal-frequency meteorological data and to compare the different conditions on either side of the divide to study the precipitation processes that can lead to catastrophic flooding in the region. Details on field sites, instrumentation used, and collection methods are discussed. Data from the study are publicly accessible from the Federated Research Data Repository at https://doi.org/10.20383/101.0221 (Thériault et al., 2020). This dataset will be used to study atmospheric conditions associated with precipitation events documented simultaneously on either side of a continental divide. This paper also provides a sample of the data gathered during a precipitation event.

Abstract. The interior of western Canada, like many similar cold mid- to high-latitude regions worldwide, is undergoing extensive and rapid climate and environmental change, which may accelerate in the coming decades. Understanding and predicting changes in coupled climate–land–hydrological systems are crucial to society yet limited by lack of understanding of changes in cold-region process responses and interactions, along with their representation in most current-generation land-surface and hydrological models. It is essential to consider the underlying processes and base predictive models on the proper physics, especially under conditions of non-stationarity where the past is no longer a reliable guide to the future and system trajectories can be unexpected. These challenges were forefront in the recently completed Changing Cold Regions Network (CCRN), which assembled and focused a wide range of multi-disciplinary expertise to improve the understanding, diagnosis, and prediction of change over the cold interior of western Canada. CCRN advanced knowledge of fundamental cold-region ecological and hydrological processes through observation and experimentation across a network of highly instrumented research basins and other sites. Significant efforts were made to improve the functionality and process representation, based on this improved understanding, within the fine-scale Cold Regions Hydrological Modelling (CRHM) platform and the large-scale Modélisation Environmentale Communautaire (MEC) – Surface and Hydrology (MESH) model. These models were, and continue to be, applied under past and projected future climates and under current and expected future land and vegetation cover configurations to diagnose historical change and predict possible future hydrological responses. This second of two articles synthesizes the nature and understanding of cold-region processes and Earth system responses to future climate, as advanced by CCRN. These include changing precipitation and moisture feedbacks to the atmosphere; altered snow regimes, changing balance of snowfall and rainfall, and glacier loss; vegetation responses to climate and the loss of ecosystem resilience to wildfire and disturbance; thawing permafrost and its influence on landscapes and hydrology; groundwater storage and cycling and its connections to surface water; and stream and river discharge as influenced by the various drivers of hydrological change. Collective insights, expert elicitation, and model application are used to provide a synthesis of this change over the CCRN region for the late 21st century.

In the future, the intensity, phases, and frequency of precipitation are expected to change due to global warming, in particular during colder seasons when temperatures are near 0°C. To investigate the impacts of warmer atmospheric conditions on the microphysical processes that lead to several precipitation types, the extreme 1998 Ice Storm was simulated using the Weather Research and Forecasting (WRF) model, with and without a pseudo-global warming. The pseudo-global warming approach simulates similar large-scale conditions but in warmer conditions, which allows for the assessment of thermodynamic feedback from cloud and precipitation microphysics. For both simulations, WRF was coupled with the Predicted Particle Properties (P3) bulk microphysics scheme that predicts the liquid fraction of mixed-phase particles. Results of the pseudo-global warming simulation show an increase of ∼828 m in the upper 0°C level and a northeastward migration (∼60 km) of the rain-snow transition region. The results also show a 20% decrease in domain-averaged freezing rain amounts, but with an increased maximum amount of 50%. The horizontal distance associated with a melting aloft and a refreezing layer near the surface is 105 km longer in southern Quebec due to the combined effects of the pseudo-warming and the presence of the Appalachian Mountains. The microphysical processes that lead to precipitation are impacted as well; the increased ice mass and riming conditions aloft in warmer temperatures result in higher liquid precipitation rates. This study contributes to our understanding of the changes in the fine-scale processes of an extreme storm, simulated with pseudo-global warming conditions.

Abstract A prognostic equation for the liquid fraction of mixed-phase particles has been recently added to the Predicted Particle Properties (P3) bulk microphysics scheme. Mixed-phase particles are necessary to simulate key microphysical processes leading to various winter precipitation types, such as ice pellets and freezing rain. To illustrate the impacts of predicting the bulk liquid fraction, the 1998 North American Ice Storm is simulated using the Weather Research and Forecasting (WRF) Model with the modified P3 scheme. It is found that simulating partial melting by predicting the bulk liquid fraction produces higher mass and number mixing ratios of rain. This leads to smaller rain sizes reaching the refreezing layer as well as a decrease in the freezing rain accumulation at the surface by up to 30% in some locations compared to when no liquid fraction is predicted. The increase in fall speed and density and decrease of particle diameter during partial melting combined with an improved representation of the refreezing process in the modified P3 leads to generally higher total solid surface precipitation rates than using the original P3 scheme. There is also an increase of solid precipitation in regions of ice pellet accumulation. Overall, the simulation of mixed-phase particles notably impacts the vertical and spatial distributions of precipitation properties.

Abstract. The 0 ∘C temperature threshold is critical for many meteorological and hydrological processes driven by melting and freezing in the atmosphere, surface, and sub-surface and by the associated precipitation varying between rain, freezing rain, wet snow, and snow. This threshold is especially important in cold regions such as Canada, because it is linked with freeze–thaw, snowmelt, and permafrost. This study develops a Canada-wide perspective on near-0 ∘C conditions using hourly surface temperature and precipitation type observations from 92 climate stations for the period from 1981 to 2011. In addition, nine stations from various climatic regions are selected for further analysis. Near-0 ∘C conditions are defined as periods when the surface temperature is between −2 and 2 ∘C. Near-0 ∘C conditions occur often across all regions of the country, although the annual number of days and hours and the duration of these events varies dramatically. Various types of precipitation (e.g., rain, freezing rain, wet snow, and ice pellets) sometimes occur with these temperatures. Near-0 ∘C conditions and the reported precipitation type occurrences tend to be higher in Atlantic Canada, although high values also occur in other regions. Trends of most temperature-based and precipitation-based indicators show little or no change despite a systematic warming in annual surface temperatures over Canada. Over the annual cycle, near-0 ∘C temperatures and precipitation often exhibit a pattern: short durations occur around summer, driven by the diurnal cycle, and a tendency toward longer durations around winter, associated with storms. There is also a tendency for near-0 ∘C surface temperatures to occur more often than expected relative to other temperature windows at some stations due, at least in part, to diabatic cooling and heating that take place with melting and freezing, respectively, in the atmosphere and at the surface.

Abstract. The Interior of Western Canada, up to and including the Arctic, has experienced rapid change in its climate, hydrology, cryosphere and ecosystems and this is expected to continue. Although there is general consensus that warming will occur in the future, many critical issues remain. In this first of two articles, attention is placed on atmospheric-related issues that range from large scales down to individual precipitation events. Each of these is considered in terms of expected change organized by season and utilizing climate scenario information as well as thermodynamically-driven future climatic forcing simulations. Large scale atmospheric circulations affecting this region are generally projected to become stronger in each season and, coupled with warming temperatures, lead to enhancements of numerous water-related and temperature-related extremes. These include winter snowstorms, freezing rain, drought as well as atmospheric forcing of spring floods although not necessarily summer convection. Collective insights of these atmospheric findings are summarized in a consistent, connected physical framework.

Abstract. The 0 °C temperature threshold is critical to many meteorological and hydrological processes driven by melting and freezing in the atmosphere, surface and sub-surface and by the associated precipitation varying between rain, freezing rain, wet snow and snow. This threshold, linked with freeze-thaw, is especially important in cold regions such as Canada. This study develops a Canada-wide perspective on near 0 °C conditions with a particular focus on the occurrence of its associated precipitation. Since this analysis requires hourly values of surface temperature and precipitation type observations, it was limited to 92 stations over the 1981–2011 period. In addition, nine stations representative of various climatic regions are selected for further analysis. Near 0 °C conditions are defined as periods when the surface temperature is between −2 °C and 2 °C. Near 0 °C conditions occur often across all regions of the country although the annual number of days and hours and the duration of these events varies dramatically. Various forms of precipitation (including rain, freezing rain, wet snow and ice pellets) are sometimes linked with these temperatures with highest fractions tending to occur in Atlantic Canada. Trends of most temperature-based and precipitation-based indicators show little or no change despite a systematic warming in annual temperatures. Over the annual cycle, near 0 °C temperatures and precipitation often exhibit a pattern with short durations near summer driven by the diurnal cycle, while longer durations tend to occur more towards winter associated with storms. There is also a tendency for near 0 °C temperatures to occur more often than expected relative to other temperature windows; due at least in part to diabatic cooling and heating occurring with melting and freezing, respectively, in the atmosphere and at the surface.

Abstract Bulk microphysics parameterizations that are used to represent clouds and precipitation usually allow only solid and liquid hydrometeors. Predicting the bulk liquid fraction on ice allows an explicit representation of mixed-phase particles and various precipitation types, such as wet snow and ice pellets. In this paper, an approach for the representation of the bulk liquid fraction into the predicted particle properties (P3) microphysics scheme is proposed and described. Solid-phase microphysical processes, such as melting and sublimation, have been modified to account for the liquid component. New processes, such as refreezing and condensation of the liquid portion of mixed-phase particles, have been added to the parameterization. Idealized simulations using a one-dimensional framework illustrate the overall behavior of the modified scheme. The proposed approach compares well to a Lagrangian benchmark model. Temperatures required for populations of ice crystals to melt completely also agree well with previous studies. The new processes of refreezing and condensation impact both the surface precipitation type and feedback between the temperature and the phase changes. Overall, prediction of the bulk liquid fraction allows an explicit description of new precipitation types, such as wet snow and ice pellets, and improves the representation of hydrometeor properties when the temperature is near 0°C.

Abstract. The phase of precipitation and its distribution at the surface can affect water resources and the regional water cycle of a region. A field project was held in March–April 2015 on the eastern slope of the Canadian Rockies to document precipitation characteristics and associated atmospheric conditions. During the project, 60 % of the particles documented were rimed in relatively warm and dry conditions. Rain–snow transitions also occurred aloft and at the surface in sub-saturated conditions. Ice-phase precipitation falling through a saturated atmospheric layer with temperatures > 0 ∘C will start melting. In contrast, if the melting layer is sub-saturated, the ice-phase precipitation undergoes sublimation, which increases the depth of the rain–snow transition. In this context, this study investigates the role of sublimation and riming in precipitation intensity and type reaching the surface in the Kananaskis Valley, Alberta, during March–April 2015. To address this, a set of numerical simulations of an event of mixed precipitation observed at the surface was conducted. This event on 31 March 2015 was documented with a set of devices at the main observation site (Kananaskis Emergency Services, KES), including a precipitation gauge, disdrometer, and micro rain radar. Sensitivity experiments were performed to assess the impacts of temperature changes from sublimation and the role of the production of graupel (riming) aloft in the surface precipitation evolution. A warmer environment associated with no temperature changes from sublimation leads to a peak in the intensity of graupel at the surface. When the formation of graupel is not considered, the maximum snowfall rate occurred at later times. Results suggest that unrimed snow reaching the surface is formed on the western flank and is advected eastward. In contrast, graupel would form aloft in the Kananaskis Valley. The cooling from sublimation and melting by rimed particles increases the vertical shear near KES. Overall, this study illustrated that the presence of graupel influenced the surface evolution of precipitation type in the valley due to the horizontal transport of precipitation particles.

Abstract. The interior of western Canada, up to and including the Arctic, has experienced rapid change in its climate, hydrology, cryosphere, and ecosystems, and this is expected to continue. Although there is general consensus that warming will occur in the future, many critical issues remain. In this first of two articles, attention is placed on atmospheric-related issues that range from large scales down to individual precipitation events. Each of these is considered in terms of expected change organized by season and utilizing mainly “business-as-usual” climate scenario information. Large-scale atmospheric circulations affecting this region are projected to shift differently in each season, with conditions that are conducive to the development of hydroclimate extremes in the domain becoming substantially more intense and frequent after the mid-century. When coupled with warming temperatures, changes in the large-scale atmospheric drivers lead to enhancements of numerous water-related and temperature-related extremes. These include winter snowstorms, freezing rain, drought, forest fires, as well as atmospheric forcing of spring floods, although not necessarily summer convection. Collective insights of these atmospheric findings are summarized in a consistent, connected physical framework.

Freezing rain occurs in complex atmospheric conditions when the temperature is close to 0°C. To better understand how its occurrence will change in the future, there is a need to assess how well regional climate models can reproduce those conditions. The goal of the present study is to investigate the influence of horizontal resolution on the simulation of freezing rain using the fifth generation of the Canadian Regional Climate Model (CRCM5). Three CRCM5 simulations driven by the European Centre for Medium-range Weather Forecasts interim reanalysis (ERA-Interim) over eastern North America at resolutions of 0.11°, 0.22°, and 0.44° were conducted over a period of 36 years (1979–2014). Freezing rain is diagnosed using an in-line diagnostic method for precipitation partitioning. A climatological study of annual and seasonal accumulated freezing rain was conducted. In addition, the ability of the three simulations to reproduce individual freezing rain events was evaluated. Our analyses include frequency and partitioning of different precipitation types and comparisons with observations. All simulations reproduced the climatology of freezing rain sufficiently well and show similar large-scale patterns. The number of freezing rain events tends to be overestimated at higher resolution and underestimated at lower resolution. Despite the overestimation, detailed maxima associated with freezing rain are well defined and located at higher resolution, notably in regions of the St. Lawrence River Valley. Overall, this study is consistent with other added-value studies, generally showing a mix of improvement and deterioration in the precipitation fields by the higher resolution simulations.

Abstract A freezing rain event, in which the Meteorological Centre of Canada’s 2.5-km numerical weather prediction system significantly underpredicted the quantity of freezing rain, is examined. The prediction system models precipitation types explicitly, directly from the Milbrandt–Yau microphysics scheme. It was determined that the freezing rain underprediction for this case was due primarily to excessive refreezing of rain, originating from melting snow and graupel, in and under the temperature inversion of the advancing warm front ultimately depleting the supply of rain reaching the surface. The refreezing was caused from excessive collisional freezing between rain and graupel. Sensitivity experiments were conducted to examine the effects of a temperature threshold for collisional freezing and on varying the values of the collection efficiencies between rain and ice-phase hydrometeors. It was shown that by reducing the rain–graupel collection efficiency and by imposing a temperature threshold of −5°C, above which collisional freezing is not permitted, excessive rain–graupel collection and graupel formation can be controlled in the microphysics scheme, leading to an improved simulation of freezing rain at the surface.

Abstract. Transfer functions are generally used to adjust for the wind-induced undercatch of solid precipitation measurements. These functions are derived based on the variation of the collection efficiency with wind speed for a particular type of gauge, either using field experiments or based on numerical simulation. Most studies use the wind speed alone, while others also include surface air temperature and/or precipitation type to try to reduce the scatter of the residuals at a given wind speed. In this study, we propose the use of the measured precipitation intensity to improve the effectiveness of the transfer function. This is achieved by applying optimized curve fitting to field measurements from the Marshall field-test site (CO, USA). The use of a non-gradient optimization algorithm ensures optimal binning of experimental data according to the parameter under test. The results reveal that using precipitation intensity as an explanatory variable significantly reduce the scatter of the residuals. The scatter reduction as indicated by the Root Mean Square Error (RMSE) is confirmed by the analysis of the recent quality controlled data from the WMO/SPICE campaign, showing that this approach can be applied to a variety of locations and catching-type gauges. We demonstrate the physical basis of the relationship between the collection efficiency and the measured precipitation intensity, due to the correlation of large particles with high intensities, by conducting a Computational Fluid-Dynamics (CFD) simulation. We use a Reynolds Averaged Navier-Stokes SST k-ω model coupled with a Lagrangian particle-tracking model. Results validate the hypothesis of using the measured precipitation intensity as a key parameter to improve the correction of wind-induced undercatch. Findings have the potential to improve operational measurements since no additional instrument other than a wind sensor is required to apply the correction. This improves the accuracy of precipitation measurements without the additional cost of ancillary instruments such as particle counters.

Abstract. Precipitation events that bring rain and snow to the Banff–Calgary area of Alberta are a critical aspect of the region's water cycle and can lead to major flooding events such as the June 2013 event that was the second most costly natural disaster in Canadian history. Because no special atmospheric-oriented observations of these events have been made, a field experiment was conducted in March and April 2015 in Kananaskis, Alberta, to begin to fill this gap. The goal was to characterize and better understand the formation of the precipitation at the surface during spring 2015 at a specific location in the Kananaskis Valley. Within the experiment, detailed measurements of precipitation and weather conditions were obtained, a vertically pointing Doppler radar was deployed and weather balloons were released. Although 17 precipitation events occurred, this period was associated with much less precipitation than normal (−35 %) and above-normal temperatures (2.5 ∘C). Of the 133 h of observed precipitation, solid precipitation occurred 71 % of the time, mixed precipitation occurred 9 % and rain occurred 20 %. An analysis of 17 504 precipitation particles from 1181 images showed that a wide variety of crystals and aggregates occurred and approximately 63 % showed signs of riming. This was largely independent of whether flows aloft were upslope (easterly) or downslope (westerly). In the often sub-saturated surface conditions, hydrometeors containing ice occurred at temperatures as high as 9 ∘C. Radar structures aloft were highly variable with reflectivity sometimes >30 dBZe and Doppler velocity up to −1 m s−1, which indicates upward motion of particles within ascending air masses. Precipitation was formed in this region within cloud fields sometimes having variable structures and within which supercooled water at least sometimes existed to produce accreted particles massive enough to reach the surface through the relatively dry sub-cloud region.

Winter weather events with temperatures near $$0\,^\circ\mathrm{{C}}$$ are often associated with freezing rain. They can have major impacts on the society by causing power outages and disruptions to the transportation networks. Despite the catastrophic consequences of freezing rain, very few studies have investigated how their occurrences could evolve under climate change. This study aims to investigate the change of freezing rain and ice pellets over southern Québec using regional climate modeling at high resolution. The fifth-generation Canadian Regional Climate Model with climate scenario RCP 8.5 at $$0.11^\circ$$ grid mesh was used. The precipitation types such as freezing rain, ice pellets or their combination are diagnosed using five methods (Cantin and Bachand, Bourgouin, Ramer, Czys and, Baldwin). The occurrences of the diagnosed precipitation types for the recent past (1980–2009) are found to be comparable to observations. The projections for the future scenario (2070–2099) suggested a general decrease in the occurrences of mixed precipitation over southern Québec from October to April. This is mainly due to a decrease in long-duration events ( $$\ge 6\,\mathrm{{h}}$$ ). Overall, this study contributes to better understand how the distribution of freezing rain and ice pellets might change in the future using high-resolution regional climate model.

Abstract A devastating, flood-producing rainstorm occurred over southern Alberta, Canada, from 19 to 22 June 2013. The long-lived, heavy rainfall event was a result of complex interplays between topographic, synoptic, and convective processes that rendered an accurate simulation of this event a challenging task. In this study, the Weather Research and Forecasting (WRF) Model was used to simulate this event and was validated against several observation datasets. Both the timing and location of the model precipitation agree closely with the observations, indicating that the WRF Model is capable of reproducing this type of severe event. Sensitivity tests with different microphysics schemes were conducted and evaluated using equitable threat and bias frequency scores. The WRF double-moment 6-class microphysics scheme (WDM6) generally performed better when compared with other schemes. The application of a conventional convective/stratiform separation algorithm shows that convective activity was dominant during the early stages, then evolved into predominantly stratiform precipitation later in the event. The HYSPLIT back-trajectory analysis and regional water budget assessments using WRF simulation output suggest that the moisture for the precipitation was mainly from recycling antecedent soil moisture through evaporation and evapotranspiration over the Canadian Prairies and the U.S. Great Plains. This analysis also shows that a small fraction of the moisture can be traced back to the northeastern Pacific, and direct uptake from the Gulf of Mexico was not a significant source in this event.

In June 2013, excessive rainfall associated with an intense weather system triggered severe flooding in southern Alberta, which became the costliest natural disaster in Canadian history. This article provides an overview of the climatological aspects and large-scale hydrometeorological features associated with the flooding event based upon information from a variety of sources, including satellite data, upper air soundings, surface observations and operational model analyses. The results show that multiple factors combined to create this unusually severe event. The event was characterized by a slow-moving upper level low pressure system west of Alberta, blocked by an upper level ridge, while an associated well-organized surface low pressure system kept southern Alberta, especially the eastern slopes of the Rocky Mountains, in continuous precipitation for up to two days. Results from air parcel trajectory analysis show that a significant amount of the moisture originated from the central Great Plains, transported into Alberta by a southeasterly low level jet. The event was first dominated by significant thunderstorm activity, and then evolved into continuous precipitation supported by the synoptic-scale low pressure system. Both the thunderstorm activity and upslope winds associated with the low pressure system produced large rainfall amounts. A comparison with previous similar events occurring in the same region suggests that the synoptic-scale features associated with the 2013 rainfall event were not particularly intense; however, its storm environment was the most convectively unstable. The system also exhibited a relatively high freezing level, which resulted in rain, rather than snow, mainly falling over the still snow-covered mountainous areas. Melting associated with this rain-on-snow scenario likely contributed to downstream flooding. Furthermore, above-normal snowfall in the preceding spring helped to maintain snow in the high-elevation areas, which facilitated the rain-on-snow event.