S. Herrera

How well do CMIP5 Earth System Models simulate present climate conditions in Europe and Africa

The present study assesses the ability of seven Earth System Models (ESMs) from the Coupled Model Intercomparison Project Phase 5 to reproduce present climate conditions in Europe and Africa. This is done from a downscaling perspective, taking into account the requirements of both statistical and dynamical approaches. ECMWF’s ERA-Interim reanalysis is used as reference for an evaluation of circulation, temperature and humidity variables on daily timescale, which is based on distributional similarity scores. To additionally obtain an estimate of reanalysis uncertainty, ERA-Interim’s deviation from the Japanese Meteorological Agency JRA-25 reanalysis is calculated. Areas with considerable differences between both reanalyses do not allow for a proper assessment, since ESM performance is sensitive to the choice of reanalysis. For use in statistical downscaling studies, ESM performance is computed on the grid-box scale and mapped over a large spatial domain covering Europe and Africa, additionally highlighting those regions where significant distributional differences remain even for the centered/zero-mean time series. For use in dynamical downscaling studies, performance is specifically assessed along the lateral boundaries of the three CORDEX domains defined for Europe, the Mediterranean Basin and Africa.

Reassessing Statistical Downscaling Techniques for Their Robust Application under Climate Change Conditions

AbstractThe performance of statistical downscaling (SD) techniques is critically reassessed with respect to their robust applicability in climate change studies. To this end, in addition to standard accuracy measures and distributional similarity scores, the authors estimate the robustness of the methods under warming climate conditions working with anomalous warm historical periods. This validation framework is applied to intercompare the performances of 12 different SD methods (from the analog, weather typing, and regression families) for downscaling minimum and maximum temperatures in Spain. First, a calibration of these methods is performed in terms of both geographical domains and predictor sets; the results are highly dependent on the latter, with optimum predictor sets including near-surface temperature data (in particular 2-m temperature), which appropriately discriminate cold episodes related to temperature inversion in the lower troposphere.Although regression methods perform best in terms of co...

On the Use of Reanalysis Data for Downscaling

AbstractIn this study, a worldwide overview on the expected sensitivity of downscaling studies to reanalysis choice is provided. To this end, the similarity of middle-tropospheric variables—which are important for the development of both dynamical and statistical downscaling schemes—from 40-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40) and NCEP–NCAR reanalysis data on a daily time scale is assessed. For estimating the distributional similarity, two comparable scores are used: the two-sample Kolmogorov–Smirnov statistic and the probability density function (PDF) score. In addition, the similarity of the day-to-day sequences is evaluated with the Pearson correlation coefficient. As the most important results demonstrated, the PDF score is found to be inappropriate if the underlying data follow a mixed distribution. By providing global similarity maps for each variable under study, regions where reanalysis data should not assumed to be “perfect” are detected. In contrast ...

Forest fire danger projections in the Mediterranean using ENSEMBLES regional climate change scenarios

We present future fire danger scenarios for the countries bordering the Mediterranean areas of Europe and north Africa building on a multi-model ensemble of state-of-the-art regional climate projections from the EU-funded project ENSEMBLES. Fire danger is estimated using the Canadian Forest Fire Weather Index (FWI) System and a related set of indices. To overcome some of the limitations of ENSEMBLES data for their application on the FWI System—recently highlighted in a previous study by Herrera et al. (Clim Chang 118:827–840, 2013)—we used an optimal proxy variable combination. A robust assessment of future fire danger projections is undertaken by disentangling the climate change signal from the uncertainty derived from the multi-model ensemble, unveiling a positive signal of fire danger potential over large areas of the Mediterranean. The increase in the fire danger signal is accentuated towards the latest part of the transient period, thus pointing to an elevated fire potential in the region with time. The fire-climate links under present and future conditions are further discussed building upon observed climate data and burned area records along a representative climatic gradient within the study region.

Assessing the skill of precipitation and temperature seasonal forecasts in Spain: windows of opportunity related to ENSO events

The skill of state-of-the-art operational seasonal forecast models in extratropical latitudes is assessed using a multimodel ensemble from the Development of a European Multimodel Ensemble System for Seasonalto-Interannual Prediction (DEMETER) project. In particular, probabilistic forecasts of surface precipitation and maximum temperature in Spain are analyzed using a high-resolution observation gridded dataset (Spain02). To this aim, a simple statistical test based on the observed and predicted tercile anomalies is used. First, the whole period 1960–2000 is considered and it is shown that the only significant skill is found for dry events in autumn. Then, the influence of ENSO events as a potential source of conditional predictability is

Testing MOS precipitation downscaling for ENSEMBLES regional climate models over Spain

[1] Model Output Statistics (MOS) has been recently proposed as an alternative to the standard perfect prognosis statistical downscaling approach for Regional Climate Model (RCM) outputs. In this case, the model output for the variable of interest (e.g. precipitation) is directly downscaled using observations. In this paper we test the performance of a MOS implementation of the popular analog methodology (referred to as MOS analog) applied to downscale daily precipitation outputs over Spain. To this aim, we consider the state‐of‐the‐art ERA40‐driven RCMs provided by the EU‐funded ENSEMBLES project and the Spain02 gridded observations data set, using the common period 1961–2000. The MOS analog method improves the representation of the mean regimes, the annual cycle, the frequency and the extremes of precipitation for all RCMs, regardless of the region and the model reliability (including relatively low‐performing models), while preserving the daily accuracy. The good performance of the method in this complex climatic region suggests its potential transferability to other regions. Furthermore, in order to test the robustness of the method in changing climate conditions, a cross‐validation in driest or wettest years was performed. The method improves the RCM results in both cases, especially in the former.

Robust projections of Fire Weather Index in the Mediterranean using statistical downscaling

The effect of climate change on wildfires constitutes a serious concern in fire-prone regions with complex fire behavior such as the Mediterranean. The coarse resolution of future climate projections produced by General Circulation Models (GCMs) prevents their direct use in local climate change studies. Statistical downscaling techniques bridge this gap using empirical models that link the synoptic-scale variables from GCMs to the local variables of interest (using e.g. data from meteorological stations). In this paper, we investigate the application of statistical downscaling methods in the context of wildfire research, focusing in the Canadian Fire Weather Index (FWI), one of the most popular fire danger indices. We target on the Iberian Peninsula and Greece and use historical observations of the FWI meteorological drivers (temperature, humidity, wind and precipitation) in several local stations. In particular, we analyze the performance of the analog method, which is a convenient first choice for this problem since it guarantees physical and spatial consistency of the downscaled variables, regardless of their different statistical properties. First we validate the method in perfect model conditions using ERA-Interim reanalysis data. Overall, not all variables are downscaled with the same accuracy, with the poorest results (with spatially averaged daily correlations below 0.5) obtained for wind, followed by precipitation. Consequently, those FWI components mostly relying on those parameters exhibit the poorest results. However, those deficiencies are compensated in the resulting FWI values due to the overall high performance of temperature and relative humidity. Then, we check the suitability of the method to downscale control projections (20C3M scenario) from a single GCM (the ECHAM5 model) and compute the downscaled future fire danger projections for the transient A1B scenario. In order to detect problems due to non-stationarities related to climate change, we compare the results with those obtained with a Regional Climate Model (RCM) driven by the same GCM. Although both statistical and dynamical projections exhibit a similar pattern of risk increment in the first half of the 21st century, they diverge during the second half of the century. As a conclusion, we advocate caution in the use of projections for this last period, regardless of the regionalization technique applied.

Large biases and inconsistent climate change signals in ENSEMBLES regional projections

In this paper we analyze some caveats found in the state-of-the-art ENSEMBLES regional projections dataset focusing on precipitation over Spain, and highlight the need of a task-oriented validation of the GCM-driven control runs. In particular, we compare the performance of the GCM-driven control runs (20C3M scenario) with the ERA40-driven ones (“perfect” boundary conditions) in a common period (1961–2000). Large deviations between the results indicate a large uncertainty/bias for the particular RCM-GCM combinations and, hence, a small confidence for the corresponding transient simulations due to the potential nonlinear amplification of biases. Specifically, we found large biases for some RCM-GCM combinations attributable to RCM in-house problems with the particular GCM coupling. These biases are shown to distort the corresponding climate change signal, or “delta”, in the last decades of the 21st century, considering the A1B scenario. Moreover, we analyze how to best combine the available RCMs to obtain more reliable projections.

Precipitation variability and trends in Ghana: An intercomparison of observational and reanalysis products

Inter-annual variability and trends of annual/seasonal precipitation totals in Ghana are analyzed considering different gridded observational (gauge- and/or satellite-based) and reanalysis products. A quality-controlled dataset formed by fourteen gauges from the Ghana Meteorological Agency (GMet) is used as reference for the period 1961–2010. Firstly, a good agreement is found between GMet and all the observational products in terms of variability, with better results for the gauge-based products—correlations in the range of 0.7–1.0 and nearly null biases—than for the satellite-gauge merged and satellite-derived products. In contrast, reanalyses exhibit a very poor performance, with correlations below 0.4 and large biases in most of the cases. Secondly, a Mann-Kendall trend analysis is carried out. In most cases, GMet data reveal the existence of predominant decreasing (increasing) trends for the first (second) half of the period of study, 1961–1985 (1986–2010). Again, observational products are shown to reproduce well the observed trends—with worst results for purely satellite-derived data—whereas reanalyses lead in general to unrealistic stronger than observed trends, with contradictory results (opposite signs for different reanalyses) in some cases. Similar inconsistencies are also found when analyzing trends of extreme precipitation indicators. Therefore, this study provides a warning concerning the use of reanalysis data as pseudo-observations in Ghana.

Bias correction and downscaling of future RCM precipitation projections using a MOS‐Analog technique

In this study we assess the suitability of a recently introduced analog-based Model Output Statistics (MOS) downscaling method (referred to as MOS-Analog) for climate change studies and compare the results with a quantile mapping bias correction method. To this aim, we focus on Spain and consider daily precipitation output from an ensemble of Regional Climate Models provided by the ENSEMBLES project. The reanalysis-driven Regional Climate Model (RCM) data provide the historical data (with day-to-day correspondence with observations induced by the forcing boundary conditions) to conduct the analog search of the control (20C3M) and future (A1B) global climate model (GCM)-driven RCM values. First, we show that the MOS-Analog method outperforms the raw RCM output in the control 20C3M scenario (period 1971–2000) for all considered regions and precipitation indices, although for the worst-performing models the method is less effective. Second, we show that the MOS-Analog method broadly preserves the original RCM climate change signal for different future periods (2011–2040, 2041–2070, 2071–2100), except for those indices related to extreme precipitation. This could be explained by the limitation of the analog method to extrapolate unobserved precipitation records. These results suggest that the MOS-Analog is a spatially consistent alternative to standard bias correction methods, although the limitation for extreme values should be taken with caution in cases where this aspect is relevant for the problem.