In many applications, it is desirable to aggregate climate model projections by combining multiple models into a single projection that aims to leverage their collective strengths, often resulting in improved performance compared to individual models. While climate models exhibit varying levels of global average bias, their local performance often displays significantly larger biases—sometimes by an order of magnitude—with each model showcasing distinct strengths and weaknesses in different regions. Aggregating models without accounting for these spatial differences can degrade the quality of projections by diluting strong regional signals from high-performing models. While many approaches ranging in complexity have been developed, including the commonly used Multi-Model Mean (MMM) and weighted MMM, these methods typically apply a global weighting to the models, overlooking the fact that certain models may excel only in specific regions.
To date, the Graph Cut optimization method (Thao et al., 2022) stands out as one of the few techniques effectively leveraging the local capabilities of different models across multi-decadal periods to produce global projections. This method involves selecting the best performing model for each grid point while also ensuring the spatial consistency of the resulting fields. Despite its promising results, which surpass those of other ensemble combination techniques, it is restricted to optimizing for a single variable. This limitation causes inconsistent model selection across variables in multivariate scenarios. This leads to a loss of the multivariate relationships captured in the models. Furthermore, this technique was limited to multi-decadal averages, and is thus unable to capture the distributional characteristics of climate variables, including extreme and compound events.
Here, we present significant enhancements to the Graph Cut optimization method, enabling the combination of distributions of daily values. This approach preserves multivariate relationships, better capturing the complete span of climate dynamics. By employing the Hellinger distance to assess model performance, we can identify, at each grid point, the model that most accurately represents the multivariate distribution of target variables (e.g., temperature, pressure, and precipitation), minimizing the emergence of unrealistic discontinuities in the combined fields.
To demonstrate the use of our method, we combine 22 models from CMIP6 using three variables: temperature, precipitation, and sea level pressure, achieving better reproduction of ERA5 reanalysis compared to the Multi-Model Mean (MMM). Additionally, a perfect model experiment was conducted to evaluate the robustness and stability of the methodology under high climate change scenarios, such as SSP8.5, and over extended timescales reaching the end of the century. These results highlight the method's ability to maintain reliable performance and spatial consistency in challenging future conditions.
REFERENCES
Thao, S., Garvik, M., Mariethoz, G. et al. Combining global climate models using graph cuts. Clim Dyn 59, 2345–2361 (2022). https://doi.org/10.1007/s00382-022-06213-4