The effect of horizontal resolution on simulation of very extreme US precipitation events in a global atmosphere model
- First Online:
- Cite this article as:
- Wehner, M.F., Smith, R.L., Bala, G. et al. Clim Dyn (2010) 34: 241. doi:10.1007/s00382-009-0656-y
- 957 Downloads
We investigate the ability of a global atmospheric general circulation model (AGCM) to reproduce observed 20 year return values of the annual maximum daily precipitation totals over the continental United States as a function of horizontal resolution. We find that at the high resolutions enabled by contemporary supercomputers, the AGCM can produce values of comparable magnitude to high quality observations. However, at the resolutions typical of the coupled general circulation models used in the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, the precipitation return values are severely underestimated.
KeywordsExtreme precipitationClimate modelsReturn valueHigh resolution
Very extreme weather events are important due to their potential to have serious impacts on human and ecological systems. As the mean climate changes due to anthropogenic causes, such rare events are also expected to change (CCSP3.3 2008). The intensity of the extreme precipitation is expected to increase with the availability of atmospheric moisture, which follows the Clausius-Clapeyron relationship (Allen and Ingram 2002).
To make credible predictions of future changes in extreme weather events, it is reasonable to ask a model to accurately simulate the observations of the recent past.
Precipitation is generally not as well simulated as air temperature in global climate models. One reason is that precipitation is influenced by vertical motions on scales smaller than the model grid. These motions, as well as relevant cloud microphysical processes, are parameterized in typical climate models. Furthermore, precipitation is often influenced by orography, which is smoothed by the finite grid sizes used in climate models. Despite these issues, climate models can have good skill in simulating large scale patterns of mean precipitation such as the zonal mean distribution. However, models generally lack skill in accurately simulating regional distributions of precipitation and vary greatly from one model to another (Covey et al. 2000).
Very extreme precipitation events are well described by the statistical formalism of the Generalized Extreme Value (GEV) theory. The statistical character of the very extreme portion of a distribution is determined by relatively few events in samples of any reasonable size. In the case of precipitation, these events are highly localized in both space and time due to the episodic nature of strong storms. In this paper, we investigate the ability of climate models to simulate the very extreme tails of the distribution of precipitation events. By varying the horizontal resolution of a model, we quantify the effect of better representing the local nature of individual storms on the statistics of very extreme precipitation. Duffy et al. (2003) and Iorio et al. (2004) showed that high-resolution simulations using the NCAR Community Climate Model version 3 produce spatial patterns of seasonal-mean precipitation that agree more closely with observed precipitation patterns than do results from the same model at coarse resolution. Iorio et al. (2004) further find that higher-resolution simulations have more realistic daily precipitation statistics in North America.
The climate model used in this study is the Community Atmospheric Model version 2 (CAM2). We use the finite volume dynamics version of the model (fvCAM), for our simulations (Lin and Rood 1996). CAM2 is the fifth generation of the NCAR atmospheric GCM (http://www.ccsm.ucar.edu/models/atm-cam/docs/cam2.0/description/index.html). As used in this study, this global atmospheric model was configured according to the protocol dictated by the Atmospheric Model Intercomparison Project (AMIP) (Gates et al. 1999) using prescribed time varying values of the sea surface temperature and sea ice extent. Although the details of a stand alone atmospheric model simulation can differ from that of a fully coupled ocean atmosphere general circulation model simulation (Covey et al. 2003), we believe that the effect of changing horizontal atmospheric resolution is similar in the two classes of climate models. We performed fvCAM AMIP integrations over the period 1979–1994 under three different horizontal grid resolutions. The first used a 2.0° × 2.5° mesh, the second used a 1.0° × 1.25° mesh, and the third used 0.5° × 0.625° mesh. No model tuning specific to resolution was involved. The initial conditions were simply regridded to each of the meshes and the surface boundary conditions (sea surface temperature and sea ice extent) were obtained by a standard AMIP request. These correspond, respectively, to roughly 240, 120 and 60 km grid spacing at the equator. For reference, the CMIP3 database of coupled climate model simulations prepared for the IPCC AR4 report ranged from 400 to 110 km at the equator (http://www-pcmdi.llnl.gov).Vertical resolution was unaltered and kept at the default 26 levels. Daily precipitation totals were saved and extracted from each integration.
High quality observations of accumulated daily precipitation for this study were taken from the “NOAA CPC (Climate Prediction Center) .25 × .25 Daily US Unified Precipitation” dataset (Higgins et al. 2000). These observations are confined to the continental United States land areas and are aggregated from three sources of station rain gauge data gridded to a 0.25°× 0.25° grid. Between 8,000 and 13,000 stations were quality controlled and gridded to about 18000 grid points using a modified Cressman (1959) scheme. Hence, there are likely many grid points with no stations as well as many with multiple stations. The density of station data is least in the Western mountainous and desert regions. The model representation of a storm is meant as an average over a grid cell. This transformed observational dataset is the closest comparison that can be made to the models and is far superior to comparison to individual station data. Data used in this study was limited to the same time period as the simulations.
Hence, return values of annual or seasonal extrema are readily obtained by this inversion of the GEV distribution function after its three parameters have been estimated.
We find that the ability of the atmospheric general circulation model, fvCAM, to simulate the extreme precipitation over the continental United States is strongly dependent on horizontal resolution. Even when the observed daily precipitation rates are regridded to the model’s native mesh, the simulated return values are substantially lower than the observations at coarse resolutions. At the highest resolutions that we were able to integrate, approximately 60 km at the equator, the simulated 20 year return value of the annual maximum daily precipitation approaches or even exceeds the observed return value when calculated on the model grid. We interpret this as somewhat of a breakthrough resolution. On coarser meshes, the resolution itself is limiting the intensity of these extreme events, rather than any particular model parameterization defect. At the 60 km scale, resolution certainly still limits precipitation intensity, particularly for tropical cyclones. However, at least for fvCAM, this resolution has allowed simulated extreme rainfall over the continental US to be at least as intense as the observations on the same mesh. This allows the model defects to be exposed whether they are circulation errors or more localized model formulation problems.
Greater model fidelity changes the distribution of simulated daily precipitation. Part of the change is a simply a shift in the distribution, improving the simulation of longer term averages of precipitation toward observations (Duffy et al. 2003). However, another part of the change is the enabling of substantially stronger precipitation events due to better resolution of the atmospheric fluid dynamics and the surface boundary conditions. By analyzing the relative amounts of cumulus convective precipitation and large scale precipitation on the days when total precipitation is the annual maximum, we find that over the continental US, either process can be responsible for the extreme values depending on the location. The large scale precipitation parameterization produces the largest precipitation rates over the Western regions while cumulus convective parameterization does so over the Southeastern US. This is likely a reflection of the nature of the most intense storms in the different regions: large scales systems coming off the Pacific Ocean for the Western US and Atlantic tropical cyclones for the Southeastern US. Higher horizontal resolution increases the realism of both processes.
While the computational burden associated with high-resolution climate models is significantly higher than that of current generation production climate models (http://www.pcmdi.llnl.gov), recent technological trends are making such simulations feasible. Worley (2007) has measured the simulation rate of fvCAM (in a beta version 3) in the 60 km configuration to be about 3,650 times faster than real time on approximately 1,000 processors of a Cray XT4 supercomputer. This compares with the simulation rate of about 1,600 times faster than real time realized on US-based computers when the Community Climate System Model (CCSM3.0) was run for the IPCC AR4 report. In fact, we estimate that for a 0.25° × 0.375° version of fvCAM, a simulation rate of about 900 times faster than real time could be presently obtained on 2,000 processors of a Cray XT4 (Wehner et al. 2008). This corresponds to about a 30 km grid spacing at the equator and is an appealing high resolution as it is still coarse enough that statistical issues in cumulus cloud process parameterization are avoided (Arakawa, private communication 1999). Also, the hydrostatic approximation employed in many current climate models is generally valid at this grid spacing.
The benefits of increased horizontal resolution to the quality of global atmosphere simulations extend well beyond those presented here for very extreme precipitation events. For instance, tropical cyclones are simulated well enough at these high resolutions to permit a direct investigation of the effect of anthropogenic climate change on these storms (Oouchi et al. 2006). Also, larger scale processes have been demonstrated (Duffy et al. 2003) to be very sensitive to horizontal grid spacing, with significant improvements to seasonal averages obtained at ~60 km resolutions over typical current production models. Regional scale climate change projections and societal impacts assessments also require higher resolution than is provided by the current typical global models. Since temporal extremes have disproportionate societal impacts, it is convenient that models using finer spatial grids also represent these extremes better.
This work was performed under the auspices of the U.S. Department of Energy (DOE) by the Lawrence Berkeley National Laboratory (LBNL) under contract No. DE-AC03-76SF00098 (LBNL) and with support from the DOE Climate Change Prediction Program. CPC US Unified Precipitation data provided by the NOAA/OAR/ESRL PSD, Boulder, Colorado, USA, from their web site at http://www.cdc.noaa.gov/.
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.