تحلیل روند تغییرات هیدروکلیماتولوژی حوضه آبریز رودخانه آجی چای(تلخه رود)
الموضوعات :علی ساجد 1 , دکتر امیر گندمکار 2 , دکتر محسن باقری بداغ آبادی 3
1 - دانشجوی دکترا اقلیم شناسی، گروه جغرافیا، واحد نجف آباد، دانشگاه آزاد اسلامی، نجف آباد، ایران
2 - گروه جغرافیا، واحد نجف آباد، دانشگاه آزاد اسلامی، نجف آباد، ایران
3 - موسسه تحقیقات خاک و آب، سازمان تحقیقات، آموزش و ترویج کشاورزی، کرج، ایران
الکلمات المفتاحية: روند, دبی, من-کندال, حوضه آجی چای,
ملخص المقالة :
باتوجهبه نقش مهم و ضروری رودخانهها در اقتصاد، صنعت و کشاورزی، مدیریت نادرست این منابع و استفاده ناصحیح از آبهای سطحی و عدم شناخت کامل از میزان منابع موجود، ما را بر آن داشت تا تأثیر منفی اثرات تغییرات اقلیمی بر منابع آبی که باعث ایجاد و تداوم خشکی و خشکسالی و اثرات زیستمحیطی بسیاری را موجب گردیده و خواهد داشت در حوضه آبریز رودخانه آجی چای بررسی نماییم. به دلیل وسعت زیاد حوضه دریاچه ارومیه برای منطقه مطالعاتی فقط حوضه آجی چای در نظر گرفته شد. تحقیق حاضر با این نگاه به این مهم میپردازد و در این زمینه، اثرات و عوامل تأثیرگذار بر میزان بارندگی و دبی رودخانه آجی چای را بررسی میکند. دبی و بارش در حوضه رودخانه آجی چای واقع در استان آذربایجان شرقی در مقیاسهای ماهانه، فصلی و سالانه از سال 1348 تا سال 1396 موردمطالعه قرار گرفت جهت بررسی و تحلیل وجود روند و یا عدم وجود روند از آزمون ناپارامتری من - کندال استفاده گردید. همچنین میزان ضریب همبستگی و ضریب تعیین با استفاده از رگرسیون خطی و ضریب تغییرات در نرمافزار xlstat استفاده گردید. نتایج تحلیل دبی در حوضه آبریز آجی چای نشان میدهد که جریان رودخانه آجی چای در مقیاس سالانه در ایستگاه ونیار روند نزولی با شیب تند داشته است. همچنین روند نزولی معناداری در دادهای دبی در مقیاس فصلی در فصول پاییز، زمستان و بهار مشاهده شد که در آن شدیدترین روند مربوط به فصل زمستان با سطح اطمینان 95 درصد و با 86 درصد ضریب همبستگی است؛ ولی در فصل زمستان دادهها دارای جهش هستند ولی روند خاصی برای تابستان متصور نیست.
1) اسفندیاری درآباد، فریبا (1392). آشکار سازی آماری تاثیر پدیده گرمایش جهانی بر ناهنجاری های دبی حوضه رودخانه ارس، پژوهش های ژئومورفولوژی کمی، 4، 60-43.
2) ایمانی امیرآباد، سمیه، فرخ نیا، اشکان، مرید، سعید، و روزبهانی، رضا (1399). بررسی روند تغییرات زمانی-مکانی دما و بارش در حوضه آبریز طشک-بختگان. مهندسی عمران امیرکبیر (امیرکبیر)، 52(11)، 2931-2944.
3) برهانی داریان، علیرضا، و حاتمی مجومرد، شادی (1394). کاربرد مدلهای داده محور در مدیریت حوزه آبریز: دریاچه ارومیه. وزارت علوم، تحقیقات و فناوری، دانشگاه صنعتی خواجه نصیرالدین طوسی، دانشکده عمران.
4) بهره مند، عبدالرضا (1392). تحلیل روند تغییرات بلند مدت بارندگی و دبی در غرب دریاجه ارومیه، پژوهشنامه و مدیریت حوضه آبخیز،4(8)، 57-43.
5) برون، اشرف، ظهوریان پردل، منیژه، لشکری، حسن، شکیبا، علیرضا، و محمدی، زینب. (1404). شناسایی امواج گرمایی استان خوزستان و تحلیل همدیدی نقش پرفشار عربستان در ایجاد آن ها. اندیشه های نو در علوم جغرافیایی، 8(3)، 1-20.
6) جاهدی اسفنجانی، ناهیده، و قربانی، محمد علی (1394). تحلیل روند تغییرات بارندگی و دبی ایستگاه های دوست بیگلو و سامیان حوضه آبریز قره سو، جغرافیا و برنامه ریزی، 19(52)، 63-43.
7) رحیمی، لیلا، دهقانی، امیراحمد، قربانی، خلیلی، و عبدالحسینی، محمد (1393). بررسی روند تغییرات دبی کل و دبی پایه ایستگاه هیدرومتری اراز کوسه «حوضه آبخیز گرگان رود استان گلستان»، پژوهش های حفاظت آب وخاک، 21(2).
8) رضایی، علی (1388). بررسی اثر گسترش اراضی دیم بر روند تغییرات دبی سیلابی سالانه، مطالعه موردی، حوضه قزل اوزن، مهندسی و مدیریت آبخیز، 2 (1)، 11-17.
9) سازمان جغرافیایی نیروهای مسلح (۱۳۸۰). فرهنگ جغرافیایی شهرستان تبریز، چاپ اول، سازمان جغرافیایی نیروهای مسلح.
10) سوره، احسان، زنگنه، محمدصادق، و کرامت، اکرم (1404). بررسی روند خشکسالی با استفاده از تصاویر ماهواره ای (مظالعه موردی: استان آذربایجان غربی)، نشریه علمی اندیشه های نو در علوم جغرافیایی، دوره 3(7)، 19-38.
11) عزیزی، قاسم، و حنفی، علی. (1389). برآورد حداکثر بارش محتمل (PMP) حوضه آجی چای به روش سینوپتیکی. مطالعات جغرافیایی مناطق خشک، 1(2)، 55-71.
12) عساکره، حسین ( 1386 ). تغییر اقلیم، انتشارات دانشگاه زنجان.
13) علیجانی، بهلول (1388). اقلیم شناسی سینوپتیک، انتشارات سمت.
فضلی فرد، پوریا و همکاران، (1398). بررسی روند تغییرات بارش و دبی در حوضه دریاچه ارومیه در طول جهار دهه اخیر، نشریه دانش آب و خاک، جلد29، شماره4،.
14) فضلی فرد، پوریا، بردی شیخ، واحد، سعدالدین، امیر و حصاری، بهزاد (1398). بررسی روند تغییرات بارش و دبی در حوضه آبریز دریاچه اورمیه در طول چهار دهه اخیر، دانش آب و خاک، 29(4)، 27-41.
15) قویدل رحیمی، یوسف (1389). آشکارسازی آماری اثر گرمایش جهانی بر ناهنجاری های بارش سالانه جلفا با استفاده از شبکه های عصبی مصنوعی، جغرافیا و برنامه ریزی محیطی، 38(2)، 82-65.
16) کنعانی، رضا، فاخری فرد، احمد، قربانی، محمد علی و دین پژوه، یعقوب (1398). تحلیل روند تغییرات تدریجی و سریع عوامل هیدرواقلیمی حوضه آبریز (مطالعه موردی: حوضه لیقوان چای)، دانش آب و خاک، 29(1).
17) معروفی، صفر، و طبری، حسین (1390). آشکار سازی روند تغییرات دبی رودخانه مارون با استفاده از روش های پارامتری و ناپارامتری، فصلنامه تحقیقات جغرافیایی، 62(2)، 17119-17141.
18) نوری، مهراب، مرید، سعید، کریمی، نعمتالله و غلامی، هوشنگ (1400). تغییرات مکانی و زمانی روند در دما و بارش حوضه آبریز فرامرزی ارس، تحقیقات منابع آب ایران، 17(3)، 104-117.
19) Booij, M.J.,Tollenaar, D., van Beek, E., Kwadijk, J. C.J. (2011). Simulating Impacts of Climate Change on River Discharges in the Nile Basin, Physics and Chemistry of the Earth, 36(13), 696-709.
20) Birsan , Marius-Victor, Dumitrescu , Alexandru, Magdalena Micu , Dana, Cheval, Sorin, 2014, Changes in annual temperature extremes in the Carpathians since AD 1961, Nat Hazards, 74:1899–1910.
21) Chen, H.; Yang, J.; Ding, Y.; Tan, C.; He, Q.; Wang, Y.; Qin, J.; Tang, F.; Ge, Q.(2022) Variation in Extreme Temperature and Its Instability in China. Atmosphere Journal, 13(19).
22) Fischer, T, Gemmer, M, Luliu, L, Buda, S. (2010). Temperature and precipitation trends and dryness pattern in the Zhujiang River Basin, south china, 2007-1961, Quaternary International, pp 1-11.
23)Gadedjisso-Tossou, A.; Adjegan, K.I.; Kablan, A.K.M. (2021). Rainfall and Temperature Trend Analysis by Mann–Kendall Test and Significance for Rainfed Cereal Yields in Northern Togo, Sci, 3(17), 1-20.
24) Khamidov, Sardor, Li, Zhi, Nasirova, Makhliyo, Pulatov, Bakhtiyor, Pulatov, Alim (2023). Assessment of temperature and precipitation trends in Kashkadarya, Uzbekistan, E3S Web of Conferences 365, 01005.
25) Kliengchuay, W., Mingkhwan, R., Kiangkoo, N. et al. (2024). Analyzing temperature, humidity, and precipitation trends in six regions of Thailand using innovative trend analysis. Sci Journal, Rep 14, 7800.
26) Salnikov, V.; Talanov, Y.; Polyakova, S.; Assylbekova, A.; Kauazov, A.; Bultekov, N.; Musralinova, G.; Kissebayev, D.; Beldeubayev, Y (2023). An Assessment of the Present Trends in Temperature and Precipitation Extremes in Kazakhstan. Climate Journal,11(2), 33.
27) Zhen, Yu, Xilin, Li (2015). Recent trends in daily temperature extremes over northeastern China (1960–2011), Quaternary International, 380–381, 35–48.
نشريه علمي اندیشههای نو در علومجغرافیایی، دوره 3، شماره 8، تابستان 1404، صفحات: 111-131 شاپا: 1473-2981
|
تحلیل روند تغییرات هیدروکلیماتولوژی حوضه آبریز رودخانه آجی چای(تلخه رود)
علی ساجد
دانشجوی دکترا اقلیم شناسی، گروه جغرافیا، واحد نجف آباد، دانشگاه آزاد اسلامی، نجف آباد، ایران
امیر گندمکار11
گروه جغرافیا، واحد نجف آباد، دانشگاه آزاد اسلامی، نجف آباد، ایران
محسن باقری بداغ آبادی
موسسه تحقیقات خاک و آب، سازمان تحقیقات، آموزش و ترویج کشاورزی، کرج، ایران
چکيده
باتوجهبه نقش مهم و ضروری رودخانهها در اقتصاد، صنعت و کشاورزی، مدیریت نادرست این منابع و استفاده ناصحیح از آبهای سطحی و عدم شناخت کامل از میزان منابع موجود، ما را بر آن داشت تا تأثیر منفی اثرات تغییرات اقلیمی بر منابع آبی که باعث ایجاد و تداوم خشکی و خشکسالی و اثرات زیستمحیطی بسیاری را موجب گردیده و خواهد داشت در حوضه آبریز رودخانه آجی چای بررسی نماییم. به دلیل وسعت زیاد حوضه دریاچه ارومیه برای منطقه مطالعاتی فقط حوضه آجی چای در نظر گرفته شد. تحقیق حاضر با این نگاه به این مهم میپردازد و در این زمینه، اثرات و عوامل تأثیرگذار بر میزان بارندگی و دبی رودخانه آجی چای را بررسی میکند. دبی و بارش در حوضه رودخانه آجی چای واقع در استان آذربایجان شرقی در مقیاسهای ماهانه، فصلی و سالانه از سال 1348 تا سال 1396 موردمطالعه قرار گرفت جهت بررسی و تحلیل وجود روند و یا عدم وجود روند از آزمون ناپارامتری من - کندال استفاده گردید. همچنین میزان ضریب همبستگی و ضریب تعیین با استفاده از رگرسیون خطی و ضریب تغییرات در نرمافزار xlstat استفاده گردید. نتایج تحلیل دبی در حوضه آبریز آجی چای نشان میدهد که جریان رودخانه آجی چای در مقیاس سالانه در ایستگاه ونیار روند نزولی با شیب تند داشته است. همچنین روند نزولی معناداری در دادهای دبی در مقیاس فصلی در فصول پاییز، زمستان و بهار مشاهده شد که در آن شدیدترین روند مربوط به فصل زمستان با سطح اطمینان 95 درصد و با 86 درصد ضریب همبستگی است؛ ولی در فصل زمستان دادهها دارای جهش هستند ولی روند خاصی برای تابستان متصور نیست.
کلمات کلیدی: روند، دبی، من - کندال، حوضه آجی چای
مقدمه
برآورد میزان دبی به دلیل استفاده در طراحی سازههای بزرگ آبی بهخصوص سدها موردتوجه هیدرولوژیستها است. نخستین گام و ضروریترین اقدام در اغلب فعالیتهای بشری شناخت اقلیم یک منطقه است. در این میان به دلیل رشد روزافزون جمعیت و افزایش نیازهای آبی، اجرای پروژههای کوچک و بزرگ هیدرولوژیکی بهمنظور ذخیرهسازی و یا انحراف یا انتقال آب و نیز مطالعات آبخیزداری هر چه بیشتر مطرح میشود اجرای چنین طرحهایی نیاز به مطالعات دقیق و تجزیهوتحلیل دادههای هیدروکلیماتولوژی دارد (سازمان جغرافیایی نیروهای مسلح،1380). دریاچه ارومیه بزرگترین دریاچه داخل کشور است که در سالهای اخیر تراز سطح آب آن به میزان چشمگیری کاهشیافته است. تغییر اقلیم، خشکسالیهای مستمر، توسعه سطح زیر کشت، تغییر کاربری اراضی و احداث سدها و... از دلایل این مشکلات شناخته شده است. ولی مستندات و ادله علمی کمتری ارائه کردند (عزیزی و همکاران، 1389). به دلیل وسعت زیاد حوضه دریاچه ارومیه برای منطقه مطالعاتی، فقط حوضه آجی چای در نظر گرفته شد که شاید بیشترین تغییرات را از این منظر داشته باشد (عزیزی و همکاران، 1389). تحقیق حاضر با این نگاه به این مهم میپردازد و در این زمینه اثرات و عوامل تاثیرگذار بر میزان دبی رودخانه آجی چای را بررسی میکند مطالعه و بررسی حوضههای دریاچه ارومیه به این دلیل اهمیت پیدا میکند که این دریاچه شور بوده و میزان شوری در وضعیت عادی دریاچه ارومیه باید 250 گرم نمک در لیتر داشته باشد که این رقم در حال حاضر به 500 گرم در لیتر رسیده است (برهانی داریان، 1394). به دلیل اهمیت و تأثیر میزان ورودی آبهای رودخانهها بر شرایط محیطی و نیز نقش آن در برنامهریزیهای مبتنی بر دانستههای اقلیمی، بررسی تغییرات دبی رودخانهها در سالهای اخیر موردتوجه قرار گرفته است. تغییرات اقلیم و نوسان دما و بارش یکی از مسائل مهم زندگی بشر بهحساب میآید، با بررسی روند تغییرات میانگین بارش میتوان تحولات اقلیمی یک منطقه را ردیابی نمود. بر اساس شواهد موجود این تغییرات در جنبههای مختلف اقلیم درگذشته وجود داشته و اکنون نیز این مسئله متحمل است. مسئله ای که این جا مهم است سرعت و ماهیت تغییر پارامترهایی اقلیمی در سالهای اخیر متفاوت بوده و شتاب بیشتری به خود گرفته و روند آن با گذشته کمی متفاوت شده است. این مسئله یعنی سرنوشت و آینده اقلیم، ذهن بشریت را به خود مشغول نموده است، این دلمشغولی هم نسبت به نسل آینده و هم به سرنوشت زمین مربوط میشود. اثبات وقوع پدیده تغییر اقلیم در سطح جهانی بهسهولت امکانپذیر نیست و نیازمند بررسیهای جامع و طولانیمدت آمارهای ثبت شده از پارامترهای اقلیمی است، هر چند روند گرمتر شدن دمای سطح زمین و افزایش غلظت گازهای گلخانهای تقریباً امری است که در حال حاضر بیشازپیش مشهود گشته است. نتایج مطالعات گستردهای که در سطوح ملی، منطقهای و جهانی صورتگرفته، بیانگر افزایش دما در بسیاری از نقاط جهان و بهطورکلی افزایش میانگین دمای هوا است. نتایج این مطالعات نشان میدهد دمای حداقل و حداکثر در بازدههای زمانی گذشته رفتارهای متفاوتی از خود نشان دادهاند و دمای حداقل در بسیاری از نقاط جهان به طور آشکار افزایشی بوده، همچنین دمای حداکثر هم در خیلی از نقاط جهان روند افزایشی را نشان میدهد، ولی بهخاطر افزایش کمتر در مقایسه با دمای حداقل، موجب کاهش دامنه شبانهروزی دما شده است (عساکره، 1386). ایران کشوری با آبوهوای خشک تا نیمهخشک است و به علت دارابودن اشکال گوناگون ناهمواریها، نسبت به دیگر کشورهای خاورمیانه از آبوهوای متنوعی برخوردار است. البته این حالت تنها ناشی از تنوع ناهمواریها نبوده؛ بلکه بیشتر از جریانهای جوی در مقیاس جهانی و سینوپتیکی نشأت میگیرد. خصوصیات جغرافیای کشور پهناور ایران، ویژگیهای اقلیمی متفاوتی را به وجود آورده است. این ویژگیها حاصل فعالیتهای طبیعی و انسانی است. از آن جایی که این پدیدهها دورهای هستند، بررسی آنها میتواند در توجیه و پیشبینی سالانه وضعیت اقلیمی کشور کمک بزرگی نماید. بررسیها حاکی از آن است که همزمان با تغییر در میزان دمای حداقل و حداکثر کشور تغییراتی در بارش ایجاد میشود. پژوهشهای انجام شده در سطح ایران بیانگر افزایش تدریجی در دمای متوسط است. بااینحال گزارشهایی نیز مبنی بر کاهش درجهحرارت وجود دارد. این افزایش و کاهش منجر به تغییرات اقلیم و الگوهای آبوهوایی در کشور شده که بر جنبههای مختلف فعالیتهای ساکنان این مرزوبوم تأثیرگذار بوده است. در ایران دمای هوا علاوه بر پدیدههای کلان اقلیمی بهشدت تابع ارتفاع، عرض جغرافیایی و محتوای رطوبتی جو است. اثر ارتفاع بر دمای هوا بسیار چشمگیرتر و صدها برابر اثر عرض جغرافیایی است. حاکمیت پر فشار جنب حاره بهویژه در دوره گرم سال در ایران مرکزی و حاکمیت پرفشار سیبری در دوره سرد سال در شمال شرق، بر دمای ایران اثر چشمگیری دارد. در بخشهای میانی ایران فقر رطوبتی جو حاکمیت یک رژیم گرمایی قارهای را آسان میسازد. مجموعه شرایط یاد شده موجب شده است تا ایران دارای اقلیمی گرم با رژیم دمایی قارهای عرض میانه باشد. بهنحویکه میتوان گفت نقاط مختلف کشور عملاً دارای فصول متفاوتی هستند. مسئله اصلی در این جا بررسی روند تغییرات دما و بارش و به سبب آن تغییرات دبی در حوضه رودخانه آجی چای است. در حقیقت این پژوهش به دنبال پاسخگویی به این مسائل است که آیا میزان بارش در این محدوده در طی دورهای که سنجش عناصر اقلیمی در ایران آغاز شده، تغییریافته یا خیر؟ افزایش و کاهش آن در چه بازههای زمانی بوده و میزان دبی در این بازهای زمانی افزایش یا کاهش داشته است؟ افزایش و کاهش دما چه میزان بوده است؟ افزایش و کاهش در کدام پارامترهای بارش صورتگرفته است؟
برای رسیدن به جواب این سؤالات پارامترهای دبی روزانه و ماهانه رودخانه آجی چای در ایستگاه آبسنجی ونیار و ایستگاه آبسنجی آخولا، بارش حداکثر، حداقل، دمای حداکثر، حداقل، میانگین دمای سالانه و ماهانه ایستگاههای هواشناسی موجود در حوضه رودخانه آچی چای در طول دوره آماری ۵۰ساله مورد تحلیل قرار گرفته است.
پیشینه تحقیق
فیسچر وهمکاران(2010) درمطالعه روندهای دما وبارش درحوضه رودخانه ژوجیانگ در جنوب چین نشان داد در دماهای سالانه و دورههای گرم روندهای معنیدار مثبت در مقابل در دادههای بارش و دورههای سرد روندهای منفی مشاهده شده است. بوئیچ و همکاران (2011) در پژوهشی با شبیهسازی اثر تغییر اقلیم بر دبی رودخانه نیل با استفاده از مدل GCM به این نتیجه رسیدند که تغییرات اقلیمی اثرات چشمگیری بر آبدهی حوضه رود نیل داشته و آبدهی آن را طی سالهای آتی کاهش خواهد داد. بیرسان و دومیترسو (2014) تغییرات دماهای حدی کارپاتیانزا را بررسی و به این نتیجه دست یافتند که پارامترهای حدی سرد از روند کاهشی و شاخصهای حدی گرم از روند افزایشی برخوردار بودهاند. ژن و شیلین (2015) روند دماهای حدی در چین شمالی را بررسی و دریافت که روندهای افزایشی در روزهای تابستان و شبهای حارهای دماهای حداقل و حداکثر وجود داشته است. نتایج مطالعه گاددجیسو و همکاران (2021) نشاندهنده روند افزایشی دمای حداکثر در کارا، مانگو و نیامتوگوآ و روند کاهشی در داپانگ بوده است. چن و همکاران (2022) روند تغییرات دماهای حدی در چین را بررسی کردند. نتایج نشان میدهد که حداکثر سالانه دمای روزانه و مقدار حداقل سالانه دمای روزانه در چین به ترتیب از سال 1966 تا 2015 مقدار 18/0 و 52/0 درجه سانتیگراد در سال افزایشیافته است. سالنیکو و همکاران (2023) روندهای دما و بارش در قزاقستان را بررسی کردند. نتایج نشان داد که قلمرو قزاقستان بر اساس تغییرات آبوهوایی به نیمه جنوب غربی و شمال شرقی تقسیم شده است. خامیدو و همکاران (2023) روندهای دما و بارش در کشکادریا ازبکستان را بررسی کردند. نتایج نشان داد که میانگین دما از 53/14 درجه سانتیگراد طی سالهای 2006-2010 به 14.40 درجه سانتیگراد طی سالهای 2011-2015 کاهشیافته است. کلینگچوری و همکاران (2024) روند دما، رطوبت و بارش در شش منطقه تایلند با استفاده از تحلیل نوآورانه روند را بررسی کردند. نتایج روند من کندال نشان داد که هر شش منطقه یکروند صعودی در شاخص دما و رطوبت داشتند. درحالیکه رطوبت نسبی و بارش هر دو روند صعودی و نزولی را در مناطق مختلف نشان دادند. رضایی (1388) به بررسی اثر گسترش اراضی دیم بر روند تغییرات دبی سیلابی سالانه حوضه آبخیز قزلاوزن در یک دوره ۳۰ساله به روش کندال تائو تنها در مورد دو ایستگاه هواشناسی گیتو و سلامت آباد و دو ایستگاه آب سنجی سرچم و دهگلان پرداخته است. نتایج نشان میدهند روند افزایشی و معنیدار آماری تغییر دبیهای اوج و بارندگی یکروزه بیشینه نسبت به زمان وجود داشته و در بقیه موارد تغییرات آن فاقد هماهنگی معنیدار با روند زمان بهصورت مثبت یا منفی است. قویدل رحیمی (1389) بامطالعه دادههای مربوط به ناهنجاریهای دمایی کره زمین و بارش متوسط سالانه ایستگاه جلفا، در طی دوره آماری 1960-2003 را با استفاده از شبکههای عصبی مصنوعی و تعیین ضریب همبستگی پیرسون، تحلیل مؤلفه روند سریهای زمانی نشان داد که همبستگی منفی و معکوس معنیداری بین بارش سالانه شهرستان جلفا و ناهنجاریهای دمایی کرة زمین است. این بدان معنی است که غالباً با منفیشدن ناهنجاریهای دمایی کرة زمین بارش سالانه شهرستان جلفا افزایشیافته و ترسالی به وقوع میپیوندد و برعکس، با مثبتشدن ناهنجاریهای دمایی کره زمین بارش سالانه جلفا کاهشیافته، خشکسالی به وقوع میپیوندد. معروفی و طبری (1390) در حوضه رودخانه مارون تغییرات دبی رودخانه را با استفاده از روشهای پارامتری و ناپارامتری مورد با آزمونهای من - کندال و سن و همچنین تحلیل پارامتری رگرسیونی با استفاده از دادههای دبی رودخانه در پنج ایستگاه آبسنجی و طی یک دوره ۲۰ساله مورد ارزیابی قرار دادند. نتایج تحقیق نشان میدهد که دبی سالانه در همه ایستگاهها دارای روند نزولی در دو دهه اخیر بوده است. بهرهمند و همکاران (1392) با استفاده از آمار ۳۰ساله 6 ایستگاه آبسنجی در یک دوره ۳۰ساله و با استفاده از نرمافزار R و با روشهای آماری ناپارامتریک من - کندال تغییرات بلندمدت بارندگی و دبی در غرب دریاچه ارومیه را مورد بررسی قرار دادند و نتایج نشان داد که بارش در بعضی ایستگاهها بدون روند و در بعضی دیگر دارای روند کاهشی است. اما دبی در بیشتر موارد دارای روند کاهشی معنیداری است. اسفندیاری و همکاران (1392) بر اساس روش همبستگی پیرسون و رگرسیون خطی و غیرخطی در بازه زمانی ۴۰ساله شامل دادههای سالانه و فصلی گرمایش جهانی و دادههای دبی رودخانه ارس پرداختند. مطالعات نشان از ارتباط معکوس و قوی گرمایش جهانی با آبدهی حوضه ارس دارد. این فرایند بهخصوص از سال 1994 به بعد آشکار است و با افزایش گرمایش جهانی آبدهی حوضه ارس کاهش چشمگیر از خود نشان میدهد. رحیمی و همکاران (1393) با استفاده از آزمون ناپارامتری من - کندال، در یک دوره آماری ۴۰ساله در حوضه آبخیز گرگان رود استان گلستان با بررسی روند تغییرات دبی کل و دبی پایه ایستگاه هیدرومتری ارازکوسه دریافتند که کاهش دبی پایه نسبت به متوسط بلندمدت آن نشاندهنده وضعیت خشکسالی است. اما کاهش متوسط دورهای دبی پایه نشان میدهد که منطقه به سمت خشکی پیش میرود و تحلیلهای معنیداری روند در سطح 5 درصد بیانگر آن است که در مقیاس سالانه، میانگین دبی پایه و دبی کل، روند نزولی داشته؛ اما حداکثر دبی پایه و حداکثر دبی کل هیچ روندی ندارند. جاهدی و قربانی (1394) تغییرات بارندگی و دبی ایستگاههای دوست بیگلو و سامیان حوضه آبریز قرهسو در استان اردبیل را در یک دوره آماری ۳۰ساله با استفاده از آزمون ناپارامتری من - کندال مورد بررسی قرار دارند و نتایج تحلیل بارندگی و دبی نشان میدهد که جریان رودخانه قرهسو در مقیاس سالانه در هر دو ایستگاه این رودخانه روند نزولی دارد. ظهرابی و همکاران (1395) با استفاده از آزمون آماری من کندال به بررسی تغییرات زمانی و مکانی متغیرهای اقلیمی دما و بارش در حوضه آبریز کرخه واقع در غرب کشور ایران پرداخته است. نتایج آزمونهای آماری برای سریهای زمانی دما و بارش سالانه به ترتیب روند افزایشی و کاهشی با سطوح معنیداری بالای 1 و 5 درصد را برای اکثر ایستگاههای موردمطالعه نشان داده است. کنعانی و همکاران (1398) روند تغییرات تدریجی و سریع عوامل هیدرو اقلیمی حوضه آبریز لیقوان چای را بررسی کردند. نتایج نشان داد مقادیر جریان رودخانه لیقوان چای (در ایستگاه هروی) در همه مقیاسهای زمانی مورد بررسی دارای روند نزولی است. در فصل بهار، تمام متغیرهای مرتبط با بارش روند کاهشی نشان دادند. در مورد متغیرهای دمای هوا، نتایج نشان داد که در اکثر مقیاسهای زمانی روند افزایشی وجود دارد که این روند در ماههای دی تا اردیبهشت بسیار معنیدار بود. همچنین، اختلاف دمای حداکثر و حداقل نیز روند افزایشی نشان داد. فضلی فرد و همکاران (1398) بهمنظور بررسی علل افت سطح آب دریاچه ارومیه، روند دو عامل دخیل در تغذیه آن، شامل بارش و دبی رودخانهها را با استفاده از آزمون من کندال موردمطالعه قراردادند. نتایج نشان داد متغیر بارندگی سالانه در برخی نقاط روند افزایشی و در برخی روند کاهشی دارد اگر چه از نظر آماری روندهای مشاهده شده در اغلب ایستگاهها معنیدار نیستند. علیرغم عدم وجود نشانههایی از روند کاهشی یکنواخت بارندگی سالانه در بیشتر ایستگاههای هیدرومتری دبی سالانه در سطح معنیداری 5درصد روند کاهشی دارد. ایمانی و همکاران (1399) به بررسی تغییرات مکانی و زمانی بارندگی و دما، در حوضه آبریز طشک - بختگان را در یک دوره آماری ۳۰ساله را با استفاده از روشهای آماری ناپارامتریک من کندال بررسی کردند و با استفاده از تغییرات زمانی و مکانی بارندگی و دما، بهعنوان پارامترهای تأثیرگذار بر وضعیت منابع آب، برای یک دوره زمانی ۳۰ساله آماری 2010 تا 1981 با استفاده از روش آماری غیرپارامتریک من کندال بررسی شد. نتایج نشان داد بارندگی در طول دوره آماری به میزان 3/14 درصد کاهش و دما به میزان 5/3 درصد افزایش داشته است که روند تغییرات برای بارش و دما به ترتیب در سال 2004 و 1985 اتفاق افتاده است. اما روند کاهش بارندگی در منطقه برخلاف افزایش دما معنیدار نیست. نوری و همکاران (1400) بهمنظور بررسی روند تغییرات دما و بارش حوضه آبریز فرامرزی ارس از آزمون ناپارامتری من - کندال و شیب سنس استفاده کردند. نتایج نشان داد که روند تغییرات بارش در مقیاس سالانه جز در قسمتهای مرکزی حوضه بهصورت افزایشی بوده است. این تغییرات در مقیاس فصلی بیشتر خود را نشان میدهد، بهنحویکه در فصل بهار روند تغییرات بهصورت کاهشی و در فصل پاییز بهصورت افزایشی بوده است. دما نیز به طور یکسان در کل منطقه مطالعاتی سیر صعودی را نشان میدهد. برون و همکاران (1404) امواج گرمایی استان خوزستان را بررسی کردند. نتایج نشان داد که در اکثر قریب بهاتفاق امواج گرمایی در لایههای زیرین و رد سپهر زبانههای چرخند پاکستان و عربستان بهعنوان مهمترین سامانه تأثیرگذار در انطباق با الگوی توپوگرافی منطقه (ارتفاعات زاگرس جنوبی) و الگوی چرخشی سامانه سبب فرارفت گرم دمای تابشی بیابانهای گرم لوت و عربستان بر روی استان خوزستان بودهاند. سوره و همکاران (1404) روند خشکسالی استان آذربایجان غربی را بررسی کردند. نتایج نشان داد در پهنه استان آذربایجان غربی شاخص VCI و SPI حساسترین شاخص برای شناسایی خشکسالی بوده که قادرند انواع خشکسالیها را بهخوبی شناسایی کنند. همچنین شاخص NDVI باتوجهبه حساسیت نسبتاً خوبی که برای شناسایی خشکسالی دارد، قادر است کاهش پوشش گیاهی در سالهای خشکسالی را نشان دهد.
مبانی نظری
تلخهرود به زبان ترکی «آجیچای» یکی از مهمترین و اصلیترین رودخانههای داخلی جاری در استان آذربایجان شرقی است. این رود همچنین یکی از بزرگترین رود تأمینکننده آب ورودی به دریاچه ارومیه است. مساحت حوضه آبریز این رودخانه ۹۲۰۰ کیلومترمربع است و حداکثر دبی ثبت شده آن در هر ثانیه ۴۰٫۶ مترمکعب است. طول رودخانه ۲۶۵ کیلومتر است و در مناطق مرکزی استان آذربایجان شرقی در جریان است (سازمان جغرافیایی نیروهای مسلح،1380). دریاچه ارومیه از اواسط دهه ۱۳۸۰ شروع به خشکشدن کرد و امروزه در خطر خشکشدن کامل قرار دارد. بررسی تصاویر ماهوارهای نشان میدهد که در سال ۲۰۱۵ دریاچه ۸۸ درصد مساحت خود را ازدستداده، گزارشهای قبلی تنها به ازدسترفتن ۲۵ تا ۵۰ درصد مساحت دریاچه اشاره کرده بودند. دلایل بسیاری برای خشکشدن دریاچه ذکر شده است از جمله خشکسالی احداث بزرگراه بر روی دریاچه و استفاده بیرویه از منابع آب حوزه آبریز دریاچه. خشکسالی یک پدیده آبوهواشناسی پیچیده بوده که متکی به چندین عامل محیطی است (سوره و همکاران، 1404). تحقیق جدیدی توسط چند تن از محققان در امریکای شمالی نشان میدهد که خشکسالی تنها باعث کاهش ۵ درصدی بارش در حوزه آبریز دریاچه شده و عوامل انسانی شامل پروژههای جاهطلبانه توسعه اقتصادی - آبی به همراه ساخت بزرگراه ۱۵ کیلومتری بر روی دریاچه با دریچه کوچک ۱/۲ کیلومتری وضعیت دریاچه را به بحران کشانیده است. تا سال ۲۰۱۲ بیش از دویست سد بر روی رودخانههای حوزه آبریز دریاچه در مرحله آماده بهرهبرداری یا پایان مراحل طراحی بودند. بر اساس آخرین بررسیها در آبان 1394 تراز آب دریاچه ارومیه 1270.04 متر عنوان شد که نسبت به زمان مشابه سال گذشته 40 سانتیمتر کاهش نشان میدهد (برهانی داریان، 1394). در صورت خشکشدن این دریاچه هوای معتدل منطقه تبدیل به هوای گرمسیری با بادهای نمکی خواهد شد و زیست محیط منطقه را تغییر خواهد داد. علاوه بر نمک بسیاری از آلودگیها شامل فلزات سمی سنگین مورداستفاده در صنعت و مواد سمی مورداستفاده در کشاورزی به آبهای سطحی و زیر سطحی مرتبط با دریاچه نفوذ کرده و در صورت خشکشدن دریاچه بسیاری از مواد سمی هوازی شده و خطرات بیماریهای تنفسی برای زیستبوم و مردم منطقه به وجود خواهد آورد. بااینحال هنوز اقدام جدی برای نجات دریاچه صورت نگرفته کارشناسان و فعالان محیطزیست اعتقاد دارند خطر خشکیدن دریاچه ارومیه تا شعاع 500 کیلومتری این دریاچه را تهدید میکند (برهانی داریان، 1394).
روششناسی تحقیق
در این پژوهش، از دادههای بلندمدت (حداقل 30 سال) ایستگاههای هواشناسی و هیدرومتری واقع در حوضه آبریز رودخانه آجی چای استفاده شده است. میانگین بارش سالانة حوضة آبریز تلخه رود حدود 355 میلیمتر، و میانگین آبدهی سالانة تلخه رود 000 ، 000 ، 392 متر مکعب است. بر روی تلخه رود و سرشاخه های آن، هفده ایستگاه آب سنجی نصب شده که سه تای آن بر روی شاخه اصلی است. يكي از مسائل بسيار مهم آب شناختي ، محاسبه ي حداكثر بارش محتمل به منظور طراحي سدهاي بزرگ است. از سال 1950 اولين بار براي طراحي سازه هاي بزرگ مقادير بارش محاسبه گرديد. به كارگيري حداكثر بارش محتمل به منظور طراحي سازه هاي آبي با ريسك يا خسارت بالا قدمتي در حدود صد سال دارد. در دسترس نبودن داده هاي دبي قابل اعتماد در طراحي سد يا مخازن سدها سبب توجه بيشتر به داده هاي بارش و بيشينه سازي داده هاي بارش براي طراحي سدها و سازه هاي آبي مدنظر قرارگرفت. در انجام چنين مطالعه اي شناسايي حد بالاي بارش براي يك دوره زماني و منطقه جغرافيايي معین داراي اهميت مي باشد.
شکل 2: نقشه هم بارش حوضه رودخانه آجی چای
مآخذ: نگارندگان،1398
دادههای روزانه و سالانه دبی حوضه رودخانه آجی چای در ایستگاه نهند، ایستگاه خواجه و ایستگاه آخوله، دمای حداقل و حداکثر حوضه رودخانه آجی چای، بارش روزانه و ماهانه حوضه رودخانه آجی چای، شامل ایستگاه جدول (1) در طول دوره آماری ۴۳ساله از دادههای موجود استخراج گردیده است.
جدول 1: مشخصات ایستگاه های هواشناسی
ردیف | نام ايستگاه | نوع ايستگاه | سازمان | طول جغرافيايي | عرض جغرافيايي | ارتفاع از سطح دريا (متر) |
1 | سهراب | بارانسنج معمولي | وزارت نيرو - تماب | 40-47 | 57-38 | 1850 |
2 | ميركوه تاجيار | تبخير سنجي | وزارت نيرو - تماب | 31-47 | 02-38 | 1830 |
3 | سراب | كليماتولوژي | سازمان هواشناسي | 34-47 | 59-37 | 1650 |
4 | سراب | بارانسنج ثبات | وزارت نيرو - تماب | 32-47 | 57-37 | 1682 |
5 | سراب | سينوپتيك | سازمان هواشناسي | 32-47 | 56-37 | 1682 |
6 | سراب تاجيار | بارانسنج معمولي | وزارت نيرو - تماب | 32-47 | 57-37 | 1680 |
7 | هريس (هارزند) | بارانسنج معمولي | سازمان هواشناسي | 31-47 | 49-37 | 1820 |
8 | سرانسر | بارانسنج معمولي | وزارت نيرو - تماب | 27-47 | 53-37 | 1700 |
9 | قزلچه سادات | بارانسنج معمولي | سازمان هواشناسي | 17-47 | 06-38 | 1844 |
10 | قوشچي | بارانسنج معمولي | سازمان هواشناسي | 16-47 | 47-37 | 1980 |
11 | کوه سالار | بارانسنج معمولي | سازمان هواشناسي | 09-47 | 37-37 | 1670 |
12 | اهر | بارانسنج ثبات | وزارت نيرو - تماب | 03-47 | 29-38 | 1310 |
13 | هريس – برازين | تبخير سنجي | وزارت نيرو - تماب | 06-47 | 16-38 | 1962 |
14 | مهربان | بارانسنجی معمولي | وزارت نيرو - تماب | 07-47 | 04-38 | 1510 |
15 | زرنق چاي | بارانسنجی معمولي | سازمان هواشناسي | 04-47 | 05-38 | 1600 |
16 | مهربان | بارانسنجی معمولي | سازمان هواشناسي | 07-47 | 04-38 | 1510 |
17 | اهر | سينوپتيك | سازمان هواشناسي | 04-47 | 26-38 | 1390 |
18 | اسب خوران | بارانسنج | وزارت نيرو - تماب | 58-46 | 20-38 | 1940 |
19 | چرو خوران | بارانسنجی معمولي | سازمان هواشناسي | 59-46 | 45-37 | 1980 |
20 | بستانآباد | كليماتولوژي | سازمان هواشناسي | 50-46 | 51-37 | 1750 |
21 | بستانآباد | بارانسنج ثبات | وزارت نيرو - تماب | 50-46 | 51-37 | 1750 |
22 | ديزناب | بارانسنج معمولي | وزارت نيرو - تماب | 44-46 | 50-37 | 1860 |
23 | آقچه كهل زمان | بارانسنج معمولي | سازمان هواشناسي | 48-46 | 45-37 | 1900 |
24 | سعيدآباد | بارانسنج معمولي | سازمان هواشناسي | 35-46 | 57-37 | 1950 |
25 | قوري گل (اماملو) | بارانسنج معمولي | وزارت نيرو - تماب | 41-46 | 53-37 | 1950 |
26 | ورزقان | بارانسنجی معمولي | سازمان هواشناسي | 36-47 | 29-38 | 1698 |
27 | نهند (خواجه) | بارانسنجی معمولي | وزارت نيرو - تماب | 36-46 | 09-38 | 1580 |
28 | نهند (محل سد) | بارانسنجی معمولي | وزارت نيرو - تماب | 29-46 | 13-38 | 1850 |
29 | سعيدآباد | بارانسنجی معمولي | وزارت نيرو - تماب | 35-46 | 59-37 | 1850 |
30 | سعيدآباد | كليماتولوژي | سازمان هواشناسي | 34-46 | 55-37 | 1760 |
31 | هروي | بارانسنج معمولي | وزارت نيرو - تماب | 29-46 | 55-37 | 1920 |
رودخانه آجی چای یکی از منابع آبی مهم در منطقه آذربایجان شرقی، نقش حیاتی در تأمین آب شرب، کشاورزی و صنعت ایفا میکند. تغییرات اقلیمی و فعالیتهای انسانی، تأثیر قابلتوجهی بر منابع آبی این حوضه داشته است. این مقاله به بررسی روند تغییرات مؤلفههای هیدروکلیماتولوژی حوضه آبریز رودخانه آجی چای، شامل بارش، دما، رواناب میپردازد. با استفاده از دادههای آماری بلندمدت (حداقل 30 سال) و روشهای آماری مناسب مانند آزمون من - کندال، روند تغییرات این مؤلفهها مورد بررسی قرار گرفته و عوامل مؤثر بر این تغییرات شناسایی شدند. بررسی رابطه بین عوامل و متغیرهای هیدرولوژیکی از خروجی مرحله کتابخانهای و میدانی تحقیق و با استفاده از نرمافزارهای آماری و گرافیکی بوده است. بهطوریکه این مرحله تحقیق با استفاده از تکنیکهای GIS و تکنیک آماری من - کندال بهکارگیری و با استفاده از بستههای نرمافزاری ARC-GIS قابلاجرا در محیط GIS،Excel و بسته نرمافزاری SPSS جهت تجزیهوتحلیل دادهها به روش همبستگی از طریق ایجاد روابط رگرسیونی ساده و چندمتغیره برای دستیابی به رابطه بین بارش، دما و میزان آبدهی حوضه رودخانه آجی چای استفاده شده است. یکی از روشهای متداول جهت تحلیل سریهای زمانی هیدرومتئورولوژیکی، بررسی وجود و یا عدم وجود روند در آنها با استفاده از آزمونهای آماری است. اصولاً وجود روند در سریهای زمانی دما و بارش ممکن است تغییرات تدریجی طبیعی و تغییر اقلیمی اثر فعالیتهای انسانی باشد. تغيير اقليم جهاني و افزايش دما و کاهش بارش و یا تغییر زمانی بارش از مسائل مهم زیستمحیطی است كه در سالهای اخير مطالعات زيادي بر روي آن صورتگرفته است. اين موضوع به دليل اثرات زیستمحیطی، اقتصادي و اجتماعي از اهميت بالايي برخوردار است چرا كه فعالیتهای انساني مانند كشاورزي، صنايع و مانند آن بر مبناي ثبات و پايداري تغيير اقليم هستند. بهطورکلی پيامدهايي نظير خشکسالیها، سیلابهای شديد و ناگهاني، امواج هواي سردوگرم، از جمله آثار و شواهد ناهنجاریهای اقليمي است كه كره زمين را با بحرانهای مختلف مواجه كرده است. بدون شناخت و آگاهي از وضعيت اقليمي حال و آينده، مديران و برنامهریزان قادر به اجراي برنامههای مختلف نخواهند بود. آبوهوای کره زمین در طول قرن بیستم، بهویژه در دو دهه اخیر تعادل خود را ازدستداده و تمایل به افزایش دما نشان داده است. (IPCC,2001a,b). تغيير اقليم، معادل تغييرات معنیدار آماري براي متوسط وضع آبوهوا در يك دوره طولاني است. اين تغييرات میتواند در متوسط دما، بارندگي، الگوهاي آبوهوایی، باد، تابش و پارامترهاي مشابه آن باشد. اقليم میتواند گرمتر و يا سردتر شود و مقادير سالانه دما میتواند افزايش و يا كاهش يابد. در مقياس كلي، افزايش تدريجي دماي كره زمين و اقیانوسها را در اثر افزايش گازهاي گلخانهای، مهمترین عامل تغيير اقليم میدانند. پارامترهاي اقليمي در مقياس زمان و مكان به دلايل زيادي تغيير میکنند كه بايد نحوه تغييرات آنها بر اساس مشاهدات با استفاده از روشهای آماري تعيين شود. تحليل روند از جمله مهمترین روشهای آماري است كه به طور گسترده براي ارزيابي اثرات بالقوه تغيير اقليم بر روي سریهای زماني هيدرولوژيكي مانند سریهای مشاهداتي دما، بارش و جريان رودخانه در نقاط مختلف جهان استفاده شده است. اثبات وجود روند معنیداری در یک سری زمانی بارندگی بهتنهایی نمیتواند دلیل قاطع بر وقوع تغییر اقلیم در یک منطقه باشد؛ بلکه فرض صفر رخداد آن را تقویت مینماید. این ویژگی ناشی از متعدد بودن امکان کنترلکننده سامانه اقلیمی است. وجود یا عدم وجود روند و تحلیل روند و تحلیل سریهای زمانی و تغییر اقلیم ارائه گردیدهاند که این روشها در دودسته روشهای پارامتریک و ناپارامتریک تقسیمبندی میشوند.
روشهای ناپارامتریک از کاربرد نسبتاً وسیعتر و چشمگیرتری نسبت به روشهای پارامتریک برخوردارند. روشهای پارامتریک عمدتاً بر اساس روشهای رابطه رگرسیونی سری دادهها با زمان استوار هستند. آزمون T استیودنت از جمله متداولترین روشهای ناپارامتریک به شمار میرود. براي سریهایی كه توزيع آماري خاصي بر آنها قابل برازش نيست و چولگي يا كشيدگي زيادي دارند استفاده از روشهای ناپارامتريك مناسبتر است. آزمون من - کندال متداولترین و پرکاربردترین روشهای ناپارامتریک تحلیل روند سریهای زمانی هیدرومتئورولوژیکی به شمار میرود. آزمون ناپارامتری من – کندال ابتدا توسط من (1945) ارائه و سپس توسط کندال (1975) بر پایه رتبه دادهها در یک سری زمانی بسط و توسعه یافت. این روش به طور متداول و گستردهای در تحلیل روند سریهای هیدرولوژیکی و هواشناسی بکار گرفته میشود. از نقاط قوت این روش میتوان به مناسببودن کاربرد آن برای سریهای زمانی و مکانی که از توزیع آماری خاصی پیروی میکنند اشاره نمود.
برای تعيين نقاط جهش میتوان با رسم نمودار سري در برابر زمان نقطه جهش را حدس زد. اما این کار از طريق روشهای آماري نيز امکانپذیر است. يكي از اين آزمونها، آزمون من كندال دنبالهای (Sneyers, 1990) است. در حقيقت اين روش مقادير آماره را در كليه زمانهای سري (i امين مرتبه) با روش رتبه دادن من كندال محاسبه میکند، كه همين عمل بهصورت عكس انجام میشود. يعني میتوان فرض كرد انتهاي سري ابتداي آن باشد، و دنبالة را بر اساس چنين سري بيان كرد. اگر دنباله بر اساس u و'u بهصورت نمودار رسم شود در حالت معنیداری روند، نقطه شروع پديده يكديگر را قطع خواهند نمود. در حالتی كه سري ايستا باشد دو دنباله u و ' u بهصورت موازي عمل خواهند نمود و يا با چند بار برخورد بهطوریکه بـه تغيير جــهش آنها منجر نشود، در خواهند آمد. اگر a = 0.05 را در نظر بگيريم هرگاه اندازه I u I بيشتر از 1.96 باشد روند معنیدار است . u+ روند افزايش و u - روند كاهشي را نشان میدهد بهعبارتدیگر اگر منحني u از خطوط بالا و پايين معنیدار (1.96) خارج شود روند وجود دارد. اگر منحني u و'u در محدوده معنیدار تلاقي بكنند ولي خارج نشوند بدين معني است كه فقط تغيير ناگهاني در ميانگين رخداده است ولي روند وجود ندارد.
ابتدا صحت دادهها با آزمون توالي يا دنبالهها صورت گرفت. در اين روش آمار به ترتيب صعودي يا نزولي مرتب شد و ميانه آنها مشخص گرديد كه مقدار آن برابر رقمي است که در وسط سري آماري قرار دارد در سریهای زوج، ميانه برابر ميانگين دو رقم وسط است سپس هريك از ارقام سري را با ميانه بهدستآمده مقايسه نموده چنانچه از آن بزرگتر باشد با علامت A و چنانچه از آن کوچکتر باشد با علامت B مشخص میگردد. براي ارقامي كه برابر ميانه هستند علامتي منظور نمیگردد و چنين حالتي مانع از ادامه توالي اعداد نمیگردد. باتوجهبه علامتهای بهدستآمده میتوان دنبالهها را كه از يك يا چندين سال متوالي با علامت A يا B مشخص شدهاند و امكان دارد كه در بين آنها سالي بدون علامت وجود داشته باشد تعيين نمود. حد اپتيمم مجموع تعداد دنبالهها را بر اساس جدول ران تست در سطح 5 درصد به دست آورد و اگر مجموع دنبالهها بين دو عدد بهدستآمده از جدول قرار گيرد میتوان آن سري داده را همگن دانست.
كفايت دادهها عامل مهمي در بررسي بوده بهنحویکه هر اندازه طول مدت آماري بيشتر باشد نتايج حاصله، از دقت بيشتري برخوردار خواهند بود كه در اين تحقيق نيز مورد بررسي قرار گرفت. براي سنجش ميزان كفايت دادهها، از رابطههای مختلف از جمله رابطه ماكوسکه N حداقل تعداد دادههای لازم، t مقدار t استيودنت که از جدول مربوطه به دست میآید و بهازای درجه آزادي برابر (N-6) به دست میآید و R نسبت مقدار عددي متغير موردنظر است میتوان استفاده كرد كه با استفاده از آن تعداد دادههای لازم را در سطح اعتماد موردنظر میتوان به دست آورد. براي تعيين تصادفي بودن دادهها از رابطه(1) زير استفاده میشود:
رابطه:(1)
كه T آماره من - كندال، N تعداد كل سالهای آماري مورداستفاده و P مجموع تعداد رتبههای بزرگتر از رديف ni كه بعد از آن قرار میگیرند بوده و از رابطه(2) زير به دست میآید:
رابطه:(2)
اين آماره براي N> 10 به توزيع نرمال با ميانگين صفر و واريانس شبيه است. در نهايت بهمنظور سنجش معنیدار بودن آماره T از رابطه(3) زير محاسبه میشود:
رابطه:(3)
كه در آن N تعداد كل سال هاي آماري، tg برابر سطح احتمال معني دار بودن آزمون و t(T) آماره من- كندال مي باشد. با توجه به مقدار بحراني بدست آمده براي t(T)، حالات مختلف بدين شرح مشاهده خواهد شد:
اگر t(T) T >باشد روند مثبت در سري زماني غالب خواهد بود، اگر t(T) <T باشد نشاندهنده روند منفي و درصورتیکه +(T ) t > T > - (T ) t باشد هیچگونه روند مهمي در سریها مشاهده نمیشود و سریها تصادفي هستند.
در اين پژوهش براي تعيين وجود يا عدم وجود روند و تعيين نوع تغييرات و زمان آن، از روش من - كندال استفاده شد. جهت شناسايي روندهاي جزئي و کوتاهمدت، نقاط جهش و نقاط شروع روند سري زماني از نمودار سري زماني بر حسب مقادير U(t) و U'(t) استفاده میگردد. براي ترسيم نمودار سري زماني مقادير متوالي، آمارهاي U(t) و U' (t) با استفاده از آزمون من - كندال محاسبه میشود. در اين روش ضريب t آزمون من - كندال از رابطه(4) زير به دست میآید:
رابطه:(4) |
|
كه تابع توزيع آن در شرايطي كه فرض صفر حاكم باشد از لحاظ مجانبي با ميانگين و واريانس برابر است.
رابطه:(5) |
|
واريانس آن برابر است با:
رابطه:(6) |
|
در صورت وجود روند در جهت مشخص اين آزمون فقط در شكل دوطرفه آن دقيق است. ازاینرو فرض صفر براي مقادير بالاي | U(t) | رد میگردد. براي بهدستآوردن U(t) از رابطه(7) زير استفاده میشود:
رابطه:(7) |
|
وقتي مقدار U(t) معنیدار است كه روند افزايشي يا كاهشي در آن مشاهده شود و در نتيجهU(t) > 0 يا U(t)< 0 حاصل گردد.
براي تعيين زمان وقوع تغيير لازم است علاوه برU(t)، مؤلفه U'(t) نيز محاسبه میشود که ابتدا بايد اميد رياضي معکوس و واريانس معکوس را از رابطه(8) و(9) و(10) زير به دست آورد:
رابطه:(8) |
|
رابطه:(9) |
|
رابطه:(10) |
|
مقدار ti برابر با مجموع niوt'i برابر با مجموع n'iمي باشد.
در اين روش، مقادير متوالي از مقدار Ui و U'i حاصله از آزمون من - كندال بهصورت گرافيكي نمايش داده میشود كه اگر مقادير Uiو U'i از منحنیها چندين بار روي همديگر قرار بگيرند روند يا تغييري وجود نخواهد داشت ولي در جايي كه منحنیها همديگر را قطع میکنند منحنیها محل شروع روند يا تغييرات را بهصورت تقريبي به نمايش میگذارند. اگر منحنیها همديگر را در داخل محدوده قطع كنند نشانه زمان آغاز تغيير ناگهاني و درصورتیکه خارج از محدوده بحراني همديگر را قطع كنند بيانگر وجود روند در سریهای زماني است.
معرفی منطقه موردمطالعه
حوضه آجي چاي بهعنوان يكي از حوضههای آبخيز رودخانهای واقع در حوضه منطقهای درياچه اروميه و از نظر موقعيت جغرافيايي محدودۀ موردمطالعه از طول جغرافیایی 46 درجه و 13 دقیقه و 22 ثانیه تا 46 درجه و 46 دقیقه و 21 ثانیه شرقی و از عرض جغرافیایی 37 درجه و 59 دقیقه و 19 ثانیه تا 38 درجه و 16 دقیقه و 6 ثانیه شمالی گسترده شده است.
دادههای روزانه و سالانه دبی رودخانه آجی چای (از سال 1348 تا سال 1392) در طول دوره آماری 43 ساله -دمای حداقل و حداکثر حوضه رودخانه آجی چای- بارش حداقل و حداکثر حوضه رودخانه آجی چای در طول دوره آماری (از سال 1965 تا سال 2017).
شکل 1: موقیت حوضه آجی چای در نقشه ایران
مآخذ: نگارندگان،1398
حوضه آبریز رودخانه آجی چای، بهعنوان یکی از حوضههای مهم ایران، از نظر اقتصادی و اجتماعی دارای اهمیت بالایی است. تأمین آب کشاورزی، شرب و صنعت در این منطقه به طور مستقیم به وضعیت هیدرولوژیکی این حوضه وابسته است. در سالهای اخیر، تغییرات اقلیمی و فعالیتهای انسانی همچون تغییر کاربری اراضی، برداشت بیرویه آبهای زیرزمینی و ساخت سدها، تغییرات قابلتوجهی در مؤلفههای هیدروکلیماتولوژی این حوضه ایجاد کرده است. این تغییرات میتواند منجر به کاهش منابع آبی، خشکسالیهای شدید، افزایش سیلابها و در نهایت، تهدید امنیت آبی منطقه شود؛ لذا، بررسی روند تغییرات این مؤلفهها و شناسایی عوامل مؤثر بر آن، برای برنامهریزی و مدیریت پایدار منابع آبی حوضه ضروری است.
بحث و یافتهها
برای انجام تحقیق حاضر دادههای بلندمدت و سالانه دبی در ایستگاه ونیار رودخانه آجی چای در سطح شهرستان تبریز از سازمان آب منطقهای استان آذربایجان شرقی اخذ و سپس دادهها در سه مقیاس ماهانه فصلی و سالانه مورد بررسی و محاسبه قرار گرفت. با استفاده از روش آماری ناپارامتریک من – کندال برای ایستگاه ونیار در طول دوره آماری از سال 1348 تا سال 1396 برای دادههای ماهانه فصلی و سالانه روند دادهها برسی و نتایج ذیل حاصل گردید.
شکل 3: بررسی روند سالانه میانگین دبی ایستگاه ونیار در بین سال های 1348 تا سال 1396
مآخذ: نگارندگان،1398
باتوجهبه اینکه در نقطه شروع نمودارها همدیگر را قطع نمودند سریها دارای روند معنادار و همچنین دارای جهش است از سال 1348 تا سال 1380 بررخود و جهش در سری ها مشاهده می شود. سری ها دارای روند می باشند ولی از سال 1380 تا سال 1396روند منفی کاهشی با شیب زیاد مشاهده می شود. ضریب تعیین 85/0 درصد نشان دهنده وجود رابطه مستقیم و قوی در رابطه حاکم است. همچنین ضریب همبستگی 92/0 تعیین کننده ضریب همبستگی قوی بین داده ها است.نمودار خطی میانگین ۵۰ساله دبی در ایستگاه ونیار از سال 1348 تا سال 1396 نشاندهنده کاهش دبی در طول دوره آماری به زیر 5 مترمکعب در سال 1386 را نشان میدهد.
شکل 4: نمودار خطی میانگین 50 ساله دبی در ایستگاه ونیار از سال 1348 تا سال 1396
مآخذ: نگارندگان،1398
این کاهش تأییدکننده نمودار سری روند کاهشی در ایستگاه ونیار در طول دورة آماری است. بیشترین میزان خروجی برای ایستگاه ونیار در طول دوره آماری در سال 1354به میزان 385/42 مترمکعب در ثانیه است و کمترین میزان ثبت شده به میزان 719/0 مربوط به میانگین سال 1386 است. دبی ایستگاه ونیار در بین سالهای 1348 تا سال 1396 با استفاده از آزمون ناپارامتریک من کندال در طول دورة آماری و تفسیر آن در میانگین سالانه و بهصورت فصلی استخراج گردیده است. بهطوریکه در شکل شماره (5) مشاهده میشود در هر چهارفصل نمودارها و سریهای زمانی پایینتر از نرمال نسبت به سال قبل مشاهده میشود.
بهار
بررسی و تحلیل روند در میانگین دبی بهار ماه در حوضه آجی چای سریها دارای روند و جهش و تغییر ناگهانی در میانگین دادهها را نشان میدهند و منحنیها با شیب تند کاهشی و معنیداری در سطح اطمینان 95درصد است. ضریب تعیین 82 درصد است و ضریب همبستگی 90درصد است که گویای همبستگی قوی در بین دادهها است و یک رابطه مستقیم و قوی در بین دادهها حاکم است.
پاییز
با استفاده از آزمون من - کندال روند و معناداری دبی رودخانه آجی چای در سهماهه پاییز در طول دوره آماری از سال 1347 تا سال 1396در ایستگاه دبی سنجی بررسی و همان طوری از شکل پیداست منحنی ها در طول دوره آماری دارای روند با سطح اطمینان 95 درصد می باشند که از سال 1370روند دارای معکوس و کاهشی بوده است. داده ها دارای ضریب تعیین 85 درصد می باشد. در معادله نیز خط رگرسیونی نیز به وضوح این کاهش را نشان می دهد. همچنین ضریب همبستگی 92/0درصد گویا همبستگی قوی بین دادها است.
شکل 5: بررسی روند سالانه میانگین سه ماهه بهار، تابستان، پاییز
و زمستان دبی ایستگاه ونیار در سالهای 1348تا سال 1396
مآخذ: نگارندگان،1398
تابستان
با استفاده از آزمون من کندال روند و معناداری دادههای میانگین دبی در ایستگاه ونیار در طول دوره آماری از سال 1348 تا سال 1396 مورد بررسی و تحلیل قرار گرفت. همانطوری از شکل پیداست منحنیها در تابستان، در طول دوره آماری دارای جهش بودند ولی روند برای آنها متصور نیست. دادها دارای ضریب تعیین 0092/0 درصد است. همچنین ضریب همبستگی 04/0 درصد گویای همبستگی ضعیف بین دادهها است.
زمستان
با بررسی و تحلیل روند دبی زمستان در حوضه آجی چای مشخص گردید که سریها دارای روند و جهش و تغییر ناگهانی در میانگین دادهها را نشان میدهند و منحنیها با شیب تند کاهشی و معنیداری در سطح اطمینان 95 درصد است. ضریب تعیین 74 درصد است و ضریب همبستگی 86 درصد است که گویای همبستگی قوی در بین دادهها است و یک رابطه مستقیم و قوی در بین دادهها حاکم است.
پاییز
در بررسی نمودار خطی میانگین دبی رودخانة آجی چای در سهماهه پاییز در طول دوره آماری از سال 1347 تا سال 1396در ایستگاه دبی سنجی ونیار مشخص است که میزان دبی نسبت به سال های نخست اندازه گیری سیر کاهشی داشته است. به طوری که حتی در دوره های نوسان در سال 1386 میانگین دبی خروجی به صفر رسیده است. همچنین بیشترین میانگین دبی خروجی در سال 1351یعنی 71/19متر مکعب در ثانیه ثبت شده است.
شکل 6: بررسی میانگین خطی سه ماهه بهار، تابستان، پاییز و
زمستان دبی ایستگاه ونیار در بین سالهای 1348تا سال 1396
مآخذ: نگارندگان،1398
زمستان
بیشترین میزان میانگین دبی برای فصل زمستان در طول دوره آماری در ایستگاه دبیسنجی ونیار در حوضه آبخیز رودخانه آجی چای مربوط به سال 1354 به میزان 66/36 مترمکعب بر ثانیه ثبت شده است. همچنین کمترین میزان دبی مربوط به سال 1385 که 319/0 مترمکعب بر ثانیه است. منحنی نمودار نشاندهنده کاهشی محسوس دبی در طول دوره آماری است.
تابستان
در بررسی نمودار خطی میانگین دبی رودخانه آجی چای در طول تابستان در طول دوره آماری از سال 1347 تا سال 1396در ایستگاه دبی سنجی ونیار مشخص است که میزان دبی نسبت به سال های نخست اندازه گیری سیر کاهشی داشته و از نوسان بشتری برخودار بوده است به طوری که در سال 1359 میانگین خروجی تابستان به بیش از 60/9 متر مکعب در ثانیه افزایش پیدا کرده است و در بقیه سال ها میزان خروجی بین صفر تا چهار متر مکعب در ثانیه در نوسان بوده است.
بهار
بیشترین میزان دبی برای بهار در طول دوره آماری در ایستگاه دبیسنجی ونیار در حوضه آبخیز رودخانه آجی چای مربوط به سال 1354 به میزان 44/120 مترمکعب بر ثانیه ثبت شده است. همچنین کمترین میزان دبی مربوط به سال 131386 که 3/1 مترمکعب بر ثانیه است. منحنی نمودار حالت کاهشی محسوسی را نشان میدهد.
نتیجهگیری و پیشنهادها
در تحقیق حاضر میزان ضریب همبستگی و ضریب تعیین با استفاده از رگرسیون خطی و ضریب تغییرات در نرمافزار xlstat استفاده گردید. نتایج تحلیل دبی در حوضه آبریز آجی چای نشان میدهد که جریان رودخانه آجی چای در مقیاس سالانه در ایستگاه ونیار روند نزولی با شیب تند داشته است همچنین روند نزولی معناداری در دادههای دبی در مقیاس فصلی در فصول پاییز، زمستان و بهار مشاهده شد که در آن شدیدترین روند مربوط به فصل زمستان با سطح اطمینان 95 درصد و با 86 درصد ضریب همبستگی است؛ ولی در فصل زمستان دادهها دارای جهش هستند ولی روند خاصی برای تابستان متصور نیست. نتایج این تحقیق نشاندهنده ارتباط قوی و معکوس آبدهی رودخانه آجی چای با افزایش گرمایش جهانی دارد. در این مطالعه مشخص گردید که تغییرات سریهای زمانی آبدهی حوضه رودخانه آجی چای از روند نزولی آبدهی نشان میدهد که در بلندمدت همواره در جهت کاهش مستمر دبی عمل کرده است. این عملکرد بهخوبی نشان میدهد که با بیشترشدن سیر صعودی گرمایش جهانی، میزان آبدهی این رود نیز همواره زیر میانگین بلندمدت خود قرار گرفته است. همچنین اگر روند گرمایش جهانی به همین روال ادامه یابد، میانگین دبی سالانه این رودخانه در آینده به کمترین حد خود خواهد رسید و برگشتن مقدار دبی رودخانه آجی چای به روال نرمال خود بعید به نظر میرسد. باتوجهبه افزایش دما در این منطقه، بارش کاهشیافته و مقدار تبخیر نیز افزایش مییابد. اثرات تغییر اقلیم و اثرات افزایش دما بر منابع آب میتواند از چند جهت موردتوجه قرار گیرد. بخش هواشناسی شاهد اولین تأثیر این اثرات است. در این شرایط با افزایش دما و جابهجایی زمان بارندگیها از فصل سرد به فصول گرم، میشود. این موضوع باعث کاهش ذخیره برفی و جابهجایی رژیم آبدهی رودخانهها میشود. بهاینترتیب دبی پایه رودخانهها و ظرفیت تنظیمی میزان خروجی آبدهی پایین آمده و تعادل در سیستم آبی مختل میشود. این موضوع همراه با تغییر الگوی بارش باعث کاهش پوشش گیاهی شده و تشدید شرایط آبوهوایی نظیر خشکسالی و سیلابها میگردد و نهایتاً منجر به کمشدن پوشش گیاهی و تشدید تنش شرایط آبوهوایی و کمشدن منابع آبی میشود و با تغییر در منابع آبی، برنامهریزی و مدیریت آب در منطقه مورد چالش قرار میگیرد. از طرف دیگر با کاهش دبی پایه ناشی از تغییر اقلیم و افزایش غلظت آلودگیها، ظرفیت خود پالایی رودخانهها و بهتبع آن کیفیت آب رودخانهها و مخازن بهشدت کاهش مییابد. همچنین به دلیل تأمین بخش اعظمی از آب مصرفی کشاورزی بهوسیله آب جاری در رودخانه آجی چای، کاهش کیفیت آبهای سطحی تأثیر مستقیمی بر عملکرد بخش کشاورزی خواهد گذاشت. از جنبههای دیگر هم بخش کشاورزی نیز در دایره تأثیر پدیده تغییر اقلیم قرار دارد بهطوریکه ازیکطرف نیاز آبی گیاهان با افزایش دما افزایش مییابد و از طرف دیگر تغییر الگوی زمانی بارشها، عدم قطعیت در تأمین نیازهای آبی گیاهان در فصل کشت توسط بارشهای مؤثر را افزایش خواهد داد. با بروز این عدم قطعیتها، نیاز تأمین آب مصرفی بخش کشاورزی توسط منابع آب استحصال شده بالا خواهد رفت. در همین راستا به دلیل کاهش منابع آب در اثر کاهش میزان بازندگی عملکرد محصولات کاهش پیدا کرده و امنیت غذایی به مخاطره میافتد. از سوی دیگر افزایش شوری خاک در اثر تبخیر باعث تسریع فرایند بیابانزایی و گسترش مناطق لمیزرع میگردد.
منابع و مآخذ
1) اسفندیاری درآباد، فریبا (1392). آشکار سازی آماری تاثیر پدیده گرمایش جهانی بر ناهنجاری های دبی حوضه رودخانه ارس، پژوهش های ژئومورفولوژی کمی، 4، 60-43.
2) ایمانی امیرآباد، سمیه، فرخ نیا، اشکان، مرید، سعید، و روزبهانی، رضا (1399). بررسی روند تغییرات زمانی-مکانی دما و بارش در حوضه آبریز طشک-بختگان. مهندسی عمران امیرکبیر (امیرکبیر)، 52(11)، 2931-2944.
3) برهانی داریان، علیرضا، و حاتمی مجومرد، شادی (1394). کاربرد مدلهای داده محور در مدیریت حوزه آبریز: دریاچه ارومیه. وزارت علوم، تحقیقات و فناوری، دانشگاه صنعتی خواجه نصیرالدین طوسی، دانشکده عمران.
4) بهره مند، عبدالرضا (1392). تحلیل روند تغییرات بلند مدت بارندگی و دبی در غرب دریاجه ارومیه، پژوهشنامه و مدیریت حوضه آبخیز،4(8)، 57-43.
5) برون، اشرف، ظهوریان پردل، منیژه، لشکری، حسن، شکیبا، علیرضا، و محمدی، زینب. (1404). شناسایی امواج گرمایی استان خوزستان و تحلیل همدیدی نقش پرفشار عربستان در ایجاد آن ها. اندیشه های نو در علوم جغرافیایی، 8(3)، 1-20.
6) جاهدی اسفنجانی، ناهیده، و قربانی، محمد علی (1394). تحلیل روند تغییرات بارندگی و دبی ایستگاه های دوست بیگلو و سامیان حوضه آبریز قره سو، جغرافیا و برنامه ریزی، 19(52)، 63-43.
7) رحیمی، لیلا، دهقانی، امیراحمد، قربانی، خلیلی، و عبدالحسینی، محمد (1393). بررسی روند تغییرات دبی کل و دبی پایه ایستگاه هیدرومتری اراز کوسه «حوضه آبخیز گرگان رود استان گلستان»، پژوهش های حفاظت آب وخاک، 21(2).
8) رضایی، علی (1388). بررسی اثر گسترش اراضی دیم بر روند تغییرات دبی سیلابی سالانه، مطالعه موردی، حوضه قزل اوزن، مهندسی و مدیریت آبخیز، 2 (1)، 11-17.
9) سازمان جغرافیایی نیروهای مسلح (۱۳۸۰). فرهنگ جغرافیایی شهرستان تبریز، چاپ اول، سازمان جغرافیایی نیروهای مسلح.
10) سوره، احسان، زنگنه، محمدصادق، و کرامت، اکرم (1404). بررسی روند خشکسالی با استفاده از تصاویر ماهواره ای (مظالعه موردی: استان آذربایجان غربی)، نشریه علمی اندیشه های نو در علوم جغرافیایی، دوره 3(7)، 19-38.
11) عزیزی، قاسم، و حنفی، علی. (1389). برآورد حداکثر بارش محتمل (PMP) حوضه آجی چای به روش سینوپتیکی. مطالعات جغرافیایی مناطق خشک، 1(2)، 55-71.
12) عساکره، حسین ( 1386 ). تغییر اقلیم، انتشارات دانشگاه زنجان.
13) علیجانی، بهلول (1388). اقلیم شناسی سینوپتیک، انتشارات سمت.
فضلی فرد، پوریا و همکاران، (1398). بررسی روند تغییرات بارش و دبی در حوضه دریاچه ارومیه در طول جهار دهه اخیر، نشریه دانش آب و خاک، جلد29، شماره4،.
14) فضلی فرد، پوریا، بردی شیخ، واحد، سعدالدین، امیر و حصاری، بهزاد (1398). بررسی روند تغییرات بارش و دبی در حوضه آبریز دریاچه اورمیه در طول چهار دهه اخیر، دانش آب و خاک، 29(4)، 27-41.
15) قویدل رحیمی، یوسف (1389). آشکارسازی آماری اثر گرمایش جهانی بر ناهنجاری های بارش سالانه جلفا با استفاده از شبکه های عصبی مصنوعی، جغرافیا و برنامه ریزی محیطی، 38(2)، 82-65.
16) کنعانی، رضا، فاخری فرد، احمد، قربانی، محمد علی و دین پژوه، یعقوب (1398). تحلیل روند تغییرات تدریجی و سریع عوامل هیدرواقلیمی حوضه آبریز (مطالعه موردی: حوضه لیقوان چای)، دانش آب و خاک، 29(1).
17) معروفی، صفر، و طبری، حسین (1390). آشکار سازی روند تغییرات دبی رودخانه مارون با استفاده از روش های پارامتری و ناپارامتری، فصلنامه تحقیقات جغرافیایی، 62(2)، 17119-17141.
18) نوری، مهراب، مرید، سعید، کریمی، نعمتالله و غلامی، هوشنگ (1400). تغییرات مکانی و زمانی روند در دما و بارش حوضه آبریز فرامرزی ارس، تحقیقات منابع آب ایران، 17(3)، 104-117.
19) Booij, M.J.,Tollenaar, D., van Beek, E., Kwadijk, J. C.J. (2011). Simulating Impacts of Climate Change on River Discharges in the Nile Basin, Physics and Chemistry of the Earth, 36(13), 696-709.
20) Birsan , Marius-Victor, Dumitrescu , Alexandru, Magdalena Micu , Dana, Cheval, Sorin, 2014, Changes in annual temperature extremes in the Carpathians since AD 1961, Nat Hazards, 74:1899–1910.
21) Chen, H.; Yang, J.; Ding, Y.; Tan, C.; He, Q.; Wang, Y.; Qin, J.; Tang, F.; Ge, Q.(2022) Variation in Extreme Temperature and Its Instability in China. Atmosphere Journal, 13(19).
22) Fischer, T, Gemmer, M, Luliu, L, Buda, S. (2010). Temperature and precipitation trends and dryness pattern in the Zhujiang River Basin, south china, 2007-1961, Quaternary International, pp 1-11.
23)Gadedjisso-Tossou, A.; Adjegan, K.I.; Kablan, A.K.M. (2021). Rainfall and Temperature Trend Analysis by Mann–Kendall Test and Significance for Rainfed Cereal Yields in Northern Togo, Sci, 3(17), 1-20.
24) Khamidov, Sardor, Li, Zhi, Nasirova, Makhliyo, Pulatov, Bakhtiyor, Pulatov, Alim (2023). Assessment of temperature and precipitation trends in Kashkadarya, Uzbekistan, E3S Web of Conferences 365, 01005.
25) Kliengchuay, W., Mingkhwan, R., Kiangkoo, N. et al. (2024). Analyzing temperature, humidity, and precipitation trends in six regions of Thailand using innovative trend analysis. Sci Journal, Rep 14, 7800.
26) Salnikov, V.; Talanov, Y.; Polyakova, S.; Assylbekova, A.; Kauazov, A.; Bultekov, N.; Musralinova, G.; Kissebayev, D.; Beldeubayev, Y (2023). An Assessment of the Present Trends in Temperature and Precipitation Extremes in Kazakhstan. Climate Journal,11(2), 33.
27) Zhen, Yu, Xilin, Li (2015). Recent trends in daily temperature extremes over northeastern China (1960–2011), Quaternary International, 380–381, 35–48.