تعداد نشریات | 44 |
تعداد شمارهها | 1,321 |
تعداد مقالات | 16,203 |
تعداد مشاهده مقاله | 52,830,935 |
تعداد دریافت فایل اصل مقاله | 15,507,544 |
تحلیل تغییرات بارش حوضه ی آبریز هراز با استفاده از روش های آماری و تکنیک تحلیل طیفی | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
هیدروژئومورفولوژی | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مقاله 4، دوره 3، شماره 7، شهریور 1395، صفحه 59-86 اصل مقاله (998.69 K) | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
نویسندگان | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
عباسعلی داداشی رودباری1؛ غلام عباس فلاح قالهری2؛ مختار کرمی3؛ محمد باعقیده4 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1دانشجوی کارشناسی ارشد آب و هواشناسی دانشگاه حکیم سبزواری | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
2استادیار اقلیم شناسی دانشگاه حکیم سبزواری | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3استادیار اقلیمشناسی دانشگاه حکیم سبزواری | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
4استادیار اقلیمشناسی دانشگاه حکیم سبزواری. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
چکیده | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
شناسایی رفتار بارش از مهمترین اصول برنامهریزیهای محیطی به شمار میرود. چرا که رفتار بارش در بلندمدت راه را برای برنامهریزیهای محلی و منطقهای هموار مینماید. پژوهش پیشرو با هدف پیشآگاهی مدیریت سیلاب در حوضهی آبریز هراز واقع در جنوب استان مازندران انجام شده است. بدین منظور از تحلیل روند و تکنیک تحلیل طیفی استفاده شده است. دادههای بارش مورد نیاز در بازهی زمانی 2007-1951 از پایگاه دادهی بارش آفرودیت اخذ شده است. نتایج نشان داد که ضریب تغییرات بارش در سطح حوضه در ماههای گرم سال نسبت به دورهی سرد سال بیشتر بوده است. همچنین با توجه به آمارههای پراکندگی به دست آمده، بزرگترین خوشههای بارش در حوضهی آبریز هراز در فصل زمستان ایجاد میشوند که نشاندهندهی نظم نسبی بارش در منطقهی مورد مطالعه میباشد. تحلیل طیفی مشاهدات بارش، حاکی از وجود چرخههای سینوسی معنیدار 2-3، 5-3، 5-11 و 11 سال به بالا در منطقهی مورد مطالعه است. وجود این چرخهها را میتوان به تأثیر عوامل کلانمقیاس جوی- اقیانوسی نسبت داد. شیب خط روند به دست آمده برای بارش سالانهی حوضهی آبریز هراز 67/7 میلیمتر بر دهه به دست آمد. مهمترین ساز و کاری که میتوان بر اساس آن با وجود روند افزایشی دما، مقادیر مثبت بارش را در حوضهی آبریز هراز توجیه نمود، ساز و کار «مرطوبترشدن گرمتر» است که بر اساس آن هر چه دمای سطحی منطقه بالاتر باشد، مقدار بارندگی بیشتر است. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
کلیدواژهها | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
آزمون من کندال؛ آزمون سنس؛ تحلیل طیفی؛ روند؛ همساز؛ حوضهی آبریز هراز | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
اصل مقاله | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مقدمه تغییر اقلیم از طریق ردیابی مشخصاتی نظیر کمیت و الگوی زمانی– کانی عناصر اقلیمی قابل مطالعه است. در این میان مطالعهی بارش از رواج زیادی برخوردار است. بارش به عنوان یکی از عناصر بنیادی اقلیم به ویژه طی قرن بیستم در عرضهای میانه و بالای نیمکره شمالی در حدود 1/0 درصد افزایش یافته است. از طرف دیگر در خشکیهای جنب حاره (عرضهای 30-10 درجه شمالی) روند بارش در حدود 3/0 درصد کاهش یافته و در نواحی حاره، افزایش نسبی بارش در حدود 2/0 تا 3/0 درصد بوده است (موسمان و همکاران[1]، 2004: 44). به تعبیری دیگر میتوان اذعان داشت که یکی از رفتارهای بارندگی که طی دهههای گذشته نمایانتر شده است، تغییر و جابجایی رژیم آن است که بـه صورت تغییر در مقدار بارنـدگی، زمان و تداوم رویداد بارش نمایان میشود (شیفته سومی و همکاران[2]، 2012: 9؛ دوهان و پاندی[3]، 2013: 138). عناصر آب و هوایی با تغییرات کمی و کیفی خود اوضاع محیطی را تحت تأثیر قرار میدهند. این اثرات هر چند که ممکن است بطئی و تدریجی باشند، لیکن میتوانند پیامدهای زیادی به دنبال داشته باشند. با شناخت ویژگیهای آب و هوائی میتوان قابلیتهای سازگاری انسان در محیط را ارزیابی کرد و متناسب با شرایط اقلیمی اشکال مطلوب ساختهای فضائی را با ابزار متناسب به وجود آورد (رزمی، 1389: 14). در یک نگاه منطقی، اقلیم و تغییرات آن بخصوص تغییرپذیری بارش میتواند سبب آسیبپذیری بیشتری به حیات بشر و همچنین بخشهای آب و کشاورزی باشد (مریانجی، 1391: 45). وجود دوره یا سالهایی که طی آن بارشی کمتر یا بیشتر از میانگین ثبت میشود، یک خصیصهی شناخته شده در اقلیم اکثر مناطق کرهی زمین میباشد. از آنجا که شرایط آب و هوایی هر منطقهای تابع متغیرهای زیادی است و این متغیرها هر ساله از شرایط و مقدار ثابتی برخوردار نیست، وجود نوسان در بارش مناطق مختلف را میتوان امری بدیهی و اجتنابناپذیر به حساب آورد. همچنین آب و هوای کرهی زمین در طول قرن بیستم به ویژه در دو دههی اخیر تعادل خود را از دست داده و تمایل به افزایش دما نشان داده است که به طور بالقوه میتواند منجر به تغییر چرخهی آب و بارش شود و ناهنجاریهای هیدرولوژیک همچون خشکیها و سیل را سبب شود (رزمی، 1389: 27). با توجه به مطرح شدن بحث تغییر اقلیم و تغییرات زمانی و مکانی عنصر بارش تا کنون از سوی محققان پژوهشهای متعددی با استفاده از روشهای پارامتری و نا پارامتری انجام شـده است که از جمله مـیتوان بـه پژوهشهای (کاهیا و کالایکی[4]، 2004؛ زانگ و همکاران[5]، 2006؛ لی یانگ و همکاران[6]، 2010؛ تانگ و همکاران[7]، 2011؛ شیفته سومئی و همکاران[8]، 2012؛ دوهان و پاندی[9]، 2013) اشاره کرد. توموزیو و همکاران[10] (2000) با استفاده از روش من-کندال و توابع تجربی متعامد[11]، تغییرپذیری میانگین بارندگی تابستانه را در امیلیا روماگنا[12] در ایتالیا بررسی کردند. مطالعات آنها نشان داد که همهی ایستگاهها دارای روند کاهشی در بارندگی تابستانه میباشند. گمر و همکاران[13] (2004) روند تغییرات بارش ماهانهی چین را مورد مطالعه قرار دادند و نشان دادند که تعداد، توزیع و کاهش روندها از ماهی به ماه دیگر متفاوت است. مارنگو و همکاران[14] (2004) الگوهای بلندمدت و دههای بارش در حوضهی آمازون را تجزیه و تحلیل کردند و به این نتیجه رسیدند که کل بارندگی در این منطقه دارای روندی منفی میباشد، اما این روند از توزیع مکانی ناهمگون برخوردار بوده است، برای مثال بخش شمالی حوضه دارای روندی منفی و بخش جنوبی دارای روندی مثبت میباشد. توشیج و انکاشویچ[15] (2005) با استفاده از روش تجزیه و تحلیل طیفی، تغییرات سریهای زمانی بارش ماهانه، فصلی و سالانه ایستگاه بلگراد را مورد ارزیابی قراردادند. نتایج نشان داد که تغییرات بارش با دورهی بازگشت دو ساله در ایستگاه بلگراد وجود دارد. در ایران نیز رزمی (1389) در پایاننامهی کارشناسی ارشد خود به ارزیابی تغییر رژیم بارندگی آذربایجان ایران پرداخت، روش اصلی به کار گرفته شده در پژوهش وی، روش تحلیل همسازها بوده است. بر اساس تحلیل همسازها مشخص گردید که دو سامانهی همدید- محلی تا بیش از 90 درصد بارش آذربایجان را توجیه میکنند. کمترین سهم بارش حاصل از سامانههای همدید در شمال شرق رخ میدهد. مریانجی (1391) در رساله دکتری خود تغییرپذیری بارش ایران را مورد بررسی قرار داد. نتایج پژوهش وی نشان داد که روند سالانه بارش در 54 درصد از مساحت کشور دارای روند مثبت با میانگین 72/1 میلیمتر بوده است. در نهایت نتایج تحلیل خوشهای نشان داد که تغییرپذیری بارش کشور در 5 گروه عمده قابل طبقهبندی میباشد. فرسادنیا و همکاران (1391) به بررسی تحلیل روند بارندگی در استان مازندران با استفاده از روش من-کندال منطقهای پرداختند. نتایج پژوهش آنان نشان داد که بررسی وجود روند در سری زمانی میانگین بارش سالانه در کل استان مازندران دارای روند صعودی در سطح معنیدار 1% است. در پژوهش حاضر ضمن بررسی ویژگیهای عمومی بارش حوضهی آبریز هراز سعی شده است روند دادههای بارش و دورهی بازگشت بارش در مقیاسهای ماهیانه، سالانه و فصلی نیز مورد ارزیابی قرار گیرد تا بتوان بر اساس دورههای بازگشت و سامانههای مؤثر در بارش منطقه، پیشآگاهی مناسبی به جهت مدیریت سیلاب ارائه گردد. مواد و روشها منطقهی مورد مطالعه در پژوهش حاضر، حوضهی آبریز هراز در جنوب استان مازندران میباشد. این منطقه از نظر تقسیمات سیاسی جزو شهرستان آمل میباشد. وسعت حوضهی آبریز هراز، 200051 هکتار است. محیط حوضهی هراز بـرابـر با 5/296062 کیلومتر، حداقل ارتفاع حوضه 200 متر و حداکثر ارتفاع آن 5600 متر میباشد. تـمامی رودخانههای حوضه به رودخانهی هراز منتهی مـیگردد. از مهمترین رودخانههای حوضه میتوان بـه رودخانههای آخنسر، شیرکلارود و نمارستاق اشاره کرد. از ارتفاعات مهم حوضهی کوههای دماوند که مرتفعترین منطقه در حوضه به شمار میآید و نیز شیمکوه و کوه امامزاده قاسم اشاره نمود. مراکز مهم مسکونی حوضهی آبریز هراز عبارتاند از پلور، نشل، تیران، رینه، کندوا، آب اسک، گزنک، بایجان، بلقلم و... اشاره کرد. شکل (1) موقعیت مطالعاتی هراز را نسبت به استان مازندران و کل کشور نشان میدهد.
شکل (1) موقعیت منطقهی مورد مطالعه جنوب استان مازندران سرزمینی ناهموار است و همین عامل به همراه چندین عامل دیگر موجب تغییرپذیری بارش در استان مازندران و بخصوص مناطق جنوبی آن میشود. با توجه به ویژگیهای اقلیمی خاص استان مازندران و اهمیت بارش در مدیریت منابع آب نیاز به مطالعات در مورد بارش و رفتارهای آن به شدت احساس میشود. اما در سطح استان مازندران دسترسی به دادههای بلندمدت بارش با گسترهی وسیع مسیر نیست. هم اکنون 15 ایستگاه همدید در استان مازندران وجود دارد که کمتر از یک سوم آنها آمار بلندمدت (30 ساله و یا بیشتر) متغیرهای هواشناسی را در اختیار دارند که با توجه به گسترهی وسیع استان، این تعداد ایستگاه برای بررسی تغییرات بارش که خود سرکشترین عنصر اقلیمی به شمار میرود و همچنین دارای نوسان مکانی - زمانی زیاد است؛ کافی نیست. افزون بر مشکلاتی که جهت دسترسی به دادههای بلند مدت اشاره شد؛ بخش اعظم منطقهی مورد مطالعه را مناطق مرتفع تشکیل میدهد (قرارگیری بخشی از مناطق شرقی قلهی دماوند و ارتفاعات شمالی و جنوبی آن) و با توجه به اینکه در نواحی با ارتفاع بیش از 2600 متر (کانونهای آبگیر کشور)، فاقد ایستگاههای سنجش و اندازهگیری بارش هستیم. همچنین نواحی مجاور که از تباین اقلیمی بالایی برخوردارند (دامنههای شمالی- جنوبی البرز) از توزیع مناسب ایستگاه برخوردار نیستند (عساکره، 1387: 25). در این پژوهش با استفاده از بروندادهای یاختهای بارش پایگاه آفرودیت[16]، به ارزیابی تغییرات زمانی- مکانی بارش ماهیانهی حوضهی آبریز هراز پرداخته شده است. دادههای این پایگاه با بازهی مکانـی خاورمیانه و جـنوب شرق آسیا و بازهی زمانی 1/1/1951 تا 31/12/2007 مـیلادی و تفکیک مکانی 25/0×25/0 و 5/0×5/0 درجـه در تـارنمای ایـنترنتی (http://www.chikyu.ac.jp/precip/index.html) موجود میباشند (یاتاگانی و همکاران[17]، 2012: 1407). در مطالعهی حاضر دادههای منطقهی خاورمیانه (APHRO_ME) از آخرین محصول پایگاه داده APHRODITE تحت عنوان v1101، با قدرت تفکیک 25/0×25/0 و فرمت ".nc" NetCDF از سایت مربوطه اخذ گردید. لازم به ذکر است کنترل کیفی بر روی دادههای فوق انجام شده است. در مرحلهی بعد با استفاده از تکنیکهای زمین آمار و روش کریجینگ، یاختهای به ابعاد 15×15 کیلومتر مربع بر منطقهی مورد مطالعه گسترانیده شد. حال به جای انجام تحلیل بر روی گرهگاهای 25/0×25/0 درجه و انتساب نتایج به دست آمده به پهنهی مورد مطالعه؛ از دادههای ریزمقیاس نمایی شده با تفکیک فضایی 15×15 کیلومتر استفاده شده است.
ارزیابی روند دادهها جهت برآورد شیب واقعی یک روند در یک سری زمانی، استفاده از روش ناپارامتریک سنس[18] میتواند یکی از روشهای مناسب در این زمینه باشد. این روش ابتدا توسط تیل در سال (1950) ارائه و سپس توسط سن در سال (1968) بسط و گسترش داده شد. روش حاضر همانند بسیاری دیگر از روشهای ناپارامتریک همچون من-کندال بر تحلیل تفاوت بین مشاهدادت سری زمانی استوار است. این روش زمانی میتواند مورد استفاده قرار گیرد که روند موجود در سری زمانی یک روند خطی باشد. این بدین معناست که در رابطهی (1) برابر است با:
در رابطهی (1) ، شیبخط روند و ، مقدار ثابت است. جهت محاسبه شیب خط روند یعنی، ابتدا بایستی شیب بین هر جفت داده مشاهداتی، با استفاده از رابطهی (2) محاسبه گردد:
در رابطهی (2) است. در این رابطه و به ترتیب دادههای مشاهداتی در زمانهای و است. با اعمال این رابطه برای هر جفت داده مشاهدهای، یک شیب به دست میآید. با قرار دادن این شیبها در کنار یکدیگر یک سری زمانی از شیبهای محاسبه شده حاصل میآید. یعنی اگر تعداد در سری زمانی وجود داشته باشد (در این تحقیق) ما به اندازهی برآورد شیب خواهیم داشت. در مرحلهی بعد سری زمانی مورد مطالعه به دست میآید. برای این کار تعداد از ها از کوچک به بزرگ مرتب میشوند و سپس با استفاده از یکی از رابطههای زیر اقدام به تعیین میانهی سری زمانی میگردد. اگر تعداد مشاهدات سری زمانی مورد مطالعه فرد باشد، از رابطهی (3) و اگر زوج باشد از رابطهی (4) استفاده میگردد (ویوکاناندان[19]، 2007: 104):
نتیجهی حاصل از این رابطه، به دست آمدن شیبخط روند است. اگر شیب خط روند مثبت باشد، حاکی از صعودی بودن روند و اگر منفی باشد، دال بر نزولی بودن روند است. در گام بعدی، آزمون نمودن شیب به دست آمده در فاصله اطمینان 95 درصد است. جهت انجام این آزمون از رابطهی زیر استفاده میگردد:
در رابطهی (5) عبارت است از آماره توزیع نرمال استاندارد در یک آزمون دوطرفه که برای سطح اطمینان 95 درصد برابر با است و نیز واریانس پارامتر است. جهت به دست آوردن مقدار پارامتر و همچنین مراحل زیر بایستی پیموده شود: الف) محاسبه اختلاف بین تک تک جملات سری با همدیگر و اعمال تابع و استخراج پارامتر :
در رابطهی (6) تعداد مشاهدات سری (در این تحقیق برابر با 57 سال)، و نیز به ترتیب دادههایام و ام سری میباشند. ب) محاسبهی تابع علامت[20] (ریچارد[21]، 1987: 25) که به صورت تابع 7 قابل محاسبه است:
ج) محاسبه واریانس توسط یکی از روابط زیر. اگر تعداد دادههای سری زمانی بزرگتر از 10 مورد باشد از رابطهی 8 و اگر کوچکتر از 10 باشد، از رابطهی 9 بهره گرفته میشود:
در رابـطهی 8 و 9، تعداد دادههای مشاهداتی، تعداد سریهائی است که در آنها حداقل یک دادهی تکراری وجود دارد و نیز بیانگر فراوانی دادههای با ارزش یکسان میباشد و در نهایت حدود اعتماد بالا و پایین به کمک روابط زیر محاسبه میگردد:
در رابطهی (10) تعداد شیبهائی است که از طریق رابطهی (2) به دست آمده است. حال امین و امین شیبها را از بین شیبهای محاسبه شده استخراج مینمائیم. در صورتی که عدد صفر در دامنهی بین دو شیب استخراج شده فوق قرار گیرد، فرض صفر پذیرفته شده و عدم وجود روند در سری دادهها تأیید میگردد. در غیر این صورت فرض صفر رد شده و وجود روند در سطح اطمینان مورد آزمون پذیرفته میگردد. در نهایت جهت به دست آوردن مقداردر رابطهی 1، تعداد تفاضل محاسبه میشود. سپس میانه این مقادیر، برآورد را در اختیار میگذارد. در پژوهش حاضر برای ارزیابی روند دادههای بارش حوضهی آبریز هراز، از آزمون ناپارامتریک من-کندال[22] نیز استفاده شده است. این آزمون ابتدا توسط من (1945) ارائه و سپس توسط کندال (1975) تکامل یافت. این روش نیز همانند آزمون شیب خط سنس در مباحث علوم محیطی به طور گستردهای کاربرد دارد. در آزمون من-کندال نمره از رابطهی 11، محاسبه میگردد:
در رابطه (11) نشانهی تفاوت مقادیر با یکدیگر و واریانس است:
که تعداد مشاهدات سری (در این تحقیق برابر با 57 سال)، و نیز به ترتیب دادههایام و ام سری میباشند. تابع علامت است که شرح آن در رابطهی (7) آمده است. عاملی مربوط به تصحیح پراش است. در صورتیکه دادههای تکراری در اطلاعات وجود داشته باشد از رابطهی (14) استفاده میشود:
در رابطهی (14) تعداد دادههای مشاهده شده و معرف تعداد سریهایی است که در آنها حداقل یک دادهی تکراری وجود دارد. آمارههای پراکندگی در پژوهش حاضر برای بررسی ساختار مکانی یاختههای بارش حوضهی آبریز هراز، از آمارههای پراکندگی استفاده شده است. آمارههای پراکندگی استفاده شده عبارتاند از: شاخص پراکندگی[23]، شاخص تجمع، شاخص گرین، شاخص فرکانس خوشه، شاخص تراکم میانگین، شاخص لکهها و شاخص موریسیتا. شاخص پراکندگی نسبت واریانس به میانگین را بررسی میکند. اگر دادههای مورد مطالعه از توزیع تصادفی پیروی کنند انتظار میرود ID برابر یک باشد. این آماره از توزیعبا درجه آزادی پیروی میکند (رابطهی 15 را ملاحظه نمایید). شاخص تجمع[24] (دیوید و مور[25]، 1954: 50) تابع مستقیم از شاخص پراکندگی است. اگر دادههای مورد مطالعه از توزیع تصادفی پیروی کنند، انتظار میرود این شاخص برابر صفر باشد. مقادیر ارزشی مثبت در این آماره حاکی از خوشهای بودن دادهها و مقادیر ارزشی منفی بیانگر ناخوشه بودن دادهها میباشد (رابطهی 16). آماره گرین[26] (گرین[27]، 1966) در حقیقت اصلاح شده شاخص تجمع میباشد که مستقل از N است. در این آماره ارزش عددی صفر برای توزیعهای تصادفی و مقادیر ارزشی یک برای توزیعهای با حداکثر خوشه است (رابطهی 17):
شاخص فرکانس خوشه[28] (داگلاس[29]، 1975: 403) معیاری برای اندازهگیری خوشه بر اساس تابع K با استفاده از توزیع دوجملهای منفی است (رابطهی 18). شاخص تراکم میانگین[30] (لوید[31]، 1967: 9) متوسط تعداد نقاط موجود در منطقهی مورد مطالعه است که به طور تصادفی از یک نقطه سنجش میشود (رابطهی 19):
شاخص لکهها[32] (لوید، 1967: 14) شدت الگوی تأثیرپذیر را مورد مطالعه قرار میدهد (رابطهی20). شاخص موریسیتا[33] (موریسیتا[34]، 1959: 13) احتمال مقیاسبندی خوشههای دو یاخته از کل یاختههای مورد مطالعه را مورد ارزیابی قرار میدهد (رابطهی 21):
علاوه بر موارد فوق، ضریب چولگی (g1) و ضریب درجهی اوج (g2) نیز برای هر یک از ماههای مورد مطالعه به طور جداگانه با استفاده از رابطهی 22 محاسبه شده است:
تکنیک تحلیل طیفی (تحلیل سریهای فوریه[35]) برای بـررسی چرخههای بارش حوضهی آبـریز هراز از روش تحلیل طیفی استفاده شد. در تـکنیک تحلیل طیفی ابتدا سریهای زمانی به تابع فرکانسی (به صورت تابعی دورهای[36] با دامنه[37] و فراوانی[38]) تبدیل شد. روشهای پرشمار و متنوعی برای تبدیل سریهای بسامد وجود دارد (عساکره و رزمی، 1391: 152). یکی از روشهای ساده و کارآ، روش همسازهاست. در این روش سری زمانی به طول ، به فرم یک مدل فوریه و به شکل رابطهی (23) نوشته میشود (ماسلن و روکمور[39]، 1997: 200).
که در آن و ضرایب فوریه خوانده میشوند و از طریق روش کمترین مربعات خطا قابل برآورد هستند. و فراوانی تکرار مشاهدات (عکس دوره بازگشت) است و با نشان داده میشود. واریانس برای فرکانس و طول دوره آماری از رابطهی (24) به دست میآید:
نموداری که مقادیر را در مقابل نشان میدهد، دوره نگار[40] نامیده میشد. از آنجا که دوره نگار نسبتی از توزیع واریانس در فرکانسهای بنیادی و برآوردی پر افت و خیز از طیف است، برای کاهش واریانس، اریبی و ایستاسازی این برآورد، باید آن را در مجاورت فرکانس مورد نظر و به وسیلهی صافیهای مناسب هموار ساخت. به عبارت دیگر دوره نگار هموار شده، برآوردی از طیف جامعه است که حاصل میانگین موزن مقدار راست و چپ فرکانس مورد نظر میباشد (گلدسته، 1377: 245-246). یکی از روشهای ساده و مناسب برای هموارسازی دادههای سالانه به وسیلهی میتچل و همکاران[41] (1966: 34) معرفی شد. در این رویه، دوره نگار با یک میانگین متحرک موزون سهجملهای هموار میشود.
بحث و نتایج مشخصات آماری بارش حوضهی آبریز هراز جدول (1) بعضی از مشخصات آماری بارش را طی ماههای مختلف نشان میدهد. مطابق این جدول بیشترین میانـگین بارش با 99/54 مـیلیمتر در ماه مارس رخ داده است. کـمترین مـقدار بـارش بـه میزان 46/16 میلیمتری در ژوئیه اتفاق افتاده است. این در حالی است که بیشترین ضریب تغییرات بارش نیز با 66/102 در ماه ژوئیه اتفاق افتاده است. چولگی و توزیع ضریب چولگی (g1) بارش در همهی ماهها مثبت بوده است. سامانههای بارشزا با توجه به شرایط دینامیک و ترمودینامیک و بسته به موقعیت جغرافیایی خود، در برخورد با شرایط متفاوت محلی میتوانند نواحی بارشی متفاوتی ایجاد کنند (غیور و همکاران، 1390: 19). لذا مقدار بارش دارای شاخصهای آماری متفاوتی خواهد بود. همانطور که در جدول (1) نیز مشخص میباشد، بیشترین مقدار عددی شاخص g2 (ضریب درجه اوج) مربوط به ماه ژوئیه (26/10) و کمترین مقدار آن متعلق به ماه سپتامبر (32/5) بوده است. ضریب درجهی اوج نشاندهندهی تفاوت مکانی بارش میباشد. دلیل بالاتر بودن شاخص درجه اوج در ماه ژوئیه ناشی از تفاوت مکانی بالای بارش میباشد. علت کاهش شدید شاخص با شروع فصل پاییز را میتوان با فرارفت رطوبت بیشتر همزمان با تقویت و نفوذ زبانههای پرفشار سیبری نسبت داد که با مطالعات علیجانی (1389) نیز همخوانی دارد. همرفت محلی در منطقهی مورد مطالعه در فصول انتقالی که رطوبت حاصل از سامانههای مدیترانهای هنوز در محیط وجود دارد، نیز میتواند در تغییرات مکانی بارش منطقه موثر باشد. اطلاعات کامل مشخصههای آماری هر یک از ماههای سال در دورهی آماری 57 ساله در جدول (1) آمده است. آمارههای پراکندگی بارش در حوضهی آبریز هراز جدول (2) مشخصات آمارههای پراکندگی بارش را طی ماههای مختلف نشان میدهد. بر اساس شاخص ID ماههای ژوئیه و اکتبر بالاترین مقادیر عددی را نشان دادهاند. این شاخص بهخوبی ناهنجاری مکانی بارش حوضهی مورد مطالعه را نشان داده است، به طوری که در ماههای نامبرده ایستگاه کره سنگ به نمایندگی از مناطق شمالی حوضه در مقایسه با ایستگاه پلور به نمایندگی از مناطق جنوبی حوضه در یک فاصله کمتر از 60 کیلومتری، اختلاف 10 میلیمتری دارند. همانگونه که در جدول (2) نشان داده شده است، شاخص ICS برای هر 12 ماه مورد مطالعه مثبت و بالای 4 بوده است که نشاندهندهی الگوی خوشهای بالای بارش در حوضه میباشد. شاخصهای GI تبیینکنندهی همین امر میباشد. شاخص ICF نیز اندازهی خوشه را نشان میدهد. بر اساس این شاخص در فصل زمستان بزرگترین خوشههای بارش در حوضهی آبریز هراز ایجاد میشوند که نشاندهندهی نظم نسبی بارش در منطقهی مورد مطالعه میباشد. شاخصهای IP و GI نشاندهندهی تفاوت مکانی خوشههای بارش میباشند که فصل تابستان بالاترین مقادیر تفاوت مکانی خوشهای را به خود اختصاص داده است. شاخص MI نیز مقدار تعدیلشده شاخص IPمیباشد که مبین تفاوت ساختار مکانی بارش میباشد. جدول (1) مشخصات آماری بارش برای هر یک ماههای سال در حوضهی آبریز هراز
جدول (2) مشخصات آمارههای پراکندگی بارش حوضهی آبریز هراز
ارزیابی روند دادهها از آنجا که مبانی زندگی انسان با بارش هماهنگ شده است؛ وجود روند کاهشی یا افزایشی در بارش به معنای تغییر شکل زندگی خواهد بود. لذا جهت شناسایی روندهای کاهشی یا افزایشی در بارشهای حوضهی آبریز هراز از دو روش ناپارامتریک برآوردکنندهی شیب سنس و من-کندال استفاده گردید. نتایج حاصل از این روشها نشان داد هر چند نشانههایی از روند کاهشی (سه ماه فوریه، مارس و آوریل در آزمون من-کندال و شیب سنس) و افزایشی (9 ماه دیگر سال) در بارشهای حوضه مشاهده میشود، اما از بین این روندهای کاهشی و افزایشی تنها دو ماه سپتامبر و اکتبر از نظر آماری در سطح معنیدار 5 درصد معنیدار بودهاند. اطلاعات هر یک از روندها در سطوح 5 و 1 درصد در جدول (3) آمده است.
جدول (3) نتایج حاصل از آزمون سنس و من-کندال در سطوح احتمالاتی 5 و 1 درصد برای بارش ماهانه حوضهی آبریز هراز
نتایج به دست آمده نشان میدهد یاختههایی که روند کاهشی در آنها مشاهده شده است، همگی در بخشی جنوبی حوضه و اطراف قله دماوند واقع شدهاند که میتواند نشانههایی از تغییر شرایط بارش در این قسمت از ایران باشد. هر چند اثبات وجود روند معنیدار در یک سری زمانی بارندگی نمیتواند به تنهایی دلیلی قاطع بر وقوع تغییر اقلیم در یک منطقه باشد، اما فرض رخداد آن را میتواند تقویت کند. لذا همین تغییرات رخ داده میتواند حائز اهمیت باشد. به طور کلی بارش حوضهی آبریز هراز در مقیاس سالیانه بر اساس نتایج آزمون شیب سنس افزایش 67/7 میلیمتری در هر دهه داشته است. بر این اساس در طول دورهی 57 ساله به رقم باور عمومی مبنی بر کاهش بارش در سطح کشور و منطقه، بارش سالانه حوضهی آبریز هراز افزایش 719/43 میلیمتری داشته است. بر این اساس میتوان برای روند تغییرات بارش در سطح حوضهی آبریز هراز رابطهای به شکل ارایه داد. شکل (1) برازش نمودار خط سنس بر سری زمانی دادههای میانگین یاختهای بارش حوضهی آبریز هراز را در دورهی زمانی مورد مطالعه نشان میدهد.
شکل (1) برازش نمودار شیب سنس بر سری زمانی دادههای میانگین یاختهای بارش حوضهی آبریز هراز
در نگاه اول به نظر میرسد نتایج حاضر دور از واقعیت است، چرا که با افزایش دمای متوسط جهانی باور عمومی مبنی بر کاهش مقدار نزولات جوی است، اما نتایج پژوهش حاضر نشان داده است که مقدار بارش در 57 سال اخیر در حوضهی آبریز هراز افزایش یافته است. بررسیها نشان داده است که این افزایش بارش بیشتر به شکل رگباری بوده و بارشهای آخر فصل سرد و فصول گرم سال در حال افزایش است. چنین افزایشی نه تنها سودمند نیست، بلکه موجب وقوع رخداد وقایع فرین اقلیمی همچون سیل شده است. روند مشاهده شده در حوضهی آبریز هراز به خوبی با گزارشهای سازمانهای مختلف حاکی از تغییرات جدی در اقلیم همخوانی دارد؛ به طوری که بر اساس گزارشهای ارایه شده از تغییرات آب و هوایی (IPCC[42],2013: 18)، از نیمهی دوم قرن بیستم روند بارش سنگین (بارشهای رگباری) در نـقاط مختلف جهان رو به افزایش است. نتایج سایر مطالعات در سطح جهانی تأییدی بر نتایج حاضر است، زیرا بررسی دادههای بارندگی بیش از ۸ هزار ایستگاه هواشناسی در سطح جهان و در فاصلهی زمانی بین سالهای ۱۹۰۰ تا ۲۰۰۹ نشان داده است که در این دورهی زمانی گرمای جو زمین به اندازهی یک درجهی سانتیگراد زیاد شده و در نتیجه آن تعداد بارندگیهای شدید نیز افزایش یافته است (وسترا و همکاران[43]، 2013: 3904). ساز و کار این امر را نیز میتوان به گونهای تشریح کرد که با افزایش گازهای گلخانهای و افزایش گرما، مقدار رطوبت بیشتری در جو میماند و بارشها شدیدتر میشود. به ازای هر یک درجه گرمایش جو زمین، شدت باران ۷ درصد افزایش مییابد و انتظار میرود دمای کرهی زمین در انتهای قرن جاری ۳ تا ۵ درجه سانتیگراد افزایش یابد و تغییرات آب و هوا شدت بارندگی انتهای قرن را نیز به میزان قابل توجهی افزایش دهد (همان، 3916). مهمترین ساز و کاری که میتوان بر اساس آن به تبع وجود روند افزایشی دما، مقادیر مثبت بارش را در حوضهی آبریز هراز توجیه نمود، مکانیسم «مرطوبترشدن گرمتر[44]» (تن و همکاران[45]، 2015: 451) است که بر اساس آن هر چه دمای سطحی منطقه بالاتر باشد، مقدار بارندگی بیشتر است. در دهههای اخیر روند چشمگیری در بارش اکثر ماههای سال مشاهده نشده است و تنها شاهد تغییر و تنوع سالانه بارندگیها هستیم و حتی در برخی از موارد با گرمتر بودن سالها بارندگی نیز بیشتر صورت گرفته است. در توجیه این تغییر و تنوع سالانه میتوان به نقش بارز پدیدهی النینو اشاره کرد. تلفیق پدیدهی گرمایش زمین و پدیدهی النینو و خروجی آن دو تحت مدلهای عددی در سطح جهان نشان داده است که، بارندگی در مناطق مرطوب، بیشتر و در مناطق خشک، کمتر میشود (وسترا و همکاران[46]، 2013: 3916). این تأییدی بر مکانیسم مرطوبتر شدن مرطوب[47] است که در حوضهی آبریز هراز با توجه به مرطوب بودن منطقه قابل تعمیم است. ادغام ساز و کار «مرطوبترشدن گرمتر» و «مرطوبتر شدن مرطوب» حدود ۸۰ درصد تغییرات بارندگی ناشی از گرمایش زمین را توضیح میدهد که در حوضهی آبریز هراز نیز مکانیسم مرطوبتر شدن گرمتر و مرطوبتر شدن مرطوب مشاهده گردید. تحلیل همسازهای بارش حوضهی آبریز هراز برای تحلیل دقیقتر تغییرات بارش حوضهی آبریز هراز چرخههای حاکم در هر ماه و فصل با استفاده از تحلیل طیفی مورد بررسی قرار گرفت. اساس روشهای تحلیل طیفی، تقسیمبندی تغییرپذیری سریهای زمانی به اجزاء یا بخشهای حاوی دورههای مختلف یا فرکانس است که به عنوان یکی از روشهای استخراج و تحلیل نوسانات آشکار و نهان طول موجهای مختلف، مطرح میباشد. جداول (4 تا 9) دوره نگار، طیف و فاصله اطمینان را برای سریهای زمانی بارش حوضهی آبریز هراز در دورههای مختلف ارائه دادهاند.
جدول (4) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 99 درصد در مقیاس ماهانه
جدول (5) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 99 درصد در مقیاس سالانه و فصلی
تغییرات چرخههای بارش نشان میدهد که در فصول تابستان و پاییز در سطح اطمینان 99 درصد چرخههای 2 ساله بر بارش حوضهی مورد مطالعه حاکم بوده است این در حالی است که در فصول بهار و زمستان بارش در حوضهی آبریز هراز هیچگونه چرخهای نداشته است. بررسی چرخههای بارش برای متوسط سالیانه نیز حاکی از چرخه کوتاهمدت 2 ساله بوده است. تحلیل همسازهای بارش در حوضهی آبریز هراز در سطوح اطمینان 95 و 90 درصد نشان داد که علاوه بر چرخههای 2 ساله، چرخههای بلندمدت 5، 8، 14، 25 و 28 ساله هم بر بارش این منطقه از کشور حاکم بوده است. در مجموع میتوان گفت که بارش هراز در شمارهی چرخههای 5 و 8 با احتمال 49/0 دارای دورهی بازگشت کوتاهمدت 2 ساله است. جدول (6) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 95 درصد در ماههای مختلف
جدول (7) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 95 درصد در مقیاس فصلی و سالانه
در بین نتایج ارایه شده آنچه بیش از همه خودنمایی مینماید، چرخههای 2-3 ساله هستند. بررسیهای یاختههای مورد نظر نشان میدهد که غالباً این چرخهها در ارتفاعات دماوند و دامنههای مرتفع آن قابل رؤیت هستند. بیشتر دانشمندان این چرخهها را به ال نینو- نوسانات جنوبی (ENSO[48])، تغییرات دوسالانه (QBO[49]) الگوی بزرگمقیاس گردش عمومی جو، جریانات مداری و سایر فرایندهای جوی- اقیانوسی نسبت دادهاند. برای مثال، کانه و تگزیرا[50] (1991: 89) چرخههای 2-3 ساله در بارش ماساچوست، هارتمن و همکاران[51] (2002: 161) چرخههای 2-3 ساله را در چین و لانا و همکاران[52] (2005: 193) چرخههای 6/4 و 1/2 ساله بارش در ایستگاه فابرا[53] در شمال شرق اسپانیا را به تغییرات دوسالانه (QBO) و چرخههای 2/9 و 5/5 ساله را به نوسانات اطلس شمالی (NAO[54]) و چرخههای 8/11 ساله را به فعالیت لکههای خورشیدی نسبت دادهاند.
جدول (8) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح90درصد در ماههای مختلف
جدول (9) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 90 درصد در مقیاس فصلی و سالانه
تورنس و وبستر[55] (1999: 2679) چرخههای 2-7 ساله را در بارش موسمی هند، گارسیا[56] و همکاران (2002: 78) چرخههای 7/2 ساله را در بارش شبهجزیره ایبری و آزاد و همکاران (2010: 2290) چرخههای 3-5 ساله را در بارش موسمی هند ناشی از تأثیر پدیده انسو دانستهاند. سلوام و جلوشی[57] (1995: 613) چرخههای 7-2 ساله و چرخههای 3-2 ساله را در دمای سطح زمین کشف کردند. ایشان به ترتیب این چرخهها را به انسو و تغییرات دوسالانه نسبت دادهاند. همانطور که مشاهده شد، چرخههای کوتاهمدت 2 تا 4 ساله و 3-5 ساله بیشترین حاکمیت را بر بارش حوضهی آبریز هراز داشتهاند که این نوع چرخهها را میتوان به رخداد النینو نسبت داد. بیشتر پژوهشگران چرخههای 2 تا 4 ساله را به النینو- نوسانات جنوبی (ENSO) (هولرینگ و کومار، 2003؛ اسچیبورت و همکاران[58]، 2009) و تغییرات دوسالانه (QBO) الگوی بزرگمقیاس گردش عمومی جو و جریانات مداری و سایر فرایندهای جوی- اقیانوسی نسبت دادهاند (گارسیا[59] و همکاران، 2002: 77-95، هارتمن[60] و همکاران، 2008:156 ، کوک و همکاران[61]، 2011a,b). بررسی چرخههای بارش طی دورههای مورد مطالعه (2007-1951 میلادی) در سطوح اطمینان سهگانهی 90، 95 و 99 درصد نشاندهندهی چندین الگوی مختلف بارش میباشد که یکی از آنها، چرخههای 11-5 ساله میباشد. بر این اساس در سطح 99 درصد ماه می با دوره بازگشت 5 ساله؛ در سطح 95 درصد ماه ژانویه با دورهی بازگشت 9 ساله، اکتبر با دوره بازگشت 6 ساله و نوامبر با دورهی بازگشت 8 ساله و در سطح 90 درصد ماه آوریل دارای دورههای بازگشت 8 و 5 ساله، اکتبر دارای دورهی بازگشت 6 ساله و نوامبر دارای دورهی بازگشت 8 ساله میباشد. تغییرپذیری درون دههای (دهه به دهه) فقط برای ماهها و یا فصولی که چرخههای 11 سال به بالا دارند، قابلمشاهده است (مانند ماههای ژولای، اوت و سپتامبر در سطح اطمینان 99 درصد). مـطالعهای که توسط جهانـبخش و عدالتدوست (1387) در شـمال غرب ایران انـجام شده است نشان داده است که چرخههای نامبرده مبین تأثیر فعالیت لکه (کلف)های خورشیدی است که با نتایج پژوهش حاضر نیز همخوانی دارد. نتایج نشان میدهند نوسانات اطلس شمالی تابعی از عملکرد لکههای خورشیدی است. همچنین بررسیها مؤید این امر است که چرخههای غیرسینوسی (روند) نیز در سری زمانی بارش منطقهی مورد مطالعه با توجه به مشاهده چرخههای درون دههای 20 تا 40 ساله در پهنهی مورد مطالعه قابل استنباط است. ارزیابی چرخههای بارش حوضهی آبریز هراز نشان داد که مقدار بارش تحت تأثیر نوسانات کوتاهمدت است؛ همچنین مشخص شد که تغییرات غیرخطی نوسانات بارش به عنوان تابعی از جریان جو، میزان بارش را تحت تأثیر قرار میدهد. تغییرات بارش حوضهی آبریز هراز همچون تغییرات بارش تمامی پهنهی سیاره از الگوی زمانی- مکانی متنوعی پیروی میکند. بررسیهای چرخهها نشان داده است پدیدههای النینو- نوسانات جنوبی و تغییرات دوسالانه بیشترین تأثیر را بر بارشهای جنوب حوضهی آبریز هراز و به طور کلی جنوب استان مازندران داشتهاند که با نتایج مریانجی (1391) همخوانی دارد. همچنین مشخصشده است ترابری رطوبت و کوههای البرز به عنوان مهمترین ساز و کار بارش و تغییرات آن در منطقه است. نتیجهگیری به طور کلی، تحلیلهای آماری و مکانی- زمانی رویدادهای بارش ابزار مناسبی هستند که میتوانند شناخت خوبی از شرایط بارش منطقه ارائه کنند. برخی از یافتههای اصلی پژوهش حاضر بدین شرح است: علیرغم وجود متغیرهای بارش متفاوت در نواحی مختلف، توزیع چگالی احتمال فراوانی رویدادهای بارش در همه ماههای سال در منطقهی مورد مطالعه، دارای چولگی مثبت و کشیدگی بیشتر از نرمال است. به عبارت دیگر، تراکم نسبی رویدادها پیرامون مرکز توزیع، بیشتر است. نتایج نشان داد که وضعیت ضریب تغییرات بارش در سطح حوضه در ماههای گرم سال نسبت به دورهی سرد سال بیشتر بوده است، به طوری که ژوئیه در فصل تابستان بیشترین ضریب تغییرات بارش را با 66/102 درصد داشته است. تحلیل ماهانهی فراوانی رویدادهای بارش در این منطقه نشان میدهد که شاخص پراکندگی بارش از ماه ژانویه تا آوریل یک روند افزایشی و در ماههای دیگر (به استثنای ژوئیه و اکتبر) روندی کاهشی دارد که همین شرایط برای واریانس بارش نیز مشاهده شده است که به خوبی ناهنجاری بالای بارش در سطح حوضه را در طول سال نمایان ساخته است. نتایج شاخص ICF که اندازه خوشههای بارش را در سطح حوضه نشان میدهد که در فصل زمستان بزرگترین خوشههای بارش در حوضهی آبریز هراز ایجاد میشوند که نشاندهندهی نظم نسبی بارش در منطقهی مورد مطالعه میباشد. نتایج حاصل از روشهای ارزیابی روند با دو آزمون شیب سنس و من-کندال نشان داد که هر چند نشانههایی از روند کاهشی (سه ماه فوریه، مارس و آوریل در آزمون من-کندال و شیب سنس) و افزایشی (9 ماه دیگر سال) در بارشهای حوضه مشاهده میشود، اما از بین این روندهای کاهشی و افزایشی، تنها دو ماه سپتامبر و اکتبر از نظر آماری در سطح معنیدار 5 درصد معنادار بودهاند. به طور کلی بارش حوضهی آبریز هراز در مقیاس سالانه بر اساس نتایج آزمون شیب سنس افزایش 67/7 میلیمتری در دهه داشته است. بر این اساس در طول دورهی 57 ساله، بارش سالانه حوضهی آبریز هراز افزایش 719/43 میلیمتری داشته است. از آنجا که سه ماه فوریه، مارس و آوریل از پربارش ترین ماههای سال هستند و روند کاهشی داشتهاند، به معنی کاهش بارش فصل سرد و افزایش بارش فصل گرم در منطقهی مورد مطالعه است و همین خود دلیل قوی بر تغییر نامطلوب بارش در حوضهی مورد مطالعه است که میتواند خطرات جدی برای فرسایش، سیلخیزی و مدیریت منابع آب در پی داشته باشد. تحلیل طیفی بارش در حوضهی آبریز هراز نشان داد که بارش در این پهنه دارای چرخههای نوسانی معنیدار 3-2، 5-3، 11-5 و 11 سال به بالاست. چرخههای 11-5 ساله بیشترین گستره را در برداشته است. ورود این چرخهها را به تأثیر عوامل کلانمقیاس جوی- اقیانوسی نسبت دادهاند. به عنوان مثال همانطور که در متن نیز بدان اشاره شده است، چرخههای 3-2 ساله را عمدتاً به تغییرات دو سالانه (QBO) الگوی بزرگمقیاس گردش عمومی جو و جریانات مداری و چرخههای 5-3 ساله را به پدیدهی انسو نسبت دادهاند. تغییرات بارش منطقهی هراز واقع در جنوب استان مازندران همچون تغییرات بارش تمامی پهنهی سیاره از الگوی زمانی- مکانی متغیر و متنوعی پیروی میکند. روند غالب در پهنهی پژوهش حاضر را روند افزایشی تشکیل میدهد؛ که در توجیه علت این روند افزایشی میتوان به نقش ساز و کارهای مرطوبتر شدن گرمتر و مرطوبتر شدن مرطوب اشاره کرد. تغییرات غیرخطی نوسانات بارش به عنوان تابعی از جریان جو میزان بارش را تحت تأثیر قرار میدهد. بررسی چرخههای بارش مبین این امر است که رشته کوههای البرز به عنوان مهمترین ساز و کار بارش و تغییرات آن در منطقه عمل میکند. تقدیر و تشکر مقالهی حاضر مستخرج از پایاننامهی مقطع کارشناسی ارشد رشته آب و هواشناسی کاربردی دانشگاه حکیم سبزواری با عنوان «ارزیابی سیلخیزی با استفاده از مدل ریاضی HEC-HMS، تحلیلهای آماری و GIS در حوضهی آبخیز هراز» میباشد. جا دارد از زحمات مسئولان زحمتکش دانشکده جغرافیا و علوم محیطی و کلیه اساتید محترم بخصوص اساتید راهنمایم آقایان غلام عباس فلاح قالهری و مختار کرمی و استاد مشاورم جناب آقای دکتر محمد باعقیده کمال تقدیر و تشکر را به عمل آورم.
[1]- Mosmann et al., [2]- Shifteh Some’e et al., [3]- Duhan and Pandey [4]- Kahya Kalayci [5]- Zhang et al., [6]- Liang et al., [7]- Tang et al., [8]- Some'e et al., [9]- Duhan and Pandey [10]- Tomozeiu [11]- Emprical Orthogonal Functions(EOF) [12]- Emilia Romagna [13]- Gemmer et al., [14]- Marengo et al., [15]- Tosic´ and M. Unkasevic [16]- Aphrodite [17]- Yatagani and et al., [18]- Sen’s Slope Estimator Nonparametric Method [19]-Vivekananda [20]- Sign Function [21]- Richard [22]- Mann-Kendall [23]- Index of Dispersion (ID) [24]- Index of Cluster Size (ICS) [25]- David and Moore [26]- Green’s Index (GI) [27]- Green [28]- Index of Cluster Frequency (ICF) [29]- Douglas [30]- Index of Mean Crowding (IMC) [31]- Lloyd [32]- Index of Patchiness (IP) [33]- Morisita’s Index (IM) [34]- Morisita [35]- fourier series [36]- Periodic [37]- Amplitude [38]- Frequency [39]- Maslen and Rockmore [40]- Periodogram [41]- Mitchell et al., [42]- Intergovernmental Panel on Climate Change [43]- Westra et al., [44]- warmer-gets-wetter [45]- Tan et al., [46]- Westra etl al., [47]- wet-gets-wetter [48]- EL Nino Southern Oscillation(ENSO) [49]- Quasi Binomial Oscillation(QBO) [50]- Kane and Teixeira [51]- Hartman et al., [52]- Lana et al., [53]- Fabra [54]- North Atlantic Oscillations(NAO) [55]- Torrence and Webster [56]- Garcia et al., [57]- Selvam and Joshi [58]- Schubert [59]- Garcia et al., [60]- Hartman et al., [61]- Cook and et al., مقدمه تغییر اقلیم از طریق ردیابی مشخصاتی نظیر کمیت و الگوی زمانی– کانی عناصر اقلیمی قابل مطالعه است. در این میان مطالعهی بارش از رواج زیادی برخوردار است. بارش به عنوان یکی از عناصر بنیادی اقلیم به ویژه طی قرن بیستم در عرضهای میانه و بالای نیمکره شمالی در حدود 1/0 درصد افزایش یافته است. از طرف دیگر در خشکیهای جنب حاره (عرضهای 30-10 درجه شمالی) روند بارش در حدود 3/0 درصد کاهش یافته و در نواحی حاره، افزایش نسبی بارش در حدود 2/0 تا 3/0 درصد بوده است (موسمان و همکاران[1]، 2004: 44). به تعبیری دیگر میتوان اذعان داشت که یکی از رفتارهای بارندگی که طی دهههای گذشته نمایانتر شده است، تغییر و جابجایی رژیم آن است که بـه صورت تغییر در مقدار بارنـدگی، زمان و تداوم رویداد بارش نمایان میشود (شیفته سومی و همکاران[2]، 2012: 9؛ دوهان و پاندی[3]، 2013: 138). عناصر آب و هوایی با تغییرات کمی و کیفی خود اوضاع محیطی را تحت تأثیر قرار میدهند. این اثرات هر چند که ممکن است بطئی و تدریجی باشند، لیکن میتوانند پیامدهای زیادی به دنبال داشته باشند. با شناخت ویژگیهای آب و هوائی میتوان قابلیتهای سازگاری انسان در محیط را ارزیابی کرد و متناسب با شرایط اقلیمی اشکال مطلوب ساختهای فضائی را با ابزار متناسب به وجود آورد (رزمی، 1389: 14). در یک نگاه منطقی، اقلیم و تغییرات آن بخصوص تغییرپذیری بارش میتواند سبب آسیبپذیری بیشتری به حیات بشر و همچنین بخشهای آب و کشاورزی باشد (مریانجی، 1391: 45). وجود دوره یا سالهایی که طی آن بارشی کمتر یا بیشتر از میانگین ثبت میشود، یک خصیصهی شناخته شده در اقلیم اکثر مناطق کرهی زمین میباشد. از آنجا که شرایط آب و هوایی هر منطقهای تابع متغیرهای زیادی است و این متغیرها هر ساله از شرایط و مقدار ثابتی برخوردار نیست، وجود نوسان در بارش مناطق مختلف را میتوان امری بدیهی و اجتنابناپذیر به حساب آورد. همچنین آب و هوای کرهی زمین در طول قرن بیستم به ویژه در دو دههی اخیر تعادل خود را از دست داده و تمایل به افزایش دما نشان داده است که به طور بالقوه میتواند منجر به تغییر چرخهی آب و بارش شود و ناهنجاریهای هیدرولوژیک همچون خشکیها و سیل را سبب شود (رزمی، 1389: 27). با توجه به مطرح شدن بحث تغییر اقلیم و تغییرات زمانی و مکانی عنصر بارش تا کنون از سوی محققان پژوهشهای متعددی با استفاده از روشهای پارامتری و نا پارامتری انجام شـده است که از جمله مـیتوان بـه پژوهشهای (کاهیا و کالایکی[4]، 2004؛ زانگ و همکاران[5]، 2006؛ لی یانگ و همکاران[6]، 2010؛ تانگ و همکاران[7]، 2011؛ شیفته سومئی و همکاران[8]، 2012؛ دوهان و پاندی[9]، 2013) اشاره کرد. توموزیو و همکاران[10] (2000) با استفاده از روش من-کندال و توابع تجربی متعامد[11]، تغییرپذیری میانگین بارندگی تابستانه را در امیلیا روماگنا[12] در ایتالیا بررسی کردند. مطالعات آنها نشان داد که همهی ایستگاهها دارای روند کاهشی در بارندگی تابستانه میباشند. گمر و همکاران[13] (2004) روند تغییرات بارش ماهانهی چین را مورد مطالعه قرار دادند و نشان دادند که تعداد، توزیع و کاهش روندها از ماهی به ماه دیگر متفاوت است. مارنگو و همکاران[14] (2004) الگوهای بلندمدت و دههای بارش در حوضهی آمازون را تجزیه و تحلیل کردند و به این نتیجه رسیدند که کل بارندگی در این منطقه دارای روندی منفی میباشد، اما این روند از توزیع مکانی ناهمگون برخوردار بوده است، برای مثال بخش شمالی حوضه دارای روندی منفی و بخش جنوبی دارای روندی مثبت میباشد. توشیج و انکاشویچ[15] (2005) با استفاده از روش تجزیه و تحلیل طیفی، تغییرات سریهای زمانی بارش ماهانه، فصلی و سالانه ایستگاه بلگراد را مورد ارزیابی قراردادند. نتایج نشان داد که تغییرات بارش با دورهی بازگشت دو ساله در ایستگاه بلگراد وجود دارد. در ایران نیز رزمی (1389) در پایاننامهی کارشناسی ارشد خود به ارزیابی تغییر رژیم بارندگی آذربایجان ایران پرداخت، روش اصلی به کار گرفته شده در پژوهش وی، روش تحلیل همسازها بوده است. بر اساس تحلیل همسازها مشخص گردید که دو سامانهی همدید- محلی تا بیش از 90 درصد بارش آذربایجان را توجیه میکنند. کمترین سهم بارش حاصل از سامانههای همدید در شمال شرق رخ میدهد. مریانجی (1391) در رساله دکتری خود تغییرپذیری بارش ایران را مورد بررسی قرار داد. نتایج پژوهش وی نشان داد که روند سالانه بارش در 54 درصد از مساحت کشور دارای روند مثبت با میانگین 72/1 میلیمتر بوده است. در نهایت نتایج تحلیل خوشهای نشان داد که تغییرپذیری بارش کشور در 5 گروه عمده قابل طبقهبندی میباشد. فرسادنیا و همکاران (1391) به بررسی تحلیل روند بارندگی در استان مازندران با استفاده از روش من-کندال منطقهای پرداختند. نتایج پژوهش آنان نشان داد که بررسی وجود روند در سری زمانی میانگین بارش سالانه در کل استان مازندران دارای روند صعودی در سطح معنیدار 1% است. در پژوهش حاضر ضمن بررسی ویژگیهای عمومی بارش حوضهی آبریز هراز سعی شده است روند دادههای بارش و دورهی بازگشت بارش در مقیاسهای ماهیانه، سالانه و فصلی نیز مورد ارزیابی قرار گیرد تا بتوان بر اساس دورههای بازگشت و سامانههای مؤثر در بارش منطقه، پیشآگاهی مناسبی به جهت مدیریت سیلاب ارائه گردد. مواد و روشها منطقهی مورد مطالعه در پژوهش حاضر، حوضهی آبریز هراز در جنوب استان مازندران میباشد. این منطقه از نظر تقسیمات سیاسی جزو شهرستان آمل میباشد. وسعت حوضهی آبریز هراز، 200051 هکتار است. محیط حوضهی هراز بـرابـر با 5/296062 کیلومتر، حداقل ارتفاع حوضه 200 متر و حداکثر ارتفاع آن 5600 متر میباشد. تـمامی رودخانههای حوضه به رودخانهی هراز منتهی مـیگردد. از مهمترین رودخانههای حوضه میتوان بـه رودخانههای آخنسر، شیرکلارود و نمارستاق اشاره کرد. از ارتفاعات مهم حوضهی کوههای دماوند که مرتفعترین منطقه در حوضه به شمار میآید و نیز شیمکوه و کوه امامزاده قاسم اشاره نمود. مراکز مهم مسکونی حوضهی آبریز هراز عبارتاند از پلور، نشل، تیران، رینه، کندوا، آب اسک، گزنک، بایجان، بلقلم و... اشاره کرد. شکل (1) موقعیت مطالعاتی هراز را نسبت به استان مازندران و کل کشور نشان میدهد.
شکل (1) موقعیت منطقهی مورد مطالعه جنوب استان مازندران سرزمینی ناهموار است و همین عامل به همراه چندین عامل دیگر موجب تغییرپذیری بارش در استان مازندران و بخصوص مناطق جنوبی آن میشود. با توجه به ویژگیهای اقلیمی خاص استان مازندران و اهمیت بارش در مدیریت منابع آب نیاز به مطالعات در مورد بارش و رفتارهای آن به شدت احساس میشود. اما در سطح استان مازندران دسترسی به دادههای بلندمدت بارش با گسترهی وسیع مسیر نیست. هم اکنون 15 ایستگاه همدید در استان مازندران وجود دارد که کمتر از یک سوم آنها آمار بلندمدت (30 ساله و یا بیشتر) متغیرهای هواشناسی را در اختیار دارند که با توجه به گسترهی وسیع استان، این تعداد ایستگاه برای بررسی تغییرات بارش که خود سرکشترین عنصر اقلیمی به شمار میرود و همچنین دارای نوسان مکانی - زمانی زیاد است؛ کافی نیست. افزون بر مشکلاتی که جهت دسترسی به دادههای بلند مدت اشاره شد؛ بخش اعظم منطقهی مورد مطالعه را مناطق مرتفع تشکیل میدهد (قرارگیری بخشی از مناطق شرقی قلهی دماوند و ارتفاعات شمالی و جنوبی آن) و با توجه به اینکه در نواحی با ارتفاع بیش از 2600 متر (کانونهای آبگیر کشور)، فاقد ایستگاههای سنجش و اندازهگیری بارش هستیم. همچنین نواحی مجاور که از تباین اقلیمی بالایی برخوردارند (دامنههای شمالی- جنوبی البرز) از توزیع مناسب ایستگاه برخوردار نیستند (عساکره، 1387: 25). در این پژوهش با استفاده از بروندادهای یاختهای بارش پایگاه آفرودیت[16]، به ارزیابی تغییرات زمانی- مکانی بارش ماهیانهی حوضهی آبریز هراز پرداخته شده است. دادههای این پایگاه با بازهی مکانـی خاورمیانه و جـنوب شرق آسیا و بازهی زمانی 1/1/1951 تا 31/12/2007 مـیلادی و تفکیک مکانی 25/0×25/0 و 5/0×5/0 درجـه در تـارنمای ایـنترنتی (http://www.chikyu.ac.jp/precip/index.html) موجود میباشند (یاتاگانی و همکاران[17]، 2012: 1407). در مطالعهی حاضر دادههای منطقهی خاورمیانه (APHRO_ME) از آخرین محصول پایگاه داده APHRODITE تحت عنوان v1101، با قدرت تفکیک 25/0×25/0 و فرمت ".nc" NetCDF از سایت مربوطه اخذ گردید. لازم به ذکر است کنترل کیفی بر روی دادههای فوق انجام شده است. در مرحلهی بعد با استفاده از تکنیکهای زمین آمار و روش کریجینگ، یاختهای به ابعاد 15×15 کیلومتر مربع بر منطقهی مورد مطالعه گسترانیده شد. حال به جای انجام تحلیل بر روی گرهگاهای 25/0×25/0 درجه و انتساب نتایج به دست آمده به پهنهی مورد مطالعه؛ از دادههای ریزمقیاس نمایی شده با تفکیک فضایی 15×15 کیلومتر استفاده شده است.
ارزیابی روند دادهها جهت برآورد شیب واقعی یک روند در یک سری زمانی، استفاده از روش ناپارامتریک سنس[18] میتواند یکی از روشهای مناسب در این زمینه باشد. این روش ابتدا توسط تیل در سال (1950) ارائه و سپس توسط سن در سال (1968) بسط و گسترش داده شد. روش حاضر همانند بسیاری دیگر از روشهای ناپارامتریک همچون من-کندال بر تحلیل تفاوت بین مشاهدادت سری زمانی استوار است. این روش زمانی میتواند مورد استفاده قرار گیرد که روند موجود در سری زمانی یک روند خطی باشد. این بدین معناست که در رابطهی (1) برابر است با:
در رابطهی (1) ، شیبخط روند و ، مقدار ثابت است. جهت محاسبه شیب خط روند یعنی، ابتدا بایستی شیب بین هر جفت داده مشاهداتی، با استفاده از رابطهی (2) محاسبه گردد:
در رابطهی (2) است. در این رابطه و به ترتیب دادههای مشاهداتی در زمانهای و است. با اعمال این رابطه برای هر جفت داده مشاهدهای، یک شیب به دست میآید. با قرار دادن این شیبها در کنار یکدیگر یک سری زمانی از شیبهای محاسبه شده حاصل میآید. یعنی اگر تعداد در سری زمانی وجود داشته باشد (در این تحقیق) ما به اندازهی برآورد شیب خواهیم داشت. در مرحلهی بعد سری زمانی مورد مطالعه به دست میآید. برای این کار تعداد از ها از کوچک به بزرگ مرتب میشوند و سپس با استفاده از یکی از رابطههای زیر اقدام به تعیین میانهی سری زمانی میگردد. اگر تعداد مشاهدات سری زمانی مورد مطالعه فرد باشد، از رابطهی (3) و اگر زوج باشد از رابطهی (4) استفاده میگردد (ویوکاناندان[19]، 2007: 104):
نتیجهی حاصل از این رابطه، به دست آمدن شیبخط روند است. اگر شیب خط روند مثبت باشد، حاکی از صعودی بودن روند و اگر منفی باشد، دال بر نزولی بودن روند است. در گام بعدی، آزمون نمودن شیب به دست آمده در فاصله اطمینان 95 درصد است. جهت انجام این آزمون از رابطهی زیر استفاده میگردد:
در رابطهی (5) عبارت است از آماره توزیع نرمال استاندارد در یک آزمون دوطرفه که برای سطح اطمینان 95 درصد برابر با است و نیز واریانس پارامتر است. جهت به دست آوردن مقدار پارامتر و همچنین مراحل زیر بایستی پیموده شود: الف) محاسبه اختلاف بین تک تک جملات سری با همدیگر و اعمال تابع و استخراج پارامتر :
در رابطهی (6) تعداد مشاهدات سری (در این تحقیق برابر با 57 سال)، و نیز به ترتیب دادههایام و ام سری میباشند. ب) محاسبهی تابع علامت[20] (ریچارد[21]، 1987: 25) که به صورت تابع 7 قابل محاسبه است:
ج) محاسبه واریانس توسط یکی از روابط زیر. اگر تعداد دادههای سری زمانی بزرگتر از 10 مورد باشد از رابطهی 8 و اگر کوچکتر از 10 باشد، از رابطهی 9 بهره گرفته میشود:
در رابـطهی 8 و 9، تعداد دادههای مشاهداتی، تعداد سریهائی است که در آنها حداقل یک دادهی تکراری وجود دارد و نیز بیانگر فراوانی دادههای با ارزش یکسان میباشد و در نهایت حدود اعتماد بالا و پایین به کمک روابط زیر محاسبه میگردد:
در رابطهی (10) تعداد شیبهائی است که از طریق رابطهی (2) به دست آمده است. حال امین و امین شیبها را از بین شیبهای محاسبه شده استخراج مینمائیم. در صورتی که عدد صفر در دامنهی بین دو شیب استخراج شده فوق قرار گیرد، فرض صفر پذیرفته شده و عدم وجود روند در سری دادهها تأیید میگردد. در غیر این صورت فرض صفر رد شده و وجود روند در سطح اطمینان مورد آزمون پذیرفته میگردد. در نهایت جهت به دست آوردن مقداردر رابطهی 1، تعداد تفاضل محاسبه میشود. سپس میانه این مقادیر، برآورد را در اختیار میگذارد. در پژوهش حاضر برای ارزیابی روند دادههای بارش حوضهی آبریز هراز، از آزمون ناپارامتریک من-کندال[22] نیز استفاده شده است. این آزمون ابتدا توسط من (1945) ارائه و سپس توسط کندال (1975) تکامل یافت. این روش نیز همانند آزمون شیب خط سنس در مباحث علوم محیطی به طور گستردهای کاربرد دارد. در آزمون من-کندال نمره از رابطهی 11، محاسبه میگردد:
در رابطه (11) نشانهی تفاوت مقادیر با یکدیگر و واریانس است:
که تعداد مشاهدات سری (در این تحقیق برابر با 57 سال)، و نیز به ترتیب دادههایام و ام سری میباشند. تابع علامت است که شرح آن در رابطهی (7) آمده است. عاملی مربوط به تصحیح پراش است. در صورتیکه دادههای تکراری در اطلاعات وجود داشته باشد از رابطهی (14) استفاده میشود:
در رابطهی (14) تعداد دادههای مشاهده شده و معرف تعداد سریهایی است که در آنها حداقل یک دادهی تکراری وجود دارد. آمارههای پراکندگی در پژوهش حاضر برای بررسی ساختار مکانی یاختههای بارش حوضهی آبریز هراز، از آمارههای پراکندگی استفاده شده است. آمارههای پراکندگی استفاده شده عبارتاند از: شاخص پراکندگی[23]، شاخص تجمع، شاخص گرین، شاخص فرکانس خوشه، شاخص تراکم میانگین، شاخص لکهها و شاخص موریسیتا. شاخص پراکندگی نسبت واریانس به میانگین را بررسی میکند. اگر دادههای مورد مطالعه از توزیع تصادفی پیروی کنند انتظار میرود ID برابر یک باشد. این آماره از توزیعبا درجه آزادی پیروی میکند (رابطهی 15 را ملاحظه نمایید). شاخص تجمع[24] (دیوید و مور[25]، 1954: 50) تابع مستقیم از شاخص پراکندگی است. اگر دادههای مورد مطالعه از توزیع تصادفی پیروی کنند، انتظار میرود این شاخص برابر صفر باشد. مقادیر ارزشی مثبت در این آماره حاکی از خوشهای بودن دادهها و مقادیر ارزشی منفی بیانگر ناخوشه بودن دادهها میباشد (رابطهی 16). آماره گرین[26] (گرین[27]، 1966) در حقیقت اصلاح شده شاخص تجمع میباشد که مستقل از N است. در این آماره ارزش عددی صفر برای توزیعهای تصادفی و مقادیر ارزشی یک برای توزیعهای با حداکثر خوشه است (رابطهی 17):
شاخص فرکانس خوشه[28] (داگلاس[29]، 1975: 403) معیاری برای اندازهگیری خوشه بر اساس تابع K با استفاده از توزیع دوجملهای منفی است (رابطهی 18). شاخص تراکم میانگین[30] (لوید[31]، 1967: 9) متوسط تعداد نقاط موجود در منطقهی مورد مطالعه است که به طور تصادفی از یک نقطه سنجش میشود (رابطهی 19):
شاخص لکهها[32] (لوید، 1967: 14) شدت الگوی تأثیرپذیر را مورد مطالعه قرار میدهد (رابطهی20). شاخص موریسیتا[33] (موریسیتا[34]، 1959: 13) احتمال مقیاسبندی خوشههای دو یاخته از کل یاختههای مورد مطالعه را مورد ارزیابی قرار میدهد (رابطهی 21):
علاوه بر موارد فوق، ضریب چولگی (g1) و ضریب درجهی اوج (g2) نیز برای هر یک از ماههای مورد مطالعه به طور جداگانه با استفاده از رابطهی 22 محاسبه شده است:
تکنیک تحلیل طیفی (تحلیل سریهای فوریه[35]) برای بـررسی چرخههای بارش حوضهی آبـریز هراز از روش تحلیل طیفی استفاده شد. در تـکنیک تحلیل طیفی ابتدا سریهای زمانی به تابع فرکانسی (به صورت تابعی دورهای[36] با دامنه[37] و فراوانی[38]) تبدیل شد. روشهای پرشمار و متنوعی برای تبدیل سریهای بسامد وجود دارد (عساکره و رزمی، 1391: 152). یکی از روشهای ساده و کارآ، روش همسازهاست. در این روش سری زمانی به طول ، به فرم یک مدل فوریه و به شکل رابطهی (23) نوشته میشود (ماسلن و روکمور[39]، 1997: 200).
که در آن و ضرایب فوریه خوانده میشوند و از طریق روش کمترین مربعات خطا قابل برآورد هستند. و فراوانی تکرار مشاهدات (عکس دوره بازگشت) است و با نشان داده میشود. واریانس برای فرکانس و طول دوره آماری از رابطهی (24) به دست میآید:
نموداری که مقادیر را در مقابل نشان میدهد، دوره نگار[40] نامیده میشد. از آنجا که دوره نگار نسبتی از توزیع واریانس در فرکانسهای بنیادی و برآوردی پر افت و خیز از طیف است، برای کاهش واریانس، اریبی و ایستاسازی این برآورد، باید آن را در مجاورت فرکانس مورد نظر و به وسیلهی صافیهای مناسب هموار ساخت. به عبارت دیگر دوره نگار هموار شده، برآوردی از طیف جامعه است که حاصل میانگین موزن مقدار راست و چپ فرکانس مورد نظر میباشد (گلدسته، 1377: 245-246). یکی از روشهای ساده و مناسب برای هموارسازی دادههای سالانه به وسیلهی میتچل و همکاران[41] (1966: 34) معرفی شد. در این رویه، دوره نگار با یک میانگین متحرک موزون سهجملهای هموار میشود.
بحث و نتایج مشخصات آماری بارش حوضهی آبریز هراز جدول (1) بعضی از مشخصات آماری بارش را طی ماههای مختلف نشان میدهد. مطابق این جدول بیشترین میانـگین بارش با 99/54 مـیلیمتر در ماه مارس رخ داده است. کـمترین مـقدار بـارش بـه میزان 46/16 میلیمتری در ژوئیه اتفاق افتاده است. این در حالی است که بیشترین ضریب تغییرات بارش نیز با 66/102 در ماه ژوئیه اتفاق افتاده است. چولگی و توزیع ضریب چولگی (g1) بارش در همهی ماهها مثبت بوده است. سامانههای بارشزا با توجه به شرایط دینامیک و ترمودینامیک و بسته به موقعیت جغرافیایی خود، در برخورد با شرایط متفاوت محلی میتوانند نواحی بارشی متفاوتی ایجاد کنند (غیور و همکاران، 1390: 19). لذا مقدار بارش دارای شاخصهای آماری متفاوتی خواهد بود. همانطور که در جدول (1) نیز مشخص میباشد، بیشترین مقدار عددی شاخص g2 (ضریب درجه اوج) مربوط به ماه ژوئیه (26/10) و کمترین مقدار آن متعلق به ماه سپتامبر (32/5) بوده است. ضریب درجهی اوج نشاندهندهی تفاوت مکانی بارش میباشد. دلیل بالاتر بودن شاخص درجه اوج در ماه ژوئیه ناشی از تفاوت مکانی بالای بارش میباشد. علت کاهش شدید شاخص با شروع فصل پاییز را میتوان با فرارفت رطوبت بیشتر همزمان با تقویت و نفوذ زبانههای پرفشار سیبری نسبت داد که با مطالعات علیجانی (1389) نیز همخوانی دارد. همرفت محلی در منطقهی مورد مطالعه در فصول انتقالی که رطوبت حاصل از سامانههای مدیترانهای هنوز در محیط وجود دارد، نیز میتواند در تغییرات مکانی بارش منطقه موثر باشد. اطلاعات کامل مشخصههای آماری هر یک از ماههای سال در دورهی آماری 57 ساله در جدول (1) آمده است. آمارههای پراکندگی بارش در حوضهی آبریز هراز جدول (2) مشخصات آمارههای پراکندگی بارش را طی ماههای مختلف نشان میدهد. بر اساس شاخص ID ماههای ژوئیه و اکتبر بالاترین مقادیر عددی را نشان دادهاند. این شاخص بهخوبی ناهنجاری مکانی بارش حوضهی مورد مطالعه را نشان داده است، به طوری که در ماههای نامبرده ایستگاه کره سنگ به نمایندگی از مناطق شمالی حوضه در مقایسه با ایستگاه پلور به نمایندگی از مناطق جنوبی حوضه در یک فاصله کمتر از 60 کیلومتری، اختلاف 10 میلیمتری دارند. همانگونه که در جدول (2) نشان داده شده است، شاخص ICS برای هر 12 ماه مورد مطالعه مثبت و بالای 4 بوده است که نشاندهندهی الگوی خوشهای بالای بارش در حوضه میباشد. شاخصهای GI تبیینکنندهی همین امر میباشد. شاخص ICF نیز اندازهی خوشه را نشان میدهد. بر اساس این شاخص در فصل زمستان بزرگترین خوشههای بارش در حوضهی آبریز هراز ایجاد میشوند که نشاندهندهی نظم نسبی بارش در منطقهی مورد مطالعه میباشد. شاخصهای IP و GI نشاندهندهی تفاوت مکانی خوشههای بارش میباشند که فصل تابستان بالاترین مقادیر تفاوت مکانی خوشهای را به خود اختصاص داده است. شاخص MI نیز مقدار تعدیلشده شاخص IPمیباشد که مبین تفاوت ساختار مکانی بارش میباشد. جدول (1) مشخصات آماری بارش برای هر یک ماههای سال در حوضهی آبریز هراز
جدول (2) مشخصات آمارههای پراکندگی بارش حوضهی آبریز هراز
ارزیابی روند دادهها از آنجا که مبانی زندگی انسان با بارش هماهنگ شده است؛ وجود روند کاهشی یا افزایشی در بارش به معنای تغییر شکل زندگی خواهد بود. لذا جهت شناسایی روندهای کاهشی یا افزایشی در بارشهای حوضهی آبریز هراز از دو روش ناپارامتریک برآوردکنندهی شیب سنس و من-کندال استفاده گردید. نتایج حاصل از این روشها نشان داد هر چند نشانههایی از روند کاهشی (سه ماه فوریه، مارس و آوریل در آزمون من-کندال و شیب سنس) و افزایشی (9 ماه دیگر سال) در بارشهای حوضه مشاهده میشود، اما از بین این روندهای کاهشی و افزایشی تنها دو ماه سپتامبر و اکتبر از نظر آماری در سطح معنیدار 5 درصد معنیدار بودهاند. اطلاعات هر یک از روندها در سطوح 5 و 1 درصد در جدول (3) آمده است.
جدول (3) نتایج حاصل از آزمون سنس و من-کندال در سطوح احتمالاتی 5 و 1 درصد برای بارش ماهانه حوضهی آبریز هراز
نتایج به دست آمده نشان میدهد یاختههایی که روند کاهشی در آنها مشاهده شده است، همگی در بخشی جنوبی حوضه و اطراف قله دماوند واقع شدهاند که میتواند نشانههایی از تغییر شرایط بارش در این قسمت از ایران باشد. هر چند اثبات وجود روند معنیدار در یک سری زمانی بارندگی نمیتواند به تنهایی دلیلی قاطع بر وقوع تغییر اقلیم در یک منطقه باشد، اما فرض رخداد آن را میتواند تقویت کند. لذا همین تغییرات رخ داده میتواند حائز اهمیت باشد. به طور کلی بارش حوضهی آبریز هراز در مقیاس سالیانه بر اساس نتایج آزمون شیب سنس افزایش 67/7 میلیمتری در هر دهه داشته است. بر این اساس در طول دورهی 57 ساله به رقم باور عمومی مبنی بر کاهش بارش در سطح کشور و منطقه، بارش سالانه حوضهی آبریز هراز افزایش 719/43 میلیمتری داشته است. بر این اساس میتوان برای روند تغییرات بارش در سطح حوضهی آبریز هراز رابطهای به شکل ارایه داد. شکل (1) برازش نمودار خط سنس بر سری زمانی دادههای میانگین یاختهای بارش حوضهی آبریز هراز را در دورهی زمانی مورد مطالعه نشان میدهد.
شکل (1) برازش نمودار شیب سنس بر سری زمانی دادههای میانگین یاختهای بارش حوضهی آبریز هراز
در نگاه اول به نظر میرسد نتایج حاضر دور از واقعیت است، چرا که با افزایش دمای متوسط جهانی باور عمومی مبنی بر کاهش مقدار نزولات جوی است، اما نتایج پژوهش حاضر نشان داده است که مقدار بارش در 57 سال اخیر در حوضهی آبریز هراز افزایش یافته است. بررسیها نشان داده است که این افزایش بارش بیشتر به شکل رگباری بوده و بارشهای آخر فصل سرد و فصول گرم سال در حال افزایش است. چنین افزایشی نه تنها سودمند نیست، بلکه موجب وقوع رخداد وقایع فرین اقلیمی همچون سیل شده است. روند مشاهده شده در حوضهی آبریز هراز به خوبی با گزارشهای سازمانهای مختلف حاکی از تغییرات جدی در اقلیم همخوانی دارد؛ به طوری که بر اساس گزارشهای ارایه شده از تغییرات آب و هوایی (IPCC[42],2013: 18)، از نیمهی دوم قرن بیستم روند بارش سنگین (بارشهای رگباری) در نـقاط مختلف جهان رو به افزایش است. نتایج سایر مطالعات در سطح جهانی تأییدی بر نتایج حاضر است، زیرا بررسی دادههای بارندگی بیش از ۸ هزار ایستگاه هواشناسی در سطح جهان و در فاصلهی زمانی بین سالهای ۱۹۰۰ تا ۲۰۰۹ نشان داده است که در این دورهی زمانی گرمای جو زمین به اندازهی یک درجهی سانتیگراد زیاد شده و در نتیجه آن تعداد بارندگیهای شدید نیز افزایش یافته است (وسترا و همکاران[43]، 2013: 3904). ساز و کار این امر را نیز میتوان به گونهای تشریح کرد که با افزایش گازهای گلخانهای و افزایش گرما، مقدار رطوبت بیشتری در جو میماند و بارشها شدیدتر میشود. به ازای هر یک درجه گرمایش جو زمین، شدت باران ۷ درصد افزایش مییابد و انتظار میرود دمای کرهی زمین در انتهای قرن جاری ۳ تا ۵ درجه سانتیگراد افزایش یابد و تغییرات آب و هوا شدت بارندگی انتهای قرن را نیز به میزان قابل توجهی افزایش دهد (همان، 3916). مهمترین ساز و کاری که میتوان بر اساس آن به تبع وجود روند افزایشی دما، مقادیر مثبت بارش را در حوضهی آبریز هراز توجیه نمود، مکانیسم «مرطوبترشدن گرمتر[44]» (تن و همکاران[45]، 2015: 451) است که بر اساس آن هر چه دمای سطحی منطقه بالاتر باشد، مقدار بارندگی بیشتر است. در دهههای اخیر روند چشمگیری در بارش اکثر ماههای سال مشاهده نشده است و تنها شاهد تغییر و تنوع سالانه بارندگیها هستیم و حتی در برخی از موارد با گرمتر بودن سالها بارندگی نیز بیشتر صورت گرفته است. در توجیه این تغییر و تنوع سالانه میتوان به نقش بارز پدیدهی النینو اشاره کرد. تلفیق پدیدهی گرمایش زمین و پدیدهی النینو و خروجی آن دو تحت مدلهای عددی در سطح جهان نشان داده است که، بارندگی در مناطق مرطوب، بیشتر و در مناطق خشک، کمتر میشود (وسترا و همکاران[46]، 2013: 3916). این تأییدی بر مکانیسم مرطوبتر شدن مرطوب[47] است که در حوضهی آبریز هراز با توجه به مرطوب بودن منطقه قابل تعمیم است. ادغام ساز و کار «مرطوبترشدن گرمتر» و «مرطوبتر شدن مرطوب» حدود ۸۰ درصد تغییرات بارندگی ناشی از گرمایش زمین را توضیح میدهد که در حوضهی آبریز هراز نیز مکانیسم مرطوبتر شدن گرمتر و مرطوبتر شدن مرطوب مشاهده گردید. تحلیل همسازهای بارش حوضهی آبریز هراز برای تحلیل دقیقتر تغییرات بارش حوضهی آبریز هراز چرخههای حاکم در هر ماه و فصل با استفاده از تحلیل طیفی مورد بررسی قرار گرفت. اساس روشهای تحلیل طیفی، تقسیمبندی تغییرپذیری سریهای زمانی به اجزاء یا بخشهای حاوی دورههای مختلف یا فرکانس است که به عنوان یکی از روشهای استخراج و تحلیل نوسانات آشکار و نهان طول موجهای مختلف، مطرح میباشد. جداول (4 تا 9) دوره نگار، طیف و فاصله اطمینان را برای سریهای زمانی بارش حوضهی آبریز هراز در دورههای مختلف ارائه دادهاند.
جدول (4) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 99 درصد در مقیاس ماهانه
جدول (5) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 99 درصد در مقیاس سالانه و فصلی
تغییرات چرخههای بارش نشان میدهد که در فصول تابستان و پاییز در سطح اطمینان 99 درصد چرخههای 2 ساله بر بارش حوضهی مورد مطالعه حاکم بوده است این در حالی است که در فصول بهار و زمستان بارش در حوضهی آبریز هراز هیچگونه چرخهای نداشته است. بررسی چرخههای بارش برای متوسط سالیانه نیز حاکی از چرخه کوتاهمدت 2 ساله بوده است. تحلیل همسازهای بارش در حوضهی آبریز هراز در سطوح اطمینان 95 و 90 درصد نشان داد که علاوه بر چرخههای 2 ساله، چرخههای بلندمدت 5، 8، 14، 25 و 28 ساله هم بر بارش این منطقه از کشور حاکم بوده است. در مجموع میتوان گفت که بارش هراز در شمارهی چرخههای 5 و 8 با احتمال 49/0 دارای دورهی بازگشت کوتاهمدت 2 ساله است. جدول (6) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 95 درصد در ماههای مختلف
جدول (7) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 95 درصد در مقیاس فصلی و سالانه
در بین نتایج ارایه شده آنچه بیش از همه خودنمایی مینماید، چرخههای 2-3 ساله هستند. بررسیهای یاختههای مورد نظر نشان میدهد که غالباً این چرخهها در ارتفاعات دماوند و دامنههای مرتفع آن قابل رؤیت هستند. بیشتر دانشمندان این چرخهها را به ال نینو- نوسانات جنوبی (ENSO[48])، تغییرات دوسالانه (QBO[49]) الگوی بزرگمقیاس گردش عمومی جو، جریانات مداری و سایر فرایندهای جوی- اقیانوسی نسبت دادهاند. برای مثال، کانه و تگزیرا[50] (1991: 89) چرخههای 2-3 ساله در بارش ماساچوست، هارتمن و همکاران[51] (2002: 161) چرخههای 2-3 ساله را در چین و لانا و همکاران[52] (2005: 193) چرخههای 6/4 و 1/2 ساله بارش در ایستگاه فابرا[53] در شمال شرق اسپانیا را به تغییرات دوسالانه (QBO) و چرخههای 2/9 و 5/5 ساله را به نوسانات اطلس شمالی (NAO[54]) و چرخههای 8/11 ساله را به فعالیت لکههای خورشیدی نسبت دادهاند.
جدول (8) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح90درصد در ماههای مختلف
جدول (9) مشخصات آماری چرخههای بارش طی دورهی مورد مطالعه (2007-1951 میلادی) در سطح 90 درصد در مقیاس فصلی و سالانه
تورنس و وبستر[55] (1999: 2679) چرخههای 2-7 ساله را در بارش موسمی هند، گارسیا[56] و همکاران (2002: 78) چرخههای 7/2 ساله را در بارش شبهجزیره ایبری و آزاد و همکاران (2010: 2290) چرخههای 3-5 ساله را در بارش موسمی هند ناشی از تأثیر پدیده انسو دانستهاند. سلوام و جلوشی[57] (1995: 613) چرخههای 7-2 ساله و چرخههای 3-2 ساله را در دمای سطح زمین کشف کردند. ایشان به ترتیب این چرخهها را به انسو و تغییرات دوسالانه نسبت دادهاند. همانطور که مشاهده شد، چرخههای کوتاهمدت 2 تا 4 ساله و 3-5 ساله بیشترین حاکمیت را بر بارش حوضهی آبریز هراز داشتهاند که این نوع چرخهها را میتوان به رخداد النینو نسبت داد. بیشتر پژوهشگران چرخههای 2 تا 4 ساله را به النینو- نوسانات جنوبی (ENSO) (هولرینگ و کومار، 2003؛ اسچیبورت و همکاران[58]، 2009) و تغییرات دوسالانه (QBO) الگوی بزرگمقیاس گردش عمومی جو و جریانات مداری و سایر فرایندهای جوی- اقیانوسی نسبت دادهاند (گارسیا[59] و همکاران، 2002: 77-95، هارتمن[60] و همکاران، 2008:156 ، کوک و همکاران[61]، 2011a,b). بررسی چرخههای بارش طی دورههای مورد مطالعه (2007-1951 میلادی) در سطوح اطمینان سهگانهی 90، 95 و 99 درصد نشاندهندهی چندین الگوی مختلف بارش میباشد که یکی از آنها، چرخههای 11-5 ساله میباشد. بر این اساس در سطح 99 درصد ماه می با دوره بازگشت 5 ساله؛ در سطح 95 درصد ماه ژانویه با دورهی بازگشت 9 ساله، اکتبر با دوره بازگشت 6 ساله و نوامبر با دورهی بازگشت 8 ساله و در سطح 90 درصد ماه آوریل دارای دورههای بازگشت 8 و 5 ساله، اکتبر دارای دورهی بازگشت 6 ساله و نوامبر دارای دورهی بازگشت 8 ساله میباشد. تغییرپذیری درون دههای (دهه به دهه) فقط برای ماهها و یا فصولی که چرخههای 11 سال به بالا دارند، قابلمشاهده است (مانند ماههای ژولای، اوت و سپتامبر در سطح اطمینان 99 درصد). مـطالعهای که توسط جهانـبخش و عدالتدوست (1387) در شـمال غرب ایران انـجام شده است نشان داده است که چرخههای نامبرده مبین تأثیر فعالیت لکه (کلف)های خورشیدی است که با نتایج پژوهش حاضر نیز همخوانی دارد. نتایج نشان میدهند نوسانات اطلس شمالی تابعی از عملکرد لکههای خورشیدی است. همچنین بررسیها مؤید این امر است که چرخههای غیرسینوسی (روند) نیز در سری زمانی بارش منطقهی مورد مطالعه با توجه به مشاهده چرخههای درون دههای 20 تا 40 ساله در پهنهی مورد مطالعه قابل استنباط است. ارزیابی چرخههای بارش حوضهی آبریز هراز نشان داد که مقدار بارش تحت تأثیر نوسانات کوتاهمدت است؛ همچنین مشخص شد که تغییرات غیرخطی نوسانات بارش به عنوان تابعی از جریان جو، میزان بارش را تحت تأثیر قرار میدهد. تغییرات بارش حوضهی آبریز هراز همچون تغییرات بارش تمامی پهنهی سیاره از الگوی زمانی- مکانی متنوعی پیروی میکند. بررسیهای چرخهها نشان داده است پدیدههای النینو- نوسانات جنوبی و تغییرات دوسالانه بیشترین تأثیر را بر بارشهای جنوب حوضهی آبریز هراز و به طور کلی جنوب استان مازندران داشتهاند که با نتایج مریانجی (1391) همخوانی دارد. همچنین مشخصشده است ترابری رطوبت و کوههای البرز به عنوان مهمترین ساز و کار بارش و تغییرات آن در منطقه است. نتیجهگیری به طور کلی، تحلیلهای آماری و مکانی- زمانی رویدادهای بارش ابزار مناسبی هستند که میتوانند شناخت خوبی از شرایط بارش منطقه ارائه کنند. برخی از یافتههای اصلی پژوهش حاضر بدین شرح است: علیرغم وجود متغیرهای بارش متفاوت در نواحی مختلف، توزیع چگالی احتمال فراوانی رویدادهای بارش در همه ماههای سال در منطقهی مورد مطالعه، دارای چولگی مثبت و کشیدگی بیشتر از نرمال است. به عبارت دیگر، تراکم نسبی رویدادها پیرامون مرکز توزیع، بیشتر است. نتایج نشان داد که وضعیت ضریب تغییرات بارش در سطح حوضه در ماههای گرم سال نسبت به دورهی سرد سال بیشتر بوده است، به طوری که ژوئیه در فصل تابستان بیشترین ضریب تغییرات بارش را با 66/102 درصد داشته است. تحلیل ماهانهی فراوانی رویدادهای بارش در این منطقه نشان میدهد که شاخص پراکندگی بارش از ماه ژانویه تا آوریل یک روند افزایشی و در ماههای دیگر (به استثنای ژوئیه و اکتبر) روندی کاهشی دارد که همین شرایط برای واریانس بارش نیز مشاهده شده است که به خوبی ناهنجاری بالای بارش در سطح حوضه را در طول سال نمایان ساخته است. نتایج شاخص ICF که اندازه خوشههای بارش را در سطح حوضه نشان میدهد که در فصل زمستان بزرگترین خوشههای بارش در حوضهی آبریز هراز ایجاد میشوند که نشاندهندهی نظم نسبی بارش در منطقهی مورد مطالعه میباشد. نتایج حاصل از روشهای ارزیابی روند با دو آزمون شیب سنس و من-کندال نشان داد که هر چند نشانههایی از روند کاهشی (سه ماه فوریه، مارس و آوریل در آزمون من-کندال و شیب سنس) و افزایشی (9 ماه دیگر سال) در بارشهای حوضه مشاهده میشود، اما از بین این روندهای کاهشی و افزایشی، تنها دو ماه سپتامبر و اکتبر از نظر آماری در سطح معنیدار 5 درصد معنادار بودهاند. به طور کلی بارش حوضهی آبریز هراز در مقیاس سالانه بر اساس نتایج آزمون شیب سنس افزایش 67/7 میلیمتری در دهه داشته است. بر این اساس در طول دورهی 57 ساله، بارش سالانه حوضهی آبریز هراز افزایش 719/43 میلیمتری داشته است. از آنجا که سه ماه فوریه، مارس و آوریل از پربارش ترین ماههای سال هستند و روند کاهشی داشتهاند، به معنی کاهش بارش فصل سرد و افزایش بارش فصل گرم در منطقهی مورد مطالعه است و همین خود دلیل قوی بر تغییر نامطلوب بارش در حوضهی مورد مطالعه است که میتواند خطرات جدی برای فرسایش، سیلخیزی و مدیریت منابع آب در پی داشته باشد. تحلیل طیفی بارش در حوضهی آبریز هراز نشان داد که بارش در این پهنه دارای چرخههای نوسانی معنیدار 3-2، 5-3، 11-5 و 11 سال به بالاست. چرخههای 11-5 ساله بیشترین گستره را در برداشته است. ورود این چرخهها را به تأثیر عوامل کلانمقیاس جوی- اقیانوسی نسبت دادهاند. به عنوان مثال همانطور که در متن نیز بدان اشاره شده است، چرخههای 3-2 ساله را عمدتاً به تغییرات دو سالانه (QBO) الگوی بزرگمقیاس گردش عمومی جو و جریانات مداری و چرخههای 5-3 ساله را به پدیدهی انسو نسبت دادهاند. تغییرات بارش منطقهی هراز واقع در جنوب استان مازندران همچون تغییرات بارش تمامی پهنهی سیاره از الگوی زمانی- مکانی متغیر و متنوعی پیروی میکند. روند غالب در پهنهی پژوهش حاضر را روند افزایشی تشکیل میدهد؛ که در توجیه علت این روند افزایشی میتوان به نقش ساز و کارهای مرطوبتر شدن گرمتر و مرطوبتر شدن مرطوب اشاره کرد. تغییرات غیرخطی نوسانات بارش به عنوان تابعی از جریان جو میزان بارش را تحت تأثیر قرار میدهد. بررسی چرخههای بارش مبین این امر است که رشته کوههای البرز به عنوان مهمترین ساز و کار بارش و تغییرات آن در منطقه عمل میکند. تقدیر و تشکر مقالهی حاضر مستخرج از پایاننامهی مقطع کارشناسی ارشد رشته آب و هواشناسی کاربردی دانشگاه حکیم سبزواری با عنوان «ارزیابی سیلخیزی با استفاده از مدل ریاضی HEC-HMS، تحلیلهای آماری و GIS در حوضهی آبخیز هراز» میباشد. جا دارد از زحمات مسئولان زحمتکش دانشکده جغرافیا و علوم محیطی و کلیه اساتید محترم بخصوص اساتید راهنمایم آقایان غلام عباس فلاح قالهری و مختار کرمی و استاد مشاورم جناب آقای دکتر محمد باعقیده کمال تقدیر و تشکر را به عمل آورم.
[1]- Mosmann et al., [2]- Shifteh Some’e et al., [3]- Duhan and Pandey [4]- Kahya Kalayci [5]- Zhang et al., [6]- Liang et al., [7]- Tang et al., [8]- Some'e et al., [9]- Duhan and Pandey [10]- Tomozeiu [11]- Emprical Orthogonal Functions(EOF) [12]- Emilia Romagna [13]- Gemmer et al., [14]- Marengo et al., [15]- Tosic´ and M. Unkasevic [16]- Aphrodite [17]- Yatagani and et al., [18]- Sen’s Slope Estimator Nonparametric Method [19]-Vivekananda [20]- Sign Function [21]- Richard [22]- Mann-Kendall [23]- Index of Dispersion (ID) [24]- Index of Cluster Size (ICS) [25]- David and Moore [26]- Green’s Index (GI) [27]- Green [28]- Index of Cluster Frequency (ICF) [29]- Douglas [30]- Index of Mean Crowding (IMC) [31]- Lloyd [32]- Index of Patchiness (IP) [33]- Morisita’s Index (IM) [34]- Morisita [35]- fourier series [36]- Periodic [37]- Amplitude [38]- Frequency [39]- Maslen and Rockmore [40]- Periodogram [41]- Mitchell et al., [42]- Intergovernmental Panel on Climate Change [43]- Westra et al., [44]- warmer-gets-wetter [45]- Tan et al., [46]- Westra etl al., [47]- wet-gets-wetter [48]- EL Nino Southern Oscillation(ENSO) [49]- Quasi Binomial Oscillation(QBO) [50]- Kane and Teixeira [51]- Hartman et al., [52]- Lana et al., [53]- Fabra [54]- North Atlantic Oscillations(NAO) [55]- Torrence and Webster [56]- Garcia et al., [57]- Selvam and Joshi [58]- Schubert [59]- Garcia et al., [60]- Hartman et al., [61]- Cook and et al., | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مراجع | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
- جهانبخش، سعید و معصومه عدالت دوست (1387)، «مطالعهی تغییرات اقلیمی در ایران (مطالعه موردی: نوسانات اطلس شمالی به عنوان شاخصی از تأثیرگذاری فعالیتهای خورشیدی در تغییرات بارش منطقهی آذربایجان)»، سومین کنفرانس مدیریت منابع آب ایران، دانشگاه تبریز، دانشکدهی مهندسی عمران. - عساکره، حسین (1387)، «کاربرد روش کریجینگ در میانیابی بارش»، جغرافیا و توسعه، شمارهی 12، صص 25-42. - عساکره، حسین و رباب رزمی (1391)، «تحلیل تغییرات بارش سالانهی شمال غرب ایران»، جغرافیا و برنامهریزی محیطی، سال 23، شمارهی 3، صص 147-162. - علیجانی، بهلول (1389)، «آب و هوای ایران»، انتشارات دانشگاه پیام نور، چاپ دهم، ص 221. - غیور، حسنعلی؛ مسعودیان، سیدابوالفضل؛ آزادی، مجید و حمید نوری (1390)، «تحلیل زمانی و مکانی رویدادهای بارشی سواحل جنوبی خزر»، فصلنامهی تحقیقات جغرافیایی، سال 25، شمارهی 100، صص 1-30. - فرساد نیاد، فرهاد؛ رستمی کامرود، محسن و علیرضا مقدم نیا (1391)، «تحلیل روند بارندگی در استان مازندران با استفاده از روش من-کندال منطقهای»، تحقیقات منابع آب ایران، سال هشتم، شمارهی 2، صص 60-70. - گلدسته، اکبر (1377)، «راهنمای کاربران SPSS 6.0 For Windows»، جلد اول، مرکز فرهنگی انتشاراتی حامی، تهران. - مریانجی، زهره (1391)، «تغییرپذیری رژیم بارش در ایران»، پایاننامهی دکتری رشته جغرافیای طبیعی گرایش اقلیمشناسی، دانشکده علوم جغرافیایی و برنامهریزی، دانشگاه اصفهان. -Azad, S., Vignesh, T.S., & Narasimha, R. (2010), “Periodicities in Indian Monsoon Rainfall over Spectrally Homogeneous Regions”, International Journal of Climatology, 30(15), PP. 2289-2298. -Cook, B.I., Cook, E.R., Anchukaitis, K.J., Seager, R., & Miller, R. L. (2011a), “Forced and Unforced Variability of Twentieth Century North American Droughts and Pluvials”, Climate dynamics, 37(5-6), PP. 1097-1110. -Cook, B.I., Cook, E.R., Anchukaitis, K. J., Seager, R., & Miller, R. L. (2011b), “Forced and Unforced Variability of Twentieth Century North American Droughts and Pluvials”, Climate Dynamics, 37(5-6), PP. 1097-1110. -David, F.N., & Moore, P.G. (1954), “Notes on Contagious Distributions in Plant Populations”, Annals of Botany, 18(1), PP. 47-53. -Douglas, J. B. (1975), “Clustering and Aggregation”, Sankhyā: The Indian Journal of Statistics, Series B, PP. 398-417. -Duhan, D., & Pandey, A. (2013), “Statistical Analysis of Long Term Spatial and Temporal Trends of Precipitation during 1901–2002 at Madhya Pradesh”, India. Atmospheric Research, 122, PP. 136-149. -Garcia, J.A, A.Serrano and M.Cruz Gallego. (2002), “A Spectral Analysis of Iberian Peninsula Monthly Rainfall Theory”, Appl.Climatol.71, PP. 77-95. -Gemmer, M., Becker, S., & Jiang, T. (2004), “Observed Monthly Precipitation Trends in China 1951–2002”, Theoretical and Applied Climatology, 77(1-2), PP. 39-45. -Green, R.H. (1966), “Measurement of Non-randomness in Spatial Distributions”, Researches on Population Ecology, 8(1), PP. 1-7. -Hartmann, S. Becker, and L. King, (2008), “Quasi-periodicities in Chinese Precipitation Time Series”, Theory. Appl. Climatology, 92, pp. 155–163. -Hoerling, M., & Kumar, A. (2003), “The Perfect Ocean for Drought”, Science, 299(5607), PP. 691-694. -Kahya, E., & Kalaycı, S. (2004), “Trend Analysis of Stream Flow in Turkey”, Journal of Hydrology, 289(1), PP. 128-144. -Kalaycı, S., Karabörk, M. Ç. & Kahya, E. (2004), “Analysis of El Nino signals on Turkish Stream Flow and Precipitation Patterns Using Spectral Analysis”, Fresenius Environmental Bulletin, 13(8), PP. 719-725. -Kane, R. P. and Teixeira. N R. (1991), “Power Spectrum Analysis of the Annual Rainfall Series for Massachusetts (NE. U.S.A)”, Climatic Change, 18: pp. 89-94. -Kendall, M. (1975), “Multivariate Analysis” Charles Griffin. -Lana, M.D. Martı´nez, C. Serra, and A. Burguen. (2005), “Periodicities and Irregularities of Indices Describing the Daily Pluvial Metric Regime of the Fabra Observatory (NE Spain) for the Years 1917–1999”, Theory. Appl. Climatology. 82, PP. 183–198. -Liang, L., Li, L., & Liu, Q. (2010), “Temporal Variation of Reference Evapotranspiration during 1961–2005 in the Taoer River Basin of Northeast China”, Agricultural and Forest Meteorology, 150(2), PP. 298-306. -Lloyd, M. (1967), “Mean Crowding”, The Journal of Animal Ecology, PP. 1-30. -Mann, H.B. (1945), “Nonparametric Tests Against Trend”, Econometrical: Journal of the Econometric Society, PP. 245-259. Marengo, J.A. (2004), “Interdecadal Variability and Trends of Rainfall across the Amazon Basin”, Theoretical and Applied Climatology, 78(1-3), PP. 79-96. Maslen, D., & Rockmore, D. (1997), “Separation of Variables and the Computation of Fourier Transforms on Finite Groups”, I. Journal of the American Mathematical Society, 10(1), PP. 169-214. -Mitchell Jr, J.M., Dzerdzeevskii, B., Flohn, H., Hofmeyr, W.L., Lamb, H.H., Rao, K.N., & Wallen, C.C. (1966), “Climatic Change: Technical Note No. 79”, Report of a Working Group of the Commission for Climatology; WMO No. 195 TP 100. World Meteorological Organisation, Geneva, Switzerland, 81. -Morisita, M. (1959), “Measuring of the Dispersion of Individuals and Analysis of the Distributional Patterns”, Mem. Fac. Sci. Kyushu Univ. Ser. E, 2(21), PP. 5-235. -Mosmann, A. Castro, R. Fraile, J.Dessens, J.L. Sa´nchez. (2004), “Detection of Statistically Significant Trends in the Summer Precipitation of Mainland Spain”, Atmospheric Research 70, PP. 43–53. -Richard O. Gilbert. (1987), “Statistical Methods for Environmental Pollution Monitoring”, John Wiley & Sons. -Schubert, S., Gutzler, D., Wang, H., Dai, A., Delworth, T., Deser, C., & Kirtman, B. (2009), “A US CLIVAR Project to Assess and Compare the Responses of Global Climate Models to Drought-related SST Forcing Patterns: Overview and Results”, Journal of Climate, 22(19), PP. 5251-5272. -Selwam and R.Joshi. (1995), “Universal Spectrum for Internal Variability in Codas Global Air and Sea-surface Temperature”, Interntional Journal Climatology, 15, PP. 613-623. -Sen, P. K., (1968), “Estimates of the Regression Coefficient Based on Kendall’s Tau”, Journal of the American Statistical Association, 63, PP. 1379-1389. -Shifteh Some'e, B., Ezani, A., & Tabari, H. (2012), “Spatiotemporal Trends and Change Point of Precipitation in Iran”, Atmospheric Research, No. 113, PP. 1-12. -Tan, J., Jakob, C., Rossow, W.B., & Tselioudis, G. (2015), “Increases in Tropical Rainfall Driven by Changes in Frequency of Organized Deep Convection”, Nature, 519(7544), PP. 451-454. -Tang, B., Tong, L., Kang, S., & Zhang, L. (2011), “Impacts of Climate Variability on Reference Evapotranspiration over 58 Years in the Haihe River Basin of North China”, Agricultural Water Management, 98(10), PP. 1660-1670. -Thiel, H., (1950), “A Rank-invariant Method of Linear and Polynomial Regression Analysis: Part3”, Proceeding of Koninalijke Nederland's Academies van Weinenschatpen A. 53, PP. 1397-1412. -Tomozeiu, R., Busuioc, A., Marletto, V., Zinoni, F., & Cacciamani, C. (2000), “Detection of Changes in the Summer Precipitation Time Series of the Region Emilia-Romagna, Italy”, Theoretical and Applied Climatology, 67(3-4), PP. 193-200. -Torrence and Webster. (1999), “Inter Decadal Changes in the ENSO–Monsoon System”, Climatology, 12, PP. 2679-2690. -Tošić, I., & Unkašević, M. (2005), “Analysis of Precipitation Series for Belgrade”, Theoretical and Applied Climatology, 80(1), PP. 67-77. -Vivekanandan, N. (2007, July), “Analysis of Trend in Rainfall Using Non Parametric Statistical Methods”, In International Symposium on Rainfall Rate and Radio Wave Propagtion (ISRR'07), Vol. 923, No. 1, PP. 101-113, AIP Publishing. -Westra, S., Alexander, L.V., & Zwiers, F. W. (2013), “Global Increasing Trends in Annual Maximum Daily Precipitation”, Journal of Climate, 26(11), PP. 3904-3918. -Yatagai, A., Kamiguchi, K., Arakawa, O., Hamada, A., Yasutomi, N., & Kitoh, A. (2012), “APHRODITE: Constructing a Long-term Daily Gridded Precipitation Dataset for Asia Based on a Dense Network of Rain Gauges”, Bulletin of the American Meteorological Society, 93(9), PP. 1401-1415. -Zhang, Q., Liu, C., Xu, C.Y., Xu, Y., & Jiang, T. (2006), “Observed Trends of Annual Maximum Water Level and Streamflow during Past 130 Years in the Yangtze River Basin”, China. Journal of Hydrology, 324(1), PP. 255-265. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
آمار تعداد مشاهده مقاله: 3,412 تعداد دریافت فایل اصل مقاله: 1,613 |