یعقوب دین پژوه

هدف این مطالعه بررسی روند تغییرات برخی از پارامترهای هواشناسی (میانگین سرعت باد، حداقل، حداکثر و میانگین دما، اختلاف دمای حداقل و حداکثر هوا، تبخیر از تشت، میانگین رطوبت نسبی، ساعات آفتابی، مقدار بارش، تعداد روزهای همراه با بارش بیش از 1، 5 و10 میلی متر، تعداد روزهای همراه با گرد و غبار، تعداد روزهای برفی، تعداد روزهای با دید کم تر از 2 کیلومتر و درجه روز بر مبنای دمای 18 و 21 درجه سانتی گراد) در ایستگاه تبریز در دوره آماری 2005-1970 در سه مقیاس زمانی ماهانه، فصلی و سالانه با استفاده از روش ناپارامتریک مان-کندال (MK) م ی باشد. تحلیل روند با دو روشMK مرسوم (برای سری های با ضرایب خود همبستگی غیرمعنی دار) و MK اصلاح شده (برای سایر سری ها) در سه سطح معنی داری 1%، 5% و 10% انجام شد. شیب خط روند با روش سن به دست آمد. نتایج نشان داد که در مقیاس سالانه پارامترهای تبخیر از تشت، حداقل، حداکثر و میانگین دما و درجه روز بر مبنای 21 درجه سانتی گراد دارای رون د افزایشی معنی دار در سطح 5% بوده ولی بارش، تعداد روزهای همراه با بارش بیش از 1 و 5 میلی متر و درجه روز بر مبنای 18 درجه سانتی گراد دارای روند کاهشی معنی دار (در سطح 5%) بودند. بقیه متغیرها در مقیاس سالانه روند معنی دار نشان نداده اند. شیب خط روند بارش سالانه 28/2- میلی متر در سال و شیب نظیر برای هر سه عنصرمیانگین، حداقل و حداکثر دما معادل 06/0 درجه سانتی گراد در سال بود. روند مشاهده شده در پارامترهای مورد بررسی بیانگر تغییرات اقلیمی روی داده در تبریز می باشد.

در این مطالعه پنج روش تخمین پارامترهای توزیع ویبول در مورد داده های باد ارزیابی شده است. بهترین روش با استفاده از آزمون کی دو معرفی شد. برای این منظور تعداد 6 ایستگاه سینوپتیک در استان آذربایجان شرقی با دوره مشترک آماری 23 سال (1987 تا 2009) انتخاب شدند. محاسبات بر روی داده های سه ساعته سرعت باد در ارتفاع 10 متری انجام شد و پارامترهای شکل و مقیاس برای هر ایستگاه با روشهای مختلف محاسبه گردید. با استفاده از تابع تبدیل، پارامترهای فوق در دو ارتفاع 20 و 40 متری به دست آمد. بر اساس نتایج حاصله در مقایسه با سایر روش ها، روش گشتاورها که مقدار آماره کی دو کمتری به دست می داد روش مناسبی تشخیص داده شد. سپس با استفاده از مقادیر پارامتر شکل (c) و مقیاس (k) حاصل از روش فوق، خصوصیات مربوط به انرژی و سرعت باد در دوره های بازگشت 10، 25، 50 و 100 سال در سه ارتفاع 10،20 و40 متری از سطح زمین به دست آمد. در ارتفاع 10 متری از سطح زمین بالاترین مقدار پارامتر k (22/1) و c (35/3 متر بر ثانیه) در مقیاس سالانه در ایستگاه تبریز و کمترین مقدار پارامتر k و c به ترتیب برابر با 73/0 و 5/1 متر بر ثانیه در ایستگاه میانه مشاهده شد. در بین ایستگاه های مورد مطالعه، ایستگاه جلفا از نظر پارامترهای مربوط به انرژی باد نسبت به سایر ایستگاه ها وضعیت مطلوب و پتانسیل خوبی جهت بهره برداری از انرژی باد را داراست و ایستگاه اهر رتبه دوم را جهت استفاده از انرژی باد به خود اختصاص داده است. تحلیل سرعت باد در دوره های بازگشت مختلف نشان داد که بیشترین سرعت باد در استان آذربایجان شرقی متعلق به ایستگاه جلفا می باشد.

هدف این مطالعه بررسی روند تغییرات سرعت باد سطحی در ایران در دو مقیاس ماهانه و سالانه با روش ناپارامتری من کندال با چهار ویرایش متفاوت شامل الف) روش من کندال مرسوم (MK1)، ب) من کندال پس از حذف اثر ضریب خودهمبستگی مرتبه اول معنی دار (MK2)، ج) من کندال پس از حذف اثر کامل ضرایب خودهمبستگی معنی دار (MK3) و د) من کندال با لحاظ کردن ضریب هارست (MK4) می باشد. داده های مورد استفاده در این مطالعه شامل میانگین سرعت باد در دو مقیاس ماهانه و سالانه در دوره آماری 2005-1966 می باشد. بدین منظور از اطلاعات 22 ایستگاه در سطح کشور ایران استفاده شد. برای تخمین شیب خط روند تغییرات سرعت باد، از روش ناپارامتری تخمین گر سن استفاده شد. نتایج نشان داد که گر چه ترکیبی از روندهای مثبت و منفی در ایستگاه های کشور برای میانگین سرعت باد در هر دو مقیاس سالانه و ماهانه تجربه شده است، ولی تعداد ایستگاه های با روند منفی در مقایسه با تعداد ایستگاه های با روند مثبت بیش تر است. در مقیاس سالانه شدیدترین شیب خط روند منفی متعلق به ایستگاه فسا معادل 074/0- متر بر ثانیه در سال و شدیدترین شیب خط روند مثبت متعلق به ایستگاه زابل معادل 141/0 متر بر ثانیه در سال بود. در مقیاس ماهانه شدیدترین شیب خط روند منفی متعلق به ایستگاه فسا معادل 1/0- متر بر ثانیه در سال و شدیدترین شیب خط روند صعودی معادل 18/0 متر بر ثانیه در سال برای ایستگاه زابل مشاهده شده است. نتایج نشان داد که میانه شیب های خط روند اکثر ایستگاه های ایران در تمام ماه های سال (بجز فوریه و نوامبر) منفی است. بنابراین، در حالت کلی می توان نتیجه گرفت که در اکثر ایستگاه های ایران روند سرعت باد در تمام ماه های سال (بجز ماه های فوریه و نوامبر) نزولی است.

تبخیر- تعرق پتانسیل ( ) یکی از مهمترین اجزاء چرخه هیدرولوژیکی برای برنامه ریزی سیستم های آبیاری، آماده سازی داده های ورودی برای مدل های هیدرولوژیکی بیلان آب و ارزیابی اثرات هیدرولوژیکی تغییر اقلیم است. در این مطالعه مقادیر با استفاده از روش بلانی- کریدل برای 6 ایستگاه واقع در بخش جنوبی حوضه رود ارس در دوره آماری 2008-1986 محاسبه شد. افزون بر این، روند تغییرات و برخی از مهمترین متغیرهای هواشناسی موثر بر آن در دو مقیاس فصلی و سالانه تجزیه و تحلیل شد. برای این کار از آزمون ناپارامتری مان کندال پس از حذف اثر خودهمبستگی داده ها استفاده شد. برای محاسبه شیب خط روند از روش تخمین گر Sen استفاده شد. نتایج نشان داد که روند در ایستگاه ها و مقیاس های زمانی مختلف متفاوت است. در بین ایستگاه های مورد مطالعه خوی شدیدترین روند افزایشی معنی دار (P<0.05) را در مقیاس سالانه برای تجربه کرده بود. در این ایستگاه روند تغییرات حداقل رطوبت نسبی کاهشی و شدیدتر از سایر ایستگاه ها بود. در این ایستگاه پارامترهای مربوط به درجه حرارت هوا و سرعت باد روند افزایشی نشان دادند. شیب خط روند در این ایستگاه معادل 25 میلیمتر در سال به دست آمد. روند در ایستگاه های ماکو و اهر نیز در مقیاس سالانه افزایشی (ولی غیرمعنی دار) بود. در این ایستگاه روند نزولی رطوبت نسبی هوا تجربه شده است. روند سالانه سایر ایستگاه ها کاهشی (ولی غیر معنی دار) بود. در مقیاس فصلی ترکیبی از روندهای مثبت و منفی در ایستگاه های مورد بررسی مشاهده شده است. روند ایستگاه خوی در همه فصول سال مثبت (و در پاییز معنی دار در سطح 5%) و در ایستگاه پارس آباد منفی (ولی غیر معنی دار) بوده است. در بین پارامترهای مورد بررسی، تعداد ساعات آفتابی در تمام ایستگاه ها و در دو مقیاس فصلی و سالانه دارای روند افزایشی است. حداقل، حداکثر و میانگین درجه حرارت هوا در مقیاس سالانه در همه ایستگاه ها روند صعودی دارد. بارش فصل پاییز در همه ایستگاه ها روند نزولی دارد. روند سرعت باد در همه ایستگاه ها (بجز جلفا و پارس آباد) در مقیاس سالانه صعودی است.

روند تغییرات 27 شاخص مرتبط با درجه حرارت و بارش های حدی در تبریز، طی دوره آماری 2011-1961 با استفاده از آزمون های ناپارامتریک من-کندال و تخمینگر شیب sen مورد بررسی قرار گرفت. همچنین توزیع مقادیر کرانه ای تعمیم یافته به رویدادهای حدی مشاهداتی برازش داده شد و مقادیر چندک متغیرها به ازای دوره های بازگشت مختلف تخمین زده شد. نتایج نشان داد که طی نیم سده گذشته، روند شاخص های مربوط به روزهای تابستانی، شب های حاره ای و شب های گرم افزایشی بوده و از نظر آماری معنی دار می باشد. شاخص های روزهای یخی و روزهای سرد دارای روند رو به پایین و معنی دار می باشد. کمترین مقدار سالانه درجه حرارت حداکثر و حداقل روزانه روند افزایشی معنی داری را نشان می دهد. در شاخص های تعداد روزهای با بارش مساوی یا بیش از 10 و 15 میلی متر، روزهای تر متوالی، مجموع بارش روزهای تر و مجموع بارش مواقعی که بارندگی بیش از صدک 95 و 99 است، روند نزولی معنی دار تجربه شده است. منحنی های رشد و نمودارهای Q-Q پس از برازش توزیع مقادیر کرانه ای تعمیم یافته به مقادیر سالانه دمای حداقل و حداکثر روزانه و بارش حداکثر روزانه در تبریز ترسیم گردید. بررسی نمودارها نشان می دهد که تابع توزیع مذکور قابلیت زیادی در برازش سری داده های حدی حتی در قسمت انتهایی توزیع دارد.

با توجه به اهمیت بارش برف در تأمین آب در مناطق کوهستانی، برآورد دقیق آب معادل برف و همچنین تغییرات سطح پوشش آن، در بخش های کشاورزی، انرژی، مدیریت مخزن و هشدار سیل مؤثر است. در این مطالعه، به منظور برآورد رواناب حاصل از ذوب برف در حوضه شهرچای، نخست سطح پوشش برف برای سال آبی 92 با استفاده از تصاویر روزانه ماهوارهترا -مودیس با تفکیک مکانی 1 کیلومتر در محیط نرم افزار ENVI استخراج شد. آنگاه با استفاده از الحاقیه Hec_GeoHMS در نرم افزارGIS، مشخصات فیزیوگرافی حوضه به دست آمد. در مرحله بعد، با وارد کردن داده های پوشش برف، متغیرهای هواشناختی و شاخص های لازم به مدلSRM، رواناب ناشی از ذوب برف شبیه سازی شد. نتایج نشان داد، سهم جریان رودخانه از ذوب برف در ماه های فروردین و اردیبهشت چشمگیر است، ولی با افزایش درجه حرارت هوا در خرداد ماه، سهم باران در جریان پررنگ تر است. همچنین نتایج شبیه سازی بیانگر دقت بالای این مدل است، به طوری که ضریب تعیین (R 2 )معادل 9/0 و درصد خطای حجمی آن (D V )96/1 به دست آمد.

باد جزء منابع انرژی پاک و تجدید پذیر به شمار می آید. در دهه ای که گذشت، استفاده از انرژی باد در جهان با استقبال فراوان همراه بوده است. در این پژوهش شش روش (گرافیکی، تجربی، گشتاورها، عامل الگوی انرژی، حداکثر راست نمایی و گشتاورهای وزنی احتمالاتی) برآورد پارامترهای توزیع ویبول در پنج ایستگاه سینوپتیک در استان های اردبیل و زنجان بررسی شد و یافته ها نشان داد، روش گشتاورهای مرسوم روش مناسب تری است. تحلیل پتانسیل انرژی باد در ارتفاع 10، 20 و 40 متری انجام گرفت و پارامترهای پتانسیل انرژی باد (چگالی توان باد، چگالی انرژی باد، سرعت باد دارای حداکثر انرژی و محتمل ترین سرعت باد) محاسبه شد. بالاترین مقدار شاخص شکل (بی بعد) 26/1 از ایستگاه زنجان در ماه فوریه به دست آمد و بالاترین مقدار تخمینی برای شاخص مقیاس برابر با 76/4 متر بر ثانیه بوده که در ماه فوریه ایستگاه اردبیل مشاهده شد. بررسی های مربوط به پتانسیل انرژی باد نشان داد ایستگاه اردبیل پتانسیل بالایی برای بهره برداری از انرژی باد دارد. در این ایستگاه چگالی توان باد در مقیاس سالانه و در ارتفاع 10 متری برابر با 59/285 وات بر متر مربع است که در ارتفاع 40 متری به 491 وات بر متر مربع می رسد.

هدف این مطالعه تحلیل روند کیفیت آب رودخانه های استان آذربایجان شرقی با روش های ناپارامتری می باشد. بدین منظور، از اطلاعات داده های غلظت یون هایCa ،Mg، Cl، Na، HCo3، pH، SO4، TDS، EC،SAR ، درصدسدیم، مجموع آنیون ها، مجموع کاتیون هادرده ایستگاه آب سنجی (1387-1362) استفاده شد. برای آزمون روند از روش من-کندال پس از حذف اثر ضرایب خودهمبستگی معنی دار استفاده شد. به منظور تخمین شیب خط روند از روش تخمین گر Sen استفاده شد. نتایج نشان داد که در اغلب ایستگاه های مورد بررسی غلظت یون های مثبت و EC دارای روند افزایشی بود. در ایستگاه هایی که میزان رواناب آن ها روند منفی معنی دار داشته است روند تغییرات غلظت اغلب عناصر شیمیایی موجود در آب های سطحی مثبت و معنی دار بود. بررسی کیفیت آب رودخانه ها با روش ویلکاکسن نیز نشان داد که کیفیت آب رودخانه ها نسبت به دهه های گذشته افت کرده است. دلیل این افت به کاهش دبی رودخانه ها و تخلیه آلاینده های مختلف به رودخانه ها نسبت داده شد.

با توجه به اهمیت بارش برف در تأمین آب مناطق مختلف بویژه در نواحی کوهستانی، برآورد دقیق آب معادل برف و همچنین تغییرات سطح پوشش آن می تواند در بخش های کشاورزی، انرژی، مدیریت مخزن و هشدار سیل مؤثر باشد. در این مطالعه، رواناب ناشی از ذوب برف در حوضه شهرچای تحت شرایط تغییر اقلیم، بررسی شد. برای این منظور نخست، سطح پوشش برف برای سال آبی 92-91 با استفاده از تصاویر روزانه ماهواره ترا- مودیس در محیط نرم افزار ENVI استخراج گردید. سپس در نرم افزار GIS مشخصات فیزیوگرافی حوضه بدست آمد. در مرحله بعد با وارد کردن داده های پوشش برف، متغیرهای هواشناختی و پارامترهای لازم به مدلSRM ، رواناب ناشی از ذوب برف شبیه سازی شد. سپس خروجی 6 مدل گردش عمومی جو تحت سه سناریوی A 2 , A 1 B و B 1 با مدل آماری LARS-WG ریزمقیاس شدند، سپس از طریق مقایسه خروجی مدل ها در دوره آتی با دوره پایه بصورت ماهانه توسط شاخص های آماری بهترین مدل و سناریو برای تولید داده های دمای هوا و بارش در دوره آماری 2030-2011 انتخاب شد. در نتیجه برای تولید بارش از مدل HADCM3 تحت سناریوی A 1 B و برای تولید دما نیز از مدل MPEH5 و تحت سناریوی A 2 استفاده شد. برای برآورد میزان تغییر رواناب ناشی از ذوب برف در شرایط تغییر اقلیم میزان تغییرات داده های دما و بارش ماهانه تحت مدل و سناریوهای منتخب در دوره آتی با دوره پایه در طول دوره شبیه سازی به مدل SRM وارد شد. نتایج نشان داد که رواناب ناشی از ذوب برف در اواخر فصل بهار کاهش می یابد. همچنین دبی اوج در مقایسه با دوره پایه زودتر ظاهر شده و مقدار آن نسبت به دوره پایه بیشتر پیش بینی شد.

نخستین گام در مدیریت علمی منابع آب، شناخت دقیق جنبه های مختلف بارش است. یکی از جنبه های اصلی مطالعات مربوط به بارندگی، تهیه منحنی هاف در هر منطقه، به ویژه در مناطق خشک و نیمه خشک است. این کار به تصمیم گیرندگان در مدیریت مؤثر منابع آب کمک می کند که به نوبه خود به مدیریت مؤثر منابع آب شیرین منتهی می شود. هدف پژوهش حاضر، رسم منحنی های هاف و هیتوگراف رگبار در استان خوزستان است. بدین منظور از داده های باران نگار ثبات در ایستگاه های ایذه، ایدنک، عبدالخان و اهواز استفاده شد. در کل تعداد 1811 رگبار ثبت شده در بازه زمانی 1345 تا 1395 در فصول مختلف در 5 کلاس زمانی دسته بندی شدند. برای هر دسته، در هر فصل، منحنی های هاف با استفاده از تجمیع رگبارها رسم شد. نتایج نشان داد دسته بندی در فصول و در کلاس های زمانی، معنا دار است. منظور مقایسه الگوی ریزش رگبارهای ایستگاه ها، سه شاخص S، I و Q تعریف شدند که نسبت بارش تجمعی را (به ترتیب در 25، 50 و 75 درصد مدت دوام رگبار) از منحنی احتمال 1/0 به مقدار نظیر احتمال 5/0 در نظر می گیرند. نتایج نشان داد شاخص S بیش از I و هر دو بیش از شاخص Q بوده است. بازه تغییرات شاخص S، از 63/9 (برای مدت دوام کمتر از 2 ساعت در بهار و ایستگاه عبدالخان) تا 47/1 (برای بارش های 12 تا 24 ساعت فصل بهار در ایستگاه ایدنک) بود؛ در حالی که شاخص Q از 44/1 (برای بارش های 6 تا 12 ساعت در بهار و ایستگاه عبدالخان) تا 03/1 (برای بارش های با مدت دوام کمتر از 2 ساعت در بهار و زمستان در ایستگاه اهواز و بهار در ایدنک) تغییر می کند. تغییرات مقدار شاخص I در بازه 12/2 (0 تا 2 ساعت در زمستان و در عبدالخان) تا 07/1 (2 تا 6 ساعت بهار اهواز) بود. برای هر منحنی هاف، هیتوگراف رگبار طرح با استفاده از منحنی احتمال 50% استخراج و ارائه شد.

تبخیر- تعرق گیاه مرجع یکی از مهم ترین مولفه های هیدرولوژیکی بوده و تخمین مناسب آن در پژوهش های بیلان آب، مدیریت منابع آب و بهبود مصرف آب در بخش کشاورزی کاربرد دارد. در این مطالعه به بررسی اثرات مستقیم و غیرمستقیم پارامترهای هواشناسی بر روی ET 0 در استان آذربایجان غربی با استفاده از تجزیه ی علیت پرداخته شد. بدین منظور از آمار هواشناسی ایستگاه سینوپتیک استان آذربایجان غربی با بیش از 20 سال آمار، استفاده شد. جهت برآورد ET 0 ، از فرمول FAO-PM56 استفاده گردید. در ادامه مقدار ET 0 با استفاده از مدل رگرسیون چندگانه خطی نیز تخمین زده شد. سپس به منظور ارزیابی عملکرد مدل رگرسیونی، از آماره های جذر میانگین مربعات خطا ( RMSE )، میانگین قدرمطلق خطا ( MEA )، ضریب سات- ناش کلایف ( NSH ) و ضریب تعیین ( R 2 ) استفاده گردید. آنگاه ضریب همبستگی بین ET 0 و هر یک از پارامترهای هواشناسی به دست آمد و سپس تاثیر اثرات مستقیم ( P ) و غیر مستقیم هر یک از پارامترهای هواشناسی بر روی ET 0 با استفاده از روش تجزیه علیت محاسبه شد. جذر میانگین مربعات خطا مابین 0/05و 0/15و میانگین قدر مطلق خطا مابین 0/04و 0/12، ضریب سات- ناش کلایف مابین 0/91و 0/99و ضریب تعیین مابین 0/91و 0/99، بدست آمد. همچنین مشخص شد که سرعت باد در تمامی ایستگاه های مورد مطالعه، همبستگی معنی داری (در سطح 0/01درصد) با ET 0 دارد. از نتایج تجزیه علیت مشخص شد که بیشترین مقدار P (اثر مستقیم پارامترهای هواشناسی بر ET 0 ) در همه ایستگاه ها متعلق به سرعت باد می باشد که مقدار آن برای ایستگاه های ارومیه؛ 0/85، پیرانشهر؛ 0/99،تکاب؛ 0/97، خوی؛ 0/90، سردشت؛ 06/1؛ ماکو؛ 0/64و مهاباد؛ 0/78می باشد.

یکی از مدل های تخمین ET0 که مورد قبول همه ی هیدرولوژیست ها و اقلیم شناسان می باشد، مدل فائو پنمن مانتیس (FAO56PM) است. این مدل، علیرغم تخمین دقیق ET0 در اقلیم های مختلف، محدودیت هایی دارد. از جمله این محدودیت ها نیاز به داده های متنوع هواشناسی از جمله تابش آفتاب، ساعات آفتابی، دمای هوا، رطوبت نسبی، دمای نقطه شبنم، سرعت باد و فشار بخار واقعی است. متاسفانه همه ی این پارامترهای در جمیع موارد در دسترس نیستند. بنابراین، هدف اصلی این مطالعه یافتن روشی ساده جهت تخمین ET0 با داده های کم هواشناسی و دقت کافی در حوضه ی دریاچه ی ارومیه می باشد. برای این کار از مدل های هارگریوز، کیمبرلی پنمن، پریستلی تایلور و مدل های رگرسیونی خطی و غیرخطی استفاده شد. ارزیابی مدل ها با سه معیار ضریب تعیین (R2)،جذر میانگین مربعات خطا) (RMSE و میانگین خطای مطلق ) MAE ( به انجام رسید. نتایج نشان داد که میانه مقادیر R2 مدل کیمبرلی- پنمن بیش از 986/0 به دست آمده است. در حالی که میانه های مقادیر R2 برای مدل های PT و HG بترتیب معادل 902/0 و 4/0 به دست آمد. در مدل HGمیانه مقادیر RMSE حدود 9/0 (میلیمتر بر روز) بود این رقم برای مدل های KPM و PT به ترتیب، حدود3/1 و 1/2 (میلیمتر بر روز) بود. میانه آماره MAE ایستگاه های منتخب، برای مدل KPM از 1 (میلی متر بر روز) نیز کمتر بود. این رقم در مورد مدل HG حدود 7/0 (میلی متر بر روز) و برای مدل PT بیشتر از 5/1 (میلی متر بر روز) بود. بنابراین، با توجه به معیار MAE نیز مدل HG مناسبتر از مدل های دیگر تشخیص داده شد. مدل رگرسیون خطی با 2 متغیر (Tmax, Wind) دارای میانه مقادیر R2 معادل 84/0 در بین سایر مدل ها مناسب تر و دقیق تر می باشد.

در این پژوهش اثرهای مستقیم و غیرمستقیم پارامترهای هواشناسی بر روی ET0 در استان آذربایجان شرقی با استفاده از تجزیه علیت بررسی شده است. برای برآورد ET0، از فرمول فائو- پنمن مانتیث استفاده شد. مهم ترین پارامترهای هواشناسی مؤثر بر ET0، با استفاده از رگرسیون گام به گام شناسایی، به منظور ارزیابی عملکرد مدل رگرسیونی، از آماره های MAPE، R2، RMSE، و MAE استفاده شد و اثرهای مستقیم و غیر مستقیم هر یک از پارامتر بر ET0 با استفاده از تجزیه علیت محاسبه شد. مقدار MAPE مابین 43/0 و 87/0، R2 مابین 97/0 و 99/0، RMSE مابین 042/0 و 092/0، و MAE مابین 033/0 و 057/0 به دست آمد. سرعت باد در ایستگاه های مورد مطالعه (به جز ایستگاه اهر) همبستگی معنی داری با ET0 نشان داد. با توجه به نتایج تجزیه علیت، بیشترین مقدار اثر مستقیم پارامترهای هواشناسی بر ET0 در همه ایستگاه ها به جز اهر متعلق به سرعت باد بود که مقدار آن برای تبریز 637/0، جلفا 787/0، سهند 877/0، سراب، 578/0، مراغه، 850/0، و میانه 780/0 است و در ایستگاه اهر متعلق به پارامتر Tmax معادل 398/0 بود. کمترین مقدار اثر مستقیم پارامترهای هواشناسی بر ET0 در منطقه مورد مطالعه متعلق به حداقل درجه حرارت هوا بوده است.

ه دف این مطالعه، شبیه سازی رواناب حاصل از ذوب ب رف در حوضه ی ش هرچای، با استفاده از دو مدل SRM و HEC-HMS است. بدین منظور، ابتدا با وارد کردن داده های پوشش برف، متغیرهای هواشناختی و پارامترهای لازم به مدلSRM ، رواناب ناشی از ذوب برف شبیه سازی شد. در مدل HEC-HMS نیز پس از ایجاد مدل حوضه با نرم افزار HEC-GeoHMS در محیطGIS  و ایجاد مدل هواشناختی و واردکردن پارامترهای لازم مانند تلفات، روندیابی و ذوب برف شبیه سازی انجام شد. ضریب تعیین (R2) مدل SRM معادل 9/0 و درصد خطای حجمی آن (DV) 96/1 به دست آمد. از طرفی مدل HEC-HMS نیز با دقت رضایت بخش (ولی کمتر از SRM) شبیه سازی رواناب ناشی از ذوب برف را انجام داد. به طوری که ضریب تعیین ( R2) آن 85/0 و درصد خطای حجمی آن (DV) 1/2 درصد به دست آمد. صرف نظر از دقت بالای مدل SRM یکی از نقاط قوت مدل HEC-HMS، استفاده از داده های دما و بارش است؛ در صورتی که مدل SRM علاوه بر این پارامترها نیازمند داده های مساحت سطح پوشش برف است که از تصاویر ماهواره ای به دست می آید. نتایج نشان داد که مقدار کل رواناب در منطقه ی  مورد مطالعه با مدل SRM ، معادل 106×9/129 متر مکعب بود. در حالی که این رقم با مدل HEC-HMS معادل با 106×6/129 متر مکعب به دست آمد. با مقایسه ی  این ارقام با مقدار مشاهداتی (یعنی 106×4/132 متر مکعب) می توان نتیجه گرفت که عملکرد مدل SRM در مقایسه با مدل HEC-HMS نسبتأ خوب است.