Reference evapotranspiration, abbreviated ETo, is the backbone of irrigation planning worldwide. It represents the rate at which water would be lost from an idealized, uniformly green surface with ample moisture, and it is the quantity that farmers, water managers, and hydrologists multiply by crop coefficients to estimate real irrigation demand. The internationally accepted method for computing it is the FAO Penman–Monteith equation, a physically based formula that combines temperature, humidity, wind, and solar radiation into a single estimate. Because the equation draws on so many inputs simultaneously, a long-standing puzzle has been which of those inputs matters most, and whether the answer shifts from place to place and decade to decade. That question has become urgent as climate change alters the individual variables at different rates.

Zarei and Mahmoudi assembled a dataset spanning 1967 to 2024, one of the longest windows ever applied to this problem in the Middle East. From the records of 34 synoptic stations, they extracted nine candidate drivers: maximum, minimum, and mean temperature; maximum, minimum, and mean relative humidity; wind speed; sunshine hours; and precipitation. For thirteen representative years distributed across the nearly six decades, they computed ETo at every station using the FAO Penman–Monteith equation and then interpolated the results into continuous spatial maps using the Inverse Distance Weighting method, a technique that estimates values at unsampled locations by weighting nearby observations according to their distance.

To move beyond station-by-station analysis, the researchers overlaid 2,000 randomly distributed points across Iran’s territory and extracted both climatic data and interpolated ETo values at each one. This created a vast statistical sample capturing the country’s extraordinary climatic diversity, from the humid Caspian coast to the arid central deserts and the mountainous northwest. The random-point approach meant that the subsequent modeling reflected the spatial character of the whole country rather than the arbitrary placement of weather stations, which tend to cluster in populated areas.

The heart of the study lies in its dual statistical framework, and this is where the methodological story becomes genuinely interesting for anyone who follows climate science. The first tool, the Corrplot approach, examines simple linear correlations between each climatic variable and ETo, identifying which variable shows the strongest individual association in each year. The second, the Lindeman, Merenda, and Gold method, usually abbreviated LMG, belongs to a family of techniques called variance decomposition. Rather than looking at each variable in isolation, LMG averages the contribution of each predictor across all possible orderings in which it could enter a regression model, thereby partitioning the explained variance fairly among predictors that are correlated with one another.

Why does that distinction matter so much? Climatic variables are notoriously interdependent. Maximum temperature and sunshine hours rise and fall together; humidity tracks inversely with temperature; precipitation suppresses both sunshine and temperature. In statistical language, this is multicollinearity, and it is the bane of attribution studies. A naive correlation analysis can assign importance to a variable that is merely riding on the coattails of another. The LMG method, by decomposing R-squared across all orderings of predictors, corrects for this shared variance and reveals which variable carries genuine, independent explanatory power. The contrast between the two methods in this study is a textbook demonstration of why the correction matters.

According to the Corrplot results, maximum temperature was the most frequently dominant variable, leading in 46.16 percent of the years examined, followed by sunshine hours at 23.08 percent, with precipitation and wind speed each dominant in 15.38 percent of years. That picture is broadly consistent with intuition in a warming country: hotter days, more evaporative demand. But when the LMG variance decomposition was applied, the hierarchy flipped dramatically. Sunshine hours emerged as the most frequently dominant factor in 92.31 percent of the study period, with maximum temperature leading in only 15.38 percent of cases. In other words, much of the apparent importance of temperature in the simple correlation analysis appears to reflect its correlation with sunshine duration rather than an independent physical influence of that magnitude.

The authors validated both the spatial interpolations and the statistical models rigorously. The quality of the IDW-generated maps was assessed using the Nash–Sutcliffe Efficiency coefficient, a standard hydrological goodness-of-fit measure, and Willmott’s Index of Agreement, which evaluates how closely predicted values track observed ones. Both metrics confirmed the reliability of the spatial products at the 1 percent significance level, and the LMG model’s performance was likewise statistically verified at the 1 percent level for every year investigated. These are demanding thresholds, and their satisfaction across nearly six decades of data lends considerable weight to the central finding.

The practical implications for Iran are substantial. The country is already one of the most water-stressed in the world, with shrinking aquifers, recurrent drought, and intense competition between agriculture and urban demand. If sunshine hours, a proxy for the solar radiation that supplies the energy for evaporation, are the true dominant lever on ETo, then projections of future water demand need to track cloud cover, aerosol loading, and radiation trends as closely as they track temperature. It also suggests that the well-documented global phenomenon of dimming and brightening, in which the solar radiation reaching the surface has changed over decades due to aerosols and cloudiness, may have played an underappreciated role in Iran’s evolving water balance. Irrigation scheduling tools, drought indices, and crop water models calibrated around temperature sensitivity alone could be systematically misattributing the drivers of demand.

The study also offers a methodological lesson that extends well beyond Iran. Sensitivity analyses of the Penman–Monteith equation have been published for regions around the world, from the Yangtze basin to Australia’s climatic zones to tropical savannas, and their conclusions have often varied. Some of that variation is surely real, reflecting genuinely different climates. But this study demonstrates that part of the disagreement may be an artifact of method: analyses relying on simple correlations can rank drivers differently than analyses that properly account for multicollinearity. As climate change reshapes the individual variables at unequal speeds, knowing which lever actually moves evapotranspiration will determine whether water managers in arid lands anticipate the future accurately or are surprised by it. For Iran, at least, the answer now points upward, toward the sun itself.

Subject of Research: Spatiotemporal drivers of reference evapotranspiration in Iran over 58 years

Article Title: A multi-decadal spatiotemporal evaluation of reference evapotranspiration drivers in Iran: Leveraging an integrated LMG and Corrplot statistical framework

Article References: Zarei, A. R., & Mahmoudi, M. R. (2026). A multi-decadal spatiotemporal evaluation of reference evapotranspiration drivers in Iran: Leveraging an integrated LMG and Corrplot statistical framework. Theoretical and Applied Climatology, 157(11), Article 703. https://doi.org/10.1007/s00704-026-06635-8

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06635-8

Keywords: reference evapotranspiration, Iran, FAO Penman-Monteith, sunshine hours, LMG variance decomposition, multicollinearity, Corrplot, Inverse Distance Weighting, climate variables, water resources, irrigation planning, Theoretical and Applied Climatology

Dry spells, defined as stretches of consecutive days with negligible rainfall, are among the most consequential yet underappreciated features of semi-arid hydroclimatology. Unlike droughts, which are usually assessed over months or seasons, dry spells operate on the scale of days and weeks, precisely the timescale at which crops experience water stress. For a deltaic district like Thanjavur, where paddy cultivation depends on both monsoon rainfall and canal irrigation from the Cauvery system, the duration and frequency of rainless runs determine how much supplementary water must be drawn from canals or groundwater. The researchers argue that a robust climatology of dry spell behavior is therefore a prerequisite for sustainable irrigation planning and drought vulnerability assessment in the region.

The study’s methodology unfolded in four phases. First, the team assembled a long-term climatology baseline of dry spell characteristics at each of the 13 stations. Second, they developed a descriptive, percentile-based classification scheme to sort dry spells into distinct regimes, allowing comparisons across seasons and locations. Third, they examined temporal variability and long-term trends using two widely trusted non-parametric tools: the Mann-Kendall test, which detects monotonic trends in time series without assuming any particular distribution, and Sen’s Slope estimator, which quantifies the magnitude of such trends. Finally, they turned to unsupervised and supervised machine learning to confirm that their classification reflected real structure in the data rather than arbitrary thresholds.

The machine learning confirmation relied on three complementary techniques. Principal Component Analysis (PCA) reduced the dataset to its essential dimensions, revealing that the first components captured 79.3 percent of the total variance, a strong indication that the underlying dataset has a compact and interpretable structure. K-means clustering and Agglomerative Hierarchical clustering were then applied to group stations with similar dry spell behavior. The hierarchical approach produced a dendrogram, a tree-like diagram of station similarity, that delivered one of the study’s most striking findings: the Thiruvaiyaru station stands apart from all other stations in the district, exhibiting a distinctly higher level of dry spell severity than its neighbors.

This spatial heterogeneity matters for practical water management. If dry spell risk were uniform across Thanjavur, a single district-wide irrigation schedule might suffice. Instead, the clustering results show that stations just tens of kilometers apart can belong to fundamentally different dry spell regimes, meaning that water allocation, canal scheduling, and groundwater extraction plans may need to be tailored station by station. The identification of Thiruvaiyaru as an outlier gives local authorities a concrete geographic priority for drought monitoring and for investments in storage or conveyance infrastructure that can buffer longer rainless periods.

On the temporal side, the trend analysis produced a genuinely encouraging result. Across the 37-year record, the annual dry spell duration shows a declining trend, and the regional analysis reveals that the strongest decreasing trend occurs in Thanjavur itself. In other words, over nearly four decades, the longest and most punishing rainless stretches have, on average, been shortening in this deltaic district. The authors link dry spell behavior to large-scale climate teleconnections, including the El Niño-Southern Oscillation (ENSO) and the Indian Ocean Dipole-related variability they refer to as INO, both of which are known to modulate Indian monsoon rainfall. Understanding these links opens the door to seasonal forecasting: if a developing El Niño or Indian Ocean event can be tied statistically to dry spell risk, water managers could anticipate stress months before it materializes.

The supervised learning phase of the study tested whether the dry spell classification could be predicted from hydroclimatic features alone. Two gradient- and ensemble-based algorithms, Random Forest (RF) and XGBoost, both staples of modern applied machine learning, achieved classification accuracy above 93 percent. That level of performance demonstrates that dry spell regimes in the Cauvery Delta are not stochastic noise but emerge from identifiable patterns in the hydroclimatic record. Even more valuable is the feature importance analysis: the single most powerful predictor of dry spell classification turned out to be the annual number of rainfall days. This is an intuitively satisfying result, since a landscape that receives rain on more days per year has fewer opportunities for long consecutive dry runs, but quantifying that relationship gives forecasters a simple, measurable indicator to track.

The broader scientific context makes this work timely. Across India and globally, researchers have documented shifting precipitation structures under a warming climate, with changes in the frequency and intensity of both wet and dry extremes. Studies of Indian meteorological subdivisions, the Indo-Gangetic Plains, and Mediterranean basins have all highlighted the value of classifying dry spells rather than treating rainfall as a single aggregate quantity. What distinguishes the Thanjavur study is its integration of classical climatology, rigorous non-parametric trend testing, and machine learning validation within a single framework applied at the district scale, where the results can feed directly into local decision-making rather than remaining at the level of regional generalization.

For the Cauvery Delta, where disputes over water sharing and the pressures of growing demand have made every cubic meter consequential, the study provides what the authors describe as a robust scientific basis for drought monitoring, irrigation planning, water resource management, and climate adaptation strategies. The declining trend in annual dry spell duration offers a measure of reassurance, but the pronounced spatial contrasts, exemplified by Thiruvaiyaru’s outlier status, and the demonstrated influence of ENSO and Indian Ocean variability caution against complacency. A single strong teleconnection event could still deliver a season of extended dry spells, and the percentile-based classification gives planners a common vocabulary for describing how severe such a season would be relative to the 37-year baseline.

The methodological template is arguably as important as the regional findings. By showing that PCA can compress 37 years of multi-station data into a few dominant dimensions, that clustering can expose meaningful geographic structure, and that Random Forest and XGBoost can classify dry spell regimes with over 93 percent accuracy while identifying annual rainfall days as the key predictor, the study offers a replicable workflow for any monsoon-dependent region facing similar questions. As climate variability intensifies and semi-arid agricultural zones worldwide confront longer, less predictable rainless periods, the ability to characterize, classify, and anticipate dry spells at the district scale may prove to be one of the most practical tools in the adaptation toolkit. For the farmers of Thanjavur, whose fields have depended on the rhythm of the Cauvery for generations, that rhythm has now been measured with unprecedented statistical clarity.

Subject of Research: Characterization and machine learning classification of dry spells in the Cauvery Delta district of Thanjavur, India

Article Title: The characterization of dry spell through classification and clustering analysis for cauvery river deltaic district of thanjavur: a 37-year hydroclimatic assessment

Article References: Priyanka, K., & Suribabu, C. R. (2026). The characterization of dry spell through classification and clustering analysis for cauvery river deltaic district of thanjavur: a 37-year hydroclimatic assessment. Theoretical and Applied Climatology, 157(10), Article 627. https://doi.org/10.1007/s00704-026-06554-8

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06554-8

Keywords: dry spells, Thanjavur, Cauvery Delta, drought, Mann-Kendall test, Principal Component Analysis, K-means clustering, Random Forest, XGBoost, ENSO, monsoon, irrigation planning