Publications Library

A detection and attribution analysis on the multidecadal trend in snow water equivalent (SWE) has been conducted in four river basins located in British Columbia (BC). Monthly output from a suite of 10 general circulation models (GCMs) that participated in phase 5 of the Coupled Model Intercomparison Project (CMIP5) is used, including 40 climate simulations with anthropogenic and natural forcing combined (ALL), 40 simulations with natural forcing alone (NAT), and approximately 4200 yr of preindustrial control simulations (CTL). This output was downscaled to 1/16° spatial resolution and daily temporal resolution to drive the Variable Infiltration Capacity hydrologic model (VIC). Observed (manual snow survey) and VIC-reconstructed SWE, which exhibit declines across BC, are projected onto the multimodel ensemble means of the VIC-simulated SWE based on the responses to different forcings using an optimal fingerprinting approach. Results of the detection and attribution analysis shows that these declines are attributable to the anthropogenic forcing, which is dominated by the effect of increases in greenhouse gas concentration, and that they are not caused by natural forcing due to volcanic activity and solar variability combined. Anthropogenic influence is detected in three of the four basins (Fraser, Columbia, and Campbell Rivers) based on the VIC-reconstructed SWE, and in all basins based on the manual snow survey records. The simulations underestimate the observed snowpack trends in the Columbia River basin, which has the highest mean elevation. Attribution is supported by the detection of human influence on the cold-season temperatures that drive the snowpack reductions. These results are robust to the use of different observed datasets and to the treatment of low-frequency variability effects.

Warming of the climate is now unequivocal. The water holding capacity of the atmosphere increases by about 7% per °C of warming, which in turn raises the expectation of more intense extreme rainfall events. Meeting the demand for robust projections for extreme short-duration rainfall is challenging, however, because of our poor understanding of its past and future behaviour. The characterization of past changes is severely limited by the availability of observational data. Climate models, including typical regional climate models, do not directly simulate all extreme rainfall producing processes, such as convection. Recently developed convection-permitting models better simulate extreme precipitation, but simulations are not yet widely available due to their computational cost, and they have their own uncertainties. Attention has thus been focused on precipitation–temperature relationships in the hope of obtaining more robust extreme precipitation projections that exploit higher confidence temperature projections. However, the observed precipitation–temperature scaling relationships have been established almost exclusively by linking precipitation extremes with day-to-day temperature variations. These scaling relationships do not appear to provide a reliable basis for projecting future precipitation extremes. Until better methods are available, the relationship of the atmosphere's water holding capacity with temperature provides better guidance for planners in the mid-latitudes, albeit with large uncertainties.

Presented by Norman Shippee at the American Meteorological Society's 97^th Annual Meeting in January, 2017:  CanSIPS reproduces the overall spatial pattern of cyclone track density found in ERA-Interim. Individual models of CanCM3 and CanCM4 over and under-represent the observed storm track, respectively. CanSIPS mean storm track and spread of track density is well represented, though an overall overprediction of storm track density is evident across the North Pacific, particularly from 45 - 60 deg N. This overprediction of storm track density is present in all components of the CanSIPS system. Moderate to strong correlation between CanSIPS mean storm track density and ERA-Interim track density are found across the primary storm track in the North Pacific and in the primary formation regions for Atlantic cyclones affecting North America (Gulf of Mexico and Cape Hatteras) Most CanSIPS bias in the North Pacific is centralized in the exit regions of the Pacific storm track, localized to Gulf of Alaska and coastal BC.

Arctic sea ice extent (SIE) has decreased over recent decades, with record-setting minimum events in 2007 and again in 2012. A question of interest across many disciplines concerns the extent to which such extreme events can be attributed to anthropogenic influences. First, a detection and attribution analysis is performed for trends in SIE anomalies over the observed period. The main objective of this study is an event attribution analysis for extreme minimum events in Arctic SIE. Although focus is placed on the 2012 event, the results are generalized to extreme events of other magnitudes, including both past and potential future extremes. Several ensembles of model responses are used, including two single-model large ensembles. Using several different metrics to define the events in question, it is shown that an extreme SIE minimum of the magnitude seen in 2012 is consistent with a scenario including anthropogenic influence and is extremely unlikely in a scenario excluding anthropogenic influence. Hence, the 2012 Arctic sea ice minimum provides a counterexample to the often-quoted idea that individual extreme events cannot be attributed to human influence.

The Fraser River basin (FRB) of British Columbia is one of the largest and most important watersheds in Western North America, and is home to a rich diversity of biological species and economic assets that depend implicitly upon its extensive riverine habitats. The hydrology of the FRB is dominated by snow accumulation and melt processes, leading to a prominent annual peak streamflow invariably occurring in June-July. However, while annual peak daily streamflow (APDF) during the spring freshet in the FRB is historically well correlated with basin-averaged, annual maximum snow water equivalent (SWEmax), there are numerous occurrences of anomalously large APDF in below- or near-normal SWEmax years, some of which have resulted in damaging floods in the region. An imperfect understanding of which other climatic factors contribute to these anomalously large APDFs complicates future projections of streamflow magnitude and frequency.

We employ the Variable Infiltration Capacity (VIC) process-based hydrological model driven by both observations and an ensemble of CMIP3 climate models in an attempt to discover the proximate causes of anomalous APDF events in the FRB. At several hydrometric stations representing a range of elevations, the relative importance of a set of predictors characterizing the magnitude and timing of rainfall, snowfall, and temperature is examined within a regression framework. The results indicate that next to the magnitude of SWEmax, the rate of warming subsequent to the date of SWEmax is the most influential variable for predicting APDF magnitudes in the lower FRB. Finally, the role of large-scale climate modes of variability for APDF magnitude and timing in the basin will be briefly discussed.

PDRMIP investigates the role of various drivers of climate change for mean and extreme precipitation changes, based on multiple climate model output and energy budget analyses.

As the global temperature increases with changing climate, precipitation rates and patterns are affected through a wide range of physical mechanisms. The globally averaged intensity of extreme precipitation also changes more rapidly than the globally averaged precipitation rate. While some aspects of the regional variation in precipitation predicted by climate models appear robust, there is still a large degree of inter-model differences unaccounted for. Individual drivers of climate change initially alter the energy budget of the atmosphere leading to distinct rapid adjustments involving changes in precipitation. Differences in how these rapid adjustment processes manifest themselves within models are likely to explain a large fraction of the present model spread and needs better quantifications to improve precipitation predictions. Here, we introduce the Precipitation Driver and Response Model Intercomparison Project (PDRMIP), where a set of idealized experiments designed to understand the role of different climate forcing mechanisms were performed by a large set of climate models. PDRMIP focuses on understanding how precipitation changes relating to rapid adjustments and slower responses to climate forcings are represented across models. Initial results show that rapid adjustments account for large regional differences in hydrological sensitivity across multiple drivers. The PDRMIP results are expected to dramatically improve our understanding of the causes of the present diversity in future climate projections.

Gridded statistical downscaling methods are the main means of preparing climate model data to drive distributed hydrological models. Past work on the validation of climate downscaling methods has focused on temperature and precipitation, with less attention paid to the ultimate outputs from hydrological models. Also, as attention shifts towards projections of extreme events, downscaling comparisons now commonly assess methods in terms of climate extremes, but hydrologic extremes are less well explored. Here, we test the ability of gridded downscaling models to replicate historical properties of climate and hydrologic extremes, as measured in terms of temporal sequencing (i.e. correlation tests) and distributional properties (i.e. tests for equality of probability distributions). Outputs from seven downscaling methods – bias correction constructed analogues (BCCA), double BCCA (DBCCA), BCCA with quantile mapping reordering (BCCAQ), bias correction spatial disaggregation (BCSD), BCSD using minimum/maximum temperature (BCSDX), the climate imprint delta method (CI), and bias corrected CI (BCCI) – are used to drive the Variable Infiltration Capacity (VIC) model over the snow-dominated Peace River basin, British Columbia. Outputs are tested using split-sample validation on 26 climate extremes indices (ClimDEX) and two hydrologic extremes indices (3-day peak flow and 7-day peak flow). To characterize observational uncertainty, four atmospheric reanalyses are used as climate model surrogates and two gridded observational data sets are used as downscaling target data. The skill of the downscaling methods generally depended on reanalysis and gridded observational data set. However, CI failed to reproduce the distribution and BCSD and BCSDX the timing of winter 7-day low-flow events, regardless of reanalysis or observational data set. Overall, DBCCA passed the greatest number of tests for the ClimDEX indices, while BCCAQ, which is designed to more accurately resolve event-scale spatial gradients, passed the greatest number of tests for hydrologic extremes. Non-stationarity in the observational/reanalysis data sets complicated the evaluation of downscaling performance. Comparing temporal homogeneity and trends in climate indices and hydrological model outputs calculated from downscaled reanalyses and gridded observations was useful for diagnosing the reliability of the various historical data sets. We recommend that such analyses be conducted before such data are used to construct future hydro-climatic change scenarios.

Regional climate models (RCMs) are the primary source of high-resolution climate projections, and it is of crucial importance to evaluate their ability to simulate extreme events under current climate conditions. Many extreme events are influenced by circulation features that occur outside, or on the edges of, RCM domains. Thus, it is of interest to know whether such dynamically controlled aspects of extremes are well represented by RCMs. This study assesses the relationship between upstream blocking and cold temperature extremes over North America in observations, reanalysis products (ERA-Interim and NARR), and RCMs (CanRCM4, CRCM5, HIRHAM5, and RCA4). Generalized extreme value distributions were fitted to winter minimum temperature (TNn) incorporating blocking frequency (BF) as a covariate, which is shown to have a significant influence on TNn. The magnitude of blocking influence in the RCMs is consistent with observations, but the spatial extent varies. CRCM5 and HIRHAM5 reproduce the pattern of influence best compared to observations. CanRCM4 and RCA4 capture the influence of blocking in British Columbia and the northeastern United States, but the extension of influence that is seen in observations and reanalysis into the southern United States is not evident. The difference in the 20-yr return value (20RV) of TNn between high and low BF in the Pacific Ocean indicates that blocking is associated with a decrease of up to 15°C in the 20RV over the majority of the United States and in western Canada. In northern North America the difference in the 20RV is positive as blocking is associated with warmer extreme cold temperatures. The 20RVs are generally simulated well by the RCMs.

The Indo-Pacific warm pool (IPWP) has warmed and grown substantially during the past century. The IPWP is Earth’s largest region of warm sea surface temperatures (SSTs), has the highest rainfall, and is fundamental to global atmospheric circulation and hydrological cycle. The region has also experienced the world’s highest rates of sea-level rise in recent decades, indicating large increases in ocean heat content and leading to substantial impacts on small island states in the region. Previous studies have considered mechanisms for the basin-scale ocean warming, but not the causes of the observed IPWP expansion, where expansion in the Indian Ocean has far exceeded that in the Pacific Ocean. We identify human and natural contributions to the observed IPWP changes since the 1950s by comparing observations with climate model simulations using an optimal fingerprinting technique. Greenhouse gas forcing is found to be the dominant cause of the observed increases in IPWP intensity and size, whereas natural fluctuations associated with the Pacific Decadal Oscillation have played a smaller yet significant role. Further, we show that the shape and impact of human-induced IPWP growth could be asymmetric between the Indian and Pacific basins, the causes of which remain uncertain. Human-induced changes in the IPWP have important implications for understanding and projecting related changes in monsoonal rainfall, and frequency or intensity of tropical storms, which have profound socioeconomic consequences.