The Canadian Atmospheric Model version 5 (CanAM5.0.3)

Abstract

The Canadian Atmospheric Model version 5 (CanAM5) is the component of Canadian Earth System Model version 5 (CanESM5) which models atmospheric processes and coupling of the atmosphere with land and lake models. Described in this paper are the main features of CanAM5, with a focus on changes relative to the last major scientific version of the model (CanAM4). These changes are mostly related to improvements in radiative transfer, clouds, and aerosol parameterizations, as well as a major upgrade of the land surface and land carbon cycle models and addition of a small lake model. In addition to changes to parameterizations and models, changes in the adjustable parameters between CanAM4 and CanAM5 are documented. Finally, the mean climatology simulated by CanAM5 for the present day is evaluated against observations and compared with that simulated by CanAM4. Although many of the aspects of the simulated climate are similar between CanAM4 and CanAM5, there is a reduction in precipitation and temperature biases over the Amazonian basin, global cloud fraction biases, and solar and thermal cloud radiative effects, all of which are improvements relative to observations.

1 Introduction

The fifth version of the Canadian Atmospheric Model (CanAM5; for a list of acronyms, initialisms, and abbreviations, see Table ) is a major component model in the Canadian Earth System Model version 5 (CanESM5) , modelling atmospheric processes and coupling of the atmosphere with land and lake models. Both CanAM5 and CanESM5 are models developed by the Canadian Centre for Climate Modelling and Analysis (CCCma) to simulate climate to improve understanding and make predictions and projections of future climate. CanAM5 is the result of several years of development on its last major scientific version, CanAM4 , which was the atmospheric component of CanESM2 used for the fifth Coupled Model Intercomparison Project (CMIP5) . Between CanAM4 and CanAM5, there were numerous changes to radiative transfer, cloud, and aerosol parameterizations, in addition to a major upgrade of the land surface model and addition of a model of unresolved subgrid-scale lakes.

CanESM5 and CanAM5 were the basis for CCCma's contribution to CMIP6 , which included a number of experiments to better understand and characterize cloud feedbacks and radiative forcings. The mean state and response of CanESM5 to external forcing in fully coupled simulations are documented in . In this paper, the focus is on the ability of CanAM5 to simulate the historical climate for simulations in which sea surface temperature and sea ice concentration are prescribed from observations. Since there are several changes in CanAM5 relative to CanAM4, most of the evaluation with observations is performed with both CanAM5 and CanAM4 to highlight changes in the simulated climate.

In Sects.  and , changes in atmospheric and surface processes between CanAM4 and CanAM5 are summarized. The process used to tune CanAM5 is discussed in Sect. , and the values used for adjustable parameters are presented. The details of experiments used to evaluate CanAM5 are presented in Sect. , and Sect.  presents an analysis of climatological features of CanAM5. Finally, in Sect. , we conclude with a brief summary and discussion of the main results of this study.

2 Atmospheric processes

This section summarizes atmospheric parameterizations in CanAM5. As most parameterizations are described and documented in detail for CanAM4 we focus on the major changes between CanAM5 and CanAM4.

The horizontal resolution of CanAM5 is identical to CanAM4 and is defined by triangular truncation at a total wavenumber of 63 (i.e., T63). The model employs a double spectral transform allowing the physical tendencies to be evaluated on a reduced “linear” T63 Gaussian grid of dimensions $\mathrm{128}\times \mathrm{64}$, which corresponds to $\sim \mathrm{2.8}$${}^{\circ }$. The number of vertical levels in CanAM5 has increased from 35 to 49. The 49 levels are used with layer thicknesses that increase monotonically from approximately 100 m at the surface to 2 km at $\sim \mathrm{1}$ hPa, which is the upper bound of the vertical domain. The additional 14 layers in CanAM5 have been added to the upper troposphere and lower stratosphere to match those employed by the Canadian Middle Atmosphere Model .

2.1 Radiation

Radiative transfer in CanAM5 includes the new specification of optical properties for the surface, cloud, and aerosol, in addition to the computation of radiative fluxes accounting for subgrid-scale surface variability.

The parameterized absorption by gases uses a correlated $k$-distribution model that is mostly unchanged from CanAM4, using the same wavelength intervals and quadrature points . A significant modification is the addition of a solar water vapour continuum , which resulted in improved simulation of absorption at solar wavelengths . The single-scattering properties of the ice clouds in CanAM4 were parameterized for thermal and solar wavelengths under the assumption that all the ice particles are hexagonal prisms . In CanAM5 the optical properties of ice particles use the parameterization of , assuming a mixture of ice habits that is based on spaceborne observations and assuming a moderately rough surface, which is found to improve retrievals . The single-scattering properties of pure liquid clouds remain the same but can be perturbed to account for internally mixed black carbon , allowing simulation of the semi-indirect effect.

Aerosol optical properties in CanAM5 use updated single-scattering properties as well as an improved method to mix aerosol optics. The single-scattering properties for organic carbon use the refractive index from HITRAN 2012 and the properties of black carbon from . Instead of externally mixing aerosols, it is assumed that sulfate and the hydrophilic components of black carbon and organic carbon are internally mixed . The refractive index of the internally mixed aerosol is computed based on the fraction, effective radius, and effective variance of each component aerosol, as well as relative humidity, which is used to compute hydrophilic growth.

The ocean optical properties are also changed in CanAM5. In CanAM4, the whitecap albedo was wavelength-invariant with a value of 0.3. In CanAM5, each wavelength interval in the solar radiative transfer model uses a different albedo (0.216, 0.134, 0.044, 0.005) based on . The parameterization of ocean albedo is similar to that in CanAM4, but the contents of the lookup table have been updated and now include a dependence on the solar zenith angle and partitioning of the incident downwelling solar radiation into direct and diffuse components . This partitioning is estimated using the vertically integrated aerosol and cloud optical depth. In CanAM4, the ocean albedo was computed using as input optical thickness and solar zenith angles from the last radiative transfer time step, which is 1 h earlier. To improve the consistency of the ocean albedo calculation and the radiative transfer calculations, especially near sunrise and sunset, in CanAM5 the ocean albedo is calculated using cloud and aerosol information from the previous dynamical time step, 15 min earlier, and the solar zenith angle from the current time step.

Over land, the dry bare-soil albedo in CanAM4 was set to a global constant combined with a parameterization to account for the effect of surface wetting . The use of a globally constant bare-soil albedo resulted in regional biases for clear-sky albedo at the top of atmosphere, especially over the Sahara and Australian interior. The constant albedo was replaced with a regionally varying soil colormap and associated albedos . These new location-dependent bare-soil albedo maps greatly reduced biases in clear-sky albedo relative to observations. In addition to the bare soil, the albedo and emissivity of snow and sea ice were also updated. The albedo of snow on land and sea ice in CanAM5 is computed using a lookup table accounting for snowpack properties. These include snow water equivalent, snow grain size, and black carbon simulated by the land surface model . Snow present on land or sea ice uses a single, wavelength-invariant emissivity, which was reduced from 1.0 to 0.97 . Similarly, the emissivity of sea ice was reduced from 1.0 to 0.97 to be consistent with the sea ice model used in CanESM5 .

Within a CanAM5 grid box, there can be multiple surface types including land, lake, ocean, and sea ice (Sect. ). When coupling the atmosphere with ocean, sea ice, and land models, it is necessary to have surface radiative fluxes that are consistent with each surface type. While it is possible to partition the grid-mean radiative fluxes at the surface using the tiled albedo, emissivity, and temperature, in CanAM5 radiative flux profiles are instead computed for each surface type and then averaged to a grid mean. It is assumed that the same atmosphere is present over each tile. Although this approach causes a small increase ($<\mathrm{5}$ %) in the total time for global radiative transfer calculations in CanAM5, it maintains consistency between the surface and the atmosphere. The modest increase in computational time is possible because multiple surface tiles are only present in a portion of the CanAM5 grid boxes, e.g., there is only one surface type over sea-ice-free ocean.

2.2 Aerosols and chemistry

The types of natural and anthropogenic aerosols in CanAM5 include sulfate, black and organic carbon, sea salt, and mineral dust, similar to CanAM4 . Parameterizations for aerosol emissions and transport, gas-phase and aqueous-phase chemistry, and dry and wet deposition account for interactions with simulated meteorology. Natural aerosol species are represented in the model using prognostic emission fluxes. In particular, a particle-size-dependent emission scheme is used to account for aeolian erosion in arid and semi-arid regions . Sea salt concentrations in two size modes are parameterized as a function of the wind speed near the surface of the ocean . Dimethyl sulfide emissions are predicted using specified climatological concentrations in the surface ocean . Sulfur oxidation in the gas and aqueous phases is simulated using specified climatological oxidant concentrations from CMAM20 .

The chemistry parameterizations in CanAM5 are unchanged relative to CanAM4 with the exception of stratospheric water vapour which can be produced by methane oxidation using a parameterization based on that described in .

2.3 Clouds

The parameterization of clouds and cloud microphysics in CanAM5 is mostly the same as in CanAM4 . Like CanAM4 and other global climate models, CanAM5 continues to employ bulk cloud microphysical parameterizations which depend on mean water content and other moments of the droplet size distribution.

Autoconversion in liquid clouds is parameterized in CanAM4 using , but this has been replaced in CanAM5 with the autoconversion parameterization of , which is a modified version of the parameterization by , to simulate the collision and coalescence of cloud droplets. For convenience we reproduce here the main equations for and , since the adjustment of parameters in the equations is discussed in Sect. .

Autoconversion in CanAM4 is parameterized as

1$${\displaystyle \frac{\partial q_{\mathrm{r}}}{\partial t}}=A\stackrel{\mathrm{‾}}{q_{\mathrm{L}}^{\mathrm{2.47}}}{N}_{\mathrm{d}}^{-\mathrm{1.79}}{\mathit{\rho }}^{-\mathrm{1.47}},$$ where $A$ is a constant, $q_{\mathrm{r}}$ is the rainwater content (in kg m${}^{-\mathrm{3}}$), $q_{\mathrm{L}}$ is the liquid cloud water content (in kg m${}^{-\mathrm{3}}$), $N_{\mathrm{d}}$ is the cloud droplet number concentration (in m${}^{-\mathrm{3}}$), and $\mathit{\rho }$ is air density (in kg m${}^{-\mathrm{3}}$). In CanAM5, the autoconversion is parameterized as 2$${\displaystyle \frac{\partial q_{\mathrm{r}}}{\partial t}}=E\stackrel{\mathrm{‾}}{q_{\mathrm{L}}^{\mathrm{3}}}{N}_{\mathrm{d}}^{-\mathrm{1}}H(R_{\mathrm{6}}-R_{\mathrm{6}C}),$$ where $q_{\mathrm{r}}$ is the rainwater content (in kg m${}^{-\mathrm{3}}$), $q_{\mathrm{L}}$ is the liquid cloud water content (in kg m${}^{-\mathrm{3}}$), $N_{\mathrm{d}}$ is the droplet number concentration (in m${}^{-\mathrm{3}}$), and $H\left(\right)$ is the Heaviside function. Additionally, $E=\mathrm{1.3}\times {\mathrm{10}}^{\mathrm{9}}{\mathit{\beta }}_{\mathrm{6}}^{\mathrm{6}}$, $R_{\mathrm{6}C}=\mathrm{7.5}/\left({\stackrel{\mathrm{‾}}{q_{\mathrm{L}}}}^{\mathrm{1}/\mathrm{6}}{R}_{\mathrm{6}}^{\mathrm{1}/\mathrm{2}}\right)$, $R_{\mathrm{6}}={\mathit{\beta }}_{\mathrm{6}}{r}_{\mathrm{v}}$, and ${\mathit{\beta }}_{\mathrm{6}}=\left[\right(r_{\mathrm{v}}+\mathrm{3})/r_{\mathrm{v}}]$, where $r_{\mathrm{v}}$ is the mean volume radius (in $\mathrm{\mu }\mathrm{m}$). In order to account for the impacts of subgrid-scale variability in cloud liquid water content, the statistical cloud scheme in CanAM5 is used to determine the mean value of $q_{\mathrm{L}}^{\mathrm{3}}$, indicated by the bar in Eq. ().

Along with the new autoconversion parameterization, CanAM5 now accounts for indirect impacts of aerosols on cloud liquid water content and lifetime, i.e., the second aerosol indirect effect . This effect was not active in CanAM4, since it used a constant cloud droplet number concentration of 50 cm${}^{-\mathrm{3}}$ in Eq. () . Given the uncertainty of applying the parameterizations at high altitudes, cloud processes are limited to pressures greater than 10 hPa.

3 Surface processes

There were three substantial changes to the treatment of surface processes: a major upgrade of the land surface including tighter integration of the land carbon cycle model, a small lake model, and tiling to accommodate fractional land in a grid box.

3.1 CLASS-CTEM

The land component of CanAM5 is represented by the Canadian Land Surface Scheme (CLASS) and the Canadian Terrestrial Ecosystem Model (CTEM), which model physical and biogeochemical processes, respectively. CanAM5 uses version 3.6 of CLASS, which models the energy and moisture fluxes at the air–land surface interface .

Compared to its predecessor used in CanAM4 (CLASS 2.7; ) there are several major structural improvements in version 3.6 of CLASS. These include optional implementation of a user-specified number of soil layers rather than the previous hard-coded three layers, as well as the capability of supporting a mosaic of vegetation, soil, water, or ice tiles within grid cells. The capability of modelling fully organic soils has been introduced, with hydraulic properties assigned on the basis of the work of . The thermal conductivities of the organic and mineral soil layers are determined following and . The wet and dry albedos of the mineral soil are assigned based on a global soil reflectivity index described in and . Organic soil albedos are assigned following . The bare-soil surface evaporation efficiency parameter is calculated using a relation presented by . Empirical corrections are applied to the saturated hydraulic conductivity of each soil layer to take into account the viscosity of water at the layer temperature and the presence of ice . The field capacity of the lowest permeable soil layer and the baseflow at the bottom of the layer are obtained using relations derived from .

A new option is provided to model snowpack albedo and transmissivity in four wavelength intervals instead of two intervals. The thermal conductivity of snow is obtained from the snow density using a relationship derived by . The fresh snow density is calculated as an empirical function of the air temperature, using relations developed by and . The maximum snowpack density is calculated as a function of snow depth following . The amount of snowfall intercepted by vegetation and the unloading rate of intercepted snow are calculated following . The canopy interception capacity for snow is determined from the plant area index and the fresh snow density as described in and . Albedos of snow-covered vegetation canopies are set following . The sensible and latent heat fluxes between the vegetation, the underlying ground, and the overlying atmosphere are evaluated based on the analysis of , which incorporates an explicit treatment of the canopy air space.

Although CLASS 3.6 can represent vegetation as a mosaic, the composite approach of representing different vegetation types in a grid cell is employed in CanAM5. This implies that area-weighted grid-mean structural attributes of different vegetation types are used in energy and water balance calculations. The number of soil and bedrock layers remains three, the same as in CanAM4, with the first and second soil layers being 0.1 and 0.25 m thick. The maximum thickness of the permeable soil for the third layer is 3.75 m but varies geographically depending on the permeable soil depth specified following .

CTEM models vegetation as a dynamic component of the climate system and provides structural attributes of vegetation to CLASS for use in its physics calculations . These include leaf area index, vegetation height, rooting depth and distribution, and canopy mass. The biogeochemical component CTEM has not changed much from CanAM4 except the diagnostic calculation of wetland extent and methane emissions , none of which affects the physical land surface processes.

3.2 Canadian Small Lake Model

CanAM5 includes a parameterization for subgrid-scale lakes to improve surface fluxes of heat and moisture over land masses. The scheme is based on the Canadian Small Lake Model, CSLM . This scheme computes a nonlinear surface energy balance in a thin skin layer and then solves the heat equation based on thermal conduction and shortwave radiation extinction following Beer's law for both visible and near-infrared bands. A diurnal surface mixed layer is simulated based on the bulk turbulent kinetic energy approach, e.g., , developed by , , and for lakes. A seasonal thermocline arises naturally as a result of the daily excursions of the surface mixed layer. The equation of state follows , except that the effects of pressure and salinity are neglected.

The model allows for the formation of both black, i.e., congelation, and white ice. Black ice forms when the energy balance in a layer is sufficiently negative to cool it below 0 ${}^{\circ }$C. White ice forms when the weight of the overlying snowpack is sufficient to crack the ice and allow lake water to flood a layer of snow, which is then assumed to freeze immediately and completely. Latent heat from the freezing of the pore water is first used to warm the snow crystals in the slush layer to 0 ${}^{\circ }$C, with the remainder going into the overlying snowpack. Both white and congelation ice is assumed to be free of air bubbles and to have the same transmissivity.

Fractional ice cover, following , and fractional snow on ice are permitted, thus allowing for the simultaneous presence of open water, bare ice, and snow-covered ice. Fractional ice cover is especially important for larger lakes subject to sufficient wind stress, which can mechanically break ice to produce pressure ridges and open water leads. The presence of some open water will alter turbulent and radiative flux exchange with the atmosphere, as well as light availability at depth due to differences in roughness, albedo, and light extinction between water and ice. Snow itself is represented as in the Canadian Land Surface Scheme (CLASS; Sect. ), with the snowpack simulated as a layer thermally distinct from the underlying ice.

The properties and interaction of all lakes within a CanAM5 grid cell are modelled by one representative subgrid lake using CSLM. The properties of the representative lake in each CanAM5 grid cell are derived from the Global Lake Database version 2 (GLDv2) , which is provided at $\mathrm{1}/\mathrm{120}$${}^{\circ }$ $\times$ $\mathrm{1}/\mathrm{120}$${}^{\circ }$ resolution. The grid fraction covered by the representative lakes is derived from the aggregate area of lakes in GLDv2 that falls within each CanAM5 grid cell. This defines the unresolved lake tile (Sect. ). Lake dynamics are governed by three external geophysical parameters that must be specified: the visible light transparency, mean depth (or volume), and mean fetch. For all representative lakes, a constant transparency of 0.5 m${}^{-\mathrm{1}}$ is assumed and the mean depth and mean fetch in each CanAM5 grid cell are derived in an aggregate manner from GLDv2.

3.3 Tiling

To more easily facilitate conservative coupling, all previous versions of coupled atmosphere–ocean models developed at CCCma employed coincident grids with an identical binary land mask; e.g., CanESM2 employed CanAM4's land mask for its CMIP5 contribution. In these earlier versions, enhanced ocean resolution was achieved by prescribing multiple ocean grid cells below each CanAM atmospheric grid cell. For CanESM5, independent arbitrarily oriented grids are assumed for both the atmosphere and the ocean . This required the implementation of a fractional land mask in the atmospheric model and tiling of its underlying surface. In general, each CanAM5 grid cell can contain tiles representing land, ocean, sea ice, and unresolved lakes. The tiling approach used is a generalization of that discussed in Sect.  for the tiling of vegetation types over the land portion of atmospheric grid cells. For example, within each model grid cell, independent energy and water fluxes are derived over each underlying surface type given its unique properties, e.g., temperature and albedo. On the atmospheric side, these fluxes undergo a weighted aggregation based on the tile fraction to produce a single flux seen by the atmosphere. In fully coupled mode, if ocean and/or sea ice sits below some portion of an atmospheric grid cell, the flux of each representing each surface type is passed to the coupler, CanCPL, and is remapped and transferred to the underlying grid of the ocean and/or sea ice model.

Currently, surface tiling in CanAM5 has been implemented for the parameterization of radiative transfer, surface processes, and vertical diffusion. Aside from the radiation, the fluxes over each tile are derived from the prognostic variables in the lowest model level, e.g., temperature, specific humidity, and wind. For simplicity, the blending height at which the fluxes from each tile are aggregated is also taken to occur in the lowest atmospheric model level. For radiative transfer calculations, profiles of radiative fluxes are computed for each of the tiles and aggregated into a grid mean, while fluxes are maintained for each tile as described in Sect. .

3.4 Snow on sea ice

For snow on sea ice in CanAM5, the parameterization of snow cover fraction was updated and a parameterization of wet-snow grain growth added to improve consistency with the treatment of snow on land.

In CanAM4, different parameterizations of snow cover were used on land and on sea ice, the snow cover on land being

3$$X_{\mathrm{snow}}=\left\{\begin{array}{ll}{Z}_{\mathrm{snow}}/Z_{\mathrm{snow},\lim }& \text{if }{Z}_{\mathrm{snow}}\le Z_{\mathrm{snow},\lim }\\ \mathrm{1.0}& \text{if }{Z}_{\mathrm{snow}}>Z_{\mathrm{snow},\lim }\end{array}\right.$$ and snow cover over sea ice being 4$$X_{\mathrm{snow}}=\left\{\begin{array}{ll}\sqrt{\mathrm{SWE}/{\mathrm{SWE}}_{\lim }}& \text{if SWE}\le {\mathrm{SWE}}_{\lim }\\ \mathrm{1.0}& \text{if SWE}>{\mathrm{SWE}}_{\lim }\end{array}\right.,$$ where $X_{\mathrm{snow}}$ is the fractional area of the land or sea ice covered with snow, $Z_{\mathrm{snow}}$ is the depth of the snow (in m), and SWE is the snow water equivalent (in kg m${}^{-\mathrm{2}}$), with $Z_{\mathrm{snow},\lim }$ and SWE${}_{\lim }$ being adjustable limits for each. In CanAM5, Eq. () is used to determine snow cover over land and sea ice. Note that in Eq. (), $Z_{\mathrm{snow}}$ is initially computed using $Z_{\mathrm{snow}}=\text{SWE}/{\mathit{\rho }}_{\mathrm{snow}}$, where ${\mathit{\rho }}_{\mathrm{snow}}$ is the snow density (in kg m${}^{-\mathrm{3}}$). If $Z_{\mathrm{snow}}\le Z_{\mathrm{snow},\lim }$, then $Z_{\mathrm{snow}}=Z_{\mathrm{snow},\lim }$ and the SWE is adjusted accordingly .

The computation of sea ice albedo includes a contribution from snowpack on sea ice when it is present. To calculate the albedo of snow, it is necessary to simulate the relevant physical properties of the snow, including the snow grain size. The approach used to parameterize these properties in CanAM5 is described in . Described here is the addition of a parameterization to CanAM5 so that the wet growth of snow grains is included for snow on sea ice, where previously only the dry growth of snow grains was considered.

To calculate the wet growth of snow grains, the same expression is used as over land (Eq. 3 of ), which requires the liquid water fraction in the snowpack. This was not available in CanAM5, so we added a parameterization of snowpack liquid water fraction using : 5$$F_{\mathrm{liq}}=\left\{\begin{array}{ll}{F}_{\mathrm{liq},\min }& \text{if }{\mathit{\rho }}_{\mathrm{snow}}\ge {\mathit{\rho }}_{\mathrm{snow},\mathrm{thres}}\\ \\ F_{\mathrm{liq},\min }+\mathrm{\Delta }{F}_{\mathrm{liq}}& \text{if }{\mathit{\rho }}_{\mathrm{snow}}<{\mathit{\rho }}_{\mathrm{snow},\mathrm{thres}}\end{array}\right.,$$ where $F_{\mathrm{liq}}$ is the fraction of liquid in the snow pack, ${\mathit{\rho }}_{\mathrm{snow}}$ is the density of snow (in kg m${}^{-\mathrm{3}}$), and ${\mathit{\rho }}_{\mathrm{snow},\mathrm{thres}}$ is the snow density threshold at which $F_{\mathrm{liq},\min }$ occurs. The term $\mathrm{\Delta }{F}_{\mathrm{liq}}$ is $(F_{\mathrm{liq},\max }-F_{\mathrm{liq},\min })\times \frac{{\mathit{\rho }}_{\mathrm{snow},\mathrm{thres}}-{\mathit{\rho }}_{\mathrm{snow}}}{{\mathit{\rho }}_{\mathrm{snow},\mathrm{thres}}}$, where $F_{\mathrm{liq},\max }$ and $F_{\mathrm{liq},\min }$ are the maximum and minimum allowed values of $F_{\mathrm{liq}}$.

4 Setting adjustable parameters

After finalizing the new and updated physical parameterizations, they were no longer changed, except for a subset of adjustable parameters. These parameters were manually adjusted within a range of physically plausible values to obtain an acceptable preindustrial climate in the coupled atmosphere–ocean configuration of CanESM5 . This is the last exercise performed to finalize a model version and is often referred to as “tuning”. The subset of parameters and values in CanAM5 is provided in Table . They include parameters adjusted specifically for CanAM5 and parameters adjusted when tuning intermediate versions of CanAM between versions 4 and 5, e.g., CanAM4.1. In this section, we discuss the process used to arrive at the values, which is different from that used in CanAM4 and CanESM2.

Table 1

Adjustable parameters in CanAM and their settings in CanAM5. The values in bold were specifically tuned in CanAM5, while the others were used to tune intermediate versions of CanAM. The rightmost column indicates references that discuss the adjustable parameter or include further references about the parameter.

Scheme	Parameter	Physical description	CanAM5	Unit	Comment/reference
Cloud microphysics	facacc	Factor scaling mass accretion rate of cloud water to precipitation due to the collection of cloud droplets by raindrops	15	–
	facaut	Factor scaling efficiency coefficient in mass autoconversion rate of cloud water to precipitation due to the collision–coalescence processes of cloud droplets	1	–
	uicefac	Prefactor in power law describing ice crystal fall speed due to the influence of gravity	6000	s${}^{-\mathrm{1}}$
Moist convection	alf	Proportionality parameter relating vertically integrated convective kinetic energy with the cloud base mass flux	$\mathrm{5.0}\times {\mathrm{10}}^{\mathrm{8}}$	m${}^{\mathrm{4}}$ kg${}^{-\mathrm{1}}$
	ccu	Weight large-scale and pressure gradient force contributions to moist convection horizontal velocity (updrafts)	0.0	–
	ccd	Weight large-scale and pressure gradient force contributions to moist convection horizontal velocity (downdrafts)	0.0	–
Gravity wave	fcrit	Critical inverse Froude number	0.22	–
	gphil	Mountain sharpness number	1.0	–
Vertical diffusion	rkhmn	Minimum background vertical diffusivity for temperature	0.1	m${}^{\mathrm{2}}$ s${}^{-\mathrm{1}}$
	rkqmn	Minimum background vertical diffusivity for moisture	0.1	m${}^{\mathrm{2}}$ s${}^{-\mathrm{1}}$
Surface processes	drn	Scaling factor for soil drainage at the bottom of the soil levels	0.1	–
	cuscale	Scaling factor of the wind stress threshold for dust emissions	1.6	–
	reff0_sea	Background specific surface area of snow grains (on sea ice)	30	m${}^{\mathrm{2}}$ kg	Personal communication (Joshua King)
	reff0_land	Background specific surface area of snow grains (on land)	60	m${}^{\mathrm{2}}$ kg	Personal communication (Joshua King)
	albp	Depth of melt ponds on sea ice	20	cm

The tuning of CanESM2 was carried out mainly by adjusting the parameters of each of its components separately, including CanAM4, with a goal of minimal additional adjustments when fully coupled. For example, the parameters for CanAM4 were mostly tuned using transient prescribed sea surface temperature (SST) and sea ice simulations of the near present, consistent with simulations used regularly for CanAM development. Applying the same approach to the tuning of CanESM5 resulted in a coupled preindustrial (1850) control simulation with a climate that was too cold with excessive sea ice relative to observations. Therefore, CanAM5 was tuned in the context of fully coupled CanESM5 simulations, with a particular focus on obtaining preindustrial control conditions with global mean temperatures and sea ice within acceptable ranges. The combination of parameters that achieved this target was then evaluated to verify that other aspects of the climate remained acceptable.

Analyses to investigate the effect of parameter sets on the climate included CanAM5 simulations using prescribed SSTs and sea ice (Atmospheric Model Intercomparison Project, AMIP; ). For the most part, the mean climate simulated in AMIP mode was close to coupled CanESM5 simulations with the exception of the net radiative flux at the top of atmosphere (TOA). Adjustments required to ensure an acceptable preindustrial climate resulted in a net downward flux at TOA ($\sim \mathrm{3.1}$ W m${}^{-\mathrm{2}}$) in historical AMIP runs, which is larger than observations. Simulations in which the AMIP net downward fluxes at TOA were close to those observed resulted in a preindustrial global mean temperature up to 2K colder than the target value. Although the AMIP net downward fluxes at TOA are larger than those observed, their value in CanESM5 historical coupled simulations during the present day are very close to observations ($\sim \mathrm{1}$ W m${}^{-\mathrm{2}}$; Fig. S6). Details of the TOA radiative fluxes are discussed in Sect. . For the purposes of tuning, this particular bias in AMIP simulations was retained to get a reasonable preindustrial control climate.

Table  lists parameters that changed between CanAM4 and CanAM5, with those in bold specifically adjusted for CanAM5 and others having been changed when tuning intermediate CanAM versions between CanAM4 and CanAM5. This table does not include parameters that were adjusted in the ocean or in the sea ice model that only affected coupled simulations. The rightmost column of Table  provides sources and, where possible, references that explain the setting in CanAM. This final set of CanAM5 parameters allows us to simulate a climate that is on balance reasonable relative to observations in both coupled and AMIP mode.

The parameters related to cloud microphysics have notable effects on radiative energy budgets and coupled climate, including emergent CanESM5 properties such as climate sensitivity. Of particular importance are the two parameters scaling the efficiency of cloud droplet autoconversion and accretion to precipitation. In CanAM5, the accretion rate factor was the main parameter adjusted instead of autoconversion, which is opposite to the approach used when tuning CanAM4. Analysis of satellite observations by indicates global climate models may severely underestimate mean accretion rates when subgrid cloud–precipitation covariability is omitted. Furthermore, , , and showed that diagnostic parameterizations of rain processes, such as those employed in CanAM5, produce considerably lower accretion rates than prognostic and more comprehensive parameterizations. Consequently, the usual assumptions of an instantaneous and horizontally uniform precipitation flux in the cloudy portions of the grid cells in CanAM5 likely cause unrealistically low accretion rates. In an attempt to compensate for this, the original parameterization of accretion of is made more efficient through the considerable increase (by a factor of 15) in the tunable parameter. Autoconversion rates, on the other hand, are not scaled.

5 Control and CMIP6 simulations

Unlike CanAM4, CanAM5 has an interactive land carbon cycle (Sect. ) which necessitates starting CanAM5 transient simulations from a state with a land carbon cycle that is reasonably close to equilibrium. To achieve this, sufficiently long simulations with a stable climate are required. This is done using an approach similar to the spinup of the CanESM5 preindustrial simulation . A long control simulation of CanAM5 is performed using a repeating annual cycle of forcing and prescribed sea surface temperature (SST) and sea ice. For this CanAM5 control, we use forcing for the year 1870 (the first year in the historical SST and sea ice dataset), while the annual cycle of SST and sea ice is the mean over 1870–1879.

Figure 1

Net atmosphere–land CO${}_{\mathrm{2}}$ flux (Pg C yr${}^{-\mathrm{1}}$), soil carbon mass (Pg C), and total soil moisture (kg m${}^{-\mathrm{1}}$) for the years 300–500 of the CanAM5 1870 control simulation. The red and black lines show the results for CanAM5 using two different physics configurations of CanESM5, p1 and p2, respectively, the details of which are described in .

[Figure omitted. See PDF]

With this configuration, the CanAM5 control simulation was initialized from a coupled CanESM5 1850 control simulation and run for $\sim \mathrm{300}$ years until the physical and biogeochemical land states, including carbon in vegetation and soil, approached a new quasi-equilibrium. The simulation was then extended by an additional 200 years. Figure  shows the net atmosphere–land CO${}_{\mathrm{2}}$ flux, the amount of C in the soil carbon pool, and the total soil moisture during the 200 years. Most atmospheric variables reached a quasi-equilibrium within a few years and are therefore not shown. For CanESM5, transient coupled simulations were started from the control coupled simulation every 50 years . A similar approach was used for CanAM5 with transient simulations starting from the CanAM5 control simulation every 10 years beginning at the year 400. These are used to generate a 10-member ensemble using transient forcings, SSTs, and sea ice for the period 1870 to 2014.

Several CMIP6 experiments were performed using CanAM5 and prescribed SSTs and sea ice. The 10-member ensemble of 1870–2014 transient simulations was contributed to the Global Monsoon Model Intercomparison Project , while the period 1950 to 2014 from each simulation makes up the CanESM5 contribution to the Atmospheric Model Intercomparison Project (AMIP) experiments . AMIP simulations are the basis for several Cloud Feedback Model Intercomparison Project (CFMIP) experiments used to characterize and understand cloud feedbacks . To characterize radiative forcings in CanESM5, simulations were performed using the Radiative Forcing Model Intercomparison Project (RFMIP) protocols for time slice and transient historical forcings , which are summarized in .

6 Evaluation of CanAM5 climatology

The properties of coupled atmosphere–ocean experiments using CanESM5 are shown in . Documented in this section is the climatology of CanAM5 from the CMIP6 AMIP simulation, evaluated against observations and highlighting differences compared to CanAM4 from a CMIP5 AMIP simulation. Details of observations used for evaluation are summarized in Table , and model variables are summarized in Table . For all figures, the first ensemble member is used for each AMIP simulation, r1i1p1 for CanAM4 and r1i1p2f1 for CanAM5. Included in several figures is the global, or near-global, mean bias, root mean square error, and Pearson correlation coefficient between the time-averaged CanAM and observations. For the most part, the results using prescribed SSTs and sea ice are similar to coupled CanESM5 and CanESM2 simulations. For ease of comparison with CanESM5 coupled simulations , the figures in this section have been reproduced in the Supplement, Sect. S2, using the first member (r1i1p2f1) of the CanESM5 historical simulation ensemble.

Figure 2

Mean histograms of the cloud fraction equatorward of 60${}^{\circ }$ as a function of the cloud top pressure and cloud visible optical thickness from ISCCP-H and the biases in CanAM5 (b, c). To the side of each histogram is the mean cloud fraction, or cloud fraction bias, as a function of cloud top pressure, while below each histogram the cloud fraction is shown, or cloud fraction bias, as a function of cloud optical thickness. Means are averages for 1987–2008.

[Figure omitted. See PDF]

6.1 Clouds and precipitation

The near-global (equatorward of 60${}^{\circ }$) cloud fraction as a function of cloud optical thickness and cloud top pressure is shown in Fig. . For the purposes of comparing more consistent model output from CanAM5 and CanAM4 with observations, output from the International Satellite Cloud Climatology Project (ISCCP) simulator is compared with ISCCP observations, both ISCCP-D and ISCCP-H . The two versions of ISCCP observations are used to illustrate the uncertainty in the cloud properties, an uncertainty that only increases once other cloud observations are considered . For example, showed that there are large differences between ISCCP and MODIS (larger than between the two versions of ISCCP shown here). Therefore, it is important that such differences be considered in the evaluation of models, especially for optically thin clouds.

With these caveats in mind, the histograms of biases indicate that CanAM5 generally simulates too much cloud with moderate optical thickness of high- and low-altitude clouds and simulates too little cloud at mid-level altitudes. Summing the histograms over cloud top pressure to look at clouds as a function of visible optical thickness, we see that there are more optically thin ($\mathit{\tau }<\mathrm{23}$) clouds in CanAM5 than in CanAM4. The structure of these biases relative to ISCCP is consistent with previous studies, for example .

Figure 3

Zonal mean cloud fraction for the total cloud amount with $\mathit{\tau }>\mathrm{0.3}$ (a, b), cloud amount for low (cloud top pressure $>\mathrm{680}$ hPa) and non-low (cloud top pressure $<\mathrm{680}$ hPa) in (c) and (d), and the cloud amount for thin (cloud visible $\mathit{\tau }$ between 0.3 and 23) and thick (cloud visible $\mathit{\tau }>\mathrm{23}$) in (e) and (f). Observations for ISCCP-H and ISCCP-D are shown in (a), (c), and (e) and biases for CanAM5 and CanAM4 relative to ISCCP-H in (b), (d), and (f). Bracketed numbers in (b), (d), and (f) are, in order, mean bias, root mean square error, and Pearson correlation coefficient, computed over the period 1987–2008 and from 60${}^{\circ }$ S to 60${}^{\circ }$ N, with black font for CanAM5 and red font for CanAM4.

[Figure omitted. See PDF]

The zonal mean structure for cloud amount is presented (Fig. ), which illustrates that the near-global mean biases are the result of regional biases which are a source of biases and improvements in the cloud radiative effect (CRE) (Fig. ). As seen in the near-global means, the differences between ISCCP-D and ISCCP-H are smaller than biases between CanAM and ISCCP-H and the change in biases between CanAM4 and CanAM5. Although there remain biases in the total cloud amount, there is a systematic reduction in biases in CanAM5 by $\sim$ 3.4 % in the near-global mean compared to biases in CanAM4. Parsing the biases in CanAM5 by the altitude of cloud top pressure, the increase in the CanAM5 total cloud amount mostly is caused by increased non-low (cloud top pressure $<\mathrm{680}$ hPa) cloud amount. From CanAM4 to CanAM5, there is an increase in the amount of “thin” (cloud visible $\mathit{\tau }$ between 0.3 and 23) and reduction in “thick” (cloud visible $\mathit{\tau }>\mathrm{23}$) cloud at all latitudes, consistent with the near-global mean (Fig. ). This increase in cloud amounts is consistent with the change in CREs (Fig. ). The shift to more optically thin cloud would reduce the reflectively, doing so in a nonlinear manner, while the increase in the total cloud fraction will increase the CRE in a linear manner.

Figure 4

Zonal cloud fraction and cloud phase from CanAM5 compared with GOCCP and MODIS observations averaged over 2007–2009. Black contours in (b) and (d) are the zonal mean cloud fraction from CanAM5. Bracketed numbers in (e) and (f) are, in order, mean bias, root mean square error, and Pearson correlation coefficient, computed using data between 75${}^{\circ }$ S and 75${}^{\circ }$ N.

[Figure omitted. See PDF]

The CMIP6 protocol requested the additional diagnostic output consistent with retrievals based on lidar observations from Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) and Moderate Resolution Imaging Spectroradiometer (MODIS) imager measurements . These are used to evaluate the vertical structure of the clouds in CanAM5 and the ability of CanAM5 to simulate the cloud phase (Fig. ). Biases in the cross section of cloud amount for CanAM5 relative to the GCM-Oriented CALIPSO Cloud Product (GOCCP) (upper row Fig. ) are consistent with biases between CanAM4 and GOCCP . There is generally too much cloud simulated at higher altitudes and too little cloud simulated at lower altitudes.

The middle and lower rows of Fig.  use diagnostics of GOCCP cloud-phase profiles and MODIS cloud top phase. In the tropics, CanAM5 underestimates the fraction of cloud that is ice in the middle troposphere; however, it occurs in a range of altitudes where CanAM5 is already simulating too few clouds. The more notable bias is that CanAM5 simulates too much ice cloud poleward of $\sim \mathrm{50}$${}^{\circ }$, seen in the GOCCP and MODIS diagnostics. These biases in the high-latitude cloud phase can have important consequences on the radiation budget and cloud feedbacks in these regions .

Precipitation biases are an important feature of any climate model. Although the structure of the biases in CanAM5 is similar to that in CanAM4, there are improvements in some key regions. The most noticeable improvement is the increased precipitation rate over the Amazon in CanAM5 for most seasons, although dry biases remain (Fig. ). This change in precipitation is consistent with a reduction in temperatures that are too warm over the Amazon in CanAM5 (Fig. ), which may be due to more moist conditions suppressing the near-surface temperature.

Figure 5

Seasonal mean precipitation rate from GPCP (left column), the bias of CanAM5 relative to GPCP (middle column), and the bias of CanAM4 relative to GPCP (right column). Bracketed numbers to the upper right of difference plots are, in order, mean bias, root mean square error, and Pearson correlation coefficient. All plots use data from the years 1980–2009.

[Figure omitted. See PDF]

Figure 6

Global and time mean radiative fluxes (shortwave, SW, and longwave, LW) at the top of atmosphere (TOA) and surface, as well as the net flux divergence for the atmosphere, from AMIP simulations by CanAM5 and CanAM4 compared with CERES EBAF. Statistics are computed over the period 2003–2009. For each pair of bracketed numbers in (a)–(d), the left value is CERES and the right value is CanAM. In (e)–(h), the bracketed numbers are root mean square error and Pearson correlation coefficient of CanAM relative to CERES.

[Figure omitted. See PDF]

6.2 Radiation

Radiative fluxes through the top of atmosphere (TOA) and bottom of atmosphere are evaluated using CERES observations . Figure  summarizes the global mean climatology for the solar and thermal flux components of the radiative energy budget for CanAM5 and CanAM4. Although a relatively short common period is used from the models and CERES observations (2003–2009), the results are very similar to those using longer periods from CERES (2003–2020) and CanAM (1979–2009).

Figure 7

Annual global and zonal mean cloud radiative effects at the top of atmosphere (TOA), atmosphere (ATM), and surface (SFC) from CERES EBAF observations (a, c, e, g) and from CanAM5 and CanAM4 AMIP simulations (b, d, f, h). The means are averages over 2003–2009.

[Figure omitted. See PDF]

We focus first on fluxes at TOA, which can be most directly compared with space-based observations from CERES. The global mean thermal fluxes are effectively identical in CanAM4 and CanAM5. The change in the downward solar flux is due to the use of updated solar forcing that has a reduced total solar irradiance, which is more consistent with observations and CERES. The upward solar flux is reduced by 3.4 W m${}^{-\mathrm{2}}$. There is a small reduction in the clear-sky upward solar flux at TOA, $\sim \mathrm{0.2}$ W m${}^{-\mathrm{2}}$, so the remainder of this reduction is due to clouds (Fig. ).

Evaluated separately, the solar and thermal radiative fluxes are within the range of values from the CMIP6 simulations . When all fluxes are combined to compute the net flux imbalance at TOA, CanAM5 has a value that is larger than CERES and CanAM4 by 2.2 W m${}^{-\mathrm{2}}$. In CanAM4 there is a compensation between the upward thermal (too small) and upward solar (too large) fluxes, resulting in a net imbalance that is in line with observations, while in CanAM5 both the upward thermal and solar fluxes are smaller than observations. We note that at least one other CMIP6 model documented a similar difference between AMIP and coupled simulations . To put the CanAM5 results into context with other models participating in CMIP6, we compared the net flux at TOA averaged over 2003–2009 for 34 models which had at least one AMIP and one historical coupled simulation. Of the 34 models, 18 AMIP simulations have absolute global mean differences relative to CERES that are greater than 1 W m${}^{-\mathrm{2}}$ but only three historical coupled simulations have a difference greater than 1 W m${}^{-\mathrm{2}}$ (not shown). This indicates that the behaviour seen in CanAM5 and CanESM5 is not unique among CMIP6 models.

The TOA flux imbalance is larger than that for coupled CanESM5 simulations, $\sim \mathrm{1.1}$ W m${}^{-\mathrm{2}}$, averaged over 2003–2009 (Fig. S6). This is mainly due to solar fluxes which are larger (99.3 W m${}^{-\mathrm{2}}$) than when using observed SSTs (97.7 W m${}^{-\mathrm{2}}$), since upward thermal fluxes at TOA are similar (239.8 W m${}^{-\mathrm{2}}$ versus 239.5 W m${}^{-\mathrm{2}}$). An in-depth analysis of why this occurs is beyond the scope of this paper. Preliminary analysis using CanAM5 with combinations of sea ice and SST specification, from observations and coupled CanESM5 simulations, suggests that differences in SSTs are the main factor which may be due to local and nonlocal responses affecting the TOA radiative fluxes.

While the TOA radiative fluxes were regularly evaluated during the development of CanAM, radiative fluxes at the surface and within the atmosphere were not. For both CanAM5 and CanAM4 the biases at the surface relative to CERES are consistent: downward and upward longwave fluxes that are too small and downward and upward shortwave fluxes that are too large. Altogether this results in too little absorption of radiation at the surface by 3 to 5 W m${}^{-\mathrm{2}}$ and fairly consistent overestimation of net absorption in the atmosphere by 5 W m${}^{-\mathrm{2}}$ due mainly to too much absorption of longwave radiation.

Clouds strongly modulate radiative fluxes, so we next examine the simulated cloud radiative effects (CREs), defined as $\text{CRE}=F_{\text{clear sky}}-F_{\text{all sky}}$, where $F_{\text{clear sky}}$ is the radiative flux in the absence of clouds and $F_{\text{all sky}}$ is the radiative flux with clouds present. The annual mean cloud radiative effects are generally positive for longwave and negative for shortwave at TOA and the surface, while the longwave strongly controls the zonal atmospheric CRE (Fig. ). The global mean CREs simulated by CanAM5 are less biased relative to CERES than CanAM4, especially for shortwave CRE, while the longwave CRE at TOA is slightly more biased than CanAM4 (Fig. ). Zonal mean CREs show that the improvements seen in the CanAM5 global means are due to reduced biases at most latitudes (Fig. ). In addition to improved global mean biases, the root mean square error is decreased (Fig. S1). These improved CREs suggest improved simulation of cloud properties in CanAM5 (Sect. ).

Figure 8

Seasonal mean latitude–pressure plots of zonal wind from ERA5 (left column), the bias of CanAM5 relative to ERA5 (middle column), and the bias of CanAM4 relative to ERA5 (right column). For all plots, contours are the mean. For the ERA5 plot, shading is the mean, and in other plots the shading is the bias relative to ERA5. All plots use data from the years 1980–2009.

[Figure omitted. See PDF]

6.3 Circulation

In this subsection we document the climatological properties of the winds, temperature, and surface pressure in CanAM5 for an AMIP experiment. Seasonal climatologies of latitude–height zonal-mean zonal wind fields and anomalies are presented in Fig. . While overall biases are similar relative to CanAM4, CanAM5 displays anomalously positive rather than negative wind biases in the mid- to high-latitude Northern Hemisphere DJF lower stratosphere. This is consistent with weaker planetary wave forcing of the Northern Hemisphere stratosphere in CanAM5. This wintertime positive zonal wind bias is associated with a weakening of the orographic gravity-wave drag due to a change in parameter values between the two model versions (Sect. ). Similarly, this weakening of the gravity wave drag contributes to a larger positive anomaly of zonal-mean zonal winds in CanAM5 in the Southern Hemisphere wintertime stratosphere. The near-surface zonal wind climatology is consistent with the coupled CanESM5 simulations , with biases relative to ERA5 generally smaller in CanAM5, with the most significant reductions in midlatitudes in both hemispheres (Fig. ).

Figure 9

Seasonal mean latitude–pressure plots of temperature from ERA5 (left column), the bias of CanAM5 relative to ERA5 (middle column), and the bias of CanAM4 relative to ERA5 (right column). For all plots, contours are the mean. For the ERA5 plot, shading is the mean, and in other plots the shading is the bias relative to ERA5. All plots use data from the years 1980–2009.

[Figure omitted. See PDF]

In Fig. , seasonal climatologies of latitude–height zonal-mean temperature from ERA5, CanAM5, and CanAM4 are presented. In general, CanAM5 and CanAM4 have similar patterns of temperature bias in all seasons, including a warm tropical tropopause and cool extratropical tropopause. However, there are regional and seasonal differences; for example, temperatures between December and May poleward of 60${}^{\circ }$ N in the stratosphere are not as systematically biased warm in CanAM5 relative to CanAM4. The pattern and magnitude of the temperature biases are similar to those in coupled configurations .

Figure 10

Seasonal mean sea-level pressure from ERA5 (left column), the bias of CanAM5 relative to ERA5 (middle column), and the bias of CanAM4 relative to ERA5 (right column). Bracketed numbers to the upper right of difference plots are, in order, mean bias, root mean square error, and Pearson correlation coefficient. All plots use data from the years 1980–2009.

[Figure omitted. See PDF]

Figure 11

Seasonal mean near-surface temperature from ERA5 (left column), the bias of CanAM5 relative to ERA5 (middle column), and the bias of CanAM4 relative to ERA5 (right column). Bracketed numbers to the upper right of difference plots are, in order, mean bias, root mean square error, and Pearson correlation coefficient. All plots use data from the years 1980–2009.

[Figure omitted. See PDF]

The seasonal mean sea-level pressure is presented in Fig. . Relative to CanAM4, CanAM5 displays larger DJF biases in the Aleutian Low and North Pacific High but lower bias in the whole of the Atlantic Ocean in all seasons.

Seasonal mean biases in near-surface temperature for CanAM5 and CanAM4 are presented in Fig. . Persistent cold biases are found over the Tibetan Plateau and North Africa in both models. The Tibetan Plateau bias is negatively correlated with snow cover bias (too much snow cover and temperatures that are too cold), a feature found in other CMIP6 models . The source of a snow cover that is too large is complex and is present to differing degrees in CanAM4 and CanAM5. That said, it does seem to be a robust feature of CanAM, given that the land model was significantly changed between CanAM4 and CanAM5, including the parameterization of snow albedo (Sect. ). In North Africa, cold biases are thought to be due to the change from a globally constant albedo for bare soil to a more realistic distribution based on local soil conditions .

Warm biases are apparent over the Brazil basin, although somewhat reduced in CanAM5, consistent with biases related to too little precipitation (Fig. ). Over central North America in JJA, the warm bias persists in CanAM5 and is more extensive than in CanAM4. This is a common warm bias among CMIP5 models during JJA for which the cause is thought to be a complex interplay between land–atmosphere coupling, radiation, and clouds that rapidly develop in climate models .

7 Conclusions

CanAM5 is the latest atmospheric model from the Canadian Centre for Climate Modelling and Analysis. In this study, we have presented the main model differences between CanAM5 and its predecessor CanAM4. In particular, these differences are primarily related to radiation, clouds, and aerosols; a major update of the land surface model; and the addition of a parameterization of freshwater lakes. Generally, mean climatologies from CanAM5 for near-present conditions, and using observed SSTs and sea ice, are similar to those from CanAM4, with some notable improvements, including reduced precipitation and temperature biases over the Amazonian basin, reduced cloud fraction biases, and a reduction in solar and thermal CREs. Some biases persist from CanAM4 to CanAM5, e.g., cold biases over the Tibetan Plateau, and new biases are present in CanAM5 when using prescribed SSTs and sea ice, e.g., a bias in net downward flux at TOA. As noted, the bias in the net downward flux at TOA is the result of tuning to have a coupled 1850 control simulation with CanESM5 that is close to target global mean temperature and sea ice area .

Why it was necessary to tune the net downward flux at TOA higher than observations when using observed SSTs and sea ice remains a question for further research. Additional simulations with CanAM5, using combinations of observed SSTs and sea ice with SSTs and sea ice from coupled CanESM5 simulations, suggest that this is due to the SSTs in CanESM5. This was not the case for CanESM2 and CanAM4, which could be largely tuned for coupling using observed SSTs and sea ice. Further analysis, including the use of Green's functions to link regional differences in SSTs to global mean fluxes at TOA, should help inform future tuning of CanESM and CanAM. Another question considered is why CanESM5 has a significantly larger climate sensitivity than CanESM2 and nearly all CMIP6 models . At present, this is thought to be mostly due to changes in cloud feedbacks . This suggests that improved mean climatologies of clouds and radiation in CanAM5 and CanESM5 do not necessarily result in improved cloud feedbacks and climate sensitivity. A better understanding of both these questions will provide guidance for the ongoing development of CanAM.

Appendix A Table A1

Observational data used for model evaluation.

Data	Description	Version	Reference
TOA radiative fluxes	CERES EBAF-TOA	4.1	,
Surface radiative fluxes	CERES EBAF	4.1	,
Lidar-based cloud amount	GOCCP (3D_CloudFraction)	3.1.2
Lidar-based cloud phase	GOCCP (3D_CloudFraction_phase)	3.1.2	,
Cloud amount histogram	ISCCP H (HGG)	v01r00	, , ,
Cloud amount histogram	ISCCP D	2	,
Cloud top phase	MODIS	6
			http://dx.doi.org/10.5067/MODIS/MCD06COSP_D3_MODIS.061
Atmospheric and surface data	ERA	5

Table A2

CMIP6 and CMIP5 data used for model evaluation.

Figure number	CMIP6/CMIP5 variable	Description
Figure	clisccp	Histogram of cloud amount by cloud top pressure
		and cloud visible optical thickness
Figure	clisccp	Histogram of cloud amount by cloud top pressure
		and cloud visible optical thickness, consistent with ISCCP
Figure	clcalipso, clcalipsoliq, clcalipsoice,	Cloud profile consistent with CALIPSO and
	clwmodis, climodis, cltmodis	cloud fraction from MODIS
Figure	pr	Precipitation rate
Figure	rsdt, rsut, rlut, rsds, rsus, rlds, rlus	Radiative fluxes at top of atmosphere and surface
Figure	rsdt, rsut, rsutcs, rlut, rlutcs, rsds	Radiative fluxes at top of atmosphere and surface
	rsus, rsdscs, rsuscs, rlds, rlus, rldscs
Figure	ua	Zonal wind
Figure	ta	Temperature
Figure	psl	Sea-level pressure
Figure	tas	Near-surface temperature

Figure A1

Annual mean near-surface zonal wind, nominally 10 m above the surface from ERA5 (upper left) and CanAM5 and CanAM4 biases. For all plots, contours are the mean. While for the ERA5 plot shading is also the mean, in the other plots the shading is the bias relative to ERA5. Bracketed numbers to the upper right of difference plots are, in order, root mean square error and Pearson correlation coefficient. All plots use data from the years 1980–2009.

[Figure omitted. See PDF]

Table A3

Acronyms, initialisms, and abbreviations used in the paper.

AMIP	Atmospheric Model Intercomparison Project
CCCma	Canadian Centre for Climate Modelling and Analysis
CanAM	Canadian Atmospheric Model
CanESM	Canadian Earth System Model
CMAM	Canadian Middle Atmosphere Model
CERES	Clouds and the Earth's Radiant Energy System
CMIP	Coupled Model Intercomparison Project
CFMIP	Cloud Feedback Model Intercomparison Project
CLASS	Canadian Land Surface Scheme
CSLM	Canadian Small Lake Model
CRE	Cloud radiative effect
CALIPSO	Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observation
DECK	Diagnostic, Evaluation and Characterization of Klima
DJF	December–January–February
ERA	ECMWF (European Centre for Medium-Range Weather Forecasts) Reanalysis
GMMIP	Global Monsoons Model Intercomparison Project
GLD	Global Lake Database
GOCCP	GCM-Oriented Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) Cloud Product
HITRAN	High-resolution transmission
ISCCP	International Satellite Cloud Climatology Project
JJA	June–July–August
MAM	March–April–May
MODIS	Moderate Resolution Imaging Spectroradiometer
RFMIP	Radiative Forcing Model Intercomparison Project
SON	September–October–November
TOA	Top of atmosphere

Code and data availability

The full CanESM5 source code is publicly available at https://gitlab.com/cccma/canesm (last access: 14 September 2023) and includes CanAM5 as a submodule. The version of the code which can be used to produce all simulations submitted to CMIP6 and described in this paper is tagged as v5.0.3 and has the following associated DOI: 10.5281/zenodo.3251114 . The scripts used to produce all the figures are available at 10.5281/zenodo.7579680 . All CanESM5–CanAM5 and CanESM2–CanAM4 simulations conducted for CMIP6 and CMIP5, respectively, including those described in this paper, are publicly available via the Earth System Grid Federation (ESGF). All observational data used are publicly available and are listed in Table .

The supplement related to this article is available online at: https://doi.org/10.5194/gmd-16-5427-2023-supplement.

Author contributions

JNSC drafted the manuscript and created the figures, contributed to development of CanAM5, and performed simulations with CanAM4 and CanAM5 used in the paper. KvS led the development of CanAM5 and wrote the “Aerosols and chemistry” section. JL developed the radiative transfer code and wrote the Radiation section. JS wrote the Tiling, Circulation, and “Clouds and precipitation” sections and contributed to CanAM5 development. DP contributed to CanAM5 development. VA developed CLASS-CTEM and wrote the section that describes it. NM and ML contributed to CanAM5 development. MM developed the Canadian Small Lake Model and wrote the section describing it. DV developed CLASS. All authors contributed to writing the manuscript.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

We thank Carsten Abraham, Barbara Winter, and two anonymous reviewers for many helpful comments that improved the manuscript. We also thank Barbara Winter for preparing forcing datasets used by CanAM5 for CMIP6 experiments shown in this paper.

Review statement

This paper was edited by Xiaohong Liu and reviewed by two anonymous referees.