The technical backbone of the study is the self-calibrating Palmer Drought Severity Index, or scPDSI, a water-balance metric that compares the precipitation a location receives against the atmospheric demand for water, known as potential evapotranspiration, while accounting for the water-holding capacity of the local soil. Rather than relying on sparse station records alone, the team fed the index with TerraClimate precipitation and potential evapotranspiration data, a high-resolution gridded dataset that blends observations with interpolated climate surfaces. Soil water-holding properties were drawn from the FAO Harmonized World Soil Database, allowing the index to reflect the fact that sandy soils in the desert interior and deeper soils in the Zagros Mountains respond very differently to the same rainfall deficit. National meteorological observations were used independently to characterize the long-term evolution of temperature and precipitation, providing a check on the gridded reconstruction.

The headline finding is temporal: drought conditions in Iran became more widespread and more severe during several periods after the late 1990s, though with substantial spatial variability across the country’s remarkably diverse climate zones, from the humid Caspian coast to the arid Lut Desert. Trend detection relied on the Mann-Kendall test, a non-parametric method that identifies monotonic changes without assuming the data follow a normal distribution, paired with Sen’s slope estimator, which provides a robust measure of the magnitude of change. This combination is standard practice in hydroclimatology precisely because drought indices are noisy, autocorrelated and frequently non-Gaussian, conditions under which ordinary least-squares trend fitting can be badly misleading.

With the drought record in hand, the researchers turned to the teleconnection question. They examined six major climate oscillations: the El Niño-Southern Oscillation, represented by the Southern Oscillation Index and the Niño3.4 sea surface temperature index; the Pacific Decadal Oscillation; the Atlantic Multidecadal Oscillation; and the Arctic and North Atlantic Oscillations. At the national scale, Pearson correlations between the monthly scPDSI and these indices revealed significant associations with the Atlantic Multidecadal Oscillation, with a correlation coefficient of minus 0.448 and a p-value below 0.0001, as well as with the Southern Oscillation Index at minus 0.307, the Niño3.4 index at 0.315, and the Pacific Decadal Oscillation at 0.242. Crucially, the p-values were computed using an effective sample size correction based on the lag-1 autocorrelation of the paired monthly series, a step that guards against the inflated significance that arises when slowly varying climate series are treated as independent observations.

Just as revealing were the indices that failed to show a national-scale signal. The Arctic Oscillation and the North Atlantic Oscillation produced correlations of only minus 0.041 and 0.036 respectively, with p-values around 0.4, statistically indistinguishable from zero at the country level. For a country whose winter precipitation is strongly influenced by the position and strength of the mid-latitude storm track, this null result is noteworthy and suggests that whatever influence these northern annular modes exert on Iranian rainfall may be regionally confined, seasonally specific, or mediated by other circulation features in ways that a single national-average correlation cannot capture.

The study’s most sophisticated contribution lies in its treatment of space and frequency. Rather than stopping at national averages, the team produced spatial correlation maps at the grid-cell level, identifying statistically significant cells after applying the Benjamini-Hochberg procedure to control the false discovery rate at q below 0.05. This multiple-testing correction matters enormously in climate science: when thousands of grid cells are tested simultaneously, a naive significance threshold guarantees a flood of false positives, a problem that has led to overstated and overinterpreted maps in the published literature. By controlling the expected proportion of false discoveries, the researchers ensure that the geographic patterns they report, which show considerable heterogeneity in how strongly each oscillation is associated with drought in different regions of Iran, are unlikely to be statistical mirages.

To probe the timing dimension, the authors deployed two wavelet-based tools: the cross-wavelet transform and wavelet coherence. Unlike classical correlation, which collapses an entire time series into a single number and assumes any relationship is stationary, wavelet methods decompose both signals into time-frequency space, revealing when particular periodicities were coherent and when they were not. The cross-wavelet transform highlights regions of high common power between two series, while wavelet coherence measures the local correlation between them across both time and frequency, even where common power is weak. Applied to the scPDSI and the teleconnection indices, these analyses indicated recurrent associations between Iranian drought and the ENSO-related indices at interannual to multi-year timescales, consistent with the well-documented influence of El Niño and La Niña events on precipitation across the Middle East and Southwest Asia.

At longer periods, a different picture emerged. The Pacific Decadal Oscillation and the Atlantic Multidecadal Oscillation showed associations extending toward multi-year to decadal periodicities, hinting that the slow background state of the Pacific and Atlantic oceans may modulate the severity of Iranian drought episodes over decades, potentially by shifting large-scale circulation patterns and sea surface temperature gradients that govern storm tracks and moisture fluxes into the region. This separation of labor, with ENSO dominating the interannual band and the Pacific and Atlantic decadal modes operating at lower frequencies, mirrors findings from other parts of the world, including North America, where multidecadal drought frequency has been linked to the combined state of the two ocean basins.

The authors are careful, and rightly so, about what these results do and do not establish. The identified relationships are statistical associations, not demonstrations of causality, and they do not by themselves translate into predictive skill for drought early warning. A significant correlation between the Atlantic Multidecadal Oscillation and Iranian drought does not reveal the physical pathway, whether through shifts in the subtropical jet, changes in the North Atlantic storm track, or altered moisture transport from the Mediterranean and Arabian seas, that would connect ocean temperature anomalies to rainfall deficits over the Iranian plateau. The researchers explicitly call for further synoptic-dynamic and atmospheric moisture-transport analyses to clarify the mechanisms underlying the teleconnections they detect, and such work would be a prerequisite for any operational forecasting application.

Even with those caveats, the study arrives at a moment when Iran’s water crisis has made drought science a matter of national urgency. Decades of groundwater over-extraction, recurring multi-year dry episodes and rising evaporative demand under a warming climate have strained the country’s water resources to a breaking point, and understanding the natural variability superimposed on that long-term stress is essential for planning. By showing that drought-teleconnection relationships in Iran vary across space, time and frequency, the analysis argues against one-size-fits-all drought monitoring and instead points toward regionally tailored early-warning systems that incorporate large-scale climate indices alongside local observations. Whether the oceans’ slow rhythms can be harnessed to give Iranian farmers and water managers a few seasons of advance notice remains an open question, but this study provides the statistical groundwork on which that answer will have to be built.

Subject of Research: Spatiotemporal variability of drought in Iran and its statistical associations with large-scale climate teleconnection patterns from 1958 to 2022

Article Title: Spatiotemporal drought variability in Iran and the role of teleconnections

Article References: Mahdizadeh, J., Ghavidel, Y., & Farajzadeh, M. (2026). Spatiotemporal drought variability in Iran and the role of teleconnections. Theoretical and Applied Climatology, 157(11), Article 709. https://doi.org/10.1007/s00704-026-06638-5

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06638-5

Keywords: drought, Iran, teleconnections, scPDSI, ENSO, Pacific Decadal Oscillation, Atlantic Multidecadal Oscillation, TerraClimate, wavelet coherence, climate variability, drought monitoring, Theoretical and Applied Climatology

The study, conducted by Amaury de Souza of the Federal University of Mato Grosso do Sul and José Francisco de Oliveira Junior of Fluminense Federal University and the Federal University of Alagoas, focused on the capital cities of Midwest Brazil: Brasília, Goiânia, Cuiabá, and Campo Grande. The team drew on the TerraClimate dataset, a high-resolution global climate product with roughly four-kilometer spatial resolution, to extract monthly precipitation and actual evapotranspiration from 1958 to 2024. Crucially, the dataset contained no missing values, eliminating the need for gap-filling procedures that can introduce artificial structure into a time series. From these two variables, the researchers constructed what they call the Hydroclimatic Balance Index, or HBI — the difference between precipitation and actual evapotranspiration, accumulated over twelve months and standardized as a z-score.

The choice of actual evapotranspiration rather than potential evapotranspiration is what sets the HBI apart from its better-known cousin, the Standardized Precipitation Evapotranspiration Index, or SPEI. Potential evapotranspiration measures how much water the atmosphere would demand if supply were unlimited; actual evapotranspiration measures how much water is genuinely exchanged between the land surface and the atmosphere, constrained by soil moisture, vegetation dynamics, and land-surface feedbacks. In tropical continental environments where evapotranspiration is a major component of the hydrological cycle, that distinction matters. The authors stress that the HBI is not intended to replace established drought indices, but to complement them by capturing realized hydroclimatic conditions rather than atmospheric demand alone. As a benchmark, the team also computed the classic Standardized Precipitation Index at the twelve-month scale, SPI-12, which relies exclusively on rainfall.

Before turning to fractal analysis, the researchers established a rigorous statistical baseline. Linear regression of annual precipitation, evapotranspiration, and HBI series revealed only weak slopes with p-values all exceeding the five percent significance level and coefficients of determination that were vanishingly small. The non-parametric Mann-Kendall test, which operates on data ranks and is robust to the outliers that plague hydroclimatic records, told the same story: Kendall’s tau coefficients hovered near zero, confirming the absence of significant monotonic trends. The Pettitt change-point test, designed to detect abrupt structural shifts without assumptions of normality, found significant breakpoints only in the evapotranspiration series of Goiânia in 2001 and Cuiabá in 2009. By every conventional measure, the region’s climate looked stable — which is precisely what made the next set of results so striking.

That next step was Multifractal Detrended Fluctuation Analysis, or MFDFA, a technique developed by Kantelhardt and colleagues in 2002 that extends earlier methods for detecting long-range correlations in nonstationary signals. The procedure integrates the original series, divides it into segments of varying lengths, removes local polynomial trends from each segment, and then examines how fluctuation magnitudes scale with segment size. For a monofractal process, the generalized Hurst exponent h(q) remains constant regardless of the moment order q. For a multifractal process, h(q) varies systematically with q, revealing that small and large fluctuations obey different scaling rules — a hallmark of complex, intermittently organized systems. Rainfall, temperature, wind, and river flow have all been shown to behave this way, but multifractal assessments of drought indices in tropical continental settings have remained scarce.

The results were unambiguous. Across all four capitals, the generalized Hurst exponent declined systematically with increasing moment order, confirming genuine multifractal scaling in both SPI-12 and HBI-12 series. The Hurst exponent at q equals two ranged from approximately 0.87 to 1.11, values that indicate strong long-range dependence — a form of hydroclimatic memory in which wet and dry conditions cluster together rather than arriving randomly. Persistent systems of this kind are prone to drought clustering and prolonged anomalies, a property with direct implications for agriculture and hydropower in a region that anchors Brazil’s agricultural production and generates a substantial share of its electricity. The width of the singularity spectrum, denoted delta-alpha, served as the study’s measure of multifractal strength: values below 0.5 indicate weak multifractality, values between 0.5 and 1.0 moderate, and values above 1.0 strong.

The spatial patterns proved revealing. Brasília exhibited the weakest multifractality, with SPI-12 producing a narrow spectrum of just 0.194, though even there the HBI-12 spectrum widened to 0.517 when evapotranspiration entered the calculation. Goiânia showed strong multifractality for both indices, at 1.333 and 1.382 respectively, suggesting that rainfall variability alone already accounts for much of its complexity. The most dramatic contrasts emerged in Cuiabá and Campo Grande. Cuiabá jumped from weak multifractality under SPI-12 at 0.254 to strong multifractality under HBI-12 at 1.008, while Campo Grande surged from a moderate 0.755 to 1.488 — the strongest multifractal signature of any capital. In these locations, evapotranspiration processes, vegetation dynamics, soil moisture availability, and land-surface feedbacks amplify nonlinear fluctuations far beyond what precipitation alone can explain.

Perhaps the most consequential finding came from splitting the record into two sub-periods, 1958 to 1990 and 1991 to 2024. The recent period showed broader singularity spectra in several capitals, with Campo Grande reaching a delta-alpha of 1.938 and Cuiabá 1.341, indicating intensified intermittency and hydroclimatic heterogeneity in recent decades. Brasília, by contrast, showed a modest reduction, underscoring the spatially uneven character of the change. Spectral asymmetry values predominantly exceeded one, meaning that large-magnitude fluctuations — the extreme wet and dry events — contribute disproportionately to the multifractal organization. Yet throughout, generalized Hurst exponents remained above unity, confirming persistent long-range dependence across all cities and both periods. The intensification of variability, in other words, is happening within a persistently memory-laden system, a combination that favors clustering of extremes.

The contrast between the classical and fractal diagnoses carries a pointed message for drought monitoring. Linear regression, Mann-Kendall tests, and Pettitt analysis detected almost nothing, yet the multifractal framework revealed substantial structural variability, growing intermittency, and shifts in scaling behavior that no trend statistic could capture. The authors argue that hydroclimatic change in Central Brazil should not be judged solely by the presence or absence of monotonic trends; increasing persistence, enhanced intermittency, and greater scaling heterogeneity may be equally important indicators of environmental change. This resonates with a broader body of climate-complexity research showing that scaling diagnostics can expose structural variability invisible to conventional methods, particularly in regions governed by monsoon dynamics and land-atmosphere interactions.

The study’s implications extend beyond Brazil. The measured spectrum widths, ranging from roughly 0.19 to 1.94, are comparable to or exceed those reported for Mediterranean rainfall systems, Indian monsoon basins, African semi-arid regions, and East Asian monsoon domains — but the gap between precipitation-only and water-balance multifractality appears unusually large in Central Brazil, underscoring the outsized role of evapotranspiration in tropical continental hydroclimates. The authors acknowledge limitations: the analysis rests entirely on gridded TerraClimate data, which may smooth local extremes; only the twelve-month accumulation scale was evaluated; and causal attribution to land-use change, though plausible given the Cerrado’s rapid agricultural transformation, would require integration with remote sensing and land-cover datasets. Still, the core conclusion stands. Drought monitoring frameworks built on precipitation alone may systematically underestimate hydroclimatic complexity in monsoon-influenced tropical regions, and as warming and land-use change continue, complexity-based diagnostics may become essential tools for anticipating when quiet variability turns into visible crisis.

Subject of Research: Multifractal assessment of hydroclimatic variability and drought indices in Midwest Brazil

Article Title: Multifractal assessment of hydroclimatic variability in Midwest Brazil using the hydroclimatic balance index (HBI): a comparison with the standardized precipitation index (SPI)

Article References: de Souza, A., & de Oliveira Junior, J. F. (2026). Multifractal assessment of hydroclimatic variability in Midwest Brazil using the hydroclimatic balance index (HBI): a comparison with the standardized precipitation index (SPI). Theoretical and Applied Climatology, 157(10), Article 682. https://doi.org/10.1007/s00704-026-06467-6

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06467-6

Keywords: multifractal analysis, MFDFA, drought, hydroclimatic variability, Brazil, Cerrado, evapotranspiration, Standardized Precipitation Index, Hurst exponent, TerraClimate, South American Monsoon, water balance