دانشگاه محقق اردبیلی مدل سازی و مدیریت آب و خاک 2783-2546 6 1 2026 03 21 Non-stationary modeling of the meteorological drought index SPIt using generalized additive models for location, scale and shape شبیه سازی ناایستای خشکسالی هواشناسی بر مبنای شاخص SPIt با استفاده از مدل تعمیم یافته جمعی پارامترهای مکان، مقیاس و شکل 37 53 3880 10.22098/mmws.2025.17463.1602 FA محمدرضا شریفی دانشکده مهندسی آب و محیط زیست، دانشگاه شهید چمران اهواز، اهواز، ایران. کامران قیصری موزرمی دانش‌آموخته کارشناسی ارشد منابع آب، دانشکده مهندسی آب و محیط زیست، دانشگاه شهید چمران اهواز، اهواز، ایران. حیدر زارعی دانشکده مهندسی آب و محیط زیست، دانشگاه شهید چمران اهواز، اهواز، ایران 0000-0003-3414-1816 مهرداد تقیان سازمان آب و برق خوزستان، گروه برنامه ریزی منابع آب کارون بزرگ، اهواز، ایران Journal Article 2025 05 18 Introduction<br /><br />Traditional hydrological time series analyses often assume stationarity, particularly in the estimation of drought indices such as the Standardized Precipitation Index (SPI). However, increasing climate variability and anthropogenic influences have introduced significant non-stationarity into hydrological processes. This challenges the reliability of stationary-based assessments and highlights the need for models that can adapt to changing conditions. Generalized Additive Models for Location, Scale, and Shape (GAMLSS) offer a flexible framework for modeling such dynamics by allowing distribution parameters to vary over time or in relation to covariates. Recent studies suggest that non-stationary modeling improves drought characterization, particularly at longer time scales. Yet, findings remain mixed: while some report better accuracy with non-stationary approaches, others find stationary models still perform adequately, depending on regional and climatic factors. Given these variations, regional validation becomes essential. This study evaluates the performance of a non-stationary SPI-based index, referred to as SPIt, in comparison with the traditional stationary SPI. The case study is the Karkheh river basin in western Iran—a region with distinct climatic conditions compared to where SPIt was originally developed. Drought monitoring is conducted using monthly precipitation data from five stations, focusing on seasonal periods outside the dry summer months due to limited rainfall during that time. By comparing the two indices, the study aims to assess whether accounting for non-stationarity leads to more accurate drought representation in semi-arid climates.<br /><br />Materials and Methods<br /><br />The study focuses on five meteorological stations: Pol-Zal, Bostan, Pol-Kohneh, Noorabad and Halilan, spanning a historical data period ranging from 1971 to 2022 (C. E.) depending on rain gauge station availability. These stations were selected due to their diverse altitudes, geographical spread within the Karkheh basin (The latitude range is 47 to 48 degrees north and longitude 31 to 34 degrees east), and sufficiently long precipitation records, especially for the winter season (December to February), which accounts for the bulk of annual rainfall in the region. Precipitation data were analyzed for stationarity using the non-parametric Mann-Kendall trend test. Stations exhibiting significant trends were modeled using non-stationary GAMLSS, where the shape parameter of the gamma distribution was allowed to vary with time using polynomial functions optimized through the Akaike Information Criterion (AIC). The study employed a two-parameter gamma distribution to model winter precipitation in both stationary and non-stationary conditions. For drought assessment, two indices were used. SPI that Assumes stationary gamma-distributed precipitation, standardized to a normal distribution and SPIt that extends SPI by allowing the shape parameter of the gamma distribution to vary over time, thereby accommodating non-stationarity.<br /><br />Results and Discussion<br /><br />The Mann-Kendall test revealed significant decreasing trends in precipitation at Polzal and Bostan stations at 95% confidence level, and a similar albeit weaker trend at Norabad. No significant trends were detected at Pol-Kohneh and Holilan. The goodness-of-fit tests (Kolmogorov-Smirnov and Chi-square) confirmed that the gamma distribution was appropriate for all stations. GAMLSS modeling showed that non-stationary models outperformed stationary ones at stations with evident trends. For instance, AIC values were reduced by 5, 6, and 1 unit at Pol-Zal, Bostan, and Norabad, respectively, indicating better fit under non-stationary conditions. Time series analysis of the shape parameter in GAMLSS revealed temporal variability at all stations, supporting the hypothesis of non-stationarity. Worm plots for residual analysis confirmed model adequacy in both conditions, but improvements in model residuals under non-stationarity were evident at trend-affected stations. A comparison of SPI and SPIt indices indicated substantial differences in drought classification over time. At Polzal, years such as 1352 and 1354 showed no drought under SPI but were classified as moderate drought (D1) under SPIt. Similarly, years with similar rainfall amounts (e.g., 127 mm in 1971 vs. 125 mm in 2010) were categorized differently in SPIt, highlighting the model's sensitivity to underlying non-stationarity. At Holilan station, where no significant trend was observed, SPI and SPIt provided nearly identical results, reaffirming the utility of SPIt in trend-sensitive environments. A station-wise drought frequency comparison between SPI and SPIt further revealed that non-stationary modeling generally results in higher estimated drought frequencies at trend-affected stations. For example, the frequency of droughts at Polzal increased from 48% (SPI) to 52% (SPIt). Similar increases were noted at Bostan and Noorabad. Conversely, at Holilan and Polkohneh, where no significant trends were detected, the drought frequency remained the same or slightly decreased under SPIt. Moreover, the frequency of severe droughts (D4) decreased under the non-stationary model, with D4 events dropping from 2% to 0% at Pol-Zal, from 3% to 0% at Bostan, and from 6% to 4% at Noorabad. This suggests that the SPI may overestimate drought severity when stationarity is incorrectly assumed.<br /><br />Conclusion<br /><br />Long-term drought monitoring at various rain gauge stations highlights the importance of considering changes in precipitation when making decisions and setting policies in watersheds. When significant trends are present, drought analysis can be performed under either stationary or non-stationary assumptions, depending on the objective. If the primary concern is drought frequency, non-stationary analysis is strongly recommended. Results showed that at stations with trends such as Polzal, Bostan, and Norabad the frequency of droughts was underestimated under stationary analysis compared to non-stationary models. For example, frequencies increased from 48, 48, and 43 (stationary) to 52, 52, and 47 (non-stationary), respectively. However, in non-trending stations, stationary models may still provide reliable results for frequency estimation. In contrast, if the focus is on severe droughts, stationary models may outperform non-stationary ones at trend-affected stations. Non-stationary analysis yielded zero severe drought events, while stationary models identified 2, 3, and 6 cases in Palzal, Bostan, and Norabad, respectively. In non-trend stations like Polkohneh and Holilan, non-stationary analysis was more effective in detecting severe events. These findings align with previous research suggesting that while non-stationary models, such as those using GAMLSS, offer better parameter estimation, stationary models may sometimes better reflect reality in future projections. Therefore, although non-stationary modeling is essential under climate variability, the choice of model should depend on the monitoring goal. It is also recommended to incorporate time-varying variance and alternative probability distributions for better drought characterization under potential extreme rainfall events. اگرچه در نظر گرفتن خصوصیت ناایستایی سری‌های زمانی در مقایسه با ایستا فرض نمودن آن، منجر به ارتقای برآورد پارامترهای توزیع احتمالاتی متغیر می‌شود، برتری شبیه سازی ناایستا نسبت به ایستا و به‌دنبال آن پیش بینی وضعیت در آینده، تحت تأثیر محل و در نتیجه عوامل موثر بر ناایستایی، نتایج متفاوتی به‌دست داده است. از این‌رو در مطالعة حاضر با هدف ارزیابی نتایج حاصل از شرایط ناایستایی بر پایش خشک‌سالی، اقدام به مدل‌سازی ناایستای خشک‌سالی هواشناسی، با شاخص خشک‌سالی SPIt، با استفاده از مدل تعمیم یافتة جمعی پارامترهای مکان، مقیاس و شکل GAMLSS در پنج ایستگاه باران سنجی پل‌زال، بستان، پل‌کهنه، نورآباد و هلیلان، در حوزة آبریز کرخه و مقایسة آن با مدل ایستای شاخص خشک‌سالی SPI، شد. دورة زمانی مورد مطالعه، از حداقل 31 ساله (1396-1365) در ایستگاه نورآباد تا حداکثر 52 ساله (1401-1350) در ایستگاه پل‌زال، شامل بارش تجمعی فصل زمستان است. نتایج نشان داد تخمین نا ایستای پارامترها، در ایستگاه‌های دارای روند، در مقایسه با تخمین ایستایی در آن‌ها، دارای دقت بیش‌تری است. به‌طوری‌که تخمین ناایستا، سبب 5، 6 و 1 واحد کاهش مقدار آکاییک به‌ترتیب، در ایستگاهای پل‌زال، بستان و نورآباد، نسبت به شرایط ایستا، شد. هم‌چنین در نظر گرفتن شرایط ناایستایی در ایستگاه‌های دارای روند، سبب افزایش فراوانی خشک‌سالی شد. به‌طوری‌که درصد فراوانی خشک‌سالی ایستگاه‌های پل زال، بستان و نورآباد به‌ترتیب از 48، 48 و 43 در حالت ایستا به 52، 52 و 47 درصد در حالت ناایستا، افزایش یافت. این در حالی است که درصد فراوانی خشک‌سالی‌های شدید (کلاس D4) در حالت ناایستا در مقایسه با ایستا، کاهش نشان داد. به‌طوری‌که در ایستگاه‌های پل‌زال، بستان و نورآباد، درصد فراوانی خشک‌سالی‌های شدید (کلاس D4)، در شرایط نایستا به‌ترتیب صفر، صفر و 4 درصد و در شرایط ایستا به‌ترتیب، 2، 3 و 6 درصد به‌دست آمد. از این‌رو بسته به هدف پایش خشک‌سالی در مدیریت حوزه‌های آبریز، مبنی بر این‌که مسأله بحرانی ناشی از خشک‌سالی، فراوانی وقوع یا فراوانی شدیدترین خشک‌سالی مد نظر باشد، در علیرغم وجود ناایستایی، به‌ترتیب تحلیل ناایستا و ایستا، نتایج توام با ریسک کم‌تری به‌دست خواهد داد. https://mmws.uma.ac.ir/article_3880_7c8da9a6a27af69a5a31ef973ef4a3bb.pdf