The study focuses on HadGEM3-GC3.1-LL, the low-resolution configuration of the UK’s CMIP6 climate model, which simulates the atmosphere, ocean, land, and sea ice as a fully coupled system. The atmosphere component, the Met Office Unified Model, runs on a grid of roughly 135 kilometers, which means that much of Antarctica’s rugged terrain is too small to be resolved directly. Instead, the model relies on a parameterization of orographic gravity wave drag, or OGWD, to represent the forces that sub-grid-scale mountains impose on the flow. This parameterization splits the total drag into two components: flow-blocking drag, which arises when air is too stable to climb over an obstacle and is deflected around it, and gravity wave drag, which arises when air does surmount the terrain and generates waves that break aloft and deposit momentum.

Which of these two regimes dominates depends on a delicate balance captured by the non-dimensional parameter Nh/U, where U is wind speed, N is atmospheric stability, and h is obstacle height. When Nh/U greatly exceeds one, blocking prevails; when it is near or below one, the flow passes over the mountain. In the Unified Model, the effective sub-grid mountain height is set by a tuning constant called n-sigma, with a default value of 2.5 that is among the least constrained parameters in the entire model. The researchers exploited this lever to run two sets of sensitivity experiments over Antarctica: one, called F_BLOCK, in which n-sigma was doubled to 5, making virtual mountains effectively taller and forcing the flow into a blocking regime, and another, called F_OVER, in which n-sigma was reduced to 1.5, nearly halving the barrier height and letting winds sweep across the terrain while generating strong gravity waves. Crucially, the modifications were applied only over Antarctica, so any circulation change in the simulations could be attributed unambiguously to the southernmost continent.

The results were striking. In the baseline historical simulation, pan-Antarctic sea ice area begins a sustained decline around the early 1980s, driven largely by the Weddell Sea, one of the continent’s main sea ice factories. In the flow-blocking experiments, that decline is delayed by roughly ten to twenty years, with the ensemble mean holding steady until just before 1995. In the flow-over experiment, the decline begins about ten years earlier still, around 1972. Because the delayed decline appears consistently across all four members of the flow-blocking ensemble, the team could rule out natural variability as the explanation: the differences are a systematic response of the coupled climate system to how orographic drag is represented.

The immediate physical mechanism operates through the atmosphere’s vertical structure. In the flow-over configuration, larger-amplitude gravity waves break in the lower and mid-troposphere, decelerating the background westerly flow and weakening the tropospheric polar vortex. In the flow-blocking configuration, weaker wave amplitudes mean less drag and stronger zonal winds aloft. A stronger polar vortex insulates the Antarctic continent, exchanging less heat with the mid-latitudes, so surface air temperatures drop in F_BLOCK and rise in F_OVER. Colder conditions strengthen the katabatic winds, the rivers of cold, dense air that drain from the high interior toward the coast, which push sea ice offshore in the Weddell Sea and promote the formation of fresh ice. Within the first one to two years of simulation, sea ice area rises by about one million square kilometers in F_BLOCK and falls by a similar amount in F_OVER relative to the historical run.

The long-term story, however, unfolds in the Weddell Sea, and it hinges on the Antarctic Peninsula acting as a formidable barrier to the prevailing westerlies. In the flow-blocking experiment, the incoming flow cannot clear the ridge, so westerly winds both upstream over the Amundsen and Bellingshausen Seas and downstream over the Weddell Sea are weakened. Mean sea level pressure anomalies generated by the blocking propagate vast distances around the continent, echoing earlier work by Sandu and colleagues showing that orographic blocking leaves fingerprints far beyond the mountains that create it. June-to-November westerly winds at 850 hectopascals across the Weddell Sea are reduced by about half a meter per second, roughly sixteen percent, in F_BLOCK, and strengthened by a comparable margin in F_OVER.

Those wind changes translate directly into the ocean. The total surface stress exerted on the Weddell Sea, the combined push of atmosphere and sea ice on the water, is much weaker under flow-blocking conditions and much stronger under flow-over conditions, partly because stress scales with the square of wind speed and partly because thinner or reduced ice allows more efficient momentum transfer. Stronger winds spin up the cyclonic Weddell gyre, whose depth-integrated circulation, measured by the barotropic streamfunction, strengthens steadily through the historical run. A stronger gyre makes the surface flow more divergent, and by continuity, relatively warm, salty water is pumped upward from depth, with Ekman pumping in the historical simulation rising from roughly 55 meters per year in 1980 to about 85 meters per year by the end of the run.

Near Maud Rise, at about 66 degrees south and 3 degrees east, weak stratification and favorable topography set the stage for open-ocean deep convection, the vertical mixing that draws heat from the deep ocean to the surface. In the model, these events are diagnosed when the ocean mixed layer exceeds 2000 meters in depth, and polynya-scale events are flagged when the convecting area surpasses 80,000 square kilometers, the size of the real Weddell Sea polynya of 2016 and 2017. The correspondence with sea ice is dramatic: every cluster of deep convection years coincides with a sharp drop in sea ice area. In the historical simulation, convection recurs almost every year from 1989 onward, and in F_OVER it appears even earlier and more often. In F_BLOCK, by contrast, deep convection is nearly absent, with only a handful of small events across sixty-five years, and the ice survives largely intact.

The team then turned the correlation into a diagnostic tool. Plotting annual sea ice area against ocean surface stress reveals a strong negative relationship, with a correlation coefficient of minus 0.75, and the historical simulation migrates over time from a low-stress, high-ice state toward a high-stress, low-ice state that mirrors the flow-over experiment. Using the observational range of sea ice area from the HadISST1 dataset as an emergent constraint, the researchers estimate that plausible surface stress values lie between 0.042 and 0.066 newtons per square meter, a range that the flow-blocking simulation satisfies best. That result suggests the default model may underestimate atmospheric blocking by the Peninsula, and that a calibrated tuning of the sub-grid orography there would be advisable.

The broader implications reach beyond model tuning. The study shows that a more positive Southern Annular Mode does not necessarily mean stronger winds over the Weddell Sea, because the Peninsula’s barrier effect can override the large-scale annular signal. It also suggests that the springtime strengthening of westerlies over the Weddell Sea seen in the historical run, which is consistent with ERA5 reanalysis, may partly reflect a long-term shift in the Peninsula’s flow regime toward weaker blocking, driven by declining static stability, rather than greenhouse and ozone forcing alone. And because the spurious deep convection is fueled by wind-driven Ekman pumping, a large-scale process the model resolves explicitly, simply tweaking vertical mixing parameterizations will not fix it. The authors point instead to revising how ocean surface stress is represented, potentially through variable drag coefficients that respond to the evolving state of the ice and atmosphere, an approach already under development for the next version of HadGEM3. In the end, the message is humbling and exhilarating in equal measure: a single tuning constant, buried in the treatment of invisible mountains, can decide whether a flagship climate model loses its Antarctic sea ice a decade early, a decade late, or not at all.

Subject of Research: The influence of orographic gravity wave drag and flow regimes over the Antarctic Peninsula on Weddell Sea sea ice and deep convection in a coupled climate model

Article Title: Southern Hemisphere sea ice response to different flow regimes over the Antarctic Peninsula

Article References: Guarino, M.-V., Ridley, J. K., Farneti, R., Kucharski, F., & Tompkins, A. M. (2026). Southern Hemisphere sea ice response to different flow regimes over the Antarctic Peninsula. Climate Dynamics, 64(9), Article 410. https://doi.org/10.1007/s00382-026-08350-6

Image Credits: AI Generated

DOI: 10.1007/s00382-026-08350-6

Keywords: Antarctic sea ice, Antarctic Peninsula, orographic gravity wave drag, Weddell Sea, deep convection, polynyas, HadGEM3, CMIP6, wind stress, Ekman pumping, Southern Annular Mode, climate modeling

The research was carried out by Ayşe Özlem Tursun of Malatya Turgut Özal University, who combined experimental plant physiology with species distribution modeling, a pairing that remains surprisingly uncommon in climate-impact studies. Conventional distribution models treat a species as a passive follower of climate envelopes, predicting where conditions might be tolerable but saying little about whether the organism will perform better or worse once it arrives. By integrating experimentally derived physiological responses, Tursun’s framework attempts to capture how photosynthesis, growth and secondary metabolite production in lemon balm shift under the novel thermal and atmospheric regimes of the coming decades, and then translates those responses into a scenario-based sensitivity analysis of range change.

The modeling backbone of the study is the Maximum Entropy algorithm, or MaxEnt, one of the most widely used tools for mapping potential species distributions from presence-only records. Tursun assembled 117 spatially independent occurrence records for lemon balm, drawing on a formally documented download from the Global Biodiversity Information Facility, and paired them with bioclimatic variables from the WorldClim version 2.1 database. To keep the projections honest, model complexity was tuned using spatial block cross-validation, a technique that separates training and testing data geographically so the model cannot simply memorize local conditions. The resulting statistics were strong: a training area under the curve of 0.825, a test AUC of 0.878, and a True Skill Statistic of 0.620, values that indicate the model discriminates suitable from unsuitable habitat considerably better than chance.

Two climatic variables emerged as the decisive levers shaping where lemon balm can live. The minimum temperature of the coldest month, known among modelers as Bio6, contributed most by variable importance, while the annual mean temperature, Bio1, ranked highest by permutation importance. In plain terms, lemon balm’s Turkish distribution is governed above all by how cold winters get and how warm the year runs on average. That makes intuitive sense for a herb of Mediterranean and Irano-Turanian affinity, and it also explains the geography of the projected expansion: as winters soften and annual temperatures climb, the cold walls that currently fence the species into parts of western and southern Türkiye begin to dissolve.

The spatial pattern of the forecast is striking. Under all scenarios and both future time horizons, 2050 and 2070, the model projects progressive habitat expansion rather than contraction, with newly suitable territory emerging along the Black Sea coastline and across the northern Central Anatolian plateau. The expansion is pronouncedly poleward and upward in elevation, a signature consistent with the general expectation that species will track their thermal niches toward higher latitudes and altitudes. For a country whose medicinal and aromatic plant sector is economically and culturally significant, the map carries practical weight: it identifies where lemon balm cultivation might plausibly be established in the coming decades and where current strongholds could face shifting competitive and climatic conditions.

What lifts the study beyond a standard MaxEnt exercise is its treatment of carbon dioxide. Rising atmospheric CO2 does not merely warm the planet; it also fertilizes photosynthesis in many plant species, a physiological effect that pure correlative models routinely ignore. Drawing on earlier experimental work in which lemon balm was grown under different temperature and CO2 concentrations, Tursun ran a scenario-based sensitivity analysis to estimate how much of the projected expansion could be attributed to this fertilization effect. The answer was roughly 25.1 percentage points of the total expansion, a substantial share, though the analysis concluded that rising temperatures remain the dominant driver of the range shift. The result is a rare quantitative partitioning of the climatic and physiological components of a forecast range change.

Robustness was addressed through an ensemble assessment across three CMIP6 general circulation models. The coefficient of variation across the ensemble came out at just 5.6 percent, a low figure indicating that the expansion signal is not an artifact of one particular climate model’s quirks. The projections were run under two Shared Socioeconomic Pathways, SSP2-4.5 representing a moderate emissions trajectory and SSP5-8.5 representing a fossil-fuel-intensive future, and the direction of change held under both, with magnitude scaling with emissions. In a field where projections often swing wildly between climate models, that degree of agreement is notable and lends credibility to the headline number of a doubling of suitable area by mid-to-late century.

Yet the study’s own framing makes clear that a larger map is not an unalloyed good. Climate change poses what the author calls dual challenges to medicinal and aromatic plants: the contraction of suitable habitats on one hand, and the alteration of secondary metabolite biosynthesis under new thermal and atmospheric conditions on the other. For lemon balm, the commercial and medicinal value lies in its essential oils and phenolic compounds, and the scientific literature has long documented that environmental stressors, temperature and CO2 among them, can reshape the chemistry of such plants. A habitat that is climatically suitable may therefore produce raw material of different quality and potency than today’s harvests, a dimension that distribution maps alone cannot capture.

The methodological choices also deserve attention from practitioners. The occurrence data were spatially thinned to reduce sampling bias, the model was evaluated with metrics appropriate to presence-only modeling, and the threshold selection for converting continuous suitability scores into binary habitat maps followed established procedures for such data. The R script underlying the analysis was released as supplementary material, and the occurrence dataset carries a persistent citable DOI, making the workflow reproducible end to end. In an era when species distribution models increasingly inform conservation planning and agricultural policy, that transparency matters as much as the headline projection.

For Türkiye, the practical message is one of climate-smart agricultural planning. The identification of the Black Sea coast and northern Central Anatolia as zones of emerging suitability offers a forward-looking guide for growers, cooperatives and policymakers in the medicinal plant sector, suggesting where cultivation trials, land-use investments and conservation measures might be directed before the climate arrives. At the same time, the study is a reminder that even the winners of climate change live in a transformed world: a lemon balm field in 2070 will grow under different skies, breathe different air, and likely synthesize a different bouquet of compounds than its predecessors. The expansion of its habitat, dramatic as it is, is only the first chapter of the story.

Subject of Research: Projected habitat expansion of lemon balm in Türkiye under CMIP6 climate scenarios using integrated physiological and distribution modeling

Article Title: Integrating experimental physiological responses with species distribution modeling to forecast the future habitat expansion of Melissa officinalis L. in Türkiye under CMIP6 scenarios

Article References: Integrating experimental physiological responses with species distribution modeling to forecast the future habitat expansion of Melissa officinalis L. in Türkiye under CMIP6 scenarios. (n.d.). https://doi.org/10.1007/s00704-026-06613-0

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06613-0

Keywords: lemon balm, Melissa officinalis, species distribution modeling, MaxEnt, CMIP6, climate change, Türkiye, CO2 fertilization, medicinal plants, habitat expansion, bioclimatic variables, SSP scenarios

A new study published in the journal Climate Dynamics offers a fresh path forward. A team of researchers led by Shu Gui, Dayong Wen, and Jie Cao of Yunnan University, together with colleagues including Liang Duan, has shown that the atmospheric heat source over the Southeast Asian low-latitude highlands, a region abbreviated as SEALLH, can serve as a powerful early warning signal for the timing of the Meiyu onset. By combining this physical precursor with machine learning models trained in an unconventional way, the team achieved hindcast skill that outperformed more direct approaches, and their results highlight a subtle but crucial lesson for anyone applying artificial intelligence to climate prediction: a model that predicts well but predicts for the wrong physical reasons is a model that cannot be fully trusted.

The Southeast Asian low-latitude highlands encompass the elevated terrain of the Indochinese Peninsula and adjacent regions, sitting south and southwest of the better-known Tibetan Plateau. Like the plateau itself, these highlands act as an elevated heat source in the boreal spring and early summer. As the sun climbs northward, the land surface warms, convection intensifies, and enormous quantities of latent and sensible heat are released into the middle and upper troposphere. This diabatic heating, the net heating of an air parcel by processes other than direct radiation, including condensation of water vapor in deep convective clouds and turbulent transfer from the surface, helps drive the large-scale circulation transitions that eventually pull the East Asian summer monsoon northward and set up the Meiyu front over the Yangtze and Huaihe river valleys.

Previous work by members of the same team had already established that anomalies in the atmospheric heat source over these highlands in late spring are closely tied to year-to-year swings in the Meiyu onset date. When the heating is anomalously strong or weak, it alters the position and intensity of the western North Pacific subtropical high, the great anticyclonic circulation that steers moisture-laden air toward eastern China. It also modulates moisture convergence over East Asia, effectively setting the stage, months in advance, for whether the Meiyu rain band will arrive early or late. What remained unclear was whether this promising physical linkage could actually be converted into quantitative prediction skill, and if so, how best to do it.

The researchers turned to the interannual increment approach, a strategy that has gained traction in seasonal climate prediction. Rather than training models to predict the onset date itself, which varies within a fairly narrow window and carries strong decadal trends that can mislead statistical learning algorithms, the method predicts the year-to-year increment, the difference between the onset date in one year and the previous year. This differencing strips away slowly varying background signals and focuses the learning problem on the genuine interannual variability that forecasters care about, often yielding sharper and more robust prediction models.

To build their prediction system, the team trained five categories of machine learning models using historical simulations from the Coupled Model Intercomparison Project Phase 6, the international archive of state-of-the-art climate model experiments. Training on climate model output rather than observations alone gives the algorithms a much larger sample of physically consistent atmospheric behavior to learn from, helping to guard against overfitting to the short observational record. The trained models were then validated against ERA5, the high-resolution reanalysis product from the Copernicus Climate Change Service that blends observations with a numerical weather prediction model to provide a best estimate of the historical atmospheric state. The predictors were the diabatic heating anomalies over the Southeast Asian low-latitude highlands, and the target was the Meiyu onset date increment.

The team designed two distinct hindcasting strategies and pitted them against each other. In the first group, models directly predicted the interannual increment of the Meiyu onset date in a single step. In the second group, the models instead predicted the onset dates of two consecutive years separately, and the increment was then derived by taking the difference between the two hindcasts. This seemingly small architectural difference turned out to matter enormously, and the reason lies in the physics that the models implicitly learned.

When the researchers examined what the best-performing models were actually doing, they found that the top models in both groups better captured the real-world chain of cause and effect linking the highland heating to large-scale circulation anomalies, moisture convergence, and ultimately the onset date, compared with the bottom-performing models. But the top models in the direct-prediction group harbored a hidden flaw: their predicted 500-hectopascal zonal wind anomalies, the east-west winds in the mid-troposphere that are central to the Meiyu circulation, were physically inconsistent with the observed relationships. In other words, these models could reproduce the onset date reasonably well while violating the very circulation physics that makes the prediction meaningful. Such inconsistency is a red flag for any operational forecast system, because a model that gets the right answer through the wrong mechanism is unlikely to remain reliable when conditions shift outside its training experience.

Remarkably, the two-step strategy largely corrected this defect. By hindcasting two consecutive onset dates and differencing them, the second group of models produced predictions whose associated circulation anomalies were far more physically coherent, and this coherence translated directly into skill. The top-performing models in the second group achieved a temporal correlation coefficient of 0.62 between hindcast and observed onset date increments, compared with 0.52 for the best models in the direct group. In the world of seasonal climate prediction, where correlation coefficients above 0.5 are often considered useful, this improvement is substantial. The study also found that the Extra Trees model, an ensemble method that builds many randomized decision trees and averages their outputs, delivered robust performance not only for predicting the current year’s onset but also for hindcasts made one year ahead, a lead time at which most dynamical seasonal models offer little guidance for the East Asian monsoon.

The implications extend well beyond the Meiyu. The work demonstrates that diabatic heating over the Southeast Asian low-latitude highlands deserves a place among the standard precursors used for East Asian summer monsoon prediction, complementing the traditional emphasis on sea surface temperature anomalies in the tropical Indian and Pacific oceans. Just as importantly, it offers a template for how machine learning should be deployed in climate science. Rather than judging models solely by their prediction scores, the researchers evaluated whether the models preserved known physical linkages among heating, circulation, and rainfall, and they showed that enforcing this physical consistency through clever problem design can itself be the key to better forecasts. As machine learning floods into weather and climate prediction, from nowcasting to decadal projection, the lesson from the Yunnan team is clear: the most skillful artificial intelligence is the kind that understands, or at least respects, the physics of the atmosphere it is trying to predict. For the farmers, dam operators, and millions of residents of the Yangtze valley, that combination of highland heat and disciplined algorithms may soon mean a more reliable heads-up before the rains arrive.

Subject of Research: Seasonal prediction of the Meiyu onset date using diabatic heating anomalies and machine learning

Article Title: Skillful prediction of Meiyu onset using diabatic heating over Southeast Asian low-latitude highlands

Article References: Gui, S., Wen, D., Cao, J., Dong, Z., Cai, L., Yang, R., & Duan, L. (2026). Skillful prediction of Meiyu onset using diabatic heating over Southeast Asian low-latitude highlands. Climate Dynamics, 64(10), Article 412. https://doi.org/10.1007/s00382-026-08366-y

Image Credits: AI Generated

DOI: 10.1007/s00382-026-08366-y

Keywords: Meiyu onset, diabatic heating, Southeast Asian highlands, machine learning, East Asian monsoon, Climate Dynamics, interannual increments, CMIP6, ERA5, Extra Trees, seasonal prediction, Yangtze River valley

The research team, led by corresponding author Hao Liang of Hohai University together with Xiaoyue Wu, Ruidong Chen and Songjuan Gao, assembled one of the most comprehensive field datasets ever compiled for this crop. The analysis drew on 572 individual biomass measurements collected at 111 monitoring sites spread across 13 provinces of southern China. Rather than relying on simple correlations, the researchers combined a Random Forest machine learning model with SHAP, an interpretable artificial intelligence technique that reveals how much each input variable contributes to a prediction and in which direction. This pairing allowed the team to move beyond black-box predictions and identify the specific climatic, geographic and soil conditions under which the crop thrives or falters.

The baseline picture is striking. The average dry biomass of Chinese milk vetch across the surveyed sites was 3.23 metric tons per hectare, a figure with direct agronomic consequences because biomass determines how much biologically fixed nitrogen and organic carbon is returned to the paddy soil before the next rice crop. The highest biomass was concentrated in the middle and lower reaches of the Yangtze River, particularly in Hunan, Hubei and Jiangxi, where mild, moist winters create near-ideal growing conditions. Lower biomass values appeared in parts of southern and southwestern China, hinting that the crop’s productivity is far from uniform across its cultivated range.

The machine learning model explained 68 percent of the observed spatial variation in biomass, a substantial share for a field-scale ecological dataset. When the contributions of different variable groups were separated, climatic factors emerged as the dominant force, accounting for 40.5 percent of the explained variation. Geographic factors contributed 31.7 percent and soil properties 27.8 percent. In other words, while local conditions and soil management matter, the weather that a milk vetch crop experiences during its winter growing season is the single most important determinant of how much nitrogen and carbon it will ultimately deliver to the rice system.

Perhaps the most consequential finding of the study is that these climatic effects are strongly nonlinear. Biomass did not simply rise or fall with temperature and rainfall; instead, the analysis uncovered clear thresholds. Growing-season precipitation between approximately 533 and 877 millimeters was associated with favorable biomass accumulation, while rainfall below or above that window was linked to reduced growth, reflecting the twin hazards of winter drought and waterlogging in paddy fields. Mean growing-season temperatures of roughly 10.7 to 13.7 degrees Celsius formed a broad thermal buffer within which the crop performed well. Above 13.7 degrees Celsius, however, the relationship shifted, with warmer conditions increasingly associated with heat stress and declining biomass.

These thresholds matter because they can be tested against the future. The team coupled its biomass model with projections from three CMIP6 climate models run under four Shared Socioeconomic Pathway scenarios, the standard framework used in international climate assessments to explore futures ranging from low to high greenhouse gas emissions. The result was a spatially explicit forecast of how Chinese milk vetch productivity might evolve through the end of the century, with projections extending to 2098.

Across southern China as a whole, the projected decline in milk vetch biomass was moderate, on the order of roughly 2 to 4 percent by 2098. But the aggregate number conceals a deeply uneven regional picture. The Huang Huai Hai single-cropping rice region was projected to suffer some of the largest losses, with biomass reductions reaching about 13 to 14 percent under higher-emission scenarios. In sharp contrast, the middle and lower Yangtze River double-cropping region, already the crop’s productivity heartland, remained comparatively stable and could even see biomass increases of approximately 1.9 to 5.9 percent under some scenarios. The same climate change that stresses the crop at the northern edge of its range may, within limits, extend favorable conditions in its core zone.

The practical implication, the authors argue, is that a single management strategy will not work everywhere. In regions facing the steepest projected losses, adaptation measures become urgent. The study proposes region-specific approaches, including adjusting sowing dates so that the growing season avoids the most stressful temperature and moisture conditions, developing stress-tolerant milk vetch varieties for vulnerable areas, conserving soil moisture through mulching and water management, and optimizing the integration of the green manure with rice straw return and nitrogen fertilization. Each of these levers interacts with the thresholds identified by the model, giving agronomists a quantitative basis for deciding where and how to intervene.

As corresponding author Hao Liang emphasized, Chinese milk vetch is more than a winter cover crop, because its biomass directly determines how much biologically fixed nitrogen and organic carbon can be returned to rice fields. The study’s results show that climate does not affect this crop in a simple linear way, and that the clear temperature and precipitation ranges within which milk vetch performs best can guide more precise regional management under a changing climate. That framing turns what might have been a purely descriptive mapping exercise into a decision-support tool for one of China’s most important low-input rice systems.

Beyond its immediate agronomic value, the work delivers a set of field-based benchmark data that could support crop modeling and remote sensing studies aimed at improving green manure management across southern China. The 572 measurements and the quantified climate thresholds provide calibration points for simulation models, and the spatial patterns documented by the team offer ground truth for satellite-based estimates of winter cover crop biomass. As climate pressures intensify through the coming decades, the study suggests that the future of Chinese milk vetch will be decided region by region, at the precise intersection of temperature, rainfall and management that the new analysis has now made visible.

Subject of Research: Climate-driven spatial variation and future projections of Chinese milk vetch biomass in southern China's rice paddies

Article Title: Climate change could reshape the future of Chinese milk vetch in southern rice paddies

Article References: Climate change could reshape the future of Chinese milk vetch in southern rice paddies. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: Chinese milk vetch, green manure, rice paddies, climate change, machine learning, Random Forest, SHAP, CMIP6, Shared Socioeconomic Pathways, biomass thresholds, nitrogen fixation, southern China

The observational record is striking. Drawing on measurements from 150 meteorological stations operated by the National Meteorological Information Center of the China Meteorological Administration, the team found that annual mean near-surface wind speed over the Tibetan Plateau declined at a rate of 0.159 meters per second per decade between 1979 and 2014. The decline was not spread evenly through the year. Spring showed the steepest weakening, with wind speeds falling by 0.235 meters per second per decade, a seasonal signal that matters because spring winds play a central role in the plateau’s exchange of heat and moisture with the surrounding atmosphere and in shaping the Asian monsoon system downstream.

To test whether climate models could reproduce this stilling, the researchers turned to the Coupled Model Intercomparison Project Phase 6, specifically its High-Resolution Model Intercomparison Project, known as HighResMIP. This experiment was designed in part to answer a long-standing question in climate science: would simply cranking up the horizontal resolution of global models improve their simulations? The team assembled 22 historical simulations from 11 model pairs, each pair consisting of a higher-resolution and a lower-resolution version of the same model, and compared them against the station observations over the 1979 to 2014 period.

The verdict on resolution was sobering. Both the higher-resolution and lower-resolution ensemble means captured the broad spatial pattern of wind speed across the plateau, but both badly underestimated the magnitude of the decline. The higher-resolution ensemble produced an annual trend of only 0.024 meters per second per decade, while the lower-resolution ensemble managed just 0.020 meters per second per decade. In other words, the models reproduced only about one-seventh of the observed stilling. Perhaps more telling, the difference between the two resolution tiers was small, suggesting that increasing horizontal resolution alone offers limited improvement in capturing the historical wind decline over this complex, high-altitude terrain.

This widespread stilling bias is not merely an academic curiosity. Near-surface wind speed governs a host of processes with real-world consequences: the evaporation of water from land surfaces, the dispersal of air pollutants, the mixing of the atmospheric boundary layer, and the output of wind turbines. The Tibetan Plateau is also the source region of major Asian rivers and a critical driver of monsoon circulation, so systematic errors in simulating its winds can propagate into broader failures in regional climate prediction. Previous studies have linked wind declines across China to factors ranging from urbanization and increased surface roughness to weakening pressure gradients and changes in large-scale circulation, but representing these processes faithfully in global models has remained elusive.

To dig into why the models fall short, the team applied a random forest machine learning analysis, a technique that can rank the relative importance of many candidate variables in explaining an outcome. The analysis pointed to biases in key dynamical and thermodynamical processes, with snow variability and surface sensible heat flux emerging as particularly influential factors. This makes physical sense. Snow cover on the plateau alters surface albedo and the partitioning of energy between sensible and latent heat, which in turn modulates the temperature contrasts that drive near-surface winds. If a model misrepresents how snow accumulates and melts across the plateau’s rugged topography, the resulting errors in surface heating can distort the very circulation features that generate wind.

Having identified which simulations performed best against observations using multiple statistical metrics, the researchers selected seven optimal models and constructed an optimal ensemble mean for future projections. Under the high-emission Shared Socioeconomic Pathway 5-8.5, the projections indicate that the stilling of the Tibetan Plateau will continue through the coming decades. Between 2015 and 2049, the ensemble projects an annual wind speed decline of approximately 0.02 meters per second per decade, unfolding alongside continued regional warming. While the projected rate is gentler than the historical observed trend, the direction of change is consistent, lending robustness to the conclusion that winds over the plateau will keep weakening in a warming world.

The juxtaposition of a robust future decline with a poorly simulated past decline creates something of a paradox for climate scientists. If models underestimate how much wind has already slowed, should we trust their projections of future slowing? The authors argue that the consistency of the projected decline across the selected optimal models, combined with the physical plausibility of continued warming-driven circulation changes, supports confidence in the sign of the projection even if the magnitude remains uncertain. At the same time, the study makes clear that closing the gap between simulated and observed trends requires more than computational horsepower.

That message may be the study’s most consequential takeaway for the modeling community. HighResMIP was conceived on the premise that finer grids would better resolve steep orography, coastlines, and mesoscale processes, and for some variables that promise has held. But for near-surface wind over the Tibetan Plateau, the near-identical performance of high- and low-resolution model pairs indicates that the bottleneck lies elsewhere, in the representation of physical processes such as snow-atmosphere interactions, surface energy budgets, and boundary-layer dynamics. Improving these process representations, the study suggests, may matter more than simply increasing spatial resolution for reducing persistent model biases in wind simulation.

For the people and ecosystems of the plateau, the implications extend well beyond model evaluation. Weakening winds influence glacier mass balance through altered turbulent heat exchange, affect the timing and intensity of the monsoon onset that billions of people depend on, and reshape the viability of wind energy development in one of the world’s most promising but fragile high-altitude environments. As the Third Pole continues to warm faster than the global average, the slow, steady quieting of its winds serves as another reminder that climate change reshapes even the most fundamental features of the atmosphere, and that our models, for all their sophistication, are still learning to listen.

Subject of Research: Evaluation and projection of near-surface wind speed trends over the Tibetan Plateau using CMIP6 HighResMIP simulations and station observations

Article Title: Near-surface wind speed over the Tibetan Plateau: evaluation, trend and projection

Article References: You, Q., Wu, F., Wu, T., Cai, Z., Sun, G., Jin, Z., Jiang, Z., Ullah, S., & Li, M. (2026). Near-surface wind speed over the Tibetan Plateau: evaluation, trend and projection. Climate Dynamics, 64(10), Article 411. https://doi.org/10.1007/s00382-026-08342-6

Image Credits: AI Generated

DOI: 10.1007/s00382-026-08342-6

Keywords: Tibetan Plateau, near-surface wind speed, global stilling, CMIP6, HighResMIP, climate modeling, Climate Dynamics, SSP5-8.5, random forest, snow cover, sensible heat flux, monsoon

The monsoon–desert mechanism works through a chain of dynamical processes that begins over South Asia. During the summer monsoon, intense deep convection releases enormous quantities of latent heat into the mid-troposphere. This heating excites a wave-like response in the upper-level circulation, including a strengthening of the Asian jet, and ultimately promotes large-scale descending motion, or subsidence, over the Mediterranean region. Sinking air is dry and adiabatically warming, which suppresses cloud formation and precipitation. In this way, the monsoon’s heat engine contributes directly to the hot, dry character of Mediterranean summers. The mechanism has been documented in observations and reproduced in models, and it helps explain why strong monsoon years often coincide with drier conditions over the Mediterranean.

What has remained uncertain is whether this teleconnection would survive global warming. To answer that question, a research team led by Professor Zhou Tianjun of the Institute of Atmospheric Physics at the Chinese Academy of Sciences combined real-world observations, dozens of climate model simulations, and idealized numerical experiments designed to isolate the monsoon’s remote influence. The team examined how the monsoon–Mediterranean link behaves under a future high-emissions scenario, comparing the strength of the coupling in the present climate with its strength in the second half of this century. The study’s lead author is Dr. Yu Hanzhao, also of the Institute of Atmospheric Physics.

The signal that emerged was strikingly consistent across the model ensemble. In 97.5 percent of the simulations performed with the Community Earth System Model version 1, known as CESM1, the South Asian monsoon exerted a weaker influence on atmospheric circulation over the Mediterranean in the warmer climate. The statistical relationship between monsoon heating and mid-tropospheric subsidence over the central and eastern Mediterranean, which stands at a correlation of roughly 0.4 in the present-day climate, is projected to decline to essentially zero by the second half of the century. In practical terms, the monsoon’s grip on Mediterranean summer circulation, a relationship that has helped define the region’s climate for millennia, may largely dissolve within a few decades.

The researchers traced the weakening to two distinct but interacting processes, one located over South Asia and the other over the Mediterranean itself. The first concerns the vertical structure of monsoon convection. As the climate warms, deep monsoon convection is expected to shift to higher levels of the troposphere. In the CESM1 simulations, the characteristic level of monsoon convection rises from about 452 hectopascals to about 417 hectopascals, and a similar upward shift appears across models participating in the sixth phase of the Coupled Model Intercomparison Project, or CMIP6. Because the heating source moves higher in the atmosphere, the warm response it produces spreads farther to the west. This westward extension reduces the east–west temperature contrast between South Asia and the Mediterranean, and it is precisely this contrast that drives the subsidence over the Mediterranean through atmospheric dynamics. Weaken the contrast, and the sinking motion weakens with it.

The second process unfolds locally over the Mediterranean basin. The monsoon-related sinking motion weakens most strongly in the middle and upper troposphere, and through an intrinsic atmospheric relationship known as Sverdrup balance, this change also weakens the northerly wind response at mid-to-lower levels. Those northerlies, which normally accompany and reinforce the descending motion, become weaker, and their weakening in turn further reduces the subsidence. The result is a local feedback that amplifies the remote effect of the changing monsoon heating, so that the Mediterranean circulation response declines even faster than the monsoon’s direct influence alone would suggest.

These circulation changes carry direct implications for Mediterranean rainfall. In the current climate, stronger South Asian monsoon heating tends to generate stronger subsidence over the Mediterranean and therefore less summer rainfall over the region’s land areas. Under future warming, the study finds, this inverse relationship is projected to largely disappear. The fraction of summer rainfall variability over Mediterranean land that can be statistically explained by the monsoon falls from about 14.2 percent to about 5.1 percent. In other words, the monsoon will account for a much smaller share of the year-to-year ups and downs of Mediterranean summer precipitation.

Importantly, the researchers emphasize that this does not mean Mediterranean rainfall itself will become less variable. Rather, it means that the South Asian monsoon will explain far less of that variability, implying a fundamental shift in the factors that control Mediterranean summers as the monsoon–desert coupling weakens. Other drivers of regional variability, whether local sea surface temperatures, Atlantic influences, or internal atmospheric dynamics, will presumably fill the gap left by the retreating monsoon, and identifying which of them dominates will be a central task for future research on the region.

“Global warming does not only change the mean state of temperature and rainfall. It can also reorganize the dynamical links between distant parts of the climate system,” said Professor Zhou, the corresponding author of the study. According to Zhou, the changing monsoon–desert coupling across Eurasia provides a new dynamical perspective on how large-scale modes of climate variability and teleconnections may evolve in a warmer world. It also means, he noted, that future projections of Mediterranean summer rainfall need to account for changes in the factors that control its year-to-year variability, rather than assuming that today’s statistical relationships will hold.

Dr. Yu, the study’s lead author, pointed out that atmospheric connections that appear robust today may not remain so under global warming. “Understanding how such teleconnections change will be important for projecting future regional climate variability, especially in climate change hotspots such as the Mediterranean,” he said. The Mediterranean is widely regarded as one of the most responsive regions on Earth to climate change, with warming rates that exceed the global average and acute water stress across southern Europe, North Africa, and the Middle East. If the monsoon’s moderating dry-season influence on regional circulation weakens, seasonal forecasters and climate adaptation planners alike will need to recalibrate the statistical tools they use to anticipate summer conditions. The study, published in Nature Geoscience under the title reporting reduced Asian monsoon influence on Mediterranean summers in a warmer climate, adds to a growing body of evidence that climate change is not merely shifting averages but rewiring the architecture of the global atmosphere itself, severing connections between distant regions that scientists have long relied upon to understand and predict the climate.

Subject of Research: Weakening of the monsoon–desert teleconnection between South Asian monsoon heating and Mediterranean summer subsidence under global warming

Article Title: Warming weakens Asian monsoon’s reach into Mediterranean

Article References: Warming weakens Asian monsoon’s reach into Mediterranean. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: South Asian monsoon, Mediterranean climate, monsoon–desert mechanism, teleconnection, global warming, subsidence, climate models, CESM1, CMIP6, Sverdrup balance, summer rainfall, Nature Geoscience

The Indian Ocean Dipole describes a seesaw in sea surface temperatures across the tropical Indian Ocean. In its positive phase, the western basin turns unusually warm while the southeastern basin turns unusually cool. This gradient reorganizes atmospheric circulation, shifting convection westward and altering rainfall patterns across continents bordering the ocean. For decades, forecasters and climate scientists have treated positive Dipole events as a single category, but recent research has revealed that they are not all alike. In some events, the western warming dominates the temperature anomaly pattern; in others, the eastern cooling is the stronger signal; and in a third group, the two poles contribute roughly equally.

Recognizing this diversity, the Fudan University team led a systematic assessment of how well the latest generation of global climate models captures these distinctions. The researchers analyzed 40 state-of-the-art models from the Coupled Model Intercomparison Project Phase 6, or CMIP6, the same modeling framework that underpins the most recent assessments of the Intergovernmental Panel on Climate Change. Based on which pole of the Dipole is stronger, they sorted positive events into three types: west-dominated events, designated Type-W; east-dominated events, designated Type-E; and comparable events, designated Type-C, in which neither pole clearly dominates.

The observational record for 1950 through 2023 contains 14 positive Dipole events, distributed across all three groups. The models told a strikingly different story. Over the historical period from 1930 to 2014, the 40-model average produced about 12 east-dominated events per century, compared with roughly seven in observations, and simulated them as nearly twice as strong as those actually observed. The imbalance ran in both directions: barely 11 of the 40 models managed to reproduce a single west-dominated event at all. In effect, the models systematically overproduced one flavor of the Dipole while starving the others.

The researchers traced this bias to a specific and consequential error in how the models represent the eastern Indian Ocean. In the simulations, this region runs colder than it does in reality. That cold bias matters because the three types of positive Dipole events are not independent of one another — they compete. The analysis showed a clear trade-off: the more east-dominated events a model produces, the fewer west-dominated and comparable events it generates. A model whose eastern Indian Ocean is too cold is primed to develop cold anomalies there, tipping the balance toward Type-E events at the expense of the other varieties.

This systematic bias is more than a technical curiosity. Climate risk assessments, infrastructure planning, and adaptation strategies increasingly rely on projections from models like those in CMIP6. If the models misrepresent the relative frequencies and strengths of different Dipole event types in the historical climate, then projections of future change built on those same models may be similarly skewed. The study’s authors warn that assessments of future Dipole-related hazards could be misled unless this bias is recognized and accounted for, since the regional fingerprints of the event types differ markedly from one another.

Looking forward, the team projected how each type of positive Dipole event would evolve through the remainder of the century under the highest-emission scenario, SSP5-8.5. The results point to a pronounced intensification of the east-dominated variety. The frequency of Type-E events rises from about 4 per 31-year period at the start of this century to about 7.2 per 31-year period by 2100 — an increase of 80 percent. Importantly, this is not the projection of a single outlier model: 30 of the 40 models agree on the direction and rough magnitude of the change, lending considerable robustness to the result.

The other two event types follow different trajectories. Comparable events, in which the western warming and eastern cooling contribute roughly equally, also become more frequent during the first half of the century before tapering off after mid-century. West-dominated events, by contrast, show little change over the projection period. The emerging picture is therefore not simply one of more positive Dipole events overall, but of a gradual reshuffling of the Dipole’s character — with the east-dominated flavor, and its distinctive pattern of regional impacts, increasingly taking center stage.

The physical driver behind this shift is an uneven warming of the Indian Ocean basin. As greenhouse gas concentrations rise, the western basin heats up faster than the southeastern basin, producing what the researchers describe as a west-fast, east-slow warming pattern. This asymmetric warming lowers the barrier for cold anomalies to develop in the eastern Indian Ocean, making it easier for east-dominated positive Dipole events to form. In other words, the background state of the ocean is being nudged in a direction that favors one particular mode of natural variability, effectively loading the dice in favor of the event type with the widest-reaching consequences.

Zhen-Qiang Zhou, corresponding author of the study and Associate Professor at the Department of Atmospheric and Oceanic Sciences at Fudan University, emphasized the practical stakes of distinguishing among the event types. The regional impacts of east-dominated events are particularly far-reaching, he noted: they invigorate tropical cyclone activity over the Northwest Pacific, intensify extreme rainfall over East Africa, and raise drought risk in Indonesia and Australia. Just as important, he argued, future climate risk assessments should not simply count positive Dipole events without regard to type, because the regional impacts differ markedly depending on which pole of the Dipole dominates. The research team’s next step is to examine how the different types of positive Dipole events affect climate in different regions, and to test whether current climate models can reproduce those impacts — a necessary foundation, they argue, for making future climate risk assessments more reliable as the Indian Ocean continues to warm.

Subject of Research: Projected changes in the frequency of east-dominated positive Indian Ocean Dipole events under global warming

Article Title: Global warming could make "east-dominated" Indian Ocean Dipole events 80% more frequent

Article References: Global warming could make "east-dominated" Indian Ocean Dipole events 80% more frequent. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: Indian Ocean Dipole, CMIP6, climate models, global warming, sea surface temperature, tropical cyclones, East Africa rainfall, drought, SSP5-8.5, ocean warming, climate risk, Fudan University

Chuyue Xu and Yanling Wu of Hohai University set out to answer a question that operational forecasting centers have largely been unable to tackle: how has the predictability of the boreal summer IOBM changed across different climate regimes? The obstacle has always been data. State-of-the-art dynamical prediction systems, such as the North American Multi-Model Ensemble, typically only produce hindcasts beginning around 1980, because running coupled ocean-atmosphere models backward through decades of history is computationally prohibitive. That short record makes it nearly impossible to detect whether forecast skill itself rises and falls from decade to decade. To get around this bottleneck, the researchers turned to a technique called the Model-based Analog Forecast, or MAF, method, which borrows its power from a library of long climate model simulations rather than from expensive new computations.

The logic behind the MAF method is elegantly simple. It rests on the assumption that sufficiently similar climate states tend to follow similar evolutionary pathways. For any observed initial state of sea surface temperature anomalies across the Indo-Pacific, the method searches a vast library of coupled general circulation model simulations for the twenty closest analogs, and then simply follows what happened next in those simulated worlds. The library in this study was assembled from the last 500 years of pre-industrial control simulations from 29 CMIP6 models, concatenated into a 14,500-year multimodel time series. Because the forecasts unfold entirely within the model’s own phase space, the method sidesteps the initial shock that plagues assimilation-initialized dynamical forecasts, in which observations are abruptly injected into a model and generate spurious transient adjustments.

Sensitivity tests shaped the final design of the forecasting system. The researchers found that matching analogs over the full Indo-Pacific domain, spanning 20 degrees south to 20 degrees north and 40 degrees east to 90 degrees west, produced substantially better and more stable skill than restricting the search to the Indian Ocean alone, a clear indication that Pacific ENSO signals are essential for initializing the Indian Ocean’s subsequent evolution. Ensemble size also mattered: skill improved rapidly as members were added and saturated at twenty analogs, beyond which additional members offered no statistically significant gain. The forecasts were verified against the Hadley Centre’s HadISST observational dataset, with wind and precipitation benchmarks drawn from ERA5, and all fields were interpolated onto a common two-degree grid with linear trends removed, so that the reported skill reflects only internal variability rather than external forcing.

The validation results were striking. Across all calendar months, the MAF method outperformed the majority of NMME dynamical models, and for the boreal summer IOBM index specifically, it achieved anomaly correlation coefficients above 0.6 at all lead times from one to twelve months, comparable to or even exceeding individual state-of-the-art dynamical models. Skill remained highest in the southwestern Indian Ocean, where downwelling Rossby waves deepen the thermocline and preserve persistent subsurface heat anomalies, and in the North Indian Ocean, where strong air-sea coupling anchored by monsoon-related feedbacks extends predictability. With the method’s credibility established, the team could finally ask the question that had been out of reach: does the predictability of the summer IOBM itself change over time?

The answer is an emphatic yes. Using a 21-year sliding window to filter out interannual noise while preserving low-frequency signals, and applying a change-point detection algorithm, the researchers identified statistically significant transitions in forecast skill at 1978 and 2000. Three distinct regimes emerged. Before 1980, useful forecasts extended only about four months ahead. From 1980 to 2000, skill surged: anomaly correlations consistently exceeded 0.6, and normalized errors stayed below one standard deviation of the observed index for lead times up to a full year, with correlations above 0.8 at most lead times. After 2002, skill collapsed back to levels resembling the pre-1980 era. Crucially, a simple persistence forecast, though far less skillful overall, showed the same decadal rise and fall, confirming that the regime shifts reflect genuine changes in the IOBM’s intrinsic predictability rather than artifacts of the analog method.

What drives these swings? The strongest clue came from the strength of the teleconnection between ENSO and the IOBM, measured as the correlation between the summer basin index and the preceding winter’s Niño3.4 index. The correspondence was remarkable: the time series of forecast skill correlated at 0.92 with the strength of this ENSO-IOBM link. During the high-skill epoch of 1980 to 2000, El Niño events were strong enough to leave a deep imprint on the Indian Ocean. The winter El Niño forced anticyclonic wind anomalies over the southeastern Indian Ocean, exciting westward-propagating downwelling Rossby waves that warmed the southwestern Indian Ocean substantially. That warming then activated the wind-evaporation-SST feedback, a self-reinforcing loop in which weakened surface winds suppress evaporative cooling, allowing the basin-wide warming to persist through the summer monsoon onset and into boreal summer, delivering a strong, predictable signal to forecasters.

During the low-skill periods, by contrast, the chain of events broke down early. Weaker ENSO forcing failed to generate significant southwestern Indian Ocean warming or a well-defined wind-evaporation-SST feedback. The characteristic C-shaped wind anomaly over the tropical Indian Ocean never materialized in spring, the basin-wide warming terminated prematurely, and no significant anticyclonic anomaly formed over the Northwest Pacific in the following summer. The recent low-skill period from 2002 to 2022 followed essentially the same script as the 1958 to 1978 era. The researchers also examined alternative explanations, including decadal changes in ENSO intensity, the depth of the southwestern Indian Ocean thermocline, and the amplitude of the IOBM itself. While the high-skill epoch coincided with stronger ENSO events and a shallower thermocline, neither factor correlated significantly with the skill transitions across the full 75-year record, and the Atlantic capacitor effect, in which the North Tropical Atlantic stores and relays ENSO’s influence, proved negligible once the ENSO signal was removed.

The implications reach well beyond the Indian Ocean. Because the summer IOBM anchors atmospheric anomalies that govern monsoon rainfall and heat extremes across densely populated South and East Asia, the finding that its predictability is nonstationary carries a direct warning for forecast users: the current background climate appears to have reverted to a regime in which summer IOBM signals are intrinsically harder to capture at long lead times, and overconfidence in contemporary seasonal forecasts could prove costly. The study also delivers a constructive message for model developers. Accurately representing the ENSO-IOBM teleconnection, particularly the wind-evaporation-SST feedback that sustains Indian Ocean warming into summer, is essential for advancing forecast skill, and climate models that mishandle this interbasin link will systematically misjudge their own reliability. The MAF framework itself, computationally cheap and free of initialization shock, offers forecasting centers a practical complement to dynamical ensembles and a unique window into how predictability has waxed and waned across three-quarters of a century of climate history, though its authors note it has so far been applied only in hindcast mode and not yet implemented for real-time operational forecasting.

Subject of Research: Decadal variability in seasonal predictability of the boreal summer Indian Ocean Basin Mode and its link to the ENSO teleconnection

Article Title: Decadal change in seasonal prediction skills of the Indian Ocean Basin Mode during boreal summer

Article References: Xu, C., & Wu, Y. (2026). Decadal change in seasonal prediction skills of the Indian Ocean Basin Mode during boreal summer. Climate Dynamics, 64(10), Article 441. https://doi.org/10.1007/s00382-026-08397-5

Image Credits: AI Generated

DOI: 10.1007/s00382-026-08397-5

Keywords: Indian Ocean Basin Mode, ENSO, seasonal prediction, predictability, analog forecasting, CMIP6, wind-evaporation-SST feedback, air-sea interaction, Asian summer monsoon, Indian Ocean, climate dynamics, decadal variability

The study focuses on the Sokoto Rima Basin, a Sudano-Sahelian region where millions of people depend on rain-fed agriculture and outdoor labor for their livelihoods. Unlike dry heatwaves, which are dangerous primarily through temperature alone, moist heatwaves combine high temperatures with elevated humidity, dramatically reducing the human body’s ability to shed excess heat through sweating. This distinction matters enormously for health outcomes, because the wet-bulb component of heat stress determines how quickly the body’s core temperature can rise toward lethal levels. In a region where irrigation, floodplain farming, and dense rural populations coincide, the combination of heat and humidity poses risks that simple air-temperature metrics systematically underestimate. The researchers argue that this is precisely why moist heatwaves in the basin have remained poorly characterized until now, despite mounting evidence that humidity-driven heat is among the fastest-growing climate hazards in West Africa.

To quantify the hazard, the team constructed daily Wet Bulb Globe Temperature, or WBGT, data from a multi-model ensemble of four global climate models participating in the Coupled Model Intercomparison Project Phase 6, known as CMIP6. WBGT is widely regarded as one of the most policy-relevant heat stress indices because it integrates temperature, humidity, radiation, and wind into a single measure that maps directly onto occupational health thresholds. Rather than relying on a single model, which can be skewed by individual model biases, the ensemble approach averages across models to produce a more robust central estimate of future conditions. The researchers then applied a hybrid heatwave definition that combines a relative threshold, the 90th percentile of local WBGT, with an absolute threshold of 30 degrees Celsius. This dual criterion is technically significant: the relative component ensures that heatwaves are identified as unusual for the local climate, while the absolute floor guarantees that only events genuinely hazardous to human physiology are counted, avoiding the trap of flagging mild anomalies in cooler periods as heatwaves.

Defining a heatwave is notoriously contentious in climate science, and the choice of threshold can materially change projected trends. A purely relative definition would identify the hottest 10 percent of days in any climate, even a cooling one, while a purely absolute definition might miss dangerous events in regions where the local baseline is already high. By anchoring their definition to both, the authors sidestepped two well-documented pitfalls. The approach also responds to growing calls in the literature for heatwave metrics that prioritize health impacts rather than purely statistical rarity, a shift that has gained momentum as studies of deadly heat events worldwide have shown that mortality often correlates better with humidity-adjusted indices than with raw temperature records.

The statistical machinery behind the projections is as notable as the climate modeling itself. The team employed a hurdle model, a composite framework that separates the question of whether a heatwave occurs at all from the question of how severe it is when it does. Linear regression captured trends in the mean characteristics of events, logistic regression estimated the annual probability of at least one heatwave occurring, Zero-truncated Poisson regression modeled the count of events in years when heatwaves do happen, and quantile regression traced how the entire distribution of heatwave intensity is shifting, not just its average. This layered design matters because heatwave data are count-heavy and zero-inflated: in the historical record, most years contained no events at all, which violates the assumptions of ordinary regression. The hurdle structure handles this sparsity explicitly, while quantile regression, a technique introduced by Koenker and Bassett in 1978, allows the researchers to ask whether the most extreme events are changing faster than typical ones.

The answer to that last question is one of the study’s most consequential findings. Under both emission scenarios, the slopes of the quantile regression at the 99th percentile were consistently steeper than at the median, meaning that the most extreme heatwave events are intensifying disproportionately faster than ordinary ones. In practical terms, the tail of the distribution is outpacing the middle: the worst heat stress episodes of the future will not simply be slightly worse than today’s worst, but dramatically so. This pattern echoes a broader theme in climate extremes research, where changes in variability and in the tails of distributions often carry greater societal risk than changes in the mean, because infrastructure, agriculture, and human physiology are all calibrated to historical worst-case conditions rather than average ones.

The scenario comparison adds a further layer of nuance. Under the moderate SSP 2-4.5 pathway, the basin is projected to experience 6.71 heatwave events per year by 2100, while under the high-emission SSP 5-8.5 pathway the figure is 6.07 events. At first glance, the higher-emission scenario appears less severe, but the authors caution against that reading. SSP 5-8.5 projects fewer events that are significantly longer and more intense than those under SSP 2-4.5. In other words, the difference between the two futures is not whether dangerous heat arrives, but whether it arrives as a series of discrete episodes or as prolonged, punishing spells that push cooling capacity, water supplies, and human endurance to their limits. Both pathways, the study concludes, point to an irreversible shift toward a permanent heat-stress regime by mid-century, a threshold at which adaptation can no longer be deferred without substantial loss of life and livelihood.

The implications for the Sokoto Rima Basin are stark. The region’s economy is dominated by informal and agricultural sectors, where workers have little protection from outdoor heat and where labor laws currently contain few, if any, provisions for heat stress. The authors argue that their findings make a case for urgently reforming labor regulations to mandate rest periods, hydration, and shade for outdoor workers during high-WBGT conditions. They also call for revising building codes to incorporate heat-reflective materials and natural ventilation, measures that can lower indoor temperatures in a region where air conditioning remains inaccessible to most households. A third recommendation centers on strengthening heat-health early warning systems, which remain rudimentary across much of the Sahel despite evidence that such systems save lives when heat events are forecast and communicated in advance.

The study also situates itself within a rapidly expanding body of West African heat research. Recent work has documented intensifying heat stress in Nigerian cities such as Kano, characterized heatwave dynamics across the Sahel using multiple thermal indices, and projected widespread increases in heatwave severity over West Africa using CMIP6 ensembles. Global analyses have further warned that combinations of heat and humidity approaching the limits of human tolerance are already emerging in parts of the world, with the Persian Gulf, South Asia, and now increasingly the Sahel identified as hotspots. What the new study adds is a basin-scale, health-anchored projection that extends to the end of the century and explicitly quantifies the probability structure of future events, from the likelihood of any heatwave occurring in a given year to the behavior of the most extreme tail of the intensity distribution.

For a region where the historical baseline was fewer than one heatwave every seven years, the projected trajectory toward near-annual, and eventually multiple, events represents a transformation with few precedents in the observational record. The authors emphasize that the window for cost-effective adaptation is narrowing: once the probability of dangerous heat approaches certainty, as their models indicate it will by 2050, reactive responses give way to structural ones, and the costs of inaction compound across health systems, agricultural output, and labor productivity. Whether the basin’s worst-case future features many moderate episodes or fewer but far more brutal ones depends on global emission choices made far beyond its borders. What the research makes clear is that in either case, the era of treating moist heatwaves in the Sokoto Rima Basin as rare emergencies is ending, and the institutions charged with protecting its people will need to treat chronic heat stress as the defining environmental challenge of the coming decades.

Subject of Research: Projected moist heatwave frequency, duration, and intensity in Nigeria's Sokoto Rima Basin using CMIP6 ensemble simulations and hybrid statistical modeling from 1981 to 2100

Article Title: A hybrid approach analysis of cmip6 multi-model ensemble simulations of moist heatwaves in Sokoto Rima Basin, Nigeria (1981–2100)

Article References: A hybrid approach analysis of cmip6 multi-model ensemble simulations of moist heatwaves in Sokoto Rima Basin, Nigeria (1981–2100). (n.d.). https://doi.org/10.1007/s00704-026-06609-w

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06609-w

Keywords: moist heatwaves, Sokoto Rima Basin, CMIP6, Wet Bulb Globe Temperature, heat stress, Nigeria, climate projections, SSP scenarios, quantile regression, hurdle model, West Africa, climate adaptation

The FWI is the world’s most widely adopted benchmark for assessing atmospheric conditions that promote wildfire ignition, spread, and intensity. Developed by the Canadian Forest Service, it integrates four meteorological variables: air temperature, relative humidity, precipitation, and wind speed. Through a chain of intermediate components, including the Fine Fuel Moisture Code, the Duff Moisture Code, the Drought Code, and the Initial Spread Index, the index captures how readily surface fuels dry out and how fast flames could travel once ignited. Because these components interact nonlinearly, small concurrent shifts in temperature or moisture can produce disproportionate jumps in fire danger during prolonged dry spells.

A central challenge for regional fire-risk studies is that global climate models are too coarse and too biased to represent local fire weather reliably. To overcome this, the researchers turned to the NASA Earth Exchange Global Daily Downscaled Projections dataset, known as NEX-GDDP-CMIP6, which applies statistical bias correction and spatial disaggregation to CMIP6 simulations, producing daily data on a 0.25-degree grid. The team evaluated 19 models from this ensemble against ERA5-Land reanalysis data for the historical period 1980 to 2014, using mean bias, root mean square error, and Pearson’s correlation to judge how faithfully each model reproduced the four variables feeding the FWI.

The validation results were revealing. Temperature emerged as the best-simulated variable, with correlations between roughly 0.92 and 0.94 across individual models, reflecting the fact that surface temperature is governed mainly by the surface energy balance and large-scale circulation, processes that general circulation models handle well. Wind speed and relative humidity performed moderately, while precipitation proved the weakest link, with correlations of only about 0.52 to 0.58. This is unsurprising: tropical precipitation depends on deep convection, moisture transport, cloud microphysics, and interactions with the complex Andean topography, all of which must be approximated through parameterizations that vary from model to model.

Crucially, no single model outperformed the others across all variables and metrics. Instead, the multi-model ensemble mean delivered the lowest errors and the highest correlations for every variable, reaching 0.95 for temperature, 0.79 for wind speed, 0.73 for relative humidity, and 0.63 for precipitation. By averaging out the idiosyncratic errors of individual simulations, the ensemble provided the most robust foundation for the fire danger projections, though the authors caution that uncertainty in precipitation and moisture-related variables still propagates through the nonlinear structure of the FWI and should temper interpretation of the results.

The historical baseline confirmed a pronounced seasonal rhythm in fire weather. During December through February, the dry season, Moderate and High hazard categories covered 48 percent and 31 percent of the study area respectively, with the highest index values concentrated over La Guajira, the Colombian Caribbean coast, and northern Venezuela. By June through August, the Low category dominated 96 percent of the domain as the wet season peaked, and September through November remained similarly benign. This cycle is orchestrated by the seasonal migration of the Intertropical Convergence Zone and the northeasterly trade winds, which together control when fuels dry and when moisture returns.

The projections for 2025 to 2100 show that climate change amplifies this seasonal cycle rather than erasing it. Under the low-emission SSP1-2.6 scenario, the Low category still covers 87 percent of the region on average, but under the high-emission SSP5-8.5 pathway it shrinks to 64 percent, while the Moderate category expands to 35 percent and High hazard appears for the first time in the mean conditions. The contrast is even starker for extreme events: the Very High category dominates 66 percent of the study area during extreme weather under SSP1-2.6, rising to 74 percent under SSP2-4.5 and 83 percent under SSP5-8.5, a clear signal that radiative forcing expands the territory exposed to the most dangerous fire weather.

Seasonally, the changes concentrate almost entirely in the December-February dry season, with smaller shifts in March through May. Under SSP5-8.5, High-risk conditions during the dry season reach 49 percent of the region and Very High risk emerges at 17 percent, while the wet-season months remain overwhelmingly Low-risk across all scenarios. The authors attribute this pattern to a reorganization of hydroclimatic processes over the tropical Atlantic: anthropogenic warming enhances tropospheric stability and suppresses deep convection, while an intensifying Caribbean Low-Level Jet and an expanding North Atlantic Subtropical High promote the advection of dry, subsiding air. The result is a growing propensity for flash droughts, in which abrupt increases in vapor pressure deficit strip moisture from fine fuels to critical flammability thresholds even without extreme rainfall deficits.

The study’s implications reach well beyond climate modeling. Because the FWI does not account for vegetation continuity, topography, or ignition sources, which in the Neotropics are predominantly human-caused, the projections should be read as an intensification of fire-favorable meteorology rather than forecasts of individual fires. Even so, the persistent identification of La Guajira, the lower Colombian Caribbean basin, and the northern arc of Venezuela as danger hotspots offers a quantitative basis for shifting from reactive firefighting toward anticipatory risk management, including early-warning infrastructure, preventive prescribed burns, and pre-season mobilization of resources. The urgency is underscored by recent history: carbon emissions from fires in Venezuela during the 2024-2025 seasons were roughly 50 percent above the historical average, and about 15 percent of the region’s páramo habitat burned between 1985 and 2022.

This is the first comprehensive FWI assessment for northern South America built on NEX-GDDP-CMIP6 projections, and it delivers a quantitative baseline of atmospheric fire hazard independent of the stochastic variability of human ignitions. The message for the region is sobering but actionable: the climate signal for wildfire danger emerges early, meaning even moderate warming can expand fire-favorable conditions before any dramatic shift in annual rainfall. As the authors note, integrating these meteorological projections with dynamic vegetation models, fuel loads, land-cover change, and socioeconomic ignition patterns remains the next frontier, but the framework now exists to put climate-informed fire danger at the center of land-use planning and cross-border adaptation strategy in one of South America’s most ecologically diverse corners.

Subject of Research: Projection of forest fire hazard in northern South America using the Fire Weather Index driven by NEX-GDDP-CMIP6 climate models

Article Title: Projection of Forest Fire Hazard in Northern South America Using the Fire Weather Index Driven by NEX-GDDP-CMIP6 Climate Projections

Article References: Guerra, Y., Arregocés, H. A., & Rojano, R. (2026). Projection of Forest Fire Hazard in Northern South America Using the Fire Weather Index Driven by NEX-GDDP-CMIP6 Climate Projections. Environmental Challenges, 25, Article 101659. https://doi.org/10.1016/j.envc.2026.101659

Image Credits: AI Generated

DOI: 10.1016/j.envc.2026.101659

Keywords: wildfire, Fire Weather Index, CMIP6, NEX-GDDP, climate change, Colombia, Venezuela, La Guajira, SSP scenarios, precipitation, drought, fire danger