Multi-model Ensemble Research Articles

Seasonal forecasts have sufficient skill to inform some agricultural decisions

Abstract Seasonal forecasts, which look several months into the future, are currently underutilized in active decision-making, particularly for agricultural and natural resource management. This underutilization can be attributed to the absence of forecasts for decision-relevant variables at the required spatiotemporal resolution and at the time when the decisions are made and a perception of poor skill by decision-makers. Addressing these constraints, we quantified the skill of seasonal forecasts in informing two agricultural decisions with differing decision timeframes and influencer variables: (a) whether to apply fertilizer in fall or wait until spring based on expected winter temperatures, and (b) drought response, such as whether to lease water based on expectations of drought. We also looked into how early the forecast can be provided without significant degradation in skill. Currently, drought response decisions are typically formulated in April, utilizing drought forecasts issued in the same month, while fall fertilization decisions are generally made between August and September. There is growing interest among stakeholders in the availability of earlier forecasts to inform these critical choices. We utilized the North American Multi-Model Ensemble (NMME) hindcasts for the time period 1982 - 2020 over the Pacific Northwest US to obtain meteorological variables. Runoff was estimated via simulations of the coupled crop-hydrology VIC-CropSyst model. The skill assessment with the Heidke Skill Score (HSS) yielded promising outcomes in both decisions for the entire Pacific Northwest region. Notably, NMME’s positive skill (median HSS of 30%) in predicting warmer winters identifies years when fertilizer application should be avoided to prevent fertilizer loss through mineralization (and associated costs). Similarly, there is skill in forecasting drought conditions in most irrigated watersheds for up to two months in advance of April, the current decision time. In conclusion, our findings affirm that contrary to the perception of low skill and resulting underutilization, current seasonal forecasts hold the potential to inform at least some key agricultural decisions.

Tropical Interbasin Interaction as Effective Predictors of Late-Spring Precipitation Variability in the Southern Great Plains

Abstract The southern Great Plains experience fluctuating precipitation extremes that significantly impact agriculture and water management. Despite ongoing efforts to enhance forecast accuracy, the underlying causes of these climatic phenomena remain inadequately understood. This study elucidates the relative influence of the tropical Pacific and Atlantic basins on April–May–June precipitation variability in this region. Our partial ocean assimilation experiments using the Community Earth System Model unveil the prominent role of interbasin interaction, with the Pacific and Atlantic contributing approximately 70% and 30%, respectively, to these interbasin contrasts. Our statistical analyses suggest that these tropical interbasin contrasts could serve as a more reliable indicator for late-spring precipitation anomalies than El Niño–Southern Oscillation. The conclusions are reinforced by analyses of seven climate forecasting systems within the North American Multi-Model Ensemble, offering an optimistic outlook for enhancing real-time forecasting of late-spring precipitation in the southern plains. However, the current predictive skills of the interbasin contrasts across the prediction systems are hindered by the lower predictability of the tropical Atlantic Ocean, pointing to the need for future research to refine climate prediction models further. Significance Statement Agriculture and infrastructure in the southern plains face challenges from severe late-spring precipitation extremes. Traditional predictors like El Niño–Southern Oscillation (ENSO) lose effectiveness during the critical spring-to-summer transition, creating a forecasting gap. This study introduces the concept of tropical interbasin interactions, known to enhance seasonal predictability for late-spring precipitation in the southern plains. Novel climate model experiments highlight contributions from the tropical Pacific and Atlantic, offering a promising predictability that potentially surpasses the limitations of ENSO-based predictions. These outcomes hold the potential for developing operational forecasts of late-spring precipitation anomalies in the southern plains, enabling proactive risk management.

Impacts of climate change and deforestation on fire risk in the Amazon basin

Fires in the Amazon rainforest are a result of anthropogenic activity in the region. As the Amazon rainforest becomes vulnerable to fire through deforestation and climate change, it is imperative to understand how the fire risk will change in future years—not only for the potential of the forest to become a source of carbon to the atmosphere if burned but also for the potential loss of unique biodiversity and its critical role in regional and global climate. Therefore, we used the Eta regional climate model driven by 3 global climate models to simulate 2 distinct scenarios: one with a high CO2 emission (RCP8.5-DEF15) path and another with a high CO2 emission path and the total deforestation of the Amazon rainforest (RCP8.5-DEFTOT). Combining the multi-model ensemble of temperature and precipitation, we calculated the Keetch-Byram Drought Index (KBDI) and analyzed the results for the dry and wet seasons. In both scenarios, precipitation is reduced and maximum temperature increased. The greatest changes occur during the dry season and with the RCP8.5-DEFTOT scenario. In RCP8.5-DEF15, the KBDI is greatly affected, showing positive anomalies in both seasons. With the addition of deforestation, the fire risk is further exacerbated, thus highlighting the role of the forest in suppressing optimal conditions for fire. Additionally, more than half of the basin area is predicted to experience a high risk of fire in extreme years. Our results emphasize the importance of reducing global emissions and deforestation of the Amazon, given that fire activity could threaten the future of the Amazon rainforest.

A Forest Fire Prediction Model Based on Meteorological Factors and the Multi-Model Ensemble Method

More than half of South Korea’s land area is covered by forests, which significantly increases the potential for extensive damage in the event of a forest fire. The majority of forest fires in South Korea are caused by humans. Over the past decade, more than half of these types of fires occurred during the spring season. Although human activities are the primary cause of forest fires, the fact that they are concentrated in the spring underscores the strong association between forest fires and meteorological factors. When meteorological conditions favor the occurrence of forest fires, certain triggering factors can lead to their ignition more easily. The purpose of this study is to analyze the meteorological factors influencing forest fires and to develop a machine learning-based prediction model for forest fire occurrence, focusing on meteorological data. The study focuses on four regions within Gangwon province in South Korea, which have experienced substantial damage from forest fires. To construct the model, historical meteorological data were collected, surrogate variables were calculated, and a variable selection process was applied to identify relevant meteorological factors. Five machine learning models were then used to predict forest fire occurrence and ensemble techniques were employed to enhance the model’s performance. The performance of the developed forest fire prediction model was evaluated using evaluation metrics. The results indicate that the ensemble model outperformed the individual models, with a higher F1-score and a notable reduction in false positives compared to the individual models. This suggests that the model developed in this study, when combined with meteorological forecast data, can potentially predict forest fire occurrence and provide insights into the expected severity of fires. This information could support decision-making for forest fire management, aiding in the development of more effective fire response plans.

Global record-breaking recurrence rates indicate more widespread and intense surface air temperature and precipitation extremes.

We analyzed the evolution of extreme annual surface air temperature and rainfall on Earth, based on the recurrence rate of record-breaking events, and found the highest recurrence rates for record-high annual temperatures in the tropics, as opposed to the polar regions with the fastest warming. Both recurrence rates and the global surface area fraction with daily mean surface air temperatures exceeding 30° and 40°C provide further evidence for extremely hot years becoming more common and widespread. A similar analysis for precipitation highlighted some regions with more record-high annual total precipitation and others with record-low annual precipitation typically associated with drought. A multimodel ensemble of 306 runs with global climate models [Coupled Model Intercomparison Project phase 6 (CMIP) Shared Socioeconomic Pathways 2-45 (SSP2-45)] reproduced the statistics of record-breaking high temperatures, but there were some differences for the reanalysis precipitation record-breaking recurrence rates. The global climate model simulations suggested a slightly altered geographical pattern for record-breaking annual precipitation recurrence rates.

Assessing Alterations of Rainfall Variability Under Climate Change in Zengwen Reservoir Watershed, Southern Taiwan

This study aims to detect changes in rainfall variability caused by climate change for various scenarios in the CMIP6 (Coupled Model Intercomparison Project Phase 6) multi-model ensemble. Projected changes in rainfall unevenness in terms of different timescale indices using three categories, namely WD50 (number of wettest days for half annual rainfall), SI (seasonality index), and DWR (ratio of dry-season to wet-season rainfall) are analyzed in Zengwen Reservoir watershed, southern Taiwan over near future (2021–2040) and midterm future (2041–2060) relative to the baseline period (1995–2014) under SSP2-4.5 and SSP5-8.5 scenarios. The projected rainfall for both baseline and future periods is derived from 25 GCMs (global climate models). The results indicate that noticeably deteriorated rainfall unevenness is projected in the Zengwen Reservoir watershed over future periods, which include decreased WD50, increased SI, and decreased DWR. Though there were noticeable differences in the rainfall projections by the different GCMs, the overall consensus reveals that uncertainties in future rainfall should not be ignored. In addition, WD50 has the greatest deviated relative change in mean, which implies that the short-timescale rainfall unevenness index is easily affected by climate change in the study area. Distributional changes in rainfall unevenness determined by simultaneously considering alterations in relative changes in mean and standard deviation indicated that there was no single dominant category. However, the top two categories, with summed frequencies exceeding 0.5, characterize different properties of rainfall unevenness indices. The top two categories of WD50 and SI commonly have decreased mean and increased mean, respectively, but nearly equal frequencies of the top two categories in DWR exhibit opposite variations. The proposed rainfall unevenness change detection approach provides a better understanding of the impacts of climate change on rainfall unevenness, which is useful for preparing adaptive mitigation measures for coping with disasters induced by climate change.

Influence of Model Resolution on Wind Energy Simulations over Tibetan Plateau Using CMIP6 HighResMIP

The assessment of wind energy resources is critical for the transition from fossil fuel to renewable energy sources. Using the outputs from high-resolution global climate models (GCMs), such as the High Resolution Model Intercomparison Project (HighResMIP) of the Coupled Model Intercomparison Project Phase 6 (CMIP6), has become one of the most important tools in wind energy research. This study evaluated the reliability of the 22 GCMs available in the HighResMIP-PRIMAVERA project by simulating the wind energy climatology and variability over the Tibetan Plateau (TP) with reference to observations and investigated the differences in performance of the GCMs between high-resolution (HR) and low-resolution (LR) simulations. The results show that most models performed relatively well in simulating the probability distribution of the observed wind speed over the TP, but nearly half of the models generally underestimated the wind speed, whereas the others tended to overestimated the wind speed. Compared with the wind speed, the GCMs showed larger biases in reproducing the wind power density (WPD) and other wind energy resources, whereas the biases in multi-model ensembles were relatively smaller than those in most individual models. With respect to interannual variability, both the HR and LR models failed to capture interannual variations in WPD over the TP. Furthermore, more than half of the HR GCMs had a reduced bias relative to the corresponding LR GCMs, indicating the good performance of most HR models in simulating wind energy resources over the TP in terms of spatial pattern and temporal variability. However, the overall performance of HR GCMs varied among models, which suggests that solely improving the horizontal resolution is not sufficient to completely solve the uncertainties and deficiencies in the simulation of wind energy over complex terrain.

Response of drought to climate extremes in a semi-arid inland river basin in China

AbstractAgainst the backdrop of global warming, climate extremes and drought events have become more severe, especially in arid and semi-arid areas. This study forecasted the characteristics of climate extremes in the Xilin River Basin (a semi-arid inland river basin) of China for the period of 2021–2100 by employing a multi-model ensemble approach based on three climate Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, and SSP5-8.5) from the latest Coupled Model Intercomparison Project Phase 6 (CMIP6). Furthermore, a linear regression, a wavelet analysis, and the correlation analysis were conducted to explore the response of climate extremes to the Standardized Precipitation Evapotranspiration Index (SPEI) and Streamflow Drought Index (SDI), as well as their respective trends during the historical period from 1970 to 2020 and during the future period from 2021 to 2070. The results indicated that extreme high temperatures and extreme precipitation will further intensify under the higher forcing scenarios (SSP5-8.5>SSP2-4.5>SSP1-2.6) in the future. The SPEI trends under the SSP1-2.6, SSP2-4.5, and SSP5-8.5 scenarios were estimated as −0.003/a, −0.004/a, and −0.008/a, respectively, indicating a drier future climate. During the historical period (1970–2020), the SPEI and SDI trends were −0.003/a and −0.016/a, respectively, with significant cycles of 15 and 22 a, and abrupt changes occurring in 1995 and 1996, respectively. The next abrupt change in the SPEI was projected to occur in the 2040s. The SPEI had a significant positive correlation with both summer days (SU) and heavy precipitation days (R10mm), while the SDI was only significantly positively correlated with R10mm. Additionally, the SPEI and SDI exhibited a strong and consistent positive correlation at a cycle of 4–6 a, indicating a robust interdependence between the two indices. These findings have important implications for policy makers, enabling them to improve water resource management of inland river basins in arid and semi-arid areas under future climate uncertainty.

Great Gerbils (Rhombomys opimus) in Central Asia Are Spreading to Higher Latitudes and Altitudes.

The great gerbil (Rhombomys opimus) is a gregarious rodent in Central Asia and is one of the major pests found in desert forest and grassland areas. The distribution changes and migration routes of R. opimus in Central Asia under climate change remain unexplored. This study employed multi-model ensemble, correlation analysis, jackknife method, and minimum cumulative resistance (MCR) model to simulate the potential habitat of R. opimus under current and future (2030 and 2050) climate scenarios and estimate its possible migration routes. The results indicate that the ensemble model integrating Random Forest (RF), Gradient Boosting Machine (GBM), and Maximum Entropy Model (MaxEnt) performed best within the present climate context. The model predicted the potential distribution of R. opimus in Central Asia with an area under the curve (AUC) of 0.986 and a True Skill Statistic (TSS) of 0.899, demonstrating excellent statistical accuracy and spatial performance. Under future climate scenarios, northern Xinjiang and southeastern Kazakhstan will remain the core areas of R. opimus distribution. However, the optimal habitat region will expand relative to the current one. This expansion will increase with the rising CO2 emission levels and over time, potentially enlarging the suitable area by up to 39.49 × 104 km2. In terms of spatial distribution, the suitable habitat for R. opimus is shifting toward higher latitudes and elevations. For specific migration routes, R. opimus tends to favor paths through farmland and grassland. This study can provide guidance for managing and controlling R. opimus under future climate change scenarios.

Formally combining different lines of evidence in extreme-event attribution

Abstract. Event attribution methods are increasingly routinely used to assess the role of climate change in individual weather events. In order to draw robust conclusions about whether changes observed in the real world can be attributed to anthropogenic climate change, it is necessary to analyse trends in observations alongside those in climate models, where the factors driving changes in weather patterns are known. Here we present a quantitative statistical synthesis method, developed over 8 years of conducting rapid probabilistic event attribution studies, to combine quantitative attribution results from multi-model ensembles and other, qualitative, lines of evidence in a single framework to draw quantitative conclusions about the overarching role of human-induced climate change in individual weather events.

Groundwater stress in Europe—assessing uncertainties in future groundwater discharge alterations due to water abstractions and climate change

Groundwater sustains human well-being and ecosystems functioning. Many regions in Europe have experienced declining groundwater levels caused by decreasing groundwater recharge (GWR) or increasing groundwater abstractions (GWAs). These changes can lead to groundwater-related stress, threatening ecosystems and water supplies. Existing groundwater stress indicators estimate stress during a given period but do not address how stress changes or show the uncertainty of future stress. We propose a novel indicator of future groundwater stress (GWSI) due to changes in GWR and GWA and, thus, the alteration of long-term mean annual groundwater discharge (GWD). Groundwater stress is defined as any alteration in GWD since ecosystems are adapted to an equilibrium state. Focusing on decreasing GWD, which is generally more harmful than increasing GWD, we quantified the future GWSI in Europe by integrating scenarios of GWR and GWA in 2070–2099. GWR was evaluated using an ISIMIP2b multi-model ensemble of eight global hydrological models driven by the output of four global climate models under two greenhouse gas emission scenarios. GWA scenarios for irrigation, domestic and manufacturing sectors were combined with the GWR projections to generate an ensemble of GWSIs, simplified into three groundwater stress scenarios (high, intermediate, low). Projected GWSIs vary significantly among the scenarios. For the high-stress scenario, 58% of Europe’s land area is projected to experience a GWD decrease of at least 25% under RCP8.5 compared to 38% under RCP2.6, while the respective values are 26 and 1% for the intermediate-stress scenario. Groundwater demand management alone might not prevent GWD declines under the high-stress and intermediate scenarios, particularly under RCP8.5. Therefore, climate change mitigation might imperative for reducing the decline of GWD, especially in Eastern and Southeastern Europe, where changes in GWR are projected to be the primary cause of declining GWD (in the high abstraction scenario under RCP8.5). Under RCP2.6, reductions in GWAs by 25–75% might balance a GWD decline in parts of Spain and Italy where GWAs are high, even in the high-stress scenario. In line with the precautionary principle, we recommend adapting to the high-stress scenario to minimize harm to the beneficiaries of groundwater.

Accelerated North Atlantic surface warming reshapes the Atlantic Multidecadal Variability

Observed North Atlantic sea-surface temperatures are modulated by a recurrent alternation of anomalously warm and cold interdecadal phases known as Atlantic Multidecadal Variability. Here we use observations and a multi-model ensemble of climate simulations to demonstrate an ongoing acceleration of North Atlantic surface warming, which implies a smaller contribution of the Atlantic Multidecadal Variability to 21st century North Atlantic sea-surface temperature anomalies than previously thought. Future projections of the Atlantic Multidecadal Variability from realistic climate simulations are poorly constrained, yet a relaxation to a neutral phase by the mid-21st century emerges as the most probable evolution of the Atlantic Multidecadal Variability. In the simulations, the mitigating effects of a less likely upcoming cold phase of the Atlantic Multidecadal Variability are overwhelmed by fast North Atlantic surface warming, which is robustly projected to persist in upcoming decades independent of emission scenarios. Sustained North Atlantic surface warming is therefore expected to continue in the near future.

Investigating the Limitations of Multi‐Model Ensembling of Climate Model Outputs in Capturing Climate Extremes

ABSTRACTIn the context of climate change, the widespread practice of directly employing Multi‐Model Ensembles (MMEs) for projecting future climate extremes, without prior evaluation of MME performance in historical periods, remains underexplored. This research addresses this gap through a comprehensive analysis of ensemble means derived from CMIP6‐based models, including both simple and weighted averages of precipitation (SEMP and WEMP) and temperature (SEMT and WEMT) time series, as well as simple (SEME) and weighted (WEME) averages of extremes based on model‐by‐model analysis. The study evaluates the efficacy of MMEs in capturing mean annual values of ETCCDI indices over India for the period 1951–2014, utilising the IMD gridded data set as a reference. The results reveal that SEME and WEME consistently align closely with IMD data across various precipitation indices. At the same time, SEMP and WEMP consistently display underestimation biases ranging from 20% to 80% across all precipitation indices, except for CWD, where there is an overestimation bias. Moreover, SEMP and WEMP consistently underestimate CDD and overestimate CWD, indicating a systematic bias in these ensemble means, while WEME and SEME demonstrate satisfactory performance. SEMT and WEMT exhibit notable underestimation in temperature indices. In summary, adopting SEME and SEMT leads to a more robust assessment of precipitation and temperature extremes, respectively. These findings highlight the limitations of traditional MME methodologies in reproducing observed extreme precipitation events across various climatic zones in India, offering essential insights for refining climate models and improving the reliability of climate projections specific to the Indian subcontinent.

MEF-AlloSite: an accurate and robust Multimodel Ensemble Feature selection for the Allosteric Site identification model

A crucial mechanism for controlling the actions of proteins is allostery. Allosteric modulators have the potential to provide many benefits compared to orthosteric ligands, such as increased selectivity and saturability of their effect. The identification of new allosteric sites presents prospects for the creation of innovative medications and enhances our comprehension of fundamental biological mechanisms. Allosteric sites are increasingly found in different protein families through various techniques, such as machine learning applications, which opens up possibilities for creating completely novel medications with a diverse variety of chemical structures. Machine learning methods, such as PASSer, exhibit limited efficacy in accurately finding allosteric binding sites when relying solely on 3D structural information.Scientific ContributionPrior to conducting feature selection for allosteric binding site identification, integration of supporting amino-acid–based information to 3D structural knowledge is advantageous. This approach can enhance performance by ensuring accuracy and robustness. Therefore, we have developed an accurate and robust model called Multimodel Ensemble Feature Selection for Allosteric Site Identification (MEF-AlloSite) after collecting 9460 relevant and diverse features from the literature to characterise pockets. The model employs an accurate and robust multimodal feature selection technique for the small training set size of only 90 proteins to improve predictive performance. This state-of-the-art technique increased the performance in allosteric binding site identification by selecting promising features from 9460 features. Also, the relationship between selected features and allosteric binding sites enlightened the understanding of complex allostery for proteins by analysing selected features. MEF-AlloSite and state-of-the-art allosteric site identification methods such as PASSer2.0 and PASSerRank have been tested on three test cases 51 times with a different split of the training set. The Student’s t test and Cohen’s D value have been used to evaluate the average precision and ROC AUC score distribution. On three test cases, most of the p-values (<0.05\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$< 0.05$$\\end{document}) and the majority of Cohen’s D values (>0.5\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$> 0.5$$\\end{document}) showed that MEF-AlloSite’s 1–6% higher mean of average precision and ROC AUC than state-of-the-art allosteric site identification methods are statistically significant.

Greenland summer blocking characteristics: an evaluation of a high-resolution multi-model ensemble

AbstractAtmospheric blocking is a phenomenon that can lead to extreme weather events over a large region, yet its causes are not fully understood. Global climate models show limitations in representing Northern Hemisphere blocking, especially its frequency, and decadal variability in Greenland blocking in summer in the recent decades. In this study we evaluate the ability of high-resolution (HighResMIP) Earth System Models (ESMs) to simulate summer blocking over the Greenland area, using different but complementary methods to describe the characteristics of blocking. We find that the HighResMIP ensemble can reproduce the spatial pattern of Greenland blocking events, albeit with systematic biases, and capture the relative frequencies of the main blocking patterns: namely the wave breaking structure, North Atlantic ridge, and omega-type blocking. However, the HighResMIP ensemble fails to simulate the observed temporal variations of Greenland blocking index (GB2) and the extremely high values of daily GB2 observed in recent decades. In addition, we do not find clearly superior representation of blocking features from higher-resolution in HighResMIP models compared with lower-resolution models. We also find large sea surface temperature (SST) biases over the North Atlantic and seas surrounding Greenland, and biases in moisture transport over the North Atlantic toward Greenland, especially over the western flank of blocking areas, which might together contribute to model biases in the representation of blocking magnitude.

Uncertainty-informed selection of CMIP6 Earth system model subsets for use in multisectoral and impact models

Abstract. Earth system models (ESMs) and general circulation models (GCMs) are heavily used to provide inputs to sectoral impact and multisector dynamic models, which include representations of energy, water, land, economics, and their interactions. Therefore, representing the full range of model uncertainty, scenario uncertainty, and interannual variability that ensembles of these models capture is critical to the exploration of the future co-evolution of the integrated human–Earth system. The pre-eminent source of these ensembles has been the Coupled Model Intercomparison Project (CMIP). With more modeling centers participating in each new CMIP phase, the size of the model archive is rapidly increasing, which can be intractable for impact modelers to effectively utilize due to computational constraints and the challenges of analyzing large datasets. In this work, we present a method to select a subset of the latest phase, CMIP6, featuring models for use as inputs to a sectoral impact or multisector dynamics models, while prioritizing preservation of the range of model uncertainty, scenario uncertainty, and interannual variability in the full CMIP6 ensemble results. This method is intended to help impact modelers select climate information from the CMIP archive efficiently for use in downstream models that require global coverage of climate information. This is particularly critical for large-ensemble experiments of multisector dynamic models that may be varying additional features beyond climate inputs in a factorial design, thus putting constraints on the number of climate simulations that can be used. We focus on temperature and precipitation outputs of CMIP6 models, as these are two of the most used variables among impact models, and many other key input variables for impacts are at least correlated with one or both of temperature and precipitation (e.g., relative humidity). Besides preserving the multi-model ensemble variance characteristics, we prioritize selecting CMIP6 models in the subset that preserve the very likely distribution of equilibrium climate sensitivity values as assessed by the latest Intergovernmental Panel on Climate Change (IPCC) report. This approach could be applied to other output variables of climate models and, possibly when combined with emulators, offers a flexible framework for designing more efficient experiments on human-relevant climate impacts. It can also provide greater insight into the properties of existing CMIP6 models.

Assessment of CMIP6 models and multi-model averaging for temperature and precipitation over Iran

In this study, the performances of 40 Coupled Model Intercomparison Project Phase 6 are evaluated against observational data at synoptic stations in Iran using various evaluation criteria. The results reveal diverse model accuracy across different climate conditions and criteria, emphasizing particularly notable disparities in the nonstationarity R criterion compared to others. Although according to the ranking of the raw and bias-corrected outputs of CMIP6 GCMs for Iran, the NorESM2-MM, AWI-ESM-1-1-LR, and MPI-ESM1-2-LR models are consistently among the top six ranked models for precipitation in both raw and corrected outputs. For temperature, MPI-ESM1-2-LR, TaiESM1, INM-CM4-8, and IITM-ESM are consistently among the top six models for both the raw and bias-corrected outputs of CMIP6 GCMs. The Bias correction methods, including quantile mapping and linear scaling, integrated with Bayesian model averaging, were applied. While quantile mapping demonstrates superior performance and less disparity than linear scaling, it proves ineffective for correcting biases at stations with bias nonstationarity over time. The RMSE for monthly precipitation ranges from almost 0 to 200 mm, with a large RMSE value related to the high precipitation stations, and the monthly temperature exhibits a range of 0 to 4 °C. The use of a multi-model ensemble improves accuracy compared to individual models, resulting in a reduction in the differences between the minimum and maximum RMSE values from 178.6 to 91.0. Additionally, the range for mean absolute error decreases from 126.9 to 93.3, and the difference in the correlation coefficient narrows from 0.9 to 0.42. Averaging models after bias correction prevents significant fluctuations while maintaining higher accuracy, in contrast to the second method, which involves bias-correcting models after averaging.

Atmospheric Precursors of Skillful SST Prediction in the Northeast Pacific

Abstract Forecasts of sea surface temperature anomalies (SSTAs) provide essential information to stakeholders of marine resources in coastal ecosystems, such as the California Current Large Marine Ecosystem (CCLME), at management-relevant monthly-to-annual time scales. Diagnosing dynamical sources of predictability and the mechanisms differentiating skill among forecasts is required for verification and improvement in operational forecasting systems. Using retrospective forecasts (1982–2020) from a four-member subset of the North American Multi-Model Ensemble (NMME), we evaluate the conditional skill of SSTA forecasts in the CCLME at monthly resolution for lead times up to 10.5 months. Forecasts from ensemble members with relatively small SSTA errors at shorter lead times retain higher skill at longer lead times, with the most substantial and long-lasting increases for forecasts initialized in the fall and early spring. The “best” low-error SSTA forecasts are characterized by increased skill in the prediction of North Pacific atmospheric circulation [sea level pressure (SLP) and 200-hPa geopotential height] the month prior to the evaluation of SSTA errors in the CCLME and exhibit more realistic progressions of anomalous SLP. The Pacific meridional mode (PMM) emerges as a diagnostic of skillful North Pacific atmosphere–ocean coupling, as forecasts that correctly simulate the PMM and its associated SLP variability increase the SSTA prediction skill in the CCLME in the fall through spring. Predictable coupled ocean–atmosphere modes provide a target for enhancing predictability with early detection of the onset of a deterministic progression emerging from stochastic atmospheric variability. Significance Statement Global forecast systems provide near-term climate predictions that inform the management of marine resources, such as those of the California Current Large Marine Ecosystem. In this study, we probe the processes which lead forecasts to succeed or fail at predicting sea surface temperatures in the California Current at seasonal time scales among retrospective forecasts from the North American Multimodel Ensemble. We demonstrate that forecasts which best simulate sea surface temperatures at the earliest lead times sustain advantages in forecast skill and find that correctly simulating extratropical atmospheric circulation increases the predictive skill of sea surface temperatures in the northeast Pacific in the following lead times. Our results offer North Pacific atmospheric circulation as a target for forecast model improvement that would additionally enhance ocean forecasts.

El Niño-like tropical pacific ocean cooling pattern during the last glacial maximum

Many state-of-the-art climate models are unable to reproduce the observed 20th century surface warming pattern in the tropical Pacific Ocean, casting doubt on the robustness of future projections. Here, we examine past changes in the tropical Pacific upper ocean spatial pattern using paleoclimate reconstructions from proxies and simulations from a multi-model ensemble. The proxy results demonstrate that during the Last Glacial Maximum, when atmospheric carbon dioxide concentrations were lower than present, temperatures in the western tropical Pacific Ocean decreased more than in the east, leading to an El Niño-like cooling pattern. This result contrasts with the zonally uniform cooling pattern observed in model simulations, highlighting the common issue of models overestimating the sensitivity of eastern tropical Pacific sea surface temperatures to greenhouse forcing. Our proxy results imply that the western Pacific may warm more than the east in response to future climate change, producing a La Niña-like surface warming pattern.

Significant winter Atlantic Niño effect on ENSO and its future projection

The Atlantic Niño, a primary climatic variability mode in the equatorial Atlantic Ocean, exhibits pronounced variability not only in boreal summer but also in winter. However, the role of the winter Atlantic Niño in trans-basin interactions remains underexplored compared to its summer counterpart. Through analysis of observational reanalysis data since the mid-twentieth century, here we found that the winter Atlantic Niño significantly influences the development of El Niño–Southern Oscillation (ENSO), surpassing the impact of summer Atlantic Niño, with a longer lead time. This effect is reasonably captured in the CMIP6 Historical simulations from a multi-model ensemble perspective. Further analysis of the global warming scenario projects that the influence of winter Atlantic Niño on ENSO will persist into the future, in contrast to the reduced impact of summer Atlantic Niño. Therefore, these findings underscore the importance of further investigating the winter Atlantic Niño to gain a comprehensive understanding of trans-basin interactions and their future changes.