Informing BC Stakeholders

You are here

Publications Library

  • 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.

  • 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.

  • Source Publication: Atmosphere-Ocean, 59, 4-5, 269-284, doi:10.1080/07055900.2021.2011103 Authors: Sobie, S.R., Zwiers, F.W., and C.L. Curry Publication Date: Dec 2021

    Recent studies have identified stronger warming in the latest generation of climate model simulations globally, and the same is true for projected changes in Canada. This study examines differences for Canada and six sub-regions between simulations from the latest Sixth Coupled Model Intercomparison Project (CMIP6) and its predecessor CMIP5. Ensembles from both experiments are assessed using a set of derived indices calculated from daily precipitation and temperature, with projections compared at fixed future time intervals and fixed levels of global temperature change. For changes calculated at fixed time intervals most temperature indices display higher projected changes in CMIP6 than CMIP5 for most sub-regions, while greater precipitation changes in CMIP6 occur mainly in extreme precipitation indices. When future projections are calculated at fixed levels of global average temperature increase, the size and spread of differences for future projected changes between CMIP6 and CMIP5 are substantially reduced for most indices. Temperature scaling behaviour, or the regional response to increasing global temperatures, is similar in both ensembles, with annual temperature anomalies for Canada and its sub-regions increasing at between 1.5 and 2.5 times the rate of increase globally, depending on the region. The CMIP6 ensemble projections exhibit modestly stronger scaling behaviour for temperature anomalies in northern Canada, as well as for certain indices of moderate and extreme events. Such temperature scaling differences persist even if anomalously warm CMIP6 global climate models are omitted. Comparing the mean and variance of future projections for Canada in CMIP5 and CMIP6 simulations from the same modelling centre suggests CMIP6 models are significantly warmer in Canada than CMIP5 models at the same level of forcing, with some evidence that internal temperature variability in CMIP6 is reduced compared with CMIP5.

  • Source Publication: Geophysical Research Letters, 48, doi:10.1029/2021GL095500 Authors: Li, C., Z. Wang, F. Zwiers and X. Zhang Publication Date: Dec 2021

    The regression-based optimal fingerprinting is a key tool for quantifying human climate influence. Most studies over the past decade used Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations, limiting fingerprinting regression configuration options. The CMIP6 Detection and Attribution Model Intercomparison Project (DAMIP) provides several types of individual forcing simulations and thus greater configuration flexibility. To avoid overfitting the limited observational data, we suggest that a DAMIP-based perfect model study is first used to best configure the fingerprinting regression prior to its application to observations. We find that a regression using all-forcing, aerosol-only, and natural-only simulations is an overall best option for constraining human-induced global terrestrial warming, which differs from choices commonly made previously. Applying this configuration to observations, we estimate that of the observed terrestrial warming of ∼1.5°C between 1850–1900 and 2011–2020, anthropogenic greenhouse gases contributed 1.4 to 2.3°C, offset by aerosol cooling of 0.2 to 1.2°C.

  • Source Publication: Geophysical Research Letters, doi:10.1029/2021GL095500 Authors: Li, C., Z. Wang, F. Zwiers and X. Zhang Publication Date: Nov 2021

    The regression-based optimal fingerprinting is a key tool for quantifying human climate influence. Most studies over the past decade used Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations, limiting fingerprinting regression configuration options. The CMIP6 Detection and Attribution Model Intercomparison Project (DAMIP) provides several types of individual forcing simulations and thus greater configuration flexibility. To avoid overfitting the limited observational data, we suggest that a DAMIP-based perfect model study is first used to best configure the fingerprinting regression prior to its application to observations. We find that a regression using all-forcing, aerosol-only, and natural-only simulations is an overall best option for constraining human-induced global terrestrial warming, which differs from choices commonly made previously. Applying this configuration to observations, we estimate that of the observed terrestrial warming of ∼1.5°C between 1850–1900 and 2011–2020, anthropogenic greenhouse gases contributed 1.4 to 2.3°C, offset by aerosol cooling of 0.2 to 1.2°C.

  • Source Publication: Progress in Oceanography, 198, 102659. doi:10.1016/j.pocean.2021.102659 Authors: Heneghan, R. F. et al. (T.C. Tai is 21st coauthor) Publication Date: Oct 2021

    Climate change is warming the ocean and impacting lower trophic level (LTL) organisms. Marine ecosystem models can provide estimates of how these changes will propagate to larger animals and impact societal services such as fisheries, but at present these estimates vary widely. A better understanding of what drives this inter-model variation will improve our ability to project fisheries and other ecosystem services into the future, while also helping to identify uncertainties in process understanding. Here, we explore the mechanisms that underlie the diversity of responses to changes in temperature and LTLs in eight global marine ecosystem models from the Fisheries and Marine Ecosystem Model Intercomparison Project (FishMIP). Temperature and LTL impacts on total consumer biomass and ecosystem structure (defined as the relative change of small and large organism biomass) were isolated using a comparative experimental protocol. Total model biomass varied between −35% to +3% in response to warming, and -17% to +15% in response to LTL changes. There was little consensus about the spatial redistribution of biomass or changes in the balance between small and large organisms (ecosystem structure) in response to warming, an LTL impacts on total consumer biomass varied depending on the choice of LTL forcing terms. Overall, climate change impacts on consumer biomass and ecosystem structure are well approximated by the sum of temperature and LTL impacts, indicating an absence of nonlinear interaction between the models’ drivers. Our results highlight a lack of theoretical clarity about how to represent fundamental ecological mechanisms, most importantly how temperature impacts scale from individual to ecosystem level, and the need to better understand the two-way coupling between LTL organisms and consumers. We finish by identifying future research needs to strengthen global marine ecosystem modelling and improve projections of climate change impacts.

  • Source Publication: Science Advances, 7, eabh0895 Authors: Cheung, W. W. L., T.L. Frölicher, V.W.Y. Lam, M. Oyinlola, G. Reygondeau, U.R. Sumaila, T.C. Tai, L.C.L Teh and C.C.C. Wabnitz Publication Date: Oct 2021

    Extreme temperature events have occurred in all ocean basins in the past two decades with detrimental impacts on marine biodiversity, ecosystem functions, and services. However, global impacts of temperature extremes on fish stocks, fisheries, and dependent people have not been quantified. Using an integrated climate-biodiversity-fisheries-economic impact model, we project that, on average, when an annual high temperature extreme occurs in an exclusive economic zone, 77% of exploited fishes and invertebrates therein will decrease in biomass while maximum catch potential will drop by 6%, adding to the decadal-scale mean impacts under climate change. The net negative impacts of high temperature extremes on fish stocks are projected to cause losses in fisheries revenues and livelihoods in most maritime countries, creating shocks to fisheries social-ecological systems particularly in climate-vulnerable areas. Our study highlights the need for rapid adaptation responses to extreme temperatures in addition to carbon mitigation to support sustainable ocean development.

  • Source Publication: Weather and Climate Extremes, 34, 100388, doi:10.1016/j.wace.2021.100388. Authors: Ben Alaya, M.A., F.W. Zwiers and X. Zhang Publication Date: Oct 2021

    The uniform risk engineering practices that are increasingly being adopted for structural design require estimates of the extreme wind loads with very low annual probabilities of exceedance, corresponding to return periods of up to 3000-years in some cases. These estimates are necessarily based on observational wind data that typically spans only a few decades. The estimates are therefore affected by both large sampling uncertainty and, potentially, non-negligible biases. Design practices that aim to meet mandated structural reliability criteria take the sampling uncertainty of long period wind speed or wind pressure estimates into account, but reliability could be compromised if estimates are also biased. In many circumstances, estimates are obtained by fitting an extreme value distribution to annual maximum wind speed observed over a few decades. A key assumption implicit in doing so is that wind speed annual maxima are max-stable. Departures from max-stability can exacerbate the uncertainty of long-period return level estimates by inducing systematic estimation bias as well. Observational records, however, are generally too short to assess max-stability. We therefore use wind speed data from a large (50-member) ensemble of CanRCM4 historical simulations over North America to assess whether wind speed annual maxima are max-stable. While results are generally reassuring at the continental scale, disquieting evidence of a lack of max-stability is often found in the central and southern parts of the continent. Results show that when annual maximum wind speeds are not max-stable, long period return level extreme wind speeds tend to be underestimated, which would compromise reliability if used to design infrastructure such as tall buildings and towers.

Pages