Publications Library

Land surface models have many parameters that have a spatially variable impact on model outputs. In applying these models, sensitivity analysis (SA) is sometimes performed as an initial step to select calibration parameters. As these models are applied to large domains, performing sensitivity analysis across the domain is computationally prohibitive. Here, using a Variable Infiltration Capacity model (VIC) deployment to a large domain as an example, we show that watershed classification based on climatic attributes and vegetation land cover helps to identify the spatial pattern of parameter sensitivity within the domain at a reduced cost. We evaluate the sensitivity of 44 VIC model parameters with regard to streamflow, evapotranspiration and snow water equivalent over 25 basins with a median size of 5078 km2. Basins are clustered based on their climatic and land cover attributes. Performance in transferring parameter sensitivity between basins of the same cluster is evaluated by the F1 score. Results show that two donor basins per cluster are sufficient to correctly identify sensitive parameters in a target basin, with F1 scores ranging between 0.66 (evapotranspiration) and 1 (snow water equivalent). While climatic attributes are sufficient to identify sensitive parameters for streamflow and evapotranspiration, including the vegetation class significantly improves skill in identifying sensitive parameters for the snow water equivalent. This work reveals that there is opportunity to leverage climate and land cover attributes to greatly increase the efficiency of parameter sensitivity analysis and facilitate more rapid deployment of land surface models over large spatial domains.

This report outlines a method for selecting a subset of earth system models (ESMs) from the Sixth Coupled Model Intercomparison Project (CMIP6) that is sufficiently representative of an ensemble of 26 models from CMIP6 for Canada and its subregions. The specific objective is to obtain a subset of reasonably independent ESMs that captures the overall range of projected change in a representative set of climate extremes (ETCCDI or Climdex) indices constructed from the ESM outputs. Projections are calculated for a future epoch corresponding to a global mean temperature change of 2 ℃ relative to 1971-2000, using results from two of the CMIP6 Shared Socioeconomic Pathways (SSPs), SSP2-4.5, and SSP5-8.5. The selection procedure is described below and representative subsets are provided for Canada and five of its subregions.

Climate change adaptation and relevant policy-making need reliable projections of future climate. Methods based on multi-model ensemble are generally considered as the most efficient way to achieve the goal. However, their efficiency varies and inter-comparison is a challenging task, as they use a variety of target variables, geographic regions, time periods, or model pools. Here, we construct and use a consistent framework to evaluate the performance of five ensemble-processing methods, i.e., multi-model ensemble mean (MME), rank-based weighting (RANK), reliability ensemble averaging (REA), climate model weighting by independence and performance (ClimWIP), and Bayesian model averaging (BMA). We investigate the annual mean temperature (Tav) and total precipitation (Prcptot) changes (relative to 1995–2014) over China and its seven subregions at 1.5 and 2 °C warming levels (relative to pre-industrial). All ensemble-processing methods perform better than MME, and achieve generally consistent results in terms of median values. But they show different results in terms of inter-model spread, served as a measure of uncertainty, and signal-to-noise ratio (SNR). ClimWIP is the most optimal method with its good performance in simulating current climate and in providing credible future projections. The uncertainty, measured by the range of 10th-90th percentiles, is reduced by about 30% for Tav, and 15% for Prcptot in China, with a certain variation among subregions. Based on ClimWIP, and averaged over whole China under 1.5/2 °C global warming levels, Tav increases by about 1.1/1.8 °C (relative to 1995–2014), while Prcptot increases by about 5.4%/11.2%, respectively. Reliability of projections is found dependent on investigated regions and indices. The projection for Tav is credible across all regions, as its SNR is generally larger than 2, while the SNR is lower than 1 for Prcptot over most regions under 1.5 °C warming. The largest warming is found in northeastern China, with increase of 1.3 (0.6-1.7)/2.0 (1.4-2.6) °C(ensemble’s median and range of the 10th–90th percentiles) under 1.5/2 °C warming, followed by northern and northwestern China. The smallest but the most robust warming is in southwestern China, with values exceeding 0.9 (0.6–1.1)/1.5 (1.1–1.7) °C. The most robust projection and largest increase is achieved in northwestern China for Prcptot, with increase of 9.1%(-1.6–24.7%)/17.9% (0.5–36.4%) under 1.5/2 °C warming. Followed by northern China, where the increase is 6.0%(-2.6–17.8%)/11.8% (2.4–25.1%), respectively. The precipitation projection is of large uncertainty in southwestern China, even with uncertain sign of variation. For the additional half-degree warming, Tav increases more than 0.5 °C throughout China. Almost all regions witness an increase of Prcptot, with the largest increase in northwestern China.

This issue of the PCIC Update contains the following stories: Climate Projections for the City of Terrace Released and A Mystery Gremlin Resolved! It also contains an update on the Pacific Climate Seminar Series, staff changes at PCIC and PCIC's most recent publications. The staff profile in this issue is on Tom Kunkel.

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.

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.

This issue of the PCIC Update contains the following stories: Correcting CMIP6 Model Output for Downscaling, Bilingual Design Value Explorer Announcement, and IPCC Summary for Policy Makers on Synthesis Report. It also contains an update on the Pacific Climate Seminar Series, staff changes at PCIC and PCIC's most recent publications. The staff profile in this issue is on Nina Nichols.

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.

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.