Informing BC Stakeholders

You are here

Publications Library

  • Source Publication: Environmental Science and Technology, 57, 19, 7401–7409, doi:10.1021/acs.est.2c08243 Authors: Lao, I.R., A. Feinberg, and N. Borduas-Dedekind Publication Date: Jun 2023

    Selenium (Se) is an essential nutrient for humans and enters our food chain through bioavailable Se in soil. Atmospheric deposition is a major source of Se to soils, driving the need to investigate the sources and sinks of atmospheric Se. Here, we used Se concentrations from PM2.5 data at 82 sites from 1988 to 2010 from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network in the US to identify the sources and sinks of particulate Se. We identified 6 distinct seasonal profiles of atmospheric Se, grouped by geographical location: West, Southwest, Midwest, Southeast, Northeast, and North Northeast. Across most of the regions, coal combustion is the largest Se source, with a terrestrial source dominating in the West. We also found evidence for gas-to-particle partitioning in the wintertime in the Northeast. Wet deposition is an important sink of particulate Se, as determined by Se/PM2.5 ratios. The Se concentrations from the IMPROVE network compare well to modeled output from a global chemistry-climate model, SOCOL-AER, except in the Southeast US. Our analysis constrains the sources and sinks of atmospheric Se, thereby improving the predictions of Se distribution under climate change.

  • Source Publication: Journal of Climate, 36, 18, 6393-6407, doi:10.1175/JCLI-D-22-0713.1 Authors: Sun, Q. F.W. Zwiers, X. Zhang and Y. Tan Publication Date: Jun 2023

    El Niño–Southern Oscillation (ENSO) has a profound influence on the occurrence of extreme precipitation events at local and regional scales in the present-day climate, and thus it is important to understand how that influence may change under future global warming. We consider this question using the large-ensemble simulations of CESM2, which simulates ENSO well historically. CESM2 projects that the influence of ENSO on extreme precipitation will strengthen further under the SSP3–7.0 scenario in most regions whose extreme precipitation regimes are strongly affected by ENSO in the boreal cold season. Extreme precipitation in the boreal cold season that exceeds historical thresholds is projected to become more common throughout the ENSO cycle. The difference in the intensity of extreme precipitation events that occur under El Niño and La Niña conditions will increase, resulting in “more extreme and more variable hydroclimate extremes.” We also consider the processes that affect the future intensity of extreme precipitation and how it varies with the ENSO cycle by partitioning changes into thermodynamic and dynamic components. The thermodynamic component, which reflects increases in atmospheric moisture content, results in a relatively uniform intensification of ENSO-driven extreme precipitation variation. In contrast, the dynamic component, which reflects changes in vertical motion, produces a strong regional difference in the response to forcing. In some regions, this component amplifies the thermodynamic-induced changes, while in others, it offsets them or even results in reduction in extreme precipitation variation.

  • Authors: City of Terrace, The Pacific Climate Impacts Consortium, Pinna Sustainability Publication Date: May 2023

    The Climate Projections for the City of Terrace report provides projections and impacts analysis for the City of Terrace, BC and is intended to support decision making throughout the region and to help community partners better understand how their work may be affected by the changing climate

  • Source Publication: Journal of Hydrology: Regional Studies, 44, 101237, doi:10.1016/j.ejrh.2022.101237 Authors: Larabi, S., M. A. Schnorbus, and F. Zwiers Publication Date: Dec 2022

    Study region:
    Nechako Reservoir, British Columbia, Canada.

    Study focus:
    Hydrological regulation affect both hydrological and thermal conditions in the reservoir and downstream reach, subsequently disrupting fish habitats. This paper aims at developing an integrated model simulating physical processes that govern the quantity and quality of inflow, reservoir, and outflow water of the Nechako Reservoir. Such a model would help stakeholders understand the response of in-reservoir water temperature stratification and downstream water temperature to changes in inflow and reservoir operation under future climate change.

    New hydrological insights for the region:
    The model was calibrated against historical reservoir levels and in-reservoir and outlet water temperature field data. The integrated model simulated accurately the wide variation of reservoir levels as well as the in-reservoir water temperature at Kenney Dam and the outlet temperature. Sensitivity analysis shows that reservoir water temperature particularly the epilimnion is sensitive to changes in both meteorological and hydrological forcing. Forcing the model with different outflow scenarios shows the weak sensitivity of temperature of water released to outflow rates. Given epilimnion water releases at the spillway, the Summer Temperature Management Program could be inefficient to provide cool water in the Nechako River during the critical period of salmon migration in a warming climate. However, colder water remains available at depth at Kenney Dam to potentially mitigate and better control downstream water temperature.

  • Source Publication: International Journal of Climatology, 42, 10, doi: org/10.1002/joc.7833.8 Authors: Diaconescu, E., H. Sankare, K. Chow, T.Q. Murdock, and A.J. Cannon Publication Date: Dec 2022

    The projected increase in the frequency and intensity of extreme heat events due to climate change means an associated increase in risk of heat-related illnesses and mortality. Public health systems need to be prepared to identify and reduce the susceptibility of vulnerable populations to increased occurrence of heat-related illness and stress. To facilitate this, climate services have begun developing climate change projections for heat-stress indices based on exceedances of thresholds used operationally in meteorological heat warning systems. This task is complicated by the fact that heat-stress indices are generally computed using hourly data whereas climate model outputs are often archived at daily or longer time steps. This study focuses on Humidex, a heat-stress index used in heat alerts issued by the Meteorological Service of Canada. Several potential solutions for computing robust Humidex indices using daily data are examined, including a new approximation method. Indices obtained with the new method are compared with indices obtained using the classic method based on hourly data as well as with other two methods based on average daily values. The new approximation gives good estimations for humidex indices, while the daily-average-value methods present biases with respect to the hourly-value method.

  • Source Publication: Journal of Hydrology X, 17, 100144, doi:10.1016/j.hydroa.2022.100144 Authors: Tsuruta, K. and M. A. Schnorbus Publication Date: Dec 2022

    As glaciers across the world continue to recede, there is a concern that their loss as a fresh water reservoir within mountainous basins will have a negative impact on stream temperatures and downstream water resources. Currently, there are relatively few glacio-hydrological models (GHMs) appropriate to study such phenomena and studies that have used GHMs generally acknowledge the high uncertainty associated with their simulations. Calibration techniques present a particular issue in GHMs as available glacier observations are limited and errors in the glacierized portion of a basin can be compensated by errors in the non-glacierized portion. Using as a study site the Cheakamus Basin in British Columbia, Canada, we 1) present a new, fully-coupled GHM, 2) analyze the effects different calibration techniques have on the model’s summer streamflow projections, and 3) compare the fully-coupled GHM results to projections using a one-way GHM. The calibration techniques studied vary in terms of glacier representation (dynamic/static), and glacier constraint (mass balance/thinning rates/thinning rates and area change). We find projected future climate forcings are sufficiently strong in the Cheakamus Basin so as to generally make the sign and significance of changes to the basin’s hydrology insensitive to the calibration and projection procedures studied. However, the variation among these procedures produces significant changes in the projected magnitude of future hydrological changes and therefore should be carefully considered in studies where precision beyond the sign and significance of change is required. Based on analysis of the variation within each procedure’s set of model outputs, we conclude 1) the two-way GHM has benefits over the one-way model, 2) calibration using dynamic glaciers and a thinning rate constraint is preferable for the new GHM, and 3) there is a need for additional studies on the uncertainties associated with the calibration of glacio-hydrological models.

  • Source Publication: Canadian Journal of Statistics, 50, 4, 1355-1386, doi:10.1002/cjs.11743 Authors: Dean, C.B., A.H. El-Shaarawi, S.R. Esterby, J. Mills-Flemming, R.D. Routledge, S.W. Taylor, D.G. Woolford, J.V. Zidek and F.W. Zwiers Publication Date: Dec 2022

    This article focuses on the importance of collaboration in statistics by Canadian researchers and highlights the contributions that Canadian statisticians have made to many research areas in environmetrics. We provide a discussion about different vehicles that have been developed for collaboration by Canadians in the environmetrics context as well as specific scientific areas that are focused on environmetrics research in Canada including climate science, forestry, and fisheries, which are areas of importance for natural resources in Canada.

  • Source Publication: Earth's Future, 10, doi:10.1029/2022EF002797 Authors: Li, M., C. Li, Z. Jiang, X. Zhang and F.W. Zwiers Publication Date: Dec 2022

    Observations show that summer precipitation in China has undergone pronounced changes, resulting in an enigmatic “north-south drying-wetting” pattern in eastern China that is of great concern for socio-economic development. Scientific consensus on the mechanisms that are responsible for this pattern of change has not yet been achieved. We show that this complex pattern of summer total precipitation trends observed in China since the 1960s is overwhelmingly the result of changes in daily precipitation frequency, rather than being the result of changes in precipitation intensity or the frequency of synoptic circulation patterns favorable to precipitation. Changes in precipitation intensity, which are very likely due to anthropogenic greenhouse gas forcing, contribute a relatively homogeneous wetting tendency across the country while changes due to synoptic circulation change are weak. The changes in daily precipitation frequency that drive the observed patterns of change may be due to aerosols, but improved process understanding will be required to resolve that question and enable reliable projections of regional scale precipitation change in China and elsewhere.

  • Source Publication: Earth System Dynamics, 13, 1689–1713, doi:10.5194/esd-13-1689-2022 Authors: Philip, S. Y., S.F. Kew, G.J. van Oldenborgh, F.S. Anslow, S.I. Seneviratne, R. Vautard, D. Coumou, K.L. Ebi, J. Arrighi, R. Singh, M. van Aalst, C. Pereira Marghidan, M. Wehner, W. Yang, S. Li, D.L. Schumacher, M. Hauser, R. Bonnet, L.N. Luu, F. Lehner, Publication Date: Dec 2022

    Towards the end of June 2021, temperature records were broken by several degrees Celsius in several cities in the Pacific Northwest areas of the US and Canada, leading to spikes in sudden deaths and sharp increases in emergency calls and hospital visits for heat-related illnesses. Here we present a multi-model, multi-method attribution analysis to investigate the extent to which human-induced climate change has influenced the probability and intensity of extreme heat waves in this region. Based on observations, modelling and a classical statistical approach, the occurrence of a heat wave defined as the maximum daily temperature (TXx) observed in the area 45–52∘ N, 119–123∘ W, was found to be virtually impossible without human-caused climate change. The observed temperatures were so extreme that they lay far outside the range of historical temperature observations. This makes it hard to state with confidence how rare the event was. Using a statistical analysis that assumes that the heat wave is part of the same distribution as previous heat waves in this region led to a first-order estimation of the event frequency of the order of once in 1000 years under current climate conditions. Using this assumption and combining the results from the analysis of climate models and weather observations, we found that such a heat wave event would be at least 150 times less common without human-induced climate change. Also, this heat wave was about 2 ∘C hotter than a 1-in-1000-year heat wave would have been in 1850–1900, when global mean temperatures were 1.2 ∘C cooler than today. Looking into the future, in a world with 2 ∘C of global warming (0.8 ∘C warmer than today), a 1000-year event would be another degree hotter. Our results provide a strong warning: our rapidly warming climate is bringing us into uncharted territory with significant consequences for health, well-being and livelihoods. Adaptation and mitigation are urgently needed to prepare societies for a very different future.

  • Source Publication: Journal of Applied Meteorology and Climatology, 61, 9, 1139-1157. doi:10.1175/JAMC-D-21-0205.1 Authors: Lao, I. R., C. Abraham, E. Wiebe, and A.H. Monahan Publication Date: Dec 2022

    Nocturnal warming events (NWEs) are abrupt interruptions in the typical cooling of surface temperatures at night. Using temperature time series from the high-resolution Vancouver Island School-Based Weather Station Network (VWSN) in British Columbia, Canada, we investigate temporal and spatial characteristics of NWEs. In this coastal region, NWEs are more frequently detected in winter than in summer, with a seasonal shift from slowly warming NWEs dominating the winter months to rapidly warming NWEs dominating the summer months. Slow-warming NWEs are of relatively small amplitude and exhibit slow cooling rates after the temperature peaks. In contrast, fast-warming NWEs have a temperature increase of several kelvins with shorter-duration temperature peaks. The median behavior of these distinct NWE classes at individual stations is similar across the entire set of stations. The spatial synchronicity of NWEs across the VWSN (determined by requiring NWEs at station pairs to occur within given time windows) decreases with distance, including substantial variability at nearby stations that reflects local influences. Fast-warming NWEs are observed to occur either simultaneously across a number of stations or in isolation at one station. Spatial synchronicity values are used to construct undirected networks to investigate spatial connectivity structures of NWEs. We find that, independent of individual seasons or NWE classes, the networks are largely unstructured, with no clear spatial connectivity structures related to local topography or direction.

  • Authors: Markus Schnorbus Publication Date: Jul 2022

    This report, the first volume in the VIC Generation 2 deployment reports, provides a description of VIC-Glacier (VIC-GL) model changes and upgrades.

  • Authors: Markus Schnorbus Publication Date: Jul 2022

    This report, the second volume in the VIC Generation 2 deployment reports, provides a description of modelling glacier dynamics with the HydroConductor Model.

  • Authors: Markus Schnorbus Publication Date: Jul 2022

    This report, the third volume in the VIC Generation 2 deployment reports, provides a description of vegetation and topography parameterization.

  • Authors: Markus Schnorbus Publication Date: Jul 2022

    This, the fifth volume in the VIC Generation 2 deployment reports, provides a description of model calibration.

  • Authors: Markus Schnorbus Publication Date: Jul 2022

    This, the sixth volume in the VIC Generation 2 deployment reports, provides a description of model set-up and deployment for the Peace, Fraser, and Columbia basins.

  • Source Publication: Weather and Climate Extremes, 36, 100441, doi:10.1016/j.wace.2022.100441 Authors: Gillett, N.P., A.J. Cannon, E. Malinina, M. Schnorbus, F. Anslow, Q. Sun, M. Kirchmeier-Young, F.W. Zwiers, C. Seiler, X. Zhang, G. Flato, H. Wan, G. Li and A. Castellan Publication Date: Jun 2022

    A strong atmospheric river made landfall in southwestern British Columbia, Canada on November 14th, 2021, bringing two days of intense precipitation to the region. The resulting floods and landslides led to the loss of at least five lives, cut Vancouver off entirely from the rest of Canada by road and rail, and made this the costliest natural disaster in the province's history. Here we show that when characterised in terms of storm-averaged water vapour transport, the variable typically used to characterise the intensity of atmospheric rivers, westerly atmospheric river events of this magnitude are approximately one in ten year events in the current climate of this region, and that such events have been made at least 60% more likely by the effects of human-induced climate change. Characterised in terms of the associated two-day precipitation, the event is substantially more extreme, approximately a one in fifty to one in a hundred year event, and the probability of events at least this large has been increased by a best estimate of 45% by human-induced climate change. The effects of this precipitation on streamflow were exacerbated by already wet conditions preceding the event, and by rising temperatures during the event that led to significant snowmelt, which led to streamflow maxima exceeding estimated one in a hundred year events in several basins in the region. Based on a large ensemble of simulations with a hydrological model which integrates the effects of multiple climatic drivers, we find that the probability of such extreme streamflow events in October to December has been increased by human-induced climate change by a best estimate of 120–330%. Together these results demonstrate the substantial human influence on this compound extreme event, and help motivate efforts to increase resiliency in the face of more frequent events of this kind in the future.

  • Source Publication: Bulletin of the American Meteorological Society, 103, 3, S50-S54, doi:10.1175/BAMS-D-21-0143.1 Authors: Liu, Y., C. Li, Y. Sun, F.W. Zwiers, X. Zhang, Z. Jiang and F. Zheng Publication Date: Mar 2022

    On 6–8 January 2021, a cold air outbreak swept across eastern China, peaking over the North China Plain the night of 6 January, when 219 weather stations recorded the lowest nighttime temperature since 1961. In total, 498 stations recorded the lowest daytime or nighttime temperature since 1961 during the 3-day event. This event, together with two other cold outbreaks that affected the region on 13–15 December 2020 and 29 December 2020–1 January 2021, led to historic peak electricity demand and resumption of the operation of the only remaining coalfired generating plant in Beijing. This analysis puts the cold outbreak into historical perspective by considering changes in the likelihood of such events over 1961–2020 in the context of a climate that is being warmed by anthropogenic forcing.

  • Source Publication: Climate Dynamics volume, 58, 793–809, doi:10.1007/s00382-021-05933-3 Authors: Tan, Y., S. Yang, F.W. Zwiers, Z. Wang and Q. Sun Publication Date: Feb 2022

    We report on the characteristics of precipitation associated with three types of landfalling atmospheric rivers (ARs) over western North America in the winter season from 1980 to 2004. The ARs are classified according to three landfalling regions as southern, middle and northern types. Two main centers of precipitation are associated with the contributions by the ARs: one over Baja California linked to the southern type of the ARs, and the other over Washington State correlated with the northern and middle types of the ARs. ARs are seen to play a dominant role in the occurrences of extreme precipitation events, with a proportionately greater impact on more extreme events. Moisture flux convergence makes the dominant contribution to precipitation when ARs and extreme precipitation occur simultaneously in the studied areas. Moisture flux convergence in these cases is, in turn, dominated by the mean and transient moisture transported by the transient wind, with greater contribution from the latter, which is mainly concentrated in certain areas. The magnitude and direction of vertically integrated vapor transport (IVT) also play a role in determining the amount of precipitation received in the three regions considered. Larger IVT magnitude corresponds to more precipitation, while an IVT direction of about 220° (0° indicating east wind) is most favorable for high precipitation amount, which is especially obvious for the northern type of the ARs.

  • Source Publication: Environmental Research Communications, 4, 1, 015009, doi:10.1088/2515-7620/ac4bab Authors: Wu, L., A. Elshorbagy, and M.S. Alam Publication Date: Jan 2022

    Understanding the dynamics of water-energy-food (WEF) nexus interactions with climate change and human intervention helps inform policymaking. This study demonstrates the WEF nexus behavior under ensembles of climate change, transboundary inflows, and policy options, and evaluates the overall nexus performance using a previously developed system dynamics-based WEF nexus model—WEF-Sask. The climate scenarios include a baseline (1986–2014) and near-future climate projections (2021–2050). The approach is demonstrated through the case study of Saskatchewan, Canada. Results show that rising temperature with increased rainfall likely maintains reliable food and feed production. The climate scenarios characterized by a combination of moderate temperature increase and slightly less rainfall or higher temperature increase with slightly higher rainfall are easier to adapt to by irrigation expansion. However, such expansion uses a large amount of water resulting in reduced hydropower production. In contrast, higher temperature, combined with less rainfall, such as SSP370 (+2.4 °C, −6 mm), is difficult to adapt to by irrigation expansion. Renewable energy expansion, the most effective climate change mitigation option in Saskatchewan, leads to the best nexus performance during 2021–2050, reducing total water demand, groundwater demand, greenhouse gas (GHG) emissions, and potentially increasing water available for food&feed production. In this study, we recommend and use food&feed and power production targets and provide an approach to assessing the impacts of hydroclimate and policy options on the WEF nexus, along with suggestions for adapting the agriculture and energy sectors to climate change.

  • Source Publication: Journal of Applied Meteorology and Climatology, 61, 1, 77-95, doi: 10.1175/JAMC-D-20-0260.1 Authors: Sobie, S.R. and T.Q. Murdock Publication Date: Jan 2022

    Information about snow water equivalent in southwestern British Columbia, Canada, is used for flood management, agriculture, fisheries, and water resource planning. This study evaluates whether a process-based, energy balance snow model supplied with high-resolution statistically downscaled temperature and precipitation data can effectively simulate snow water equivalent (SWE) in the mountainous terrain of this region. Daily values of SWE from 1951 to 2018 are simulated at 1-km resolution and evaluated using a reanalysis SWE product [Snow Data Assimilation System (SNODAS)], manual snow-survey measurements at 41 sites, and automated snow pillows at six locations in the study region. Simulated SWE matches observed interannual variability well (R2 > 0.8 for annual maximum SWE), but peak SWE biases of 20%–40% occur at some sites in the study domain, and higher biases occur where observed SWE is very low. Modeled SWE displays lower bias relative to SNODAS reanalysis at most manual survey locations. Future projections for the study area are produced using 12 downscaled climate model simulations and are used to illustrate the impacts of climate change on SWE at 1°, 2°, and 3°C of warming. Model results are used to quantify spring SWE changes at different elevations of the Whistler mountain ski resort and the sensitivity of annual peak SWE in the Metropolitan Vancouver municipal watersheds to moderate temperature increases. The results both illustrate the potential utility of a process-based snow model and identify areas where the input meteorological variables could be improved.

Pages