Low (web) Quality - BALTEX
Low (web) Quality - BALTEX
Low (web) Quality - BALTEX
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
16<br />
Assessment of precipitation as simulated by a RCM and its driving data<br />
Alejandro Di Luca 1,3,4 , Ramón de Elía 2,3,4 and René Laprise 1,3,4<br />
1 Université du Québec à Montréal (UQÀM), Montréal (Québec), Canada, diluca@sca.uqam.ca; 2 Consortium Ouranos,<br />
Montréal (Québec), Canada; 3 Canadian Network for Regional Climate Modelling and Diagnostics (CRCMD), Montréal ,<br />
Canada; 4 Centre pour l’Étude et la Simulation du Climat à l’Échelle Régionale (ESCER), Montréal (Québec), Canada<br />
1. Introduction<br />
The primary and most comprehensive tools to study future<br />
climate are the Atmosphere-Ocean General Circulation<br />
Models (AOGCMs). Present horizontal grid intervals of the<br />
atmospheric component of AOGCMs are insufficient to<br />
capture the fine-scale structure of climatic. In this context,<br />
an alternative to obtain future regional climate projections is<br />
the use of high-resolution Regional Climate Models<br />
(RCMs), nested at their lateral boundaries with lowresolution<br />
AOGCMs (Giorgi and Bates, 1989; Laprise et al.,<br />
2008). Because differences between AOGCMs and RCMs<br />
are mainly their resolution, the small scales represent the<br />
main potential added value of the high resolution RCM over<br />
the AOGCM. We say “potential added value” because the<br />
RCM simulation should satisfy several conditions before<br />
this potentiality becomes effective. Among them: the one<br />
way nesting technique must be reliable in the sense that it<br />
should allow a good development of the fine scale with no<br />
amplification of driving fields errors; overall regional<br />
models errors should be smaller or equal than those from the<br />
global model; and, the considered variable must contain<br />
information of fine scales.<br />
In this work we focus on the models specific errors, and in<br />
order to do so we have studied daily precipitation as<br />
simulated by the Canadian RCM (CRCM) and the Canadian<br />
GCM (CGCM).<br />
2. Methodology and data<br />
The methodology used to investigate the presence of added<br />
value in CRCM simulations is based on the assessment of<br />
the RCM performance when compared to its driving model<br />
and observed data. We have evaluated some statistics of<br />
daily precipitation as simulated by both models in several<br />
regions across Canada (see Figure 1) during the period 1971<br />
- 1990.<br />
CGCM cumulative precipitation data was archived at 24-<br />
hour intervals forces to carry the comparison at this time<br />
scale, thus discarding shorter time interval information.<br />
As a result, the considered variable does not explicitly<br />
contain spatio-temporal fine scale information produced<br />
by the CRCM. The hypothesis of the existence of added<br />
value then lies not in the presence of fine scale<br />
information but in the assumption that the global model is<br />
expected to have little skill near its truncation limit<br />
(Laprise, 2003; Feser, 2006).<br />
The global model used in this study is the third generation<br />
of the Canadian Centre for Climate Modelling and<br />
Analysis Coupled Global Climate Model (CGCM3). The<br />
gridded output of precipitation occurs on a 96 by 48<br />
Gaussian grid (output data has a grid spacing of 3.75° in<br />
latitude and longitude).<br />
The CRCM simulations were performed at the Ouranos<br />
Consortium with horizontal grid spacing of 45 km (true at<br />
60° N) over a North American domain with a total of 201<br />
by 193 grid points. Two CRCM simulations were<br />
considered in the present investigation differing only in<br />
the lateral boundary conditions used as nesting data. One<br />
simulation is driven by the CGCM and will be designated<br />
as CRCM (CGCM). The other simulation is nested by the<br />
National Centers for Environmental Prediction (NCEP) -<br />
National Center for Atmospheric Research (NCAR)<br />
reanalyses and will be designated as CRCM (NCEP).<br />
3. Intensity frequency distributions<br />
Intensity frequency distributions are constructed with bin<br />
sizes that vary logarithmically in order to account for the<br />
reduction on the number of events with increasing<br />
intensity. The frequency of each category is calculated<br />
using the following thresholds: 1, 2, 4, 8, 16, 32, 64, 128,<br />
256 and 512 mm/day. As an example, Figure 2 shows the<br />
intensity frequency distributions as observed and<br />
simulated in the BC region for the winter season.<br />
QC<br />
BC<br />
ALTA SAS MAN<br />
Figure 1. Specification of areas of interest. Blue boxes<br />
indicate regions including one CGCM grid point.<br />
A direct comparison between a RCM, a GCM and observed<br />
values can be properly done when quantities are equivalent<br />
for the three sources of data. In this work, the evaluation is<br />
carried out at scales that are greater or equal than that of the<br />
coarser resolution data. In our case, the CGCM defines the<br />
minimum area (i.e., one CGCM grid box) at which perform<br />
the comparison and this forces as to transforms high<br />
resolution data into lower resolution. Upscaling RCM and<br />
observed results to the GCM level is simply done by<br />
computing the spatial-average of all grid-points and stations<br />
data within each GCM grid box. Similarly, the fact that the<br />
Figure 2. Intensity frequency distributions of precipitation<br />
rate in BC for wintertime.<br />
A simple score S (Perkins et al., 2007) that measures the<br />
overlap between simulated (f r mod ) and observed (f r obs )<br />
intensity frequency distributions is defined for each region<br />
r with the aim of obtaining an objective comparison,<br />
S r mod =<br />
9<br />
∑<br />
k=0<br />
( (), f obs r () k )<br />
min f r mod k