Publications Library

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.

We evaluate the capacity of a regional climate model to represent observed extreme temperature and precipitation events and also examine the impact of increased resolution, in an effort to identify added value in this respect. Two climate simulations of western Canada (WCan) were conducted with the Canadian Regional Climate Model (version 4) at 15 (CRCM15) and 45 km (CRCM45) horizontal resolution driven at the lateral boundaries by data from the European Centre for Medium-range Weather Forecasts (ECMWF) 40-year Reanalysis (ERA-40) for the period 1973–1995. The simulations were evaluated using the spline-interpolated dataset ANUSPLIN, a daily observational gridded surface temperature and precipitation product with a nominal resolution of approximately 10 km. We examine a range of climate extremes, comprising the 10th and 90th percentiles of daily maximum (TX) and minimum (TN) temperatures, the 90th percentile of daily precipitation (PR90), and the 27 core Climate Daily Extremes (CLIMDEX) indices.

Both simulations exhibit cold biases compared with observations over WCan, with the bias exacerbated at higher resolution, suggesting little added value for temperature overall. There are instances, however, of regional improvement in the spatial pattern of temperature extremes at the higher resolution of CRCM15 (e.g., the CLIMDEX index for the annual number of days when TX > 25°C). The high-resolution simulations also reveal similarly localized features in precipitation (e.g., rain shadows) that are not resolved at the 45 km resolution. With regard to precipitation extremes, although both simulations generally display wet biases, CRCM15 features a reduced bias in PR90 in all seasons except winter. This improvement occurs despite the fact that spatial and interannual variability of PR90 in CRCM15 is significantly overestimated relative to both CRCM45 and ANUSPLIN. We posit that these characteristics are the result of demonstrable differences between corresponding topographical datasets used in the gridded observations and CRCM, the resulting errors propagated to physical variables tied to elevation and the beneficial effect of subsequent spatial averaging. Because topographical input is often discordant between simulations and gridded observations, it is argued that a limited form of spatial averaging may contribute added value beyond that which has already been noted in previous studies with respect to small-scale climate variability.

This study investigates a physical basis for heterogeneity in hydrological changes, which suggests a greater detectability in wet than dry regions. Wet regions are those where atmospheric demand is less than precipitation (energy limited), and dry regions are those where atmospheric demand is greater than precipitation (water limited). Long-term streamflow trends in western North America and an analysis of Coupled Model Intercomparison Project Phase 5 (CMIP5) climate models at global scales show geographically heterogeneous detectability of hydrological changes. We apply the Budyko framework and state-of-the-art climate model data from CMIP5 to quantify the sensitivity and detectability of terrestrial hydrological changes. The Budyko framework quantifies the partitioning of precipitation into evapotranspiration and runoff components. We find that the terrestrial hydrological sensitivity is 3 times greater in regions where the hydrological cycle is energy limited rather than water limited. This additional source (the terrestrial part) contributes to 30–40% greater detectability in energy-limited regions. We also quantified the contribution of changes in the catchment efficiency parameter that oppose the effects of increasing evaporative demand in global warming scenarios. Incorporating changes to the catchment efficiency parameter in the Budyko framework reduces dry biases in global runoff change projections by 88% in the 21st century.

The need for future projections of extremes is growing, particularly as users planning to adapt to climate change continue to experience record-breaking events (Figure 1). Decision-making demands that such projections possess high spatial resolution. Downscaling has been carried out for Canada by the Pacic Climate Impacts Consortium for the newest Global Climate Model (GCM) and Regional Climate Model (RCM) projections.

For terrestrial British Columbia, precipitation averages and extremes can be simulated more accurately within individual regions by using gridded downscaling to increase the resolution of regional climate models. In locations where the difference between observations and RCMs is large, bias correction tends to inflate the magnitude of projected extremes. Differences in projections between the RCMs and downscaled simulations can be large enough to affect the PIEVC Risk Assessment process, leading to higher risk values. Future work will focus on correcting the inflation of extremes in the downscaling method and extending the analysis to additional regions.