ارزیابی چند مدل در پیش‌بینی سری‌های زمانی شاخص بارندگی استاندارد

نوع مقاله : مقاله پژوهشی

نویسندگان

1 دانشجوی دکترای هواشناسی‌کشاورزی دانشگاه تهران؛

2 عضو هیئت علمی

چکیده

خشکسالی یک رخداد هواشناسی تکرار­شونده و موقتی است که ناشی از کاهش بارندگی نسبت به میانگین بلندمدت آن می­باشد. با توجه به اهمیت پیش­بینی خشکسالی و نقش آن در مدیریت منابع آب، در این مطالعه به ارزیابی مهارت شش نوع مدل استوکاستیک (AR، MA، ARMA، ARIMA غیرفصلی، ARIMA فصلی و ARIMA مکثّر) در مدلسازی و پیش­بینی سری­های زمانی شاخص بارندگی استاندارد شده (SPI) پرداخته شد. برای دستیابی به این هدف، از داده­ های بارندگی ماهانه­ ی 10 ایستگاه سینوپتیک با اقلیم ­های خیلی مرطوب تا فراخشک (دوره­ی آماری 2007-1973) استفاده شد. ابتدا از طریق انتقال احتمال تجمعی توزیع بهینه بارندگی به توزیع احتمال تجمعی نرمال استاندارد، مقادیر SPI در مقیاس­های زمانی 3، 6 و 12 ماهه محاسبه شد. سپس، در یک فرآیند چند مرحله‌ای (شناسایی، برآورد پارامتر و آزمون کفایت مدل) توسعه­ ی مدل­ها بر روی مقادیر SPI مربوط به دوره­ی 1973 تا 2000 صورت گرفت و از بین مدل­های داوطلب، مناسب­ترین مدل استوکاستیک برای هر یک از مقیاس‌های زمانی تعیین شد. به منظور صحت­ سنجی مدل­ های منتخب، پس از انجام پیش ­بینی­ های یک تا دوازده گام به جلو مقادیر SPI برای دوره­ی 2001 تا 2007، به مقایسه­ی مقادیر و طبقات مشاهده­ شده و پیش­ بینی­ شده ­ی SPI پرداخته شد. نتایج ارزیابی دقت مدل‌ها در پیش‌بینی‌ مقادیر SPI در ایستگاه­های منتخب نشان داد در پیش ­بینی یک گام به جلو، در مقیاس 3 ماهه، مدل استوکاستیک ایستگاه بوشهر (0/70=r و 0/66=RMSE) و در مقیاس 6 ماهه و 12 ماهه، مدل استوکاستیک ایستگاه همدان نوژه (مقیاس 6 ماهه: 0/84=r و 0/41=RMSE؛ مقیاس 12 ماهه: 0/93=r و 0/30=RMSE) دارای بیشترین دقت در مقایسه با مدل استوکاستیک سایر ایستگاه­ها می­باشد. همچنین با افزایش مقیاس زمانی، خطای پیش‌بینی به میزان قابل توجهی کاهش می‌یابد و با افزایش گام زمانی دقت پیش ­بینی کاهش می­یابد. نتایج ارزیابی دقت مدل‌ها در پیش‌بینی‌ طبقات SPI بر مبنای آماره­ی کاپا (K) نشان داد حداکثر توافق طبقات مشاهده ­شده و پیش ­بینی شده­ی یک گام به جلو در مورد SPI3، SPI6 و SPI12 به ترتیب مربوط به ایستگاه بوشهر (0/46=K)، ایستگاه گرگان (0/66=K) و ایستگاه زاهدان (0/81=K) می­باشد. همچنین با افزایش مقیاس زمانی، میزان توافق طبقات مشاهده ­شده و پیش ­بینی ­شده افزایش می­یابد و با افزایش گام زمانی میزان توافق طبقات مشاهده­ شده و پیش ­بینی­ شده کاهش می­یابد.

کلیدواژه‌ها


عنوان مقاله [English]

Evaluation of Six Types of Stochastic Models Skill in Modeling and Forecasting the Standardized Precipitation Index Time Series

نویسندگان [English]

  • Somayeh Hejabi 1
  • Javad Bazrafshan 2
چکیده [English]

Drought is a temporary and recurring meteorological event, which originates from the lack of precipitation relative to its long-term average. Since drought forecasting has a critical role in water resources management, in this study, six types of stochastic models (AR, MA, ARMA, non-seasonal ARIMA, seasonal ARIMA and multiplicative ARIMA) skill in modeling and forecasting the Standardized Precipitation Index (SPI) time series was evaluated. For reaching this purpose, the monthly total precipitation data related to ten synoptic stations with hyper humid to very dry climates (1973-2007) were used. At first, by transferring the optimum distribution cumulative probability of precipitation to cumulative probability distribution of standard normal, the SPI values in three time scales of 3, 6 and 12 months were calculated. Then, the development of models was done on SPI values related to period of 1973 to 2000, over a multi stage process (Identification, Parameter estimation, Diagnostic check) and the most appropriate stochastic model was determined for each time series from the candidate models. In order to validate the chosen models, at first, one to twelve lead times ahead forecasting was done for period of 2001 to 2007.Then, the values and the classes of the observed and the forecasted SPI were compared. The results of evaluating the models accuracy in forecasting the SPI values over chosen stations showed that in one month lead time ahead drought forecasting, over 3-month time scale, the stochastic model of Bushehr station (r=0.70, RMSE= 0.66) and over 6 and 12-month time scales, the stochastic model of Hamedan Nojeh station (6-month time scale: r=0.84, RMSE=0.41; 12-month time scale: r=0.93, RMSE=0.30) have the most accuracy comparing with stochastic model of other stations. Moreover, the forecasting error decreases with increasing the time scale and the forecasting accuracy decreases with increasing the lead time. The results of evaluating the models accuracy in forecasting the SPI classes based on Kappa statistic (K), showed that the maximum agreement of the observed and the forecasted classes for one lead time ahead forecasting, about 3, 6 and 12-month time scales is related to Bushehr (K=0.46), Gorgan (K=0.66) and Zahedan (K=0.81), respectively. Moreover, the agreement of the observed and the forecasted classes increases with increasing the time scale and the agreement of the observed and the forecasted classes decreases with increasing the lead time.

کلیدواژه‌ها [English]

  • Forecasting
  • Kappa statistic
  • Standardized Precipitation Index
  • Stochastic Models