Abstract. The Canadian Earth System Model version 5 (CanESM5) is a global model developed to simulate historical climate change and variability, to make centennial-scale projections of future climate, and to produce initialized seasonal and decadal predictions. This paper describes the model components and their coupling, as well as various aspects of model development, including tuning, optimization, and a reproducibility strategy. We also document the stability of the model using a long control simulation, quantify the model's ability to reproduce large-scale features of the historical climate, and evaluate the response of the model to external forcing. CanESM5 is comprised of three-dimensional atmosphere (T63 spectral resolution equivalent roughly to 2.8∘) and ocean (nominally 1∘) general circulation models, a sea-ice model, a land surface scheme, and explicit land and ocean carbon cycle models. The model features relatively coarse resolution and high throughput, which facilitates the production of large ensembles. CanESM5 has a notably higher equilibrium climate sensitivity (5.6 K) than its predecessor, CanESM2 (3.7 K), which we briefly discuss, along with simulated changes over the historical period. CanESM5 simulations contribute to the Coupled Model Intercomparison Project phase 6 (CMIP6) and will be employed for climate science and service applications in Canada.

Peak runoff in streams and rivers of the western United States is strongly influenced by melting of accumulated mountain snowpack. A significant decline in this resource has a direct connection to streamflow, with substantial economic and societal impacts. Observations and reanalyses indicate that between the 1980s and 2000s, there was a 10–20% loss in the annual maximum amount of water contained in the region's snowpack. Here we show that this loss is consistent with results from a large ensemble of climate simulations forced with natural and anthropogenic changes, but is inconsistent with simulations forced by natural changes alone. A further loss of up to 60% is projected within the next 30 years. Uncertainties in loss estimates depend on the size and the rate of response to continued anthropogenic forcing and the magnitude and phasing of internal decadal variability. The projected losses have serious implications for the hydropower, municipal and agricultural sectors in the region.

[1] Evidence is presented that exchanges of water and energy between the vegetation and the atmosphere play an important role in east Asian and West African monsoon development and are among the most important mechanisms governing the development of the monsoon. The results were obtained by conducting simulations for five months of 1987 using a general circulation model (GCM) coupled with two different land surface parameterizations, with and without explicit vegetation representations, referred to as the GCM/vegetation and the GCM/soil, respectively. The two land surface models produced similar results at the planetary scale but substantial differences at regional scales, especially in the monsoon regions and some of the large continental areas. In the simulation with GCM/soil, the east Asian summer monsoon moisture transport and precipitation were too strong in the premonsoon season, and an important east Asian monsoon feature, the abrupt monsoon northward jump, was unclear. In the GCM/ vegetation simulation, the abrupt northward jump and other monsoon evolution processes were simulated, such as the large-scale turning of the low-level airflow during the early monsoon stage in both regions. With improved initial soil moisture and vegetation maps, the intensity and spatial distribution of the summer precipitation were also improved. The two land surface representations produced different longitudinal and latitudinal sensible heat gradients at the surface that, in turn, influenced the low-level temperature and pressure gradients, wind flow (through geostrophic balance), and moisture transport. It is suggested that the great east-west thermal gradient may contribute to the abrupt northward jump and the latitudinal heating gradient may contribute to the clockwise and counterclockwise turning of the low-level wind. The results showed that under unstable atmospheric conditions, not only low-frequency mean forcings from the land surface, such as monthly mean albedo, but also the perturbation processes of vegetation were important to the monsoon evolution, affecting its intensity, the spatial distribution of precipitation, and associated circulation at the continental scale.

A model evaluation approach is proposed in which weather and climate prediction models are analyzed along a Pacific Ocean cross section, from the stratocumulus regions off the coast of California, across the shallow convection dominated trade winds, to the deep convection regions of the ITCZ-the Global Energy and Water Cycle Experiment Cloud System Study/Working Group on Numerical Experimentation (GCSS/ WGNE) Pacific Cross-Section Intercomparison (GPCI). The main goal of GPCI is to evaluate and help understand and improve the representation of tropical and subtropical cloud processes in weather and climate prediction models. In this paper, a detailed analysis of cloud regime transitions along the cross section from the subtropics to the tropics for the season June-July-August of 1998 is presented. This GPCI study confirms many of the typical weather and climate prediction model problems in the representation of clouds: underestimation of clouds in the stratocumulus regime by most models with the corresponding consequences in terms of shortwave radiation biases; overestimation of clouds by the 40-yr ECMWF Re-Analysis (ERA-40) in the deep tropics (in particular) with the corresponding impact in the outgoing longwave radiation; large spread between the different models in terms of cloud cover, liquid water path and shortwave radiation; significant differences between the models in terms of vertical cross sections of cloud properties (in particular), vertical velocity, and relative humidity. An alternative analysis of cloud cover mean statistics is proposed where sharp gradients in cloud cover along the GPCI transect are taken into account. This analysis shows that the negative cloud bias of some models and ERA-40 in the stratocumulus regions [as compared to the first International Satellite Cloud Climatology Project (ISCCP)] is associated not only with lower values of cloud cover in these regimes, but also with a stratocumulus-to-cumulus transition that occurs too early along the trade wind Lagrangian trajectory. Histograms of cloud cover along the cross section differ significantly between models. Some models exhibit a quasi-bimodal structure with cloud cover being either very large (close to 100%) or very small, while other models show a more continuous transition. The ISCCP observations suggest that reality is in-between these two extreme examples. These different patterns reflect the diverse nature of the cloud, boundary layer, and convection parameterizations in the participating weather and climate prediction models.

A new approach of coordinated global and regional climate modeling is presented. It is applied to the Canadian Centre for Climate Modelling and Analysis Regional Climate Model (CanRCM4) and its parent global climate model CanESM2. CanRCM4 was developed specifically to downscale climate predictions and climate projections made by its parent global model. The close association of a regional climate model (RCM) with a parent global climate model (GCM) offers novel avenues of model development and application that are not typically available to independent regional climate modeling centers. For example, when CanRCM4 is driven by its parent model, driving information for all of its prognostic variables is available (including aerosols and chemical species), significantly improving the quality of their simulation. Additionally, CanRCM4 can be driven by its parent model for all downscaling applications by employing a spectral nudging procedure in CanESM2 designed to constrain its evolution to follow any large-scale driving data. Coordination offers benefit to the development of physical parameterizations and provides an objective means to evaluate the scalability of such parameterizations across a range of spatial resolutions. Finally, coordinating regional and global modeling efforts helps to highlight the importance of assessing RCMs’ value added relative to their driving global models. As a first step in this direction, a framework for identifying appreciable differences in RCM versus GCM climate change results is proposed and applied to CanRCM4 and CanESM2.