Evaluation of water and energy budgets in regional climate models applied over Europe

Abstract This study presents a model intercomparison of four regional climate models (RCMs) and one variable resolution atmospheric general circulation model (AGCM) applied over Europe with special focus on the hydrological cycle and the surface energy budget. The models simulated the 15 years from 1979 to 1993 by using quasi-observed boundary conditions derived from ECMWF re-analyses (ERA). The model intercomparison focuses on two large atchments representing two different climate conditions covering two areas of major research interest within Europe. The first is the Danube catchment which represents a continental climate dominated by advection from the surrounding land areas. It is used to analyse the common model error of a too dry and too warm simulation of the summertime climate of southeastern Europe. This summer warming and drying problem is seen in many RCMs, and to a less extent in GCMs. The second area is the Baltic Sea catchment which represents maritime climate dominated by advection from the ocean and from the Baltic Sea. This catchment is a research area of many studies within Europe and also covered by the BALTEX program. The observed data used are monthly mean surface air temperature, precipitation and river discharge. For all models, these are used to estimate mean monthly biases of all components of the hydrological cycle over land. In addition, the mean monthly deviations of the surface energy fluxes from ERA data are computed. Atmospheric moisture fluxes from ERA are compared with those of one model to provide an independent estimate of the convergence bias derived from the observed data. These help to add weight to some of the inferred estimates and explain some of the discrepancies between them. An evaluation of these biases and deviations suggests possible sources of error in each of the models. For the Danube catchment, systematic errors in the dynamics cause the prominent summer drying problem for three of the RCMs, while for the fourth RCM this is related to deficiencies in the land surface parametrization. The AGCM does not show this drying problem. For the Baltic Sea catchment, all models similarily overestimate the precipitation throughout the year except during the summer. This model deficit is probably caused by the internal model parametrizations, such as the large-scale condensation and the convection schemes.

Predictability and uncertainty in a regional climate model

The evaluation of the quality and usefulness of climate modeling systems is dependent upon an assessment of both the limited predictability of the climate system and the uncertainties stemming from model formulation. In this study a methodology is presented that is suited to assess the performance of a regional climate model (RCM), based on its ability to represent the natural interannual variability on monthly and seasonal timescales. The methodology involves carrying out multiyear ensemble simulations (to assess the predictability bounds within which the model can be evaluated against observations) and multiyear sensitivity experiments using different model formulations (to assess the model uncertainty). As an example application, experiments driven by assimilated lateral boundary conditions and sea surface temperatures from the ECMWF Reanalysis Project (ERA-15, 1979–1993) were conducted. While the ensemble experiment demonstrates that the predictability of the regional climate varies strongly between different seasons and regions, being weakest during the summer and over continental regions, important sensitivities of the modeling system to parameterization choices are uncovered. In particular, compensating mechanisms related to the long-term representation of the water cycle are revealed, in which summer dry and hot conditions at the surface, resulting from insufficient evaporation, can persist despite insufficient net solar radiation (a result of unrealistic cloud-radiative feedbacks).

Realism of rainfall in a very high-resolution regional climate model

The realistic representation of rainfall on the local scale in climate models remains a key challenge. Realism encompasses the full spatial and temporal structure of rainfall, and is a key indicator of model skill in representing the underlying processes. In particular, if rainfall is more realistic in a climate model, there is greater confidence in its projections of future change. In this study, the realism of rainfall in a very high-resolution (1.5 km) regional climate model (RCM) is compared to a coarser-resolution 12-km RCM. This is the first time a convection-permitting model has been run for an extended period (1989–2008) over a region of the United Kingdom, allowing the characteristics of rainfall to be evaluated in a climatological sense. In particular, the duration and spatial extent of hourly rainfall across the southern United Kingdom is examined, with a key focus on heavy rainfall. Rainfall in the 1.5-km RCM is found to be much more realistic than in the 12-km RCM. In the 12-km RCM, heavy rain events are not heavy enough, and tend to be too persistent and widespread. While the 1.5-km model does have a tendency for heavy rain to be too intense, it still gives a much better representation of its duration and spatial extent. Long-standing problems in climate models, such as the tendency for too much persistent light rain and errors in the diurnal cycle, are also considerably reduced in the 1.5-km RCM. Biases in the 12-km RCM appear to be linked to deficiencies in the representation of convection.

Origin of regional climate differences: role of boundary conditions and model formulation in two GCMs

Model differences in projections of extratropical regional climate change due to increasing greenhouse gases are investigated using two atmospheric general circulation models (AGCMs): ECHAM4 (Max Planck Institute, version 4) and CCM3 (National Center for Atmospheric Research Community Climate Model version 3). Sea-surface temperature (SST) fields calculated from observations and coupled versions of the two models are used to force each AGCM in experiments based on time-slice methodology. Results from the forced AGCMs are then compared to coupled model results from the Coupled Model Intercomparison Project 2 (CMIP2) database. The time-slice methodology is verified by showing that the response of each model to doubled CO2 and SST forcing from the CMIP2 experiments is consistent with the results of the coupled GCMs. The differences in the responses of the models are attributed to (1) the different tropical SST warmings in the coupled simulations and (2) the different atmospheric model responses to the same tropical SST warmings. Both are found to have important contributions to differences in implied Northern Hemisphere (NH) winter extratropical regional 500 mb height and tropical precipitation climate changes. Forced teleconnection patterns from tropical SST differences are primarily responsible for sensitivity differences in the extratropical North Pacific, but have relatively little impact on the North Atlantic. There are also significant differences in the extratropical response of the models to the same tropical SST anomalies due to differences in numerical and physical parameterizations. Differences due to parameterizations dominate in the North Atlantic. Differences in the control climates of the two coupled models from the current climate, in particular for the coupled model containing CCM3, are also demonstrated to be important in leading to differences in extratropical regional sensitivity.

Reliability of regional climate model trends

A necessary condition for a good probabilistic forecast is that the forecast system is shown to be reliable: forecast probabilities should equal observed probabilities verified over a large number of cases. As climate change trends are now emerging from the natural variability, we can apply this concept to climate predictions and compute the reliability of simulated local and regional temperature and precipitation trends (1950–2011) in a recent multi-model ensemble of climate model simulations prepared for the Intergovernmental Panel on Climate Change (IPCC) fifth assessment report (AR5). With only a single verification time, the verification is over the spatial dimension. The local temperature trends appear to be reliable. However, when the global mean climate response is factored out, the ensemble is overconfident: the observed trend is outside the range of modelled trends in many more regions than would be expected by the model estimate of natural variability and model spread. Precipitation trends are overconfident for all trend definitions. This implies that for near-term local climate forecasts the CMIP5 ensemble cannot simply be used as a reliable probabilistic forecast.

Conditioning model output statistics of regional climate model precipitation on circulation patterns

Dynamical downscaling of Global Climate Models (GCMs) through regional climate models (RCMs) potentially improves the usability of the output for hydrological impact studies. However, a further downscaling or interpolation of precipitation from RCMs is often needed to match the precipitation characteristics at the local scale. This study analysed three Model Output Statistics (MOS) techniques to adjust RCM precipitation; (1) a simple direct method (DM), (2) quantile-quantile mapping (QM) and (3) a distribution-based scaling (DBS) approach. The modelled precipitation was daily means from 16 RCMs driven by ERA40 reanalysis data over the 1961–2000 provided by the ENSEMBLES (ENSEMBLE-based Predictions of Climate Changes and their Impacts) project over a small catchment located in the Midlands, UK. All methods were conditioned on the entire time series, separate months and using an objective classification of Lamb's weather types. The performance of the MOS techniques were assessed regarding temporal and spatial characteristics of the precipitation fields, as well as modelled runoff using the HBV rainfall-runoff model. The results indicate that the DBS conditioned on classification patterns performed better than the other methods, however an ensemble approach in terms of both climate models and downscaling methods is recommended to account for uncertainties in the MOS methods.

Regional changes in wind energy potential over Europe using regional climate model ensemble projections

The impact of climate change on wind power generation potentials over Europe is investigated by considering ensemble projections from two regional climate models (RCMs) driven by a global climate model (GCM). Wind energy density and its interannual variability are estimated based on hourly near-surface wind speeds. Additionally, the possible impact of climatic changes on the energy output of a sample 2.5-MW turbine is discussed. GCM-driven RCM simulations capture the behavior and variability of current wind energy indices, even though some differences exist when compared with reanalysis-driven RCM simulations. Toward the end of the twenty-first century, projections show significant changes of energy density on annual average across Europe that are substantially stronger in seasonal terms. The emergence time of these changes varies from region to region and season to season, but some long-term trends are already statistically significant in the middle of the twenty-first century. Over northern and central Europe, the wind energy potential is projected to increase, particularly in winter and autumn. In contrast, energy potential over southern Europe may experience a decrease in all seasons except for the Aegean Sea. Changes for wind energy output follow the same patterns but are of smaller magnitude. The GCM/RCM model chains project a significant intensification of both interannual and intra-annual variability of energy density over parts of western and central Europe, thus imposing new challenges to a reliable pan-European energy supply in future decades.

Tackling regional climate change by leaf albedo bio-geoengineering

The likelihood that continuing greenhouse-gas emissions will lead to an unmanageable degree of climate change [1] has stimulated the search for planetary-scale technological solutions for reducing global warming [2] (“geoengineering”), typically characterized by the necessity for costly new infrastructures and industries [3]. We suggest that the existing global infrastructure associated with arable agriculture can help, given that crop plants exert an important influence over the climatic energy budget 4 and 5 because of differences in their albedo (solar reflectivity) compared to soils and to natural vegetation [6]. Specifically, we propose a “bio-geoengineering” approach to mitigate surface warming, in which crop varieties having specific leaf glossiness and/or canopy morphological traits are specifically chosen to maximize solar reflectivity. We quantify this by modifying the canopy albedo of vegetation in prescribed cropland areas in a global-climate model, and thereby estimate the near-term potential for bio-geoengineering to be a summertime cooling of more than 1°C throughout much of central North America and midlatitude Eurasia, equivalent to seasonally offsetting approximately one-fifth of regional warming due to doubling of atmospheric CO2[7]. Ultimately, genetic modification of plant leaf waxes or canopy structure could achieve greater temperature reductions, although better characterization of existing intraspecies variability is needed first.

Regional climate change mitigation with crops: context and assessment

The intention of this review is to place crop albedo biogeoengineering in the wider picture of climate manipulation. Crop biogeoengineering is considered within the context of the long-term modification of the land surface for agriculture over several thousand years. Biogeoengineering is also critiqued in relation to other geoengineering schemes in terms of mitigation power and adherence to social principles for geoengineering. Although its impact is small and regional, crop biogeoengineering could be a useful and inexpensive component of an ensemble of geoengineering schemes to provide temperature mitigation. The method should not detrimentally affect food security and there may even be positive impacts on crop productivity, although more laboratory and field research is required in this area to understand the underlying mechanisms.

Probabilistic parameterisation of the surface mass balance--elevation feedback in regional climate model simulations of the Greenland ice sheet

We present a new parameterisation that relates surface mass balance (SMB: the sum of surface accumulation and surface ablation) to changes in surface elevation of the Greenland ice sheet (GrIS) for the MAR (Modèle Atmosphérique Régional: Fettweis, 2007) regional climate model. The motivation is to dynamically adjust SMB as the GrIS evolves, allowing us to force ice sheet models with SMB simulated by MAR while incorporating the SMB–elevation feedback, without the substantial technical challenges of coupling ice sheet and climate models. This also allows us to assess the effect of elevation feedback uncertainty on the GrIS contribution to sea level, using multiple global climate and ice sheet models, without the need for additional, expensive MAR simulations. We estimate this relationship separately below and above the equilibrium line altitude (ELA, separating negative and positive SMB) and for regions north and south of 77� N, from a set of MAR simulations in which we alter the ice sheet surface elevation. These give four “SMB lapse rates”, gradients that relate SMB changes to elevation changes. We assess uncertainties within a Bayesian framework, estimating probability distributions for each gradient from which we present best estimates and credibility intervals (CI) that bound 95% of the probability. Below the ELA our gradient estimates are mostly positive, because SMB usually increases with elevation: 0.56 (95% CI: −0.22 to 1.33) kgm−3 a−1 for the north, and 1.91 (1.03 to 2.61) kgm−3 a−1 for the south. Above the ELA, the gradients are much smaller in magnitude: 0.09 (−0.03 to 0.23) kgm−3 a−1 in the north, and 0.07 (−0.07 to 0.59) kgm−3 a−1 in the south, because SMB can either increase or decrease in response to increased elevation. Our statistically founded approach allows us to make probabilistic assessments for the effect of elevation feedback uncertainty on sea level projections (Edwards et al., 2014).

Impacts of uplift of northern Tibetan Plateau and formation of Asian inland deserts on regional climate and environment

Based on the geological evidence that the northern Tibetan Plateau (NTP) had an uplift of a finite magnitude since the Miocene and the major Asian inland deserts formed in the early Pliocene, a regional climate model (RegCM4.1) with a horizontal resolution of 50 km was used to explore the effects of the NTP uplift and the related aridification of inland Asia on regional climate. We designed three numerical experiments including the control experiment representing the present-day condition, the high-mountain experiment representing the early Pliocene condition with uplifted NTP but absence of the Asian inland deserts, and the low-mountain experiment representing the mid-Miocene condition with reduced topography in the NTP (by as much as 2400 m) and also absence of the deserts. Our simulation results indicated that the NTP uplift caused significant reductions in annual precipitation in a broad region of inland Asia north of the Tibetan Plateau (TP) mainly due to the enhanced rain shadow effect of the mountains and changes in the regional circulations. However, four mountainous regions located in the uplift showed significant increases in precipitation, stretching from the Pamir Plateau in the west to the Qilian Mountains in the east. These mountainous areas also experienced different changes in the rainfall seasonality with the greatest increases occurring during the respective rainy seasons, predominantly resulted from the enhanced orographically forced upwind ascents. The appearance of the major deserts in the inland Asia further reduced precipitation in the region and led to increased dust emission and deposition fluxes, while the spatial patterns of dust deposition were also changed, not only in the regions of uplift-impacted topography, but also in the downwind regions. One major contribution from this study is the comparison of the simulation results with 11 existing geological records representing the moisture conditions from Miocene to Pliocene. The comparisons revealed good matches between the simulation results and the published geological records. Therefore, we conclude that the NTP uplift and the related formation of the major deserts played a controlling role in the evolution of regional climatic conditions in a broad region in inland Asia since the Miocene.

Robustness despite uncertainty: regional climate data reveal the dominant role of humans in explaining global extinctions of Late Quaternary megafauna

Debate over the late Quaternary megafaunal extinctions has focussed on whether human colonisation or climatic changes were more important drivers of extinction, with few extinctions being unambiguously attributable to either. Most analyses have been geographically or taxonomically restricted and the few quantitative global analyses have been limited by coarse temporal resolution or overly simplified climate reconstructions or proxies. We present a global analysis of the causes of these extinctions which uses high-resolution climate reconstructions and explicitly investigates the sensitivity of our results to uncertainty in the palaeological record. Our results show that human colonisation was the dominant driver of megafaunal extinction across the world but that climatic factors were also important. We identify the geographic regions where future research is likely to have the most impact, with our models reliably predicting extinctions across most of the world, with the notable exception of mainland Asia where we fail to explain the apparently low rate of extinction found in in the fossil record. Our results are highly robust to uncertainties in the palaeological record, and our main conclusions are unlikely to change qualitatively following minor improvements or changes in the dates of extinctions and human colonisation.