Climate Projections Using Bayesian Model Averaging and Space–Time Dependence

Bayesian weighting of climate models based on climate sensitivity

Using climate model ensembles containing members that exhibit very high climate sensitivities to increasing CO2 concentrations can result in biased projections. Various methods have been proposed to ameliorate this 'hot model' problem, such as model emulators or model culling. Here, we utilize Bayesian Model Averaging as a framework to address this problem without resorting to outright rejection of models from the ensemble. Taking advantage of multiple lines of evidence used to construct the best estimate of the earth's climate sensitivity, the Bayesian Model Averaging framework produces an unbiased posterior probability distribution of model weights. The updated multi-model ensemble projects end-of-century global mean surface temperature increases of 2 oC for a low emissions scenario (SSP1-2.6) and 5 oC for a high emissions scenario (SSP5-8.5). These estimates are lower than those produced using a simple multi-model mean for the CMIP6 ensemble. The results are also similar to results from a model culling approach, but retain some weight on low-probability models, allowing for consideration of the possibility that the true value could lie at the extremes of the assessed distribution. Our results showcase Bayesian Model Averaging as a path forward to project future climate change that is commensurate with the available scientific evidence.

Regional and tele-connected impacts of the Tibetan Plateau surface darkening

Despite knowledge of the presence of the Tibetan Plateau (TP) in reorganizing large-scale atmospheric circulation, it remains unclear how surface albedo darkening over TP will impact local glaciers and remote Asian monsoon systems. Here, we use a coupled land-atmosphere global climate model and a glacier model to address these questions. Under a high-emission scenario, TP surface albedo darkening will increase local temperature by 0.24 K by the end of this century. This warming will strengthen the elevated heat pump of TP, increasing South Asian monsoon precipitation while exacerbating the current "South Flood-North Drought" pattern over East Asia. The albedo darkening-induced climate change also leads to an accompanying TP glacier volume loss of 6.9%, which further increases to 25.2% at the equilibrium, with a notable loss in western TP. Our findings emphasize the importance of land-surface change responses in projecting future water resource availability, with important implications for water management policies.

Combining global climate models using graph cuts

Global Climate Models are the main tools for climate projections. Since many models exist, it is common to use Multi-Model Ensembles to reduce biases and assess uncertainties in climate projections. Several approaches have been proposed to combine individual models and extract a robust signal from an ensemble. Among them, the Multi-Model Mean (MMM) is the most commonly used. Based on the assumption that the models are centered around the truth, it consists in averaging the ensemble, with the possibility of using equal weights for all models or to adjust weights to favor some models. In this paper, we propose a new alternative to reconstruct multi-decadal means of climate variables from a Multi-Model Ensemble, where the local performance of the models is taken into account. This is in contrast with MMM where a model has the same weight for all locations. Our approach is based on a computer vision method called graph cuts and consists in selecting for each grid point the most appropriate model, while at the same time considering the overall spatial consistency of the resulting field. The performance of the graph cuts approach is assessed based on two experiments: one where the ERA5 reanalyses are considered as the reference, and another involving a perfect model experiment where each model is in turn considered as the reference. We show that the graph cuts approach generally results in lower biases than other model combination approaches such as MMM, while at the same time preserving a similar level of spatial continuity.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s00382-022-06213-4.

A Bayesian ensemble approach to combine PM2.5 estimates from statistical models using satellite imagery and numerical model simulation

Ambient fine particulate matter less than 2.5 μm in aerodynamic diameter (PM_2.5) has been linked to various adverse health outcomes. PM_2.5 arises from both natural and anthropogenic sources, and PM_2.5 concentrations can vary over space and time. However, the sparsity of existing air quality monitors greatly restricts the spatial-temporal coverage of PM_2.5 measurements, potentially limiting the accuracy of PM_2.5-related health studies. Various methods exist to address these limitations by supplementing air quality monitoring measurements with additional data. We develop a method to combine PM_2.5 estimated from satellite-retrieved aerosol optical depth (AOD) and chemical transport model (CTM) simulations using statistical models. While most previous methods utilize AOD or CTM separately, we aim to leverage advantages offered by both data sources in terms of resolution and coverage using Bayesian ensemble averaging. Our approach differs from previous ensemble approaches in its ability to not only incorporate uncertainties in PM_2.5 estimates from individual models but also to provide uncertainties for the resulting ensemble estimates. In an application of estimating daily PM_2.5 in the Southeastern US, the ensemble approach outperforms previously developed spatial-temporal statistical models that use either AOD or bias-corrected CTM simulations in cross-validation (CV) analyses. More specifically, in spatially clustered CV experiments, the ensemble approach reduced the AOD-only and CTM-only model's root mean squared error (RMSE) by at least 13%. Similar improvements were seen in R². The enhanced prediction performance that the ensemble technique provides at fine-scale spatial resolution, as well as the availability of prediction uncertainty, can be further used in health effect analyses of air pollution exposure.

A novel method to test non-exclusive hypotheses applied to Arctic ice projections from dependent models

A major conundrum in climate science is how to account for dependence between climate models. This complicates interpretation of probabilistic projections derived from such models. Here we show that this problem can be addressed using a novel method to test multiple non-exclusive hypotheses, and to make predictions under such hypotheses. We apply the method to probabilistically estimate the level of global warming needed for a September ice-free Arctic, using an ensemble of historical and representative concentration pathway 8.5 emissions scenario climate model runs. We show that not accounting for model dependence can lead to biased projections. Incorporating more constraints on models may minimize the impact of neglecting model non-exclusivity. Most likely, September Arctic sea ice will effectively disappear at between approximately 2 and 2.5 K of global warming. Yet, limiting the warming to 1.5 K under the Paris agreement may not be sufficient to prevent the ice-free Arctic.

Accounting for skill in trend, variability, and autocorrelation facilitates better multi-model projections: Application to the AMOC and temperature time series

We present a novel quasi-Bayesian method to weight multiple dynamical models by their skill at capturing both potentially non-linear trends and first-order autocorrelated variability of the underlying process, and to make weighted probabilistic projections. We validate the method using a suite of one-at-a-time cross-validation experiments involving Atlantic meridional overturning circulation (AMOC), its temperature-based index, as well as Korean summer mean maximum temperature. In these experiments the method tends to exhibit superior skill over a trend-only Bayesian model averaging weighting method in terms of weight assignment and probabilistic forecasts. Specifically, mean credible interval width, and mean absolute error of the projections tend to improve. We apply the method to a problem of projecting summer mean maximum temperature change over Korea by the end of the 21st century using a multi-model ensemble. Compared to the trend-only method, the new method appreciably sharpens the probability distribution function (pdf) and increases future most likely, median, and mean warming in Korea. The method is flexible, with a potential to improve forecasts in geosciences and other fields.