René Laprise

A Semi-Implicit Semi-Lagrangian Regional Climate Model: The Canadian RCM

A new regional climate model (RCM) is presented in this paper and its performance is investigated through a pair of 60-day simulations. This new model is based on the dynamical formulation of the Cooperative Centre for Research in Mesometeorology (CCRM) mesoscale nonhydrostatic community model and on the complete subgrid-scale physical parameterization package of the second-generation Canadian Centre for Climate modeling and analysis General Circulation Model (CCCma GCMII). The main feature of the Canadian RCM (CRCM) comes from the very efficient semi-implicit and semi-Lagrangian (SISL) numerical scheme used for the integration of the fully elastic nonhydrostatic Euler equations. The efficiency of the SISL scheme allows the use of longer time steps (at least by a factor of 5) for the integration of this model (e.g., the 45-km resolution version of the model uses a 15-min time step). A complete description of the numerical formulation of the model is presented with a review of the principal characteristics of the physical package. A pair of two-month-long winter simulations is also analyzed to investigate the behavior of the model and to evaluate the potential of the SISL integration scheme in the context of regional climate simulation. The two integrations, produced with a 45-km resolution version of the model, developed realistic small-scale details from the low-resolution GCMII fields used to initialize and drive the RCM.

The Euler Equations of Motion with Hydrostatic Pressure as an Independent Variable

Abstract A novel form of the Euler equations is developed through the use of a different vertical coordinate system. It is shown that the use of hydrostatic pressure as an independent variable has the advantage that the Euler equations then take a form that parallels very closely the form of the hydrostatic equations cast in isobaric coordinates. This similarity holds even when topography is incorporated through a further transformation into terrain-following coordinates. This leads us to suggest that hydrostatic-pressure coordinates could be used advantageously in nonhydrostatic atmospheric models based on the fully compressible equations.

A Semi-implicit Send-Lagrangian Fully Compressible Regional Forecast Model

Abstract The semi-implicit algorithm, originally developed by Robert for an economical integration of the primitive equations in large-scale models of the atmospheric, is here generalized in order to integrate the fully compressible, nonhydrostatic equations. We show that there is little computational overhead associated with the integration of the full, and hence presumably more correct, set of equations that do not invoke the hydrostatic assumption to exclude the high frequency, vertically propagating acoustic modes.

The Canadian Climate Centre spectral atmospheric general circulation model

Abstract A general description of the Canadian Climate Centre atmospheric general circulation model is presented. The model includes, either in explicit or parametric form, all of the physical processes deemed important for long‐term climate simulations. Detailed descriptions of the methods used to represent these processes are presented. Selected results from test runs with the model are presented to illustrate its sensitivity to some aspects of the subgrid‐scale vertical flux parameterizations and the gravity wave drag formulation.

Regional climate modelling

A new technique for climate simulations and projections has emerged over the past 15 years that allows to reach unprecedented spatial resolution at an affordable computing cost: regional climate modelling. The technique has matured and several groups around the world are participating in model intercomparison projects and collaborating in coordinated endeavour to produce ensemble climate-change projections at regional scales. The paper describes the formulation of a specific regional climate model (RCM), the Canadian RCM, commenting on the formulation choices and experience gained during the development and validation process of this model. Current issues in regional climate modelling are also reviewed and discussed.

Project to Intercompare Regional Climate Simulations (PIRCS): Description and initial results

The first simulation experiment and output archives of the Project to Intercompare Regional Climate Simulations (PIRCS) is described. Initial results from simulations of the summer 1988 drought over the central United States indicate that limited-area models forced by large-scale information at the lateral boundaries reproduce bulk temporal and spatial characteristics of meteorological fields. In particular, the 500 hPa height field time average and temporal variability are generally well simulated by all participating models. Model simulations of precipitation episodes vary depending on the scale of the dynamical forcing. Organized synoptic-scale precipitation systems are simulated deterministically in that precipitation occurs at close to the same time and location as observed (although amounts may vary from observations). Episodes of mesoscale and convective precipitation are represented in a more stochastic sense, with less precise agreement in temporal and spatial patterns. Simulated surface energy fluxes show broad similarity with the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) observations in their temporal evolution and time average diurnal cycle. Intermodel differences in midday Bowen ratio tend to be closely associated with precipitation differences. Differences in daily maximum temperatures also are linked to Bowen ratio differences, indicating strong local, surface influence on this field. Although some models have bias with respect to FIFE observations, all tend to reproduce the synoptic variability of observed daily maximum and minimum temperatures. Results also reveal the advantage of an intercomparison in exposing common tendencies of models despite their differences in convective and surface parameterizations and different methods of assimilating lateral boundary conditions.

Climate and climate change in western canada as simulated by the Canadian regional climate model

Abstract A þrst climate simulation performed with the novel Canadian Regional Climate Model (CRCM) is presented. The CRCM is based on fully elastic non‐hydrostatic þeld equations, which are solved with an efþcient semi‐implicit semi‐Lagrangian (SISL) marching algorithm, and on the parametrization package of subgrid‐scale physical effects of the second‐generation Canadian Global Climate Model (GCMII). Two 5‐year integrations of the CRCM nested with GCMII simulated data as lateral boundary conditions are made for conditions corresponding to current and doubled CO2 scenarios. For these simulations the CRCM used a grid size of 45 km on a polar‐stereographic projection, 20 scaled‐height levels and a time step of 15 min; the nesting GCMII has a spectral truncation of T32, 10 hybrid‐pressure levels and a time step of 20 min. These simulations serve to document: (1) the suitability of the SISL numerical scheme for regional climate modelling, (2) the use of GCMII physics at much higher resolution than in the nesti...

Potential for added value in precipitation simulated by high-resolution nested Regional Climate Models and observations

Regional Climate Models (RCMs) constitute the most often used method to perform affordable high-resolution regional climate simulations. The key issue in the evaluation of nested regional models is to determine whether RCM simulations improve the representation of climatic statistics compared to the driving data, that is, whether RCMs add value. In this study we examine a necessary condition that some climate statistics derived from the precipitation field must satisfy in order that the RCM technique can generate some added value: we focus on whether the climate statistics of interest contain some fine spatial-scale variability that would be absent on a coarser grid. The presence and magnitude of fine-scale precipitation variance required to adequately describe a given climate statistics will then be used to quantify the potential added value (PAV) of RCMs. Our results show that the PAV of RCMs is much higher for short temporal scales (e.g., 3-hourly data) than for long temporal scales (16-day average data) due to the filtering resulting from the time-averaging process. PAV is higher in warm season compared to cold season due to the higher proportion of precipitation falling from small-scale weather systems in the warm season. In regions of complex topography, the orographic forcing induces an extra component of PAV, no matter the season or the temporal scale considered. The PAV is also estimated using high-resolution datasets based on observations allowing the evaluation of the sensitivity of changing resolution in the real climate system. The results show that RCMs tend to reproduce relatively well the PAV compared to observations although showing an overestimation of the PAV in warm season and mountainous regions.

Climate simulation over CORDEX Africa domain using the fifth-generation Canadian Regional Climate Model (CRCM5)

The new fifth-generation Regional Climate Model (CRCM5) was driven by ERA reanalyses for the period 1984–2008 over the African continent following the CORDEX experimental protocol. Overall the model succeeds in reproducing the main features of the geographical distribution and seasonal cycle of temperature and precipitation, the diurnal cycle of precipitation, and the West African Monsoon (WAM). Biases in surface temperature and precipitation are discussed in relation with some circulation defects noted in the simulation. In the African regions near the equator, the model successfully reproduces the double peak of rainfall due to the double passage of the tropical rainbelt, although it better simulates the magnitude and timing of the second peak of precipitation. CRCM5 captures the timing of the monsoon onset for the Sahel region but underestimates the magnitude of precipitation. The simulated diurnal cycle is quite well simulated for all of the regions, but is always somewhat in advance for the timing of rainfall peak. In boreal summer the CRCM5 simulation exhibits a weak cold bias over the Sahara and the maximum temperature is located too far south, resulting in a southward bias in the position of the Saharan Heat Low. The region of maximum ascent in the deep meridional circulation of the Hadley cell is well located in the CRCM5 simulation, but it is somewhat too narrow. The core of the African Easterly Jet is of the right strength and almost at the right height, but it is displayed slightly southward, as a consequence of the southward bias in the position of the Saharan Heat Low and the thermal wind relationship. These biases appear to be germane to the WAM rainfall band being narrower and not moving far enough northward, resulting in a dry bias in the Sahel.

Internal Variability in Regional Climate Downscaling at the Seasonal Scale

To study the internal variability of the model and its consequences on seasonal statistics, large ensembles of twenty 3-month simulations of the Canadian Regional Climate Model (CRCM), differing only in their initial conditions, were generated over different domain sizes in eastern North America for a summer season. The degree of internal variability was measured as the spread between the individual members of the ensemble during the integration period. Results show that the CRCM internal variability depends strongly on synoptic events, as is seen by the pulsating behavior of the time evolution of variance during the period of integration. The existence of bimodal solutions for the circulation is also noted. The geographical distribution of variance depends on the variables; precipitation shows maximum variance in the southern United States, while 850-hPa geopotential height exhibits maximum variance in the northeast part of the domain. Results suggest that strong precipitation events in the southern United States may act as a triggering mechanism for the 850-hPa geopotential height spread along the storm track, which reaches its maximum toward the northeast of the domain. This study reveals that successive reductions of the domain size induce a general decrease in the internal variability of the model, but an important variation in its geographical distribution and amplitude was detected. The influence of the internal variability at the seasonal scale was evaluated by computing the variance between the individual member seasonal averages of the ensemble. Large values of internal variability for precipitation suggest possible repercussions of internal variability on seasonal statistics.