The effect of univariate bias adjustment on multivariate hazard estimates
- 1Institute for Atmospheric and Climate Science, ETH Zurich, Universitaetstrasse 16, 8092 Zurich, Switzerland
- 2Climate and Environmental Physics, University of Bern, Sidlerstrasse 5, 3012 Bern, Switzerland
- 3Oeschger Centre for Climate Change Research, University of Bern, 3012 Bern, Switzerland
- 4Potsdam Institute for Climate Impact Research (PIK), Member of the Leibniz Association, P.O. Box 60 12 03, 14412 Potsdam, Germany
Correspondence: Jakob Zscheischler (email@example.com)
Bias adjustment is often a necessity in estimating climate impacts because impact models usually rely on unbiased climate information, a requirement that climate model outputs rarely fulfil. Most currently used statistical bias-adjustment methods adjust each climate variable separately, even though impacts usually depend on multiple potentially dependent variables. Human heat stress, for instance, depends on temperature and relative humidity, two variables that are often strongly correlated. Whether univariate bias-adjustment methods effectively improve estimates of impacts that depend on multiple drivers is largely unknown, and the lack of long-term impact data prevents a direct comparison between model outputs and observations for many climate-related impacts. Here we use two hazard indicators, heat stress and a simple fire risk indicator, as proxies for more sophisticated impact models. We show that univariate bias-adjustment methods such as univariate quantile mapping often cannot effectively reduce biases in multivariate hazard estimates. In some cases, it even increases biases. These cases typically occur (i) when hazards depend equally strongly on more than one climatic driver, (ii) when models exhibit biases in the dependence structure of drivers and (iii) when univariate biases are relatively small. Using a perfect model approach, we further quantify the uncertainty in bias-adjusted hazard indicators due to internal variability and show how imperfect bias adjustment can amplify this uncertainty. Both issues can be addressed successfully with a statistical bias adjustment that corrects the multivariate dependence structure in addition to the marginal distributions of the climate drivers. Our results suggest that currently many modeled climate impacts are associated with uncertainties related to the choice of bias adjustment. We conclude that in cases where impacts depend on multiple dependent climate variables these uncertainties can be reduced using statistical bias-adjustment approaches that correct the variables' multivariate dependence structure.
With ongoing climate change, climate impact modeling has become an important pillar of climate research, informing decision makers and risk managers in many sectors that are affected by climate variability. Impact models such as hydrological models, crop models and epidemiological models usually rely on absolute thresholds in their driving climate variables such as temperature, precipitation, wind speed and humidity (Ruane et al., 2017; Winsemius et al., 2015); and thus require unbiased climate information as input. However, biases continue to persist in global (Flato et al., 2013; Wang et al., 2014) and regional climate models (Christensen et al., 2008; Kotlarski et al., 2014), rendering bias adjustment an undesired but often unavoidable data processing step for climate impact modeling (Piani et al., 2010).
The need for easily available information on climate impacts has led to a sometimes overly uncritical use of bias adjustment in large impact modeling projects, even though its current usage might, in some cases, result in ill-informed adaptation decisions (Maraun et al., 2017). It is often assumed that biases are stationary, and hence that bias adjustment can be developed in the current climate and applied in a warmer world. It is not clear, however, for which variable and to what extent this assumption is met. Moreover, while uncertainties related to the choice of climate model and impact model as well as future socioeconomic scenarios are frequently carried through the modeling chain (Frieler et al., 2017; Wilby and Dessai, 2010), uncertainties associated with the choice of the applied bias-adjustment method are rarely reported (Ehret et al., 2012), though there are exceptions (Addor and Fischer, 2015; Bosshard et al., 2013; Chen et al., 2011).
Assessing the appropriateness of bias adjustment for impact modeling is not an easy task (Papadimitriou et al., 2017). While the quality of bias adjustment on climate variables can be evaluated against observed quantities in the climate domain, the implications for modeled impacts are much harder to determine. On the one hand, limited length and a lack of homogeneity in observed climate impacts render an evaluation of modeled impacts very challenging (Cramer et al., 2014). On the other hand, many impacts rely on the complex interaction of multiple climate variables across time and space (Zscheischler et al., 2018), preventing an evaluation in the climate domain. For instance, temperature and precipitation variability affect crop yields (Semenov and Porter, 1995; Zscheischler et al., 2017) and other ecosystem services such as net carbon uptake (Humphrey et al., 2018). Flood occurrence, its intensity and associated damages depend on the temporal and spatial characteristics of precipitation, soil moisture, river flow and surge (Vorogushyn et al., 2018). Drought-related impacts depend on precipitation, evapotranspiration, and temperature; their spatial distribution; and their interaction with human activities (Van Loon et al., 2016). Fire occurrence, and strength, not only requires available fuel and an ignition source, but is also strongly dependent on relative humidity and temperature (Brando et al., 2014). Infrastructural damage is particularly large when strong winds and extreme precipitation occur jointly (Martius et al., 2016). Finally, climate-related human-health impacts are linked to a suite of climatic drivers (McMichael et al., 2006). It is currently largely unknown how widely used statistical bias-adjustment approaches affect modeled impacts that depend on multiple drivers such as the ones mentioned above.
Bias-adjustment methods are often designed to correct one variable at a time (Hempel et al., 2013; Teutschbein and Seibert, 2012). A possible bias in the dependence structure between climate variables is therefore not adjusted, even though climate models may not capture dependencies between climate drivers very well, for instance the correlation between temperature and precipitation in summer (Zscheischler and Seneviratne, 2017). Even worse, correcting a bias in marginal distributions may modify the dependence between variables. In recognition of these issues, several multivariate bias-adjustment methods have been suggested in the recent past (Cannon, 2016; Li et al., 2014; Mehrotra and Sharma, 2016; Piani and Haerter, 2012; Vrac, 2018; Vrac and Friederichs, 2015). In addition to adjusting the marginal distributions, these methods adjust, to a certain extent, the dependence structure between multiple variables. Some studies have recently suggested that multivariate approaches do not lead to a substantial improvement for certain specific regional impacts (Casanueva et al., 2018; Räty et al., 2018; Yang et al., 2015). However, this does not imply that multivariate bias adjustment is not necessary in any case. Many multivariate methods have not been systematically evaluated against modeled impacts and have been rarely used in larger-scale impact modeling frameworks, though there are notable exceptions. For instance, Cannon (2018) presents a highly flexible multivariate bias-adjustment method and demonstrates its effectiveness on a modeled multivariate hazard indicator. By adjusting biases and dependencies between temperature, precipitation, relative humidity and wind speed, the method substantially reduces biases in the five-dimensional Fire Weather Index (FWI), outperforming univariate bias-adjustment approaches.
In this paper we study the limitations of univariate bias-adjustment methods for multivariate impacts. We focus on global-scale climate outputs from global circulation models (GCMs), which are frequently used for global impact studies (Frieler et al., 2017; Ruane et al., 2017; Winsemius et al., 2015). However, our analyses are of a rather conceptual nature and thus also apply to other climate model outputs, for instance from regional climate models. First, we investigate whether traditional univariate bias-adjustment methods generally reduce biases in impacts that depend on multiple dependent drivers. Bias adjustment may also amplify uncertainties inherent to the observations (Chen et al., 2011). Hence, using a model environment with multiple models and multiple runs for a single model, we further estimate how uncertainties related to internal variability may be propagated and amplified through incomplete bias adjustment. We use multivariate hazard indicators as proxies for actual impacts to evaluate the effect of bias adjustment on impacts and to overcome the challenge of missing impact data. It can be assumed that “real” impacts are in many cases more complex and also depend on more driving variables. Therefore, our analysis provides a rather conservative estimate of the potential effects of univariate bias adjustment on modeled impacts.
Observations. We use daily temperature and relative humidity from the ERA-Interim reanalysis (Dee et al., 2011) as main reference dataset. We further use daily temperature and relative humidity from the observational dataset EWEMBI (Lange, 2016, 2018), which has been used in the Inter-Sectoral Impact Model Intercomparison Project phase 2b (ISIMIP2b; Frieler et al., 2017).
Model simulations. We use daily model output from the historical runs of the Coupled Model Intercomparison Project Phase 5 (CMIP5; Taylor et al., 2012). A main objective of the CMIP5 model experiments is to improve our understanding of the climate system, and to provide estimates of future climate change. All 29 model simulations used in this study are listed in Table A1.
All data have been bilinearly interpolated to a 2.5∘ by 2.5∘ regular latitude–longitude grid prior to analysis. The studied time period is 1981–1995. Our analysis is based on daily values during the hottest month of the year at each grid point, which results in about 450 samples. We focus on the hottest month to avoid dealing with seasonality and because, arguably, fire risk and heat stress are most relevant during this time period. The hottest month was identified based on the climatology of ERA-Interim monthly temperature data. Although the effectiveness of bias adjustment is typically evaluated outside the calibration period (Maraun, 2013), here we focus our analysis completely on the selected 15-year period. This allows for a separate assessment of the effect of bias adjustment on univariate versus multivariate impacts independently from other effects such as cross-validation error. A good performance in the calibration period is a necessary requirement for an effective bias-adjustment approach. Furthermore, cross-validation might not help to diagnose whether a bias-adjustment approach is effective (Maraun and Widmann, 2018; Switanek et al., 2017).
Hazard indicators. We use the Wet Bulb Globe Temperature (WBGT; Dunne et al., 2013) as an indicator for heat stress and the Chandler Burning Index (CBI; Chandler et al., 1983) as an indicator for fire risk. Both are relatively simple indicators that can be computed solely from daily temperature and relative humidity. WBGT and its variants have been extensively used to assess projections of heat stress under climate change (Li et al., 2017; Pal and Eltahir, 2015; Zhao et al., 2015). CBI is one of many fire risk indicators (Carlson and Burgan, 2003; Lee, 1980), mainly chosen here for its simplicity. The hazard intensities of WBGT and CBI vary along different gradients in the temperature–humidity domain (Fischer et al., 2013; Zscheischler et al., 2018). WBGT increases with hotter and more humid conditions and is equally dependent on temperature and humidity (Fig. 1). In contrast, CBI increases with hotter and drier conditions and its variability is mostly driven by humidity and much less by temperature. Using two hazard indicators that depend differently on the same climatic drivers enables us to study how the relationship between hazard direction and driver distribution changes the way bias adjustment affects modeled hazards (Zscheischler et al., 2018). While certainly more sophisticated indicators exist, both for fire risk and heat stress (Bröde et al., 2013), our goal here is rather to provide an illustrative example of the issues associated with bias adjustment and impact modeling than to provide the most reliable hazard projections. We estimate WBGT following the approach outlined in the supplement of Dunne et al. (2013). CBI can be computed as
where RH is relative humidity in % and T is temperature in degrees Celsius.
Bias adjustment. We employ two different bias-adjustment methods, the widely used univariate empirical quantile mapping (UQM; Casanueva et al., 2018; Maraun, 2013; Panofsky and Brier, 1968) and the multivariate bias adjustment in n dimensions (MBCn) developed by Cannon (2018). UQM applies separate corrections to a fixed number of quantiles to adjust a modeled empirical cumulative distribution towards the observed empirical cumulative distribution. Hence, if Xo and Xm are the observed and modeled values, respectively, then
where Fm is the empirical cumulative distribution function of Xm and is the inverse empirical distribution function (or quantile function) corresponding to Xo. Values between the predefined quantiles are approximated using linear interpolation. We apply UQM with the R package qmap (Gudmundsson, 2014) using 100 quantiles.
MBCn is a bias-adjustment method in which both the marginal distribution of each individual variable and the multivariate dependence structure are corrected at the same time (Cannon, 2018). This is achieved by an iterative approach, which first applies a random rotation R[j] to the multivariate observed and modeled data distribution
Subsequently, quantile mapping (Eq. 2) is applied to the rotated data with as a reference, yielding . Then the inverse rotation is applied
The observed data are carried forward to the next iteration unchanged . These steps are repeated until the modeled data distribution has converged to the observed distribution. One may interpret the algorithm such that the random rotations allow an information exchange between the different dimensions.
Finally, we use the bias-adjustment method (Lange, 2017) applied in ISIMIP2b (Frieler et al., 2017). This method adjusts relative humidity with parametric quantile mapping using beta distributions to model simulated and observed daily values (Frieler et al., 2017; Lange, 2018), and temperature with an additive correction of monthly mean temperature and a multiplicative correction of daily anomalies from the monthly mean temperature (Hempel et al., 2013).
We adjust daily values of temperature and relative humidity in CMIP5 during the hottest month at each grid point and evaluate the change in bias for WBGT and CBI.
Perfect model approach. Internal variability can lead to uncertainties in the bias adjustment, which may be amplified through an inadequate choice of bias adjustment. Since fully coupled ocean–atmosphere models provide different realizations of unforced internal variability, observations and models as well as different simulations of the same model are not expected to agree on a year-by-year basis. Even when using a time period of 2–3 decades, this constitutes a substantial source of uncertainty (Addor and Fischer, 2015). Using the multimodel environment of CMIP5, we study whether UQM increases uncertainties related to internal variability for multivariate hazard indicators. We conduct a perfect model approach to separate uncertainties associated with internal variability and the choice of bias adjustment (Elía et al., 2002; Griffies and Bryan, 1997; Hawkins et al., 2011). To this end we use the five available initial-condition members of the model CanESM to estimate the influence of internal variability as a source of uncertainty. First, we bias-adjust CanESM runs against each other to estimate the range of the uncertainty due to internal variability (“noise”). Bias adjusting all other model runs against all five CanESM runs then provides an estimate for the full uncertainty range (“full range”). Comparing the range of the noise with the full range, we estimate whether univariate bias adjustment amplifies uncertainties related to internal variability.
Statistical bias adjustment that is separately applied to each marginal distribution of a multivariate distribution such as UQM ensures that the bias-adjusted modeled marginal distributions are well aligned with the observed marginal distributions. We illustrate the effect of bias adjustment for one selected grid point in central Africa. UQM “squeezes” and “stretches” a modeled multivariate distribution along the marginal axes to match the observations of the marginals (Fig. 1c). For hazards that are a function of multiple drivers and that vary along a diagonal gradient such as WBGT, UQM may not be able to reduce biases. In fact, for many percentiles, it can even increase biases as illustrated in Fig. 2c (blue dots). CBI seems to be less affected, likely because its variability is mostly driven by a single variable, namely relative humidity (Figs. 1–2).
Bias adjustment can also affect the relationship between a multivariate hazard indicator and its individual contributing climate variables. For instance, looking at the raw model output of IPSL-CM5A-MR in Fig. 1b, one would infer that high heat stress (WBGT values exceeding their 95th percentile, highlighted in red) is reached at rather low relative humidity. Instead, for ERA-Interim data for the same grid point, high heat stress is associated with high relative humidity (Fig. 1a). UQM improves this mismatch to some extent but not entirely (Fig. 1c). Using MBCn for bias adjustment we would infer the correct contribution of individual variables. This is an important aspect, which needs to be taken into account when applying bias adjustment because it may lead to incorrect conclusions about which climate drivers are most relevant for extreme hazards and impacts. Consequently, we might also focus on the wrong aspects to improve in numerical climate models.
In the following, we quantify whether the case illustrated in Figs. 1 and 2 is representative across the globe. We find that an increase in root mean squared error (RMSE) between the cumulative distribution functions of ERA-Interim and model outputs before and after bias adjustment is an exception (Fig. 3a). However, even though biases decrease at most grid points when applying UQM, the reduction in bias is less than 50 % at 15±6 % (mean ± one standard deviation across models) of all grid points (Fig. 3b). Reducing biases for extreme percentiles (≥q90) may be even more challenging. For instance, for the 90th (95th) percentile of WBGT, UQM applied on temperature and relative humidity results in WBGT estimates that have larger biases than the estimate based on raw model output for about 15±6 % (18±8 %) of all grid points (Fig. 3a). If we ask for a reduction in bias by at least 50 %, 27±10 % (32±12 %) of grid points cannot reach this benchmark (Fig. 3b). Note that this means that in the majority of grid points, UQM reduces biases in WBGT. However, in many cases the reduction is not satisfactory. Understanding the conditions under which bias adjustment fails may help to design and use approaches that are more suitable for the given target. In contrast to WBGT, for CBI, UQM efficiently reduces the bias in most cases. Relative humidity alone explains most of the variance in CBI (Fig. 1); hence, UQM can correct it very well (Figs. 2–3). If instead of UQM the bias adjustment used in ISIMIP2b is applied, the numbers mostly fall into the CMIP5 range for WBGT (Fig. 3). However, for CBI, more grid points show no substantial improvement after bias adjustment compared to CMIP5. This is probably due to the fact that UQM is extremely flexible and thus almost perfectly adjusts temperature and relative humidity. This in turn leads to a very good adjustment of CBI, whose variability largely follows relative humidity (Fig. 1), whereas the bias adjustment used in ISIMIP2b is more conservative and therefore less flexible in adjusting CBI.
Regions for which biases only slightly increase or decrease include locations where biases are small to start with. To study for how many grid points this is the case, we compute the fraction of grid points for which the bias in WBGT is larger than 1 K either before or after bias adjustment. This is the case for 50 %–90 % of the grid points, depending on the model and the metric (Fig. 4a). Recomputing Fig. 3 based on this subset reduces the fraction of locations where bias adjustment does not achieve the two benchmarks by about half (Fig. 4b, c). Because these numbers strongly depend on the size of the accepted bias (1 K in our example), we continue the analysis with all grid points.
In Australia, the Sahel, and some parts in sub-Saharan Africa and South America, UQM increases biases in the 90th percentile of WBGT for a large fraction of models (Fig. 5a). In those regions, but also some other areas in the world, UQM reduces WBGT biases in the 90th percentile by less than 50 % in the majority of all models (Fig. 5b). Overall, in nearly 11 % of the land area, for more than 50 % of the models, UQM does not reduce the bias of the 90th percentile of WBGT by more than 50 %. MBCn, on the other hand, is able to fit the hazard estimates well where UQM fails (not shown for all grid points but Figs. 1d and 2c illustrate the worst case).
To improve the usage of bias-adjustment methods it would be important to know whether we can make any a priori statements as to whether UQM adjustment will lead to a reduction in biases of modeled multivariate hazards or impacts. Generally, regions where UQM fails to improve biases in WBGT coincide with regions where the observed correlations between temperature and relative humidity is outside the CMIP5 range (Fig. 6a). UQM fails when (univariate) mean biases in temperature are small (Fig. 6b), and when the correlation between the driving variables (in this case temperature and relative humidity) are not well captured (Fig. 6d). The mean bias in relative humidity looks similar for both cases (Fig. 6c). Overall this means that if mean biases in climate drivers are large, any bias adjustment will lead to a substantial reduction in biases of hazard or impact indicators (consistent with Fig. 4).
Internal variability can impact the effectiveness of bias adjustment because it introduces uncertainties that may be amplified by incomplete bias adjustment. For temperature and relative humidity individually, the perfect model approach reveals little difference between the range of the noise (uncertainty associated with internal variability) and the uncertainty range of all models (Fig. 7 illustrates the approach for the grid point of Austin, Texas, US). This is to be expected, as these driver variables can be perfectly adjusted. In the case of WBGT and CBI, however, UQM leads to much larger uncertainties for the full range than for the noise for some percentiles (Fig. 7). Overall, the full range of bias-adjusted model simulations for the 90th percentile of WBGT is a factor of 10 larger than the range expected from internal variability alone in many regions, including the Amazon, eastern North America and Indonesia (Fig. 8a). The regions with large increases in uncertainty roughly coincide with regions where the between-model variability in the correlation between temperature and relative humidity is very high compared to the variability within the CanESM runs (Fig. 8b). This is consistent with Fig. 6d, which shows that univariate bias adjustment is not very successful when the correlation structure between models and reference is not well matched. For a fixed percentile threshold, the fraction of grid points for which the full range is at least 2 times the range of the noise varies between 10 % and 40 % for CBI, and between 30 % and 70 % for WBGT (Fig. 9). If we use MBCn to adjust biases, the uncertainty associated with internal variability is largely indistinguishable from the full uncertainty range (Fig. 9, dashed lines).
We find that UQM of temperature and relative humidity often does not lead to a substantial reduction in WBGT biases. In a sizable number of cases UQM even leads to an increase in the original biases, particularly for high percentiles that are potentially most impact relevant. The fire risk indicator CBI is less affected by these issues because, although temperature is required for calculating it, it is largely dominated by variations in relative humidity. Our findings on the limited effectiveness of univariate bias adjustment are admittedly based on one bias-adjustment approach and two hazard indicators. Nevertheless, we expect that our results also hold for other multivariate hazards or impacts whose drivers are adjusted with bias-adjustment approaches that do not explicitly correct for dependencies between variables.
We show that for nearly 40 % of the land area, UQM fails to reduce biases in the 90th percentile of WBGT by at least 50 % in more than half of the models. The challenges in correcting high percentiles of WBGT suggest that a direct application of UQM to a warmer climate may lead to large errors. Multivariate bias adjustment such as MBCn offer remedies, though for a generic impact modeling project such as ISIMIP, all variables and dependencies would need to be corrected at once, requiring large amounts of good climate observations to fill the high-dimensional data space (Cannon, 2018). Furthermore, observational datasets would need to be carefully tested as to whether they represent the desired dependencies correctly (Cortés-Hernández et al., 2016; Zscheischler and Seneviratne, 2017).
The perfect model approach demonstrates that uncertainties related to internal variability and the use of an incomplete bias adjustment can lead to substantial uncertainty in multivariate hazards, particularly for more extreme percentiles. The large increase in uncertainty illustrated by the full range stems from a combination of initial uncertainty related to internal variability and the inability of UQM to adjust dependencies adequately. Therefore, uncertainties related to internal variability or other types of uncertainties may be strongly increased by incomplete bias adjustment. These uncertainties are likely sensitive to the length of the time period, as longer time periods will reduce the noise component related to internal variability due to a better sampling of the tails. These types of uncertainty are typically not communicated and accounted for in impact modeling chains. Hence, to be fully transparent, impact modeling should account for uncertainties associated with the chosen bias-adjustment approach in addition to uncertainties related to the choice of climate model and impact model.
Our results challenge the general applicability of conclusions from previous studies that have investigated whether a bias adjustment that is separately applied to all components of a hazard indicator can effectively reduce biases in the latter. For instance, Yang et al. (2015) studied the effect of component-wise bias adjustment on a fire risk indicator in Sweden and conclude that bias adjusting more drivers reduces biases in the indicator without the need to consider the dependence between the drivers. Similarly, Casanueva et al. (2018) studied the implications of component-wise bias adjustment on fire risk in Spain and find that there is little difference between separately adjusting the drivers before computing the hazard indicator and adjusting the hazard estimate directly. Räty et al. (2018) conclude, for hydrological variables, that for most impacts, univariate methods have a comparable performance to bivariate methods. In contrast, Cannon (2018) clearly shows that including the dependence structure of drivers into the bias-adjustment procedure strongly reduces biases in the Canadian Fire Weather Index. Thus, we cannot draw the general conclusion that multivariate bias adjustment is not necessary in any case from individual, typically regional, studies. Our findings suggest that whether bias-adjustment approaches lead to a substantial improvement of impact indicators ultimately depends on (i) how large the initial model biases are, (ii) how strongly the indicator truly depends on multiple variables, and (iii) how well the models simulate relevant dependencies between the climate variables (Fig. 6). Overall, it is difficult to pin down under which exact circumstances univariate bias adjustment might fail. We assume that modeled impacts that truly depend on multiple dependent climate drivers are particularly susceptible to these issues. Impacts that fall into that category include heat-related mortality (mostly depending on temperature and humidity but also other factors), agricultural yields and carbon uptake (temperature, precipitation, radiation), drought (precipitation, evapotranspiration), and the security of energy supply in a system that combines input from various renewable energy sources (Sterl et al., 2018). Impacts that might be less affected by these issues are those that predominantly depend on a single climate variable such as runoff or floods, which mostly rely on precipitation (Räty et al., 2018). In these cases, the adjustment of the spatial and temporal distribution of precipitation might by more relevant than the adjustment of dependencies between precipitation and other climate variables. A bias-adjustment method that is able to deal with very high dimensionality, for instance occurring when adjusting the covariance between many locations at the same time, was recently proposed by Vrac (2018).
While efforts to improve climate models and reduce their biases will continue (Davin et al., 2016; Kay et al., 2016; Wang et al., 2014), impact assessments need to become more transparent. In particular, uncertainties need to be well communicated to aid adaptation planning (Wilby and Dessai, 2010). Many artifacts of currently widely used bias-adjustment methods are probably unknown because of the way bias adjustment is evaluated (Addor and Fischer, 2015; Maraun et al., 2017). Evaluation of multivariate relationships and extremes need to become standard in the evaluation of climate and impact models as well as an evaluation of the appropriateness of the chosen bias-adjustment approach (Cortés-Hernández et al., 2016; Zscheischler et al., 2018). Overall, multivariate bias-adjustment methods should be favored in impact modeling to ensure that multivariate impacts are captured more realistically.
As long as major biases in climate models persist, some form of bias adjustment is unavoidable in model climate impacts. In the fortunate case where large multimodel ensembles are available and a clear target is identified, ensemble members can be selected according to how well they match certain criteria associated with the target (Herger et al., 2018; Maraun et al., 2017; Sippel et al., 2016). Process-based observational constraints (Hall and Qu, 2006; Sippel et al., 2017; Vogel et al., 2018) are one way forward to select the most promising subset of models (Maraun et al., 2017). Yet, in many cases, this information is not available and large model ensembles may be too expensive to obtain. To achieve the most reliable outcome with respect to modeled impacts, bias adjustment preferably takes into account the known limitations of the relevant climate model as well as characteristics of the target system for which the bias adjustment is applied (Maraun et al., 2017). In large-scale impact modeling projects such as ISIMIP, global flood modeling (Winsemius et al., 2015), global crop modeling (Ruane et al., 2017), and generic modeling of high-impact events (Done et al., 2015), achieving these standards is extremely challenging, if not unfeasible. In these cases, typically, a single bias-adjustment method is applied to a set of climate variables that then serves as input for a variety of impact models across regions and sectors. The usefulness of such one-size-fits-all approaches may be debated (Maraun et al., 2017), yet decision makers urgently require robust information on potential impacts of climate change. Impact modeling frameworks such as ISIMIP and many others sample the uncertainty related to the chosen climate model and impact model but do not take into account uncertainties related to the applied bias-adjustment method. To be transparent towards potential users, however, the scientific community should provide information on all uncertainties associated with modeled impacts.
Climate impact modeling is crucial to translate information from climate projections into potential impacts to aid decision-making and planning. Due to persistent biases in current climate models, bias adjustment is an integral part of most impact modeling activities. Our results demonstrate that univariate bias adjustment can increase biases of simulated hazards that depend on multiple correlated climatic drivers. Univariate bias adjustment can furthermore lead to a large increase in the uncertainty of such modeled hazards and impacts. Both aspects are particularly severe when studying extremes. More importantly, if univariate bias adjustment does not adequately adjust biases in hazards in the present-day climate, future projections of such hazards have to be interpreted very carefully. Our findings highlight that impact modeling chains need to incorporate uncertainties associated with the choice of bias adjustment into their uncertainty assessment to be transparent for decision makers.
All datasets and the code to compute the different bias adjustments are freely available from the sources mentioned in Sect. 2 (Data and methods).
JZ and EMF conceived the study. SL computed the bias adjustment following the ISIMIP2b scheme. JZ analyzed all data and produced all figures. JZ wrote the first draft, all authors commented on the draft and all revisions.
The authors declare that they have no conflict of interest.
We thank Alex Cannon for helpful discussions related to the application of the MBCn
approach. Jakob Zscheischler acknowledges financial support from the SNSF (Ambizione
grant PZ00P2_179876). Stefan Lange acknowledges funding from the European Union's
Horizon 2020 research and innovation program under grant agreement no. 641816
Edited by: Ben Kravitz
Reviewed by: Stefan Hagemann and one anonymous referee
Addor, N. and Fischer, E. M.: The influence of natural variability and interpolation errors on bias characterization in RCM simulations, J. Geophys. Res.-Atmos., 120, 10180–10195, https://doi.org/10.1002/2014JD022824, 2015. a, b, c
Bosshard, T., Carambia, M., Goergen, K., Kotlarski, S., Krahe, P., Zappa, M., and Schär, C.: Quantifying uncertainty sources in an ensemble of hydrological climate-impact projections, Water Resour. Res., 49, 1523–1536, https://doi.org/10.1029/2011WR011533, 2013. a
Brando, P. M., Balch, J. K., Nepstad, D. C., Morton, D. C., Putz, F. E., Coe, M. T., Silvério, D., Macedo, M. N., Davidson, E. A., Nóbrega, C. C., Alencar, A., and Soares-Filho, B. S.: Abrupt increases in Amazonian tree mortality due to drought–fire interactions, P. Natl. Acad. Sci. USA, 111, 6347–6352, 2014. a
Bröde, P., Blazejczyk, K., Fiala, D., Havenith, G., Holmér, I., Jendritzky, G., Kuklane, K., and Kampmann, B.: The Universal Thermal Climate Index UTCI Compared to Ergonomics Standards for Assessing the Thermal Environment, Ind. Health, 51, 16–24, https://doi.org/10.2486/indhealth.2012-0098, 2013. a
Cannon, A. J.: Multivariate Bias Correction of Climate Model Output: Matching Marginal Distributions and Intervariable Dependence Structure, J. Climate, 29, 7045–7064, https://doi.org/10.1175/jcli-d-15-0679.1, 2016. a
Cannon, A. J.: Multivariate quantile mapping bias correction: an N-dimensional probability density function transform for climate model simulations of multiple variables, Clim. Dynam., 50, 31–49, https://doi.org/10.1007/s00382-017-3580-6, 2018. a, b, c, d, e
Carlson, J. D. and Burgan, R. E.: Review of users' needs in operational fire danger estimation: The Oklahoma example, Int. J. Remote Sens., 24, 1601–1620, https://doi.org/10.1080/01431160210144651, 2003. a
Casanueva, A., Bedia, J., Herrera, S., Fernández, J., and Gutiérrez, J. M.: Direct and component-wise bias correction of multi-variate climate indices: the percentile adjustment function diagnostic tool, Climatic Change, 147, 411–425, https://doi.org/10.1007/s10584-018-2167-5, 2018. a, b, c
Chandler, C., Cheney, P., Thomas, P., Trabaud, L., and Williams, D.: Fire in forestry. Volume 1. Forest fire behavior and effects., John Wiley & Sons, Inc., New York, USA, 1983. a
Chen, C., Haerter, J. O., Hagemann, S., and Piani, C.: On the contribution of statistical bias correction to the uncertainty in the projected hydrological cycle, Geophys. Res. Lett., 38, L20 403, https://doi.org/10.1029/2011GL049318, 2011. a, b
Christensen, J. H., Boberg, F., Christensen, O. B., and Lucas-Picher, P.: On the need for bias correction of regional climate change projections of temperature and precipitation, Geophys. Res. Lett., 35, L20709, https://doi.org/10.1029/2008GL035694, 2008. a
Cortés-Hernández, V. E., Zheng, F., Evans, J., Lambert, M., Sharma, A., and Westra, S.: Evaluating regional climate models for simulating sub-daily rainfall extremes, Clim. Dynam., 47, 1613–1628, https://doi.org/10.1007/s00382-015-2923-4, 2016. a, b
Cramer, W., Yohe, G. W., Auffhammer, M., Huggel, C., Molau, U., Dias, M. A. F. S., Solow, A., Stone, D. A., and Tibig, L.: Detection and attribution of observed impacts, in: Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel of Climate Change, edited by: Field, C. B., Barros, V. R., Dokken, D. J., Mach, K. J., Mastrandrea, M. D., Bilir, T. E., Chatterjee, M., Ebi, K. L., Estrada, Y. O., Genova, R. C., Girma, B., Kissel, E. S., Levy, A. N., MacCracken, S., Mastrandrea, P. R., and White, L. L., 979–1037, Cambridge University Press, Cambridge, UK and New York, NY, USA, 2014. a
Davin, E. L., Maisonnave, E., and Seneviratne, S. I.: Is land surface processes representation a possible weak link in current Regional Climate Models?, Environ. Res. Lett., 11, 074027, https://doi.org/10.1088/1748-9326/11/7/074027, 2016. a
Dee, D., Uppala, S., Simmons, A., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, 2011. a
Done, J. M., Holland, G. J., Bruyère, C. L., Leung, L. R., and Suzuki-Parker, A.: Modeling high-impact weather and climate: lessons from a tropical cyclone perspective, Clim. Change, 129, 381–395, https://doi.org/10.1007/s10584-013-0954-6, 2015. a
Ehret, U., Zehe, E., Wulfmeyer, V., Warrach-Sagi, K., and Liebert, J.: HESS Opinions “Should we apply bias correction to global and regional climate model data?”, Hydrol. Earth Syst. Sci., 16, 3391–3404, https://doi.org/10.5194/hess-16-3391-2012, 2012. a
Elía, R. d., Laprise, R., and Denis, B.: Forecasting Skill Limits of Nested, Limited-Area Models: A Perfect-Model Approach, Mon. Weather Rev., 130, 2006–2023, https://doi.org/10.1175/1520-0493(2002)130<2006:fslonl>2.0.co;2, 2002. a
Fischer, E., Beyerle, U., and Knutti, R.: Robust spatially aggregated projections of climate extremes, Nat. Clim. Change, 3, 1033–1038, 2013. a
Flato, G., Marotzke, J., Abiodun, B., Braconnot, P., Chou, S., Collins, W., Cox, P., Driouech, F., Emori, S., Eyring, V., Forest, C., Gleckler, P., Guilyardi, E., Jakob, C., Kattsov, V., Reason, C., and Rummukainen, M.: Evaluation of Climate Models, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T., Qin, D., Plattner, G.-K., Tignor, M., Allen, S., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P., 741–866, Cambridge University Press, Cambridge, UK and New York, NY, USA, https://doi.org/10.1017/CBO9781107415324.020, 2013. a
Frieler, K., Lange, S., Piontek, F., Reyer, C. P. O., Schewe, J., Warszawski, L., Zhao, F., Chini, L., Denvil, S., Emanuel, K., Geiger, T., Halladay, K., Hurtt, G., Mengel, M., Murakami, D., Ostberg, S., Popp, A., Riva, R., Stevanovic, M., Suzuki, T., Volkholz, J., Burke, E., Ciais, P., Ebi, K., Eddy, T. D., Elliott, J., Galbraith, E., Gosling, S. N., Hattermann, F., Hickler, T., Hinkel, J., Hof, C., Huber, V., Jägermeyr, J., Krysanova, V., Marcé, R., Müller Schmied, H., Mouratiadou, I., Pierson, D., Tittensor, D. P., Vautard, R., van Vliet, M., Biber, M. F., Betts, R. A., Bodirsky, B. L., Deryng, D., Frolking, S., Jones, C. D., Lotze, H. K., Lotze-Campen, H., Sahajpal, R., Thonicke, K., Tian, H., and Yamagata, Y.: Assessing the impacts of 1.5 ∘C global warming – simulation protocol of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP2b), Geosci. Model Dev., 10, 4321–4345, https://doi.org/10.5194/gmd-10-4321-2017, 2017. a, b, c, d, e
Gudmundsson, L.: qmap: Statistical transformations for postprocessing climate model output, R package version 1.0-2, 2014. a
Hawkins, E., Robson, J., Sutton, R., Smith, D., and Keenlyside, N.: Evaluating the potential for statistical decadal predictions of sea surface temperatures with a perfect model approach, Clim. Dynam., 37, 2495–2509, https://doi.org/10.1007/s00382-011-1023-3, 2011. a
Hempel, S., Frieler, K., Warszawski, L., Schewe, J., and Piontek, F.: A trend-preserving bias correction – the ISI-MIP approach, Earth Syst. Dynam., 4, 219–236, https://doi.org/10.5194/esd-4-219-2013, 2013. a, b
Herger, N., Abramowitz, G., Knutti, R., Angélil, O., Lehmann, K., and Sanderson, B. M.: Selecting a climate model subset to optimise key ensemble properties, Earth Syst. Dynam., 9, 135–151, https://doi.org/10.5194/esd-9-135-2018, 2018. a
Humphrey, V., Zscheischler, J., Ciais, P., Gudmundsson, L., Sitch, S., and Seneviratne, S.: Sensitivity of atmospheric CO2 growth rate to observed changes in terrestrial water storage, Nature, 560, 628–631, https://doi.org/10.1038/s41586-018-0424-4, 2018. a
Kay, J. E., Wall, C., Yettella, V., Medeiros, B., Hannay, C., Caldwell, P., and Bitz, C.: Global Climate Impacts of Fixing the Southern Ocean Shortwave Radiation Bias in the Community Earth System Model (CESM), J. Climate, 29, 4617–4636, https://doi.org/10.1175/jcli-d-15-0358.1, 2016. a
Kotlarski, S., Keuler, K., Christensen, O. B., Colette, A., Déqué, M., Gobiet, A., Goergen, K., Jacob, D., Lüthi, D., van Meijgaard, E., Nikulin, G., Schär, C., Teichmann, C., Vautard, R., Warrach-Sagi, K., and Wulfmeyer, V.: Regional climate modeling on European scales: a joint standard evaluation of the EURO-CORDEX RCM ensemble, Geosci. Model Dev., 7, 1297–1333, https://doi.org/10.5194/gmd-7-1297-2014, 2014. a
Li, C., Sinha, E., Horton, D. E., Diffenbaugh, N. S., and Michalak, A. M.: Joint bias correction of temperature and precipitation in climate model simulations, J. Geophys. Res.-Atmos., 119, 13153–13162, https://doi.org/10.1002/2014JD022514, 2014. a
Li, C., Zhang, X., Zwiers, F., Fang, Y., and Michalak, A. M.: Recent Very Hot Summers in Northern Hemispheric Land Areas Measured by Wet Bulb Globe Temperature Will Be the Norm Within 20 Years, Earth's Future, 5, 1203–1216, https://doi.org/10.1002/2017EF000639, 2017. a
Maraun, D., Shepherd, T. G., Widmann, M., Zappa, G., Walton, D., Gutierrez, J. M., Hagemann, S., Richter, I., Soares, P. M. M., Hall, A., and Mearns, L. O.: Towards process-informed bias correction of climate change simulations, Nat. Clim. Change, 7, 764–773, https://doi.org/10.1038/nclimate3418, 2017. a, b, c, d, e, f
Martius, O., Pfahl, S., and Chevalier, C.: A global quantification of compound precipitation and wind extremes, Geophys. Res. Lett., 43, 7709–7717, 2016. a
McMichael, A. J., Woodruff, R. E., and Hales, S.: Climate change and human health: present and future risks, Lancet, 367, 859–869, 2006. a
Mehrotra, R. and Sharma, A.: A Multivariate Quantile-Matching Bias Correction Approach with Auto- and Cross-Dependence across Multiple Time Scales: Implications for Downscaling, J. Climate, 29, 3519–3539, https://doi.org/10.1175/JCLI-D-15-0356.1, 2016. a
Pal, J. S. and Eltahir, E. A. B.: Future temperature in southwest Asia projected to exceed a threshold for human adaptability, Nat. Clim. Change, 6, 197–200, https://doi.org/10.1038/nclimate2833, 2015. a
Panofsky, H. and Brier, G.: Some Applications of Statistics to Meteorology, The Pennsylvania State University, University Park, PA, 1968. a
Papadimitriou, L. V., Koutroulis, A. G., Grillakis, M. G., and Tsanis, I. K.: The effect of GCM biases on global runoff simulations of a land surface model, Hydrol. Earth Syst. Sci., 21, 4379–4401, https://doi.org/10.5194/hess-21-4379-2017, 2017. a
Piani, C. and Haerter, J. O.: Two dimensional bias correction of temperature and precipitation copulas in climate models, Geophys. Res. Lett., 39, l20401, https://doi.org/10.1029/2012GL053839, 2012. a
Piani, C., Haerter, J. O., and Coppola, E.: Statistical bias correction for daily precipitation in regional climate models over Europe, Theor. Appl. Climatol., 99, 187–192, https://doi.org/10.1007/s00704-009-0134-9, 2010. a
Räty, O., Räisänen, J., Bosshard, T., and Donnelly, C.: Intercomparison of Univariate and Joint Bias Correction Methods in Changing Climate From a Hydrological Perspective, Climate, 6, 33, https://doi.org/10.3390/cli6020033, 2018. a, b, c
Ruane, A. C., Rosenzweig, C., Asseng, S., Boote, K. J., Elliott, J., Ewert, F., Jones, J. W., Martre, P., McDermid, S. P., Müller, C., Snyder, A., and Thorburn, P. J.: An AgMIP framework for improved agricultural representation in integrated assessment models, Environ. Res. Lett., 12, 125003, https://doi.org/10.1088/1748-9326/aa8da6, 2017. a, b, c
Sippel, S., Otto, F. E. L., Forkel, M., Allen, M. R., Guillod, B. P., Heimann, M., Reichstein, M., Seneviratne, S. I., Thonicke, K., and Mahecha, M. D.: A novel bias correction methodology for climate impact simulations, Earth Syst. Dynam., 7, 71–88, https://doi.org/10.5194/esd-7-71-2016, 2016. a
Sippel, S., Zscheischler, J., Mahecha, M. D., Orth, R., Reichstein, M., Vogel, M., and Seneviratne, S. I.: Refining multi-model projections of temperature extremes by evaluation against land-atmosphere coupling diagnostics, Earth Syst. Dynam., 8, 387–403, https://doi.org/10.5194/esd-8-387-2017, 2017. a
Sterl, S., Liersch, S., Koch, H., van Lipzig, N. P. M., and Thiery, W.: A new approach for assessing synergies of solar and wind power: implications for West Africa, Environ. Res. Lett., 13, 094009, https://doi.org/10.1088/1748-9326/aad8f6, 2018. a
Switanek, M. B., Troch, P. A., Castro, C. L., Leuprecht, A., Chang, H.-I., Mukherjee, R., and Demaria, E. M. C.: Scaled distribution mapping: a bias correction method that preserves raw climate model projected changes, Hydrol. Earth Syst. Sci., 21, 2649–2666, https://doi.org/10.5194/hess-21-2649-2017, 2017. a
Teutschbein, C. and Seibert, J.: Bias correction of regional climate model simulations for hydrological climate-change impact studies: Review and evaluation of different methods, J. Hydrol., 456–457, 12–29, 2012. a
Van Loon, A. F., Gleeson, T., Clark, J., Van Dijk, A. I. J. M., Stahl, K., Hannaford, J., Di Baldassarre, G., Teuling, A. J., Tallaksen, L. M., Uijlenhoet, R., Hannah, D. M., Sheffield, J., Svoboda, M., Verbeiren, B., Wagener, T., Rangecroft, S., Wanders, N., and Van Lanen, H. A. J.: Drought in the Anthropocene, Nat. Geosci., 9, 89–91, https://doi.org/10.1038/ngeo2646, 2016. a
Vogel, M. M., Zscheischler, J., and Seneviratne, S. I.: Varying soil moisture-atmosphere feedbacks explain divergent temperature extremes and precipitation projections in central Europe, Earth Syst. Dynam., 9, 1107–1125, https://doi.org/10.5194/esd-9-1107-2018, 2018. a
Vorogushyn, S., Bates, P. D., de Bruijn, K., Castellarin, A., Kreibich, H., Priest, S., Schröter, K., Bagli, S., Blöschl, G., Domeneghetti, A., Gouldby, B., Klijn, F., Lammersen, R., Neal, J. C., Ridder, N., Terink, W., Viavattene, C., Viglione, A., Zanardo, S., and Merz, B.: Evolutionary leap in large-scale flood risk assessment needed, Wiley Interdisciplinary Reviews: Water, 5, e1266, https://doi.org/10.1002/wat2.1266, 2018. a
Vrac, M.: Multivariate bias adjustment of high-dimensional climate simulations: the Rank Resampling for Distributions and Dependences (R2D2) bias correction, Hydrol. Earth Syst. Sci., 22, 3175–3196, https://doi.org/10.5194/hess-22-3175-2018, 2018. a, b
Vrac, M. and Friederichs, P.: Multivariate-intervariable, spatial, and temporal-bias correction, J. Climate, 28, 218–237, 2015. a
Winsemius, H. C., Aerts, J. J. H., van Beek, L. H., Bierkens, M. P., Bouwman, A., Jongman, B., Kwadijk, J. J., Ligtvoet, W., Lucas, P., vanVuuren, D., and Ward, P.: Global drivers of future river flood risk, Nat. Clim. Change, 6, 381–385, https://doi.org/10.1038/nclimate2893, 2015. a, b, c
Yang, W., Gardelin, M., Olsson, J., and Bosshard, T.: Multi-variable bias correction: application of forest fire risk in present and future climate in Sweden, Nat. Hazards Earth Syst. Sci., 15, 2037–2057, https://doi.org/10.5194/nhess-15-2037-2015, 2015. a, b
Zhao, Y., Ducharne, A., Sultan, B., Braconnot, P., and Vautard, R.: Estimating heat stress from climate-based indicators: present-day biases and future spreads in the CMIP5 global climate model ensemble, Environ. Res. Lett., 10, 084013, https://doi.org/10.1088/1748-9326/10/8/084013, 2015. a
Zscheischler, J., Orth, R., and Seneviratne, S. I.: Bivariate return periods of temperature and precipitation explain a large fraction of European crop yields, Biogeosciences, 14, 3309–3320, https://doi.org/10.5194/bg-14-3309-2017, 2017. a
Zscheischler, J., Westra, S., van den Hurk, B. J. J., Pitman, A., Ward, P., Bresch, D. N., Leonard, M., Zhang, X., AghaKouchak, A., Wahl, T., and Seneviratne, S. I.: Future climate risk from compound events, Nat. Clim. Change, 8, 469–477, https://doi.org/10.1038/s41558-018-0156-3, 2018. a, b, c, d