Mathieu Vrac

Precipitation downscaling under climate change: Recent developments to bridge the gap between dynamical models and the end user

Precipitation downscaling improves the coarse resolution and poor representation of precipitation in global climate models and helps end users to assess the likely hydrological impacts of climate change. This paper integrates perspectives from meteorologists, climatologists, statisticians, and hydrologists to identify generic end user (in particular, impact modeler) needs and to discuss downscaling capabilities and gaps. End users need a reliable representation of precipitation intensities and temporal and spatial variability, as well as physical consistency, independent of region and season. In addition to presenting dynamical downscaling, we review perfect prognosis statistical downscaling, model output statistics, and weather generators, focusing on recent developments to improve the representation of space-time variability. Furthermore, evaluation techniques to assess downscaling skill are presented. Downscaling adds considerable value to projections from global climate models. Remaining gaps are uncertainties arising from sparse data; representation of extreme summer precipitation, subdaily precipitation, and full precipitation fields on fine scales; capturing changes in small-scale processes and their feedback on large scales; and errors inherited from the driving global climate model.

Multivariate—Intervariable, Spatial, and Temporal—Bias Correction*

AbstractStatistical methods to bias correct global or regional climate model output are now common to get data closer to observations in distribution. However, most bias correction (BC) methods work for one variable and one location at a time and basically reproduce the temporal structure of the models. The intervariable, spatial, and temporal dependencies of the corrected data are usually poor compared to observations. Here, the authors propose a novel method for multivariate BC. The empirical copula–bias correction (EC–BC) combines a one-dimensional BC with a shuffling technique that restores an empirical multidimensional copula. Several BC methods are investigated and compared to high-resolution reference data over the French Mediterranean basin: notably, (i) a 1D BC method applied independently to precipitation and temperature fields, (ii) a recent conditional correction approach developed for producing correct two-dimensional intervariable structures, and (iii) the EC–BC method.Assessments are realiz...

Are regional climate models relevant for crop yield prediction in West Africa

This study assesses the accuracy of state-of-the-art regional climate models for agriculture applications in West Africa. A set of nine regional configurations with eight regional models from the ENSEMBLES project is evaluated. Although they are all based on similar large-scale conditions, the performances of regional models in reproducing the most crucial variables for crop production are extremely variable. This therefore leads to a large dispersion in crop yield prediction when using regional models in a climate/crop modelling system. This dispersion comes from the different physics in each regional model and also the choice of parametrizations for a single regional model. Indeed, two configurations of the same regional model are sometimes more distinct than two different regional models. Promising results are obtained when applying a bias correction technique to climate model outputs. Simulated yields with bias corrected climate variables show much more realistic means and standard deviations. However, such a bias correction technique is not able to improve the reproduction of the year-to-year variations of simulated yields. This study confirms the importance of the multi-model approach for quantifying uncertainties for impact studies and also stresses the benefits of combining both regional and statistical downscaling techniques. Finally, it indicates the urgent need to address the main uncertainties in atmospheric processes controlling the monsoon system and to contribute to the evaluation and improvement of climate and weather forecasting models in that respect.

Regional climate downscaling with prior statistical correction of the global climate forcing

A novel climate downscaling methodology that attempts to correct climate simulation biases is proposed. By combining an advanced statistical bias correction method with a dynamical downscaling it constitutes a hybrid technique that yields nearly unbiased, high-resolution, physically consistent, three-dimensional fields that can be used for climate impact studies. The method is based on a prior statistical distribution correction of large-scale global climate model (GCM) 3-dimensional output fields to be taken as boundary forcing of a dynamical regional climate model (RCM). GCM fields are corrected using meteorological reanalyses. We evaluate this methodology over a decadal experiment. The improvement in terms of spatial and temporal variability is discussed against observations for a past period. The biases of the downscaled fields are much lower using this hybrid technique, up to a factor 4 for the mean temperature bias compared to the dynamical downscaling alone without prior bias correction. Precipitation biases are subsequently improved hence offering optimistic perspectives for climate impact studies.

Stochastic Model Output Statistics for Bias Correcting and Downscaling Precipitation Including Extremes

Precipitation is highly variable in space and time; hence, rain gauge time series generally exhibit additional random small-scale variability compared to area averages. Therefore, differences between daily precipitation statistics simulated by climate models and gauge observations are generally not only caused by model biases, but also by the corresponding scale gap. Classical bias correction methods, in general, cannot bridge this gap; they do not account for small-scale random variability and may produce artifacts. Here, stochastic model output statistics is proposed as a bias correction framework to explicitly account for random small-scale variability. Daily precipitation simulated by a regional climate model (RCM) is employed to predict the probability distribution of local precipitation. The pairwise correspondence between predictor and predictand required for calibration is ensured by driving the RCM with perfect boundary conditions. Wet day probabilities are described by a logistic regression, and precipitation intensities are described by a mixture model consisting of a gamma distribution for moderate precipitation and a generalized Pareto distribution for extremes. The dependence of the model parameters on simulated precipitation is modeled by a vector generalized linear model. The proposed model effectively corrects systematic biases and correctly represents local-scale random variability for most gauges. Additionally, a simplified model is considered that disregards the separate tail model. This computationally efficient model proves to be a feasible alternative for precipitation up to moderately extreme intensities. The approach sets a new framework for bias correction that combines the advantages of weather generators and RCMs.

Clustering of Maxima: Spatial Dependencies among Heavy Rainfall in France

AbstractOne of the main objectives of statistical climatology is to extract relevant information hidden in complex spatial–temporal climatological datasets. To identify spatial patterns, most well-known statistical techniques are based on the concept of intra- and intercluster variances (like the k-means algorithm or EOFs). As analyzing quantitative extremes like heavy rainfall has become more and more prevalent for climatologists and hydrologists during these last decades, finding spatial patterns with methods based on deviations from the mean (i.e., variances) may not be the most appropriate strategy in this context of studying such extremes. For practitioners, simple and fast clustering tools tailored for extremes have been lacking. A possible avenue to bridging this methodological gap resides in taking advantage of multivariate extreme value theory, a well-developed research field in probability, and to adapt it to the context of spatial clustering. In this paper, a novel algorithm based on this plan ...

Intercomparison of statistical and dynamical downscaling models under the EURO- and MED-CORDEX initiative framework: present climate evaluations

Given the coarse spatial resolution of General Circulation Models, finer scale projections of variables affected by local-scale processes such as precipitation are often needed to drive impacts models, for example in hydrology or ecology among other fields. This need for high-resolution data leads to apply projection techniques called downscaling. Downscaling can be performed according to two approaches: dynamical and statistical models. The latter approach is constituted by various statistical families conceptually different. If several studies have made some intercomparisons of existing downscaling models, none of them included all those families and approaches in a manner that all the models are equally considered. To this end, the present study conducts an intercomparison exercise under the EURO- and MED-CORDEX initiative hindcast framework. Six Statistical Downscaling Models (SDMs) and five Regional Climate Models (RCMs) are compared in terms of precipitation outputs. The downscaled simulations are driven by the ERAinterim reanalyses over the 1989–2008 period over a common area at 0.44° of resolution. The 11 models are evaluated according to four aspects of the precipitation: occurrence, intensity, as well as spatial and temporal properties. For each aspect, one or several indicators are computed to discriminate the models. The results indicate that marginal properties of rain occurrence and intensity are better modelled by stochastic and resampling-based SDMs, while spatial and temporal variability are better modelled by RCMs and resampling-based SDM. These general conclusions have to be considered with caution because they rely on the chosen indicators and could change when considering other specific criteria. The indicators suit specific purpose and therefore the model evaluation results depend on the end-users point of view and how they intend to use with model outputs. Nevertheless, building on previous intercomparison exercises, this study provides a consistent intercomparison framework, including both SDMs and RCMs, which is designed to be flexible, i.e., other models and indicators can easily be added. More generally, this framework provides a tool to select the downscaling model to be used according to the statistical properties of the local-scale climate data to drive properly specific impact models.

Comparison of GCM‐ and RCM‐simulated precipitation following stochastic postprocessing

In order to assess to what extent regional climate models (RCMs) yield better representations of climatic states than general circulation models (GCMs), the output of each is usually directly compared with observations. RCM output is often bias corrected, and in some cases correction methods can also be applied to GCMs. This leads to the question of whether bias-corrected RCMs perform better than bias-corrected GCMs. Here the first results from such a comparison are presented, followed by discussion of the value added by RCMs in this setup. Stochastic postprocessing, based on Model Output Statistics (MOS), is used to estimate daily precipitation at 465 stations across the United Kingdom between 1961 and 2000 using simulated precipitation from two RCMs (RACMO2 and CCLM) and, for the first time, a GCM (ECHAM5) as predictors. The large-scale weather states in each simulation are forced toward observations. The MOS method uses logistic regression to model precipitation occurrence and a Gamma distribution for the wet day distribution, and is cross validated based on Brier and quantile skill scores. A major outcome of the study is that the corrected GCM-simulated precipitation yields consistently higher validation scores than the corrected RCM-simulated precipitation. This seems to suggest that, in a setup with postprocessing, there is no clear added value by RCMs with respect to downscaling individual weather states. However, due to the different ways of controlling the atmospheric circulation in the RCM and the GCM simulations, such a strong conclusion cannot be drawn. Yet the study demonstrates how challenging it is to demonstrate the value added by RCMs in this setup.

Strengthening the link between climate, hydrological and species distribution modeling to assess the impacts of climate change on freshwater biodiversity.

To understand the resilience of aquatic ecosystems to environmental change, it is important to determine how multiple, related environmental factors, such as near-surface air temperature and river flow, will change during the next century. This study develops a novel methodology that combines statistical downscaling and fish species distribution modeling, to enhance the understanding of how global climate changes (modeled by global climate models at coarse-resolution) may affect local riverine fish diversity. The novelty of this work is the downscaling framework developed to provide suitable future projections of fish habitat descriptors, focusing particularly on the hydrology which has been rarely considered in previous studies. The proposed modeling framework was developed and tested in a major European system, the Adour-Garonne river basin (SW France, 116,000 km(2)), which covers distinct hydrological and thermal regions from the Pyrenees to the Atlantic coast. The simulations suggest that, by 2100, the mean annual stream flow is projected to decrease by approximately 15% and temperature to increase by approximately 1.2 °C, on average. As consequence, the majority of cool- and warm-water fish species is projected to expand their geographical range within the basin while the few cold-water species will experience a reduction in their distribution. The limitations and potential benefits of the proposed modeling approach are discussed.

Quantifying Differences in Circulation Patterns Based on Probabilistic Models: IPCC AR4 Multimodel Comparison for the North Atlantic*

The comparison of circulation patterns (CPs) obtained from reanalysis data to those from general circulation model (GCM) simulations is a frequent task for model validation, downscaling of GCM simulations, or other climate change‐related studies. Here, the authors suggest a set of measures to quantify the differences between CPs. A combination of clustering using Gaussian mixture models with a set of related difference measures allows for taking cluster size and shape information into account and thus provides more information than the Euclidean distances between cluster centroids. The characteristics of the various distance measures are illustrated with a simple simulated example. Subsequently, a five-component Gaussian mixture to define circulation patterns for the North Atlantic region from reanalysis data and GCM simulations is used. CPs are obtained independently for the NCEP‐NCAR reanalysis and the 40-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40), as well as for twentieth-century simulations from 14 GCMs of the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) database. After discussing the difference of CPs based on spherical and nonspherical clusters for the reanalysis datasets, the authors give a detailed evaluation of the cluster configuration for two GCMs relative to NCEP‐NCAR. Finally, as an illustration, the capability of reproducing the NCEP‐NCAR probability density function(pdf) definingtheGreenlandanticycloneCPisevaluated forall 14GCMs, consideringthat thesizeand shape of the underlying pdfs complement the commonly used Euclidean distance of CPs’ mean values.