The first multi-model ensemble of regional climate simulations at kilometer-scale resolution part 2: historical and future simulations of precipitation

This paper presents the first multi-model ensemble of 10-year, “convection-permitting” kilometer-scale regional climate model (RCM) scenario simulations downscaled from selected CMIP5 GCM projections for historical and end of century time slices. The technique is to first downscale the CMIP5 GCM projections to an intermediate 12–15 km resolution grid using RCMs, and then use these fields to downscale further to the kilometer scale. The aim of the paper is to provide an overview of the representation of the precipitation characteristics and their projected changes over the greater Alpine domain within a Coordinated Regional Climate Downscaling Experiment Flagship Pilot Study and the European Climate Prediction system project, tasked with investigating convective processes at the kilometer scale. An ensemble of 12 simulations performed by different research groups around Europe is analyzed. The simulations are evaluated through comparison with high resolution observations while the complementary ensemble of 12 km resolution driving models is used as a benchmark to evaluate the added value of the convection-permitting ensemble. The results show that the kilometer-scale ensemble is able to improve the representation of fine scale details of mean daily, wet-day/hour frequency, wet-day/hour intensity and heavy precipitation on a seasonal scale, reducing uncertainty over some regions. It also improves the representation of the summer diurnal cycle, showing more realistic onset and peak of convection. The kilometer-scale ensemble refines and enhances the projected patterns of change from the coarser resolution simulations and even modifies the sign of the precipitation intensity change and heavy precipitation over some regions. The convection permitting simulations also show larger changes for all indices over the diurnal cycle, also suggesting a change in the duration of convection over some regions. A larger positive change of frequency of heavy to severe precipitation is found. The results are encouraging towards the use of convection-permitting model ensembles to produce robust assessments of the local impacts of future climate change.

     

对流允许尺度区域气候模拟的研究进展

与以往的区域气候模式相比,对流允许区域气候模式不再依赖于对流参数化方案,其精细的分辨率可以显式表示深对流过程,在夏季对流降水的日变化和极端降水事件模拟等方面具有明显增值能力,是区域气候模拟的发展方向。对现有的对流允许尺度区域气候模拟研究进行了较为详细的回顾和介绍,简述了对流允许尺度区域气候模式中比较重要的物理过程及外部驱动条件的影响,总结了以往对流允许尺度区域气候模拟的研究成果以及当下所面临的挑战和对未来的展望,以期对中国及东亚区域对流允许区域气候模拟的研究提供有益参考。诸多研究表明,对流允许区域气候模拟作为一种有前景的气候模式,可提供更加可靠的区域尺度的气候信息。     

Extreme precipitation in WRF during the Newcastle East Coast Low of 2007

In the context of regional downscaling, we study the representation of extreme precipitation in the Weather Research and Forecasting (WRF) model, focusing on a major event that occurred on the 8^th of June 2007 along the coast of eastern Australia (abbreviated “Newy”). This was one of the strongest extra-tropical low-pressure systems off eastern Australia in the last 30 years and was one of several storms comprising a test bed for the WRF ensemble that underpins the regional climate change projections for eastern Australia (New South Wales/Australian Capital Territory Regional Climate Modelling Project, NARCliM). Newy provides an informative case study for examining precipitation extremes as simulated by WRF set up for regional downscaling. Here, simulations from the NARCliM physics ensemble of Newy available at ～10 km grid spacing are used. Extremes and spatio-temporal characteristics are examined using land-based daily and hourly precipitation totals, with a particular focus on hourly accumulations. Of the different physics schemes assessed, the cumulus and the boundary layer schemes cause the largest differences. Although the Betts-Miller-Janjic cumulus scheme produces better rainfall totals over the entire storm, the Kain-Fritsch cumulus scheme promotes higher and more realistic hourly extreme precipitation totals. Analysis indicates the Kain-Fritsch runs are correlated with larger resolved grid-scale vertical moisture fluxes, which are produced through the influence of parameterized convection on the larger-scale circulation and the subsequent convergence and ascent of moisture. Results show that WRF qualitatively reproduces spatial precipitation patterns during the storm, albeit with some errors in timing. This case study indicates that whilst regional climate simulations of an extreme event such as Newy in WRF may be well represented at daily scales irrespective of the physics scheme used, the representation at hourly scales is likely to be physics scheme dependent.     

Simulating the diurnal cycle of rainfall in global climate models: resolution versus parameterization

The effects of horizontal resolution and the treatment of convection on simulation of the diurnal cycle of precipitation during boreal summer are analyzed in several innovative weather and climate model integrations. The simulations include: season-long integrations of the Non-hydrostatic Icosahedral Atmospheric Model (NICAM) with explicit clouds and convection; year-long integrations of the operational Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts at three resolutions (125, 39 and 16 km); seasonal simulations of the same model at 10 km resolution; and seasonal simulations of the National Center for Atmospheric Research (NCAR) low-resolution climate model with and without an embedded two-dimensional cloud-resolving model in each grid box. NICAM with explicit convection simulates best the phase of the diurnal cycle, as well as many regional features such as rainfall triggered by advancing sea breezes or high topography. However, NICAM greatly overestimates mean rainfall and the magnitude of the diurnal cycle. Introduction of an embedded cloud model within the NCAR model significantly improves global statistics of the seasonal mean and diurnal cycle of rainfall, as well as many regional features. However, errors often remain larger than for the other higher-resolution models. Increasing resolution alone has little impact on the timing of daily rainfall in IFS with parameterized convection, yet the amplitude of the diurnal cycle does improve along with the representation of mean rainfall. Variations during the day in atmospheric prognostic fields appear quite similar among models, suggesting that the distinctive treatments of model physics account for the differences in representing the diurnal cycle of precipitation.     

Climate change scenarios of precipitation extremes in Central Europe from ENSEMBLES regional climate models

The study examines future scenarios of precipitation extremes over Central Europe in an ensemble of 12 regional climate model (RCM) simulations with the 25-km resolution, carried out within the European project ENSEMBLES. We apply the region-of-influence method as a pooling scheme when estimating distributions of extremes, which consists in incorporating data from a ‘region’ (set of gridboxes) when fitting an extreme value distribution in any single gridbox. The method reduces random variations in the estimates of parameters of the extreme value distribution that result from large spatial variability of heavy precipitation. Although spatial patterns differ among the models, most RCMs simulate increases in high quantiles of precipitation amounts when averaged over the area for the late-twenty-first century (2070–2099) climate in both winter and summer. The sign as well as the magnitude of the projected change vary only little for individual parts of the distribution of daily precipitation in winter. In summer, on the other hand, the projected changes increase with the quantile of the distribution in all RCMs, and they are negative (positive) for parts of the distribution below (above) the 98% quantile if averaged over the RCMs. The increases in precipitation extremes in summer are projected in spite of a pronounced drying in most RCMs. Although a rather general qualitative agreement of the models concerning the projected changes of precipitation extremes is found in both winter and summer, the uncertainties in climate change scenarios remain large and would likely further increase considerably if a more complete ensemble of RCM simulations driven by a larger suite of global models and with a range of possible scenarios of the radiative forcing is available.     

四川盆地夏季降水日变化的数值模拟

利用区域气候模式RegCM3对1991-2004年四川盆地夏季降水进行了数值模拟,通过模拟结果和NCEP/NCAR再分析资料的对比,评估了模式对四川盆地夏季＂夜雨＂现象的模拟能力。结果表明,RegCM3模式能较好地模拟出四川盆地夏季降水的空间分布和日变化规律,四川盆地夏季＂夜雨＂现象的形成与该地区的地形分布有密切关系。...     

Evaluation of the regional climate model ALADIN to simulate the climate over North America in the CORDEX framework

In this study, an ensemble of four multi-year climate simulations is performed with the regional climate model ALADIN to evaluate its ability to simulate the climate over North America in the CORDEX framework. The simulations differ in their driving fields (ERA-40 or ERA-Interim) and the nudging technique (with or without large-scale nudging). The validation of the simulated 2-m temperature and precipitation with observationally-based gridded data sets shows that ALADIN performs similarly to other regional climate models that are commonly used over North America. Large-scale nudging improves the temporal correlation of the atmospheric circulation between ALADIN and its driving field, and also reduces the warm and dry summer biases in central North America. The differences between the simulations driven with different reanalyses are small and are likely related to the regional climate model’s induced internal variability. In general, the impact of different driving fields on ALADIN is smaller than that of large-scale nudging. The analysis of the multi-year simulations over the prairie and the east taiga indicates that the ALADIN 2-m temperature and precipitation interannual variability is similar or larger than that observed. Finally, a comparison of the simulations with observations for the summer 1993 shows that ALADIN underestimates the flood in central North America mainly due to its systematic dry bias in this region. Overall, the results indicate that ALADIN can produce a valuable contribution to CORDEX over North America.     

Assessing future changes in seasonal climatic extremes in the Ganges river basin using an ensemble of regional climate models

Using an ensemble of four high resolution (~25 km) regional climate models, this study analyses the future (2021–2050) spatial distribution of seasonal temperature and precipitation extremes in the Ganges river basin based on the SRES A1B emissions scenario. The model validation results (1989–2008) show that the models simulate seasonality and spatial distribution of extreme temperature events better than precipitation. The models are able to capture fine topographical detail in the spatial distribution of indices based on their ability to resolve processes at a higher regional resolution. Future simulations of extreme temperature indices generally agree with expected warming in the Ganges basin, with considerable seasonal and spatial variation. Significantly warmer summers in the central part of the basin along with basin-wide increase in night temperature are expected during the summer and monsoon months. An increase in heavy precipitation indices during monsoon, coupled with extended periods without precipitation during the winter months; indicates an increase in the incidence of extreme events.     

Probabilistic Precipitation Forecasting Based on Ensemble Output Using Generalized Additive Models and Bayesian Model Averaging

A probabilistic precipitation forecasting model using generalized additive models (GAMs) and Bayesian model averaging (BMA) was proposed in this paper.GAMs were used to fit the spatial-temporal precipi...     

Climate-change effects on extreme precipitation in central Europe: uncertainties of scenarios based on regional climate models

Observations as well as most climate model simulations are generally in accord with the hypothesis that the hydrologic cycle should intensify and become highly volatile with the greenhouse-gas-induced climate change, although uncertainties of these projections as well as the spatial and seasonal variability of the changes are much larger than for temperature extremes. In this study, we examine scenarios of changes in extreme precipitation events in 24 future climate runs of ten regional climate models, focusing on a specific area of the Czech Republic (central Europe) where complex orography and an interaction of other factors governing the occurrence of heavy precipitation events result in patterns that cannot be captured by global models. The peaks-over-threshold analysis with increasing threshold censoring is applied to estimate multi-year return levels of daily rainfall amounts. Uncertainties in scenarios of changes for the late 21st century related to the inter-model and within-ensemble variability and the use of the SRES-A2 and SRES-B2 greenhouse gas emission scenarios are evaluated. The results show that heavy precipitation events are likely to increase in severity in winter and (with less agreement among models) also in summer. The inter-model and intra-model variability and related uncertainties in the pattern and magnitude of the change is large, but the scenarios tend to agree with precipitation trends recently observed in the area, which may strengthen their credibility. In most scenario runs, the projected change in extreme precipitation in summer is of the opposite sign than a change in mean seasonal totals, the latter pointing towards generally drier conditions in summer. A combination of enhanced heavy precipitation amounts and reduced water infiltration capabilities of a dry soil may severely increase peak river discharges and flood-related risks in this region.     

Role of orography in modulating space–time scales of convection over South Asia

Propagation of convective systems in the meridional direction during boreal summer is responsible for active and break phases of monsoon over south Asia. This region is unique in the world in its characteristics of monsoon variability and is in close proximity of mountains like the Himalayas. Here, using an atmospheric general circulation model, we try to understand the role of orography in determining spatial and temporal scales of these convective systems. Absence of orography (noGlOrog) decreased the simulated seasonal mean precipitation over India by 23 % due to delay in onset by about a month vis-a-vis the full-mountain case. In noGlOrog, poleward propagations were absent during the delayed period prior to onset. Post-onset, both simulations had similar patterns of poleward propagations. The spatial and temporal scales of propagating clouds bands were determined using wavelet analysis. These scales were found to be different in full-mountain and no-mountain experiments in June–July. However, after the onset of monsoon in noGlOrog, these scales become similar to that with orography. Simulations with two different sets of convection schemes confirmed this result. Further analysis shows that the absence (presence) of meridional propagations during early (late) phase of summer monsoon in noGlOrog was associated with weaker (stronger) vertical shear of zonal wind over south Asia. Our study shows that orography plays a major role in determining the time of onset over the Indian region. However, after onset, basic characteristics of propagating convective systems and therefore the monthly precipitation over India, are less sensitive to the presence of orography and are modulated by moist convective processes.     

MM5 simulations of interannual change and the diurnal cycle of southern African regional climate

Summary Two cumulus convection and two planetary boundary layer schemes are used to investigate the climate of southern Africa using the MM5 regional climate model. Both a wet (1988/89) and a dry (1991/92) summer (December–February, DJF) rainfall season are simulated and the results compared with three different observational sources: Climate Research Unit seasonal data (precipitation, 2 m surface temperature, number of rain days), satellite-derived diurnal precipitation and the Surface Radiation Budget diurnal short-wave fluxes and optical depth. Using the ETA model boundary layer in MM5 simulates too much incident short-wave radiation at the surface at 12 UTC, whereas the medium range forecast model boundary layer yields a diurnal cycle of short-wave radiation closer to the observed. The Betts-Miller convection scheme in MM5 simulates peak rainfall later in the day and less rain days than observed, whereas when using the Kain-Fritsch convection scheme a peak rainfall earlier in the day and more rain days than observed are simulated. The intensity of the hydrological cycle is therefore dependent on the choice of convection scheme, which in turn is further modified by the boundary layer scheme. Precipitation during the wet 1988/89 season is reasonably captured by most simulations, though using the Betts-Miller scheme more accurately simulates rainfall during the dry 1991/92 season. Mean DJF biases in the surface temperature and diurnal temperature range are consistent with biases in the number of rain days and the diurnal cycles of surface moisture and energy.     

Future evolution of extreme precipitation in the Mediterranean

Mediterranean basins can be impacted by severe floods caused by extreme rainfall, and there is a growing awareness about the possible increase in these heavy rainfall events due to climate change. In this study, the climate change impacts on extreme daily precipitation in 102 catchments covering the whole Mediterranean basin are investigated using nonstationary extreme value model applied to annual maximum precipitation in an ensemble of high-resolution regional climate model (RCM) simulations from the Euro-CORDEX experiment. Results indicate contrasted trends, with significant increasing trends in Northern catchments and conversely decreasing trends in Southern catchments. For most cases, the time of signal emergence for these trends is before the year 2000. The same spatial pattern is obtained under the two climate scenarios considered (RCP4.5 and RCP8.5) and in most RCM simulations, suggesting a robust climate change signal. The strongest multi-model agreement concerns the positive trends, which can exceed +?20% by the end of the twenty-first century in some simulations, impacting South France, North Italy, and the Balkans. For these areas, society-relevant strong impacts of such Mediterranean extreme precipitation changes could be expected in particular concerning flood-related damages.     

A review of cloud-resolving model studies of convective processes

Convective processes affect large-scale environments through cloud-radiation interaction, cloud micro- physical processes, and surface rainfall processes. Over the last three decades, cloud-resolving models （CRMs） have demonstrated to be capable of simulating convective-radiative responses to an imposed large-scale forcing. The CRM-produced cloud and radiative properties have been utilized to study the convective- related processes and their ensemble effects on large-scale circulations. This review the recent progress on the understanding of convective processes with the use of CRM simulations, including precipitation processes; cloud microphysical and radiative processes; dynamical processes; precipitation efficiency; diurnal variations of tropical oceanic convection; local-scale atmosphere-ocean coupling processes; and tropical convective-radiative equilibrium states. Two different ongoing applications of CRMs to general circulation models （GCMs） are discussed： replacing convection and cloud schemes for studying the interaction between cloud systems and large-scale circulation, and improving the schemes for climate simulations.     

The thermodynamic influence of subgrid orography in a global climate model

Assessments of the impacts of climate change typically require information at scales of 10 km or less. Such a resolution in global climate simulations is unlikely for at least two decades. We have developed an alternative to explicit resolution that provides a framework for meeting the needs of climate change impact assessment much sooner. We have applied to a global climate model a physically based subgrid-scale treatment of the influence of orography on temperature, clouds, precipitation, and land surface hydrology. The treatment represents subgrid variations in surface elevation in terms of fractional area distributions of discrete elevation classes. For each class it calculates the height rise/descent of air parcels traveling through the grid cell, and applies the influence of the rise/descent to the temperature and humidity profiles of the elevation class. Cloud, radiative, and surface processes are calculated separately for each elevation class using the same physical parametrizations used by the model without the subgrid orography parametrization. The simulated climate fields for each elevation class can then be distributed in post-processing according to the spatial distribution of surface elevation within each grid cell. Parallel 10-year simulations with and without the subgrid treatment have been performed. The simulated temperature, precipitation and snow water are mapped to 2.5-minute (~5 km) resolution and compared with gridded analyses of station measurements. The simulation with the subgrid scheme produces a much more realistic distribution of snow water and significantly more realistic distributions of temperature and precipitation than the simulation without the subgrid scheme. Moreover, the 250-km grid cell means of most other fields are virtually unchanged by the subgrid scheme. This suggests that the tuning of the climate model without the subgrid scheme is also applicable to the model with the scheme.     

The influence of synoptic airflow on UK daily precipitation extremes. Part I: Observed spatio-temporal relationships

We study the influence of synoptic scale atmospheric circulation on extreme daily precipitation across the United Kingdom, using observed time series from 689 rain gauges. To this end we employ a statistical model, that uses airflow strength, direction and vorticity as predictors for the generalised extreme value distribution of monthly precipitation maxima. The inferred relationships are connected with the dominant westerly flow, the orography, and the moisture supply from surrounding seas. We aggregated the results for individual rain gauges to regional scales to investigate the temporal variability of extreme precipitation. Airflow explains a significant fraction of the variability on subannual to decadal time scales. A large fraction of the especially heavy winter precipitation during the 1980s and 1990s in north Scotland can be attributed to a prevailing positive phase of the North Atlantic Oscillation. Our statistical model can be used for statistical downscaling and to validate regional climate model output.     

Performance of a multi-RCM ensemble for South Eastern South America

The ability of four regional climate models to reproduce the present-day South American climate is examined with emphasis on La Plata Basin. Models were integrated for the period 1991–2000 with initial and lateral boundary conditions from ERA-40 Reanalysis. The ensemble sea level pressure, maximum and minimum temperatures and precipitation are evaluated in terms of seasonal means and extreme indices based on a percentile approach. Dispersion among the individual models and uncertainties when comparing the ensemble mean with different climatologies are also discussed. The ensemble mean is warmer than the observations in South Eastern South America (SESA), especially for minimum winter temperatures with errors increasing in magnitude towards the tails of the distributions. The ensemble mean reproduces the broad spatial pattern of precipitation, but overestimates the convective precipitation in the tropics and the orographic precipitation along the Andes and over the Brazilian Highlands, and underestimates the precipitation near the monsoon core region. The models overestimate the number of wet days and underestimate the daily intensity of rainfall for both seasons suggesting a premature triggering of convection. The skill of models to simulate the intensity of convective precipitation in summer in SESA and the variability associated with heavy precipitation events (the upper quartile daily precipitation) is far from satisfactory. Owing to the sparseness of the observing network, ensemble and observations uncertainties in seasonal means are comparable for some regions and seasons.     

Field significance of performance measures in the context of regional climate model evaluation. Part 2: precipitation

A new approach for rigorous spatial analysis of the downscaling performance of regional climate model (RCM) simulations is introduced. It is based on a multiple comparison of the local tests at the grid cells and is also known as ‘field’ or ‘global’ significance. The block length for the local resampling tests is precisely determined to adequately account for the time series structure. New performance measures for estimating the added value of downscaled data relative to the large-scale forcing fields are developed. The methodology is exemplarily applied to a standard EURO-CORDEX hindcast simulation with the Weather Research and Forecasting (WRF) model coupled with the land surface model NOAH at 0.11 ^° grid resolution. Daily precipitation climatology for the 1990–2009 period is analysed for Germany for winter and summer in comparison with high-resolution gridded observations from the German Weather Service. The field significance test controls the proportion of falsely rejected local tests in a meaningful way and is robust to spatial dependence. Hence, the spatial patterns of the statistically significant local tests are also meaningful. We interpret them from a process-oriented perspective. While the downscaled precipitation distributions are statistically indistinguishable from the observed ones in most regions in summer, the biases of some distribution characteristics are significant over large areas in winter. WRF-NOAH generates appropriate stationary fine-scale climate features in the daily precipitation field over regions of complex topography in both seasons and appropriate transient fine-scale features almost everywhere in summer. As the added value of global climate model (GCM)-driven simulations cannot be smaller than this perfect-boundary estimate, this work demonstrates in a rigorous manner the clear additional value of dynamical downscaling over global climate simulations. The evaluation methodology has a broad spectrum of applicability as it is distribution-free, robust to spatial dependence, and accounts for time series structure.     

Dynamical downscaling of ERA-40 in complex terrain using the WRF regional climate model

Results from a first-time employment of the WRF regional climate model to climatological simulations in Europe are presented. The ERA-40 reanalysis (resolution 1°) has been downscaled to a horizontal resolution of 30 and 10?km for the period of 1961?C1990. This model setup includes the whole North Atlantic in the 30?km domain and spectral nudging is used to keep the large scales consistent with the driving ERA-40 reanalysis. The model results are compared against an extensive observational network of surface variables in complex terrain in Norway. The comparison shows that the WRF model is able to add significant detail to the representation of precipitation and 2-m temperature of the ERA-40 reanalysis. Especially the geographical distribution, wet day frequency and extreme values of precipitation are highly improved due to the better representation of the orography. Refining the resolution from 30 to 10?km further increases the skill of the model, especially in case of precipitation. Our results indicate that the use of 10-km resolution is advantageous for producing regional future climate projections. Use of a large domain and spectral nudging seems to be useful in reproducing the extreme precipitation events due to the better resolved synoptic scale features over the North Atlantic, and also helps to reduce the large regional temperature biases over Norway. This study presents a high-resolution, high-quality climatological data set useful for reference climate impact studies.     

Validation of precipitation over Japan during 1985–2004 simulated by three regional climate models and two multi-model ensemble means

We dynamically downscaled Japanese reanalysis data (JRA-25) for 60 regions of Japan using three regional climate models (RCMs): the Non-Hydrostatic Regional Climate Model (NHRCM), modified RAMS version 4.3 (NRAMS), and modified Weather Research and Forecasting model (TWRF). We validated their simulations of the precipitation climatology and interannual variations of summer and winter precipitation. We also validated precipitation for two multi-model ensemble means: the arithmetic ensemble mean (AEM) and an ensemble mean weighted according to model reliability. In the 60 regions NRAMS simulated both the winter and summer climatological precipitation better than JRA-25, and NHRCM simulated the wintertime precipitation better than JRA-25. TWRF, however, overestimated precipitation in the 60 regions in both the winter and summer, and NHRCM overestimated precipitation in the summer. The three RCMs simulated interannual variations, particularly summer precipitation, better than JRA-25. AEM simulated both climatological precipitation and interannual variations during the two seasons more realistically than JRA-25 and the three RCMs overall, but the best RCM was often superior to the AEM result. In contrast, the weighted ensemble mean skills were usually superior to those of the best RCM. Thus, both RCMs and multi-model ensemble means, especially multi-model ensemble means weighted according to model reliability, are powerful tools for simulating seasonal and interannual variability of precipitation in Japan under the current climate.