فهرست مطالب

پژوهش های اقلیم شناسی - پیاپی 37-38 (بهار و تابستان 1398)

نشریه پژوهش های اقلیم شناسی
پیاپی 37-38 (بهار و تابستان 1398)

  • تاریخ انتشار: 1398/04/01
  • تعداد عناوین: 7
|
  • آزیتا امیری*، بهلول علی جانی، زهرا حجازی زاده، ابراهیم فتاحی، مجید آزادی صفحات 1-12
    رابطه دمای هوا در ایران با توزیع تاوایی نسبی در خاورمیانه با استفاده از مقادیر ماهانه تاوایی نسبی در دو تراز 850 و 500 هکتوپاسکال و دمای هوای 97 ایستگاه سینوپتیک کشور (2010-1981) بررسی شد. بدین منظور از روش های آماری تحلیل مولفه های اصلی (PCA) و تحلیل همبستگی کانونیکال (CCA) استفاده شد. ابتدا توسط تحلیل مولفه های اصلی، مراکز فعالیت تاوایی نسبی در خاورمیانه تعیین شدند. سپس رابطه این مراکز با دمای هوای ایران توسط تحلیل همبستگی کانونیکال بررسی شد و مناطقی که تغییرات تاوایی نسبی در آن با تغییرات دمایی در ایران رابطه دارد شناسایی شدند. همبستگی کانونیکال بین مراکز فعالیت تاوایی نسبی تراز 500هکتوپاسکال در منطقه ای شامل ایران، عراق و شمال عربستان با دمای هوا در ایران برابر با 0.95 است. همبستگی کانونیکال بین این مراکز در تراز 850هکتوپاسکال نیز با دمای کشور برابر با 0.96 می باشد که بیانگر ارتباط بسیار بالایی بین تاوایی نسبی این دو  تراز در منطقه خاورمیانه و دمای هوا در ایران است. شاخص بار متقابل یا همبستگی واقعی بین تاوایی نسبی تراز 500 و 850 هکتوپاسکال با دمای ایران به ترتیب 0.93- و 0.93 می باشد. برای تفسیر نتایج، نقشه توزیع همبستگی بین مولفه کانونیکال تاوایی نسبی با مقادیر واقعی آن در خاورمیانه ترسیم شد. مقایسه این نقشه با توزیع ماهانه تاوایی نسبی، شباهت بسیار زیاد آن را با توزیع تاوایی نسبی در ماه های ژوئن، ژوئیه و اوت در هر دو تراز  نشان می دهد. بدین معنی که تاوایی نسبی و دمای هوا در فصل گرم سال بالاترین همبستگی را با هم نشان می دهند. در فصل گرم تاوایی نسبی منفی در تراز میانی، حضور یک پرارتفاع و تاوایی نسبی مثبت در تراز نزدیک زمین، حضور یک کم فشار گرمایی سطح زمین را نشان می دهد که هر دو با افزایش دمای هوا در ایران همراهند.
    کلیدواژگان: تاوایی نسبی، دما، تحلیل همبستگی کانونیکال، خاورمیانه، ایران
  • فردین صابری *، بهلول علی جانی، شهریار خالدی صفحات 13-30
    رویداد های فرین ناشی از تغییرات آب و هوایی در پیش بینی تراز سطح دریا اهمیت زیادی دارند. در منطقه مورد مطالعه، در آینده فراوانی و شدت رویدادهای فرین دما و بارش افزایش خواهند یافت. نمایه های فرین، نشان دهنده تغییر مقادیر فرین دما و بارش نسبت به دوره پایه 2010-1981 بوده است و این امر، نشان دهنده مجموع بارش و یا دمایی بیش از صدک 95 دوره پایه است. ضریب تغیرات بارش و دما برای کل حوضه آبگیر خزر مثبت است و در ناحیه جنوبی الگوی نامنظمی بر آستانه های بارش حاکم است. افزایش سطح دریا  (SLR) یک نگرانی عمده برای نواحی ساحلی است. پیش بینی دقیق سطح دریای خزر برای آینده غیر ممکن است، اما مدل های کامپیوتری می توانند پیش بینی احتمالی تغییرات آینده را ارائه دهند. این مشکل با استفاده از یک سیستم مدل یکپارچه به نام SIMCLIM با دقت مکانی 0/1درجه در0/1 درجه در خوشبینانه ترین حالت یعنی RCP4.5 و بدبینانه ترین حالت یعنی RCP8.5 مورد بررسی قرار گرفته است. پیش بینی های انجام شده نشان می دهد که افزایش سطح دریای خزر به آرامی تا سال 2100 افزایش می یابد. اگر افزایش حداکثر پیش بینی شده رخ دهد، سواحل خزر آسیب پذیرترین جابجایی ساحلی را تجربه خواهد کرد.
    کلیدواژگان: افزایش سطح دریا، سواحل جنوبی دریای خزر، سناریو های انتشار، رویدادهای فرین، ضریب تغییرات
  • محمد مرادی* صفحات 31-42
    در این پژوهش سریعترین سرعت باد سالانه ایستگاه بوشهر در دوره آماری 2015-1951 بررسی شد و با استفاده از رابطه گامبل نوع اول، دوره برگشت چند ساله برای این کمیت بدست آمد. برای محاسبه فراسنجهای مقیاس و محل رابطه گامبل، از روش های گرافیکی، و عددی شامل کمترین مربعات، گشتاورها، تقرب آماری مرتبه ای و بیشینه درشتنمائی استفاده شد. نتایج برآورد فراسنجهای مقیاس و محل برای محاسبه دوره برگشت سریعترین باد ایستگاه بوشهر به روش گامبل نوع اول نشان داد که بر پایه میانگین مربع خطاها، روش های گشتاورها، کمترین مربعات، گرافیکی، تقرب آماری مرتبه ای و بیشینه درشتنمائی به ترتیب کمترین تا بیشترین میانگین مربع خطاها را دارا می باشند. در نتیجه روش گشتاورها برای محاسبه دوره برگشت سریعترین باد ایستگاه بوشهر در دوره آماری انتخابی از دقت بیشتری برخوردار است. بر این اساس از روش گشتاورها برای محاسبه دوره برگشت سریعترین باد سالانه ایستگاه بوشهر  استفاده شد.از محاسبات دیده شد که سریعترین باد ایستگاه بوشهر با دوره برگشت25 ساله،29/7متر بر ثانیه، 50 ساله، 8/32 متر بر ثانیه،  100 ساله، 35/8متر بر ثانیه و با دوره برگشت 1000 ساله، 45/9متر بر ثانیه برآورد شده است. بعلاوه بیشترین سرعت باد روزانه ایستگاه بوشهر39/0متر بر ثانیه است که در سال 1959 رخ داده است که دوره برگشت آن 2/207 سال است. در سال 2014 نیز سریعترین سرعت باد ایستگاه بوشهر 35 متر بر ثانیه ثبت شده است که این مقدار نیز بر اساس روش گشتاورها دوره برگشت 82/9ساله دارد.
    کلیدواژگان: کمترین مربعات، بیشینه درشتنمائی، تقرب آماری مرتبه ای، روش گشتاورها
  • فهیمه خادم پور، عباس خاشعی *، مهدی امیرآبادی زاده صفحات 43-56
    انرژی تابشی خورشید سرمنشا انرژی های روی کره ی زمین است که یکی از منابع مهم انرژی پاک به ویژه در کشور ایران محسوب می شود. برآورد درست مقدار تابش خورشیدی به عنوان یکی از پارامترهای مهم در شبیه سازی رشد گیاهان و تخمین تبخیرتعرق اهمیت زیادی دارد. هدف از انجام این مطالعه، توانایی مدل شبیه سازی GeneXproTools 5.0در پیش بینی تشعشعات خورشیدی بر اساس داده های هواشناسی (بارش، ساعات آفتابی، رطوبت نسبی، حداکثر دما، حداقل دما، میانگین دما و تابش خورشیدی، هم چنین با در نظر گرفتن ارتفاع، طول جغرافیایی و عرض جغرافیایی) در مقیاس روزانه، در 31 ایستگاه در گستره ی ایران در بازه ی زمانی 2016-2007 می باشد. بهترین مدل بر اساس معیارهای ارزیابی RMSE،MAE ،ENSو R2انتخاب شد. از نتایج تجزیه و تحلیل، مشخص شد که مقادیر پیشی بینی شده سازگاری خوبی با مقادیر اندازه گیری شده در مدل برنامه ریزی بیان ژن(GEP) دارد. در مدل GEP، سناریوی eبه دلیل در نظر نگرفتن ارتفاع، طول جغرافیایی، عرض جغرافیایی، بارش، رطوبت نسبی، حداکثر دما و حداقل دما، دارای ضریب همبستگی کمتر و خطای بیشتری می باشد و کارایی کمتری دارد. در این مدل، در سناریوی cبا اضافه شدن پارامترهای میانگین دما و تابش خورشیدی، برآورد مطلوب تری از تابش خورشیدی حاصل می شود. نتایج سناریوهای cو eبسیار به هم نزدیک می باشد، اما سناریوی eدر پیش بینی تابش خورشیدی ضعیف تر عمل می کند و بهترین مدل در این مطالعه، سناریوی cمی باشد. زمانی که عرض جغرافیایی، بارش، رطوبت نسبی، حداکثر دما، حداقل دما، میانگین دما و تابش خورشیدی برای پیش بینی تابش خورشیدی در نظر گرفته شود، نتایج دقیق تری حاصل می شود. هم چنین، معیارهای ارزیابی تحت سناریوی cبا 0/72=R2 ، 3/59=RMSE ، 2/82= MAE و 0/72= NSEدر بخش آموزش از سایر سناریوها بیشتر می باشد. در مجموع روش GEPدقیق ترین نتایج را در تخمین تابش خورشیدی روزانه در گستره ی ایران دارد.
    کلیدواژگان: تابش خورشیدی، سناریو، نرم افزار GeneXproTools
  • جواد سالاروند، فرهاد قاسمی *، زهرا اسدالهی صفحات 57-72
    فرسایندگی باران به صورت قدرت تراکمی باران در بروز فرسایش تعریف می شود. اگر سایر خصوصیات موثر بر فرسایش ثابت در نظر گرفته شود، میزان هدررفت خاک مستقیما متناسب با میزان فرسایندگی باران خواهد بود. در پژوهش حاضر با هدف تهیه نقشه فرسایندگی باران در استان لرستان، در ابتدا میزان فرسایندگی باران در 53 ایستگاه هواشناسی در سطح استان با شاخص اصلاح شده فورنیه (MF) محاسبه شد. در ادامه برای نشان دادن همبستگی مکانی میان داده های فرسایندگی باران از ترسیم واریوگرام در محیط نرم افزارGS+5.1.1 استفاده شد. در نهایت با بکارگیری روش زمین آماری کریجینگ در نرم افزار ArcGIS10.3، توزیع مکانی فرسایندگی باران نقشه سازی شد و به منظور انتخاب روش مناسب تر، دو روش کریجینگ معمولی و کریجینگ ساده با یکدیگر مقایسه شدند. بعد از انتخاب مدل واریوگرام و درون یابی انجام شده توسط روش های مورد بررسی در این تحقیق، در ادامه صحت درون یابی با روش ارزیابی متقابل مورد ارزیابی قرار گرفت. براساس نتایج بدست آمده روش کریجینگ معمولی بدلیل بالاتر بودن مقدار R2 و پایین تر بودن مقادیر میانگین خطا (ME)، میانگین استاندارد شده خطا (MSE)، ریشه دوم میانگین استاندارد شده خطا (RMSE) دقت بالاتری را در مقایسه با روش کریجینگ ساده نشان داد. روش کریجینگ معمولی با  میزان میانگین خطای برآورد پایین تر (0/06) و مربع میانگین ریشه خطا (0/4) مناسب ترین روش برای درون یابی در این پژوهش ارزیابی شد. همچنین میزان همبستگی(R2) مقادیر برآوردی و مشاهده ای فرسایندگی باران با روش کریجینگ معمولی  0/68 برآورد شد. در نهایت با استفاده از واریوگرام نمایی و روش کریجینگ معمولی نقشه فرسایندگی باران در سطح استان ترسیم شد که  براساس آن مناطق واقع در امتداد شمالی-جنوبی و مرکزی استان لرستان دارای فرسایندگی زیاد و مناطق غربی و شرقی دارای فرسایندگی کمتری هستند. نقشه فرسایندگی باران تولید شده در این پژوهش می تواند به منظور برآورد میزان هدررفت خاک با مدل USLE و شناسایی مناطق با پتانسیل فرسایش زیاد مورد استفاده قرار گیرد.
    کلیدواژگان: فرسایش، کریجینگ، واریوگرام، درون یابی
  • احمد سیفی، نرگس صالح نیا، فاطمه گل زاده خرق* صفحات 73-89
    رویکرد سنتی اقتصاد بیشتر مبتنی بر استفاده بی رویه از منابع طبیعی و عدم توجه به حقوق نسل های آتی نسبت به این منابع است. این در حالی است که آثار زیست محیطی چنین رویکردی جبران ناپذیر خواهد بود و برای دست یابی به توسعه پایدار باید به تغییر آن پرداخت.با توجه به  اهمیت حفاظت از محیط زیست و همچنین لزوم توجه به رشد و توسعه ی پایدار، در این مطالعه با تمرکز بیشتر بر دو متغیر درجه ی باز بودن تجاری و کنترل فساد، سعی در بررسی عوامل اثرگذار بر آلودگی هوا از دو منظر ساختاری و نهادی داریم. لذا هدف از این مطالعه، بررسی اثر درجه ی باز بودن تجاری و کنترل فساد بر شاخص ذرات معلق هوا با استفاده از رویکرد حداقل مربعات کاملا تعدیل شده (FMOLS) در کشورهای منتخب سه گروه درآمدی و طی دوره ی 1996-2011 است. یافته های تجربی مدل نشان می دهند که اثر متغیر درجه ی باز بودن تجاری بر آلودگی هوا در هر سه گروه از کشورها، منفی و معنی دار است. به این معنا که بالا رفتن درجه باز بودن تجاری در هر سه گروه به بهبودی کیفیت هوا و کاهش انتشار ذرات معلق می انجامد. همچنین با بالا رفتن کنترل فساد می توان شاهد کاهش انتشار ذرات معلق هوا در گروه های درآمد متوسط و درآمد پایین بود؛ در حالی که ضریب این متغیر برای گروه درآمد بالا به لحاظ آماری معنادار نیست. بنابراینکاهش سطح فساد و مبارزه با فعالیت های آلوده به آن، ارتقای کیفیت هوا را برای دو گروه درآمد متوسط و درآمد پایین به ارمغان می آورد.
    کلیدواژگان: درجه ی باز بودن تجاری، روش حداقل مربعات کاملا تعدیل شده، شاخص ذرات معلق هوا (PM10)، کنترل فساد
  • فخرالدین ایران پور*، حسن زهره وندی صفحات 91-105
    یکی از مهم ترین عوامل موثر در مدیریت صحیح منابع آب، شناخت دقیق احتمالات رخ داد بارش و نحوهتوزیع روزهای متوالی خشک و بارانی است. چنین شناختی می تواند زمینه های مناسبی برای برنامه ریزان در جهت مقابله با اثرات مخرب خشکسالی ها و نوسانات شدید بارش ارائه دهد. در این پژوهش با استفاده ازداده های مربوط به بارش های روزانه تبریز در یک دوره آماری 60 ساله(1951- 2010) که از اداره کل هواشناسی استان آذربایجان شرقی دریافت گردید، با استفاده از قوانین احتمال، به صورت فرآیندهای تصادفی و با استفاده از  مدل زنجیره های مارکوف، احتمال تداوم و تواتر روزهای بارش و خشکی، احتمال وقوع بارش و احتمال مقدار بارش، دوره برگشت و احتمال تداوم خشکی های 3، 4، 5، 6 و 7 روزه برای تمام روزهای سال در قالب هفتگی محاسبه و مورد تحلیل قرار گردید. در این بررسی احتمال وقوعبارش هفتگی 22 درصد و احتمال عدم وقوع بارش 78 درصد بدست آمد. نتایج حاصل از محاسبه ماتریس تغییر وضعیت برای روزهای مختلف سال در قالب هفته ای نشان می دهد که در مجموع 17168 روز، روز خشک و4792روز بارندگی وجود داشته است. به عنوان نمونه طی 119 روز آمار موجود از اولین هفته فروردین در این دوره 60 ساله 191 روز، روزهای خشکی است که بعد از روز خشک رخ داده است. بالاترین احتمال وقوع بارش در فصل بهار و هفته با 46 درصد می باشد و خشکی های 7 روزه دارای کمترین احتمال تداوم برای تمام هفته های سال بدست آمد.
    کلیدواژگان: روزهای خشک متوالی، زنجیره مارکوف، احتمالات بارش، تبریز
|
  • Azita Amiri *, Bohlol Alijani, Ebrahim Fattahi, Mahid Azadi, Zahra Hejazizadeh Pages 1-12
    Recognizing the interactions of the large-scale components of atmospheric circulation creates an appropriate capacity to examine the regional climate variability and improve description of land-atmosphere connections. Assessment Iran vulnerability to climate fluctuations is very important. Part of this recognition will be achieved by assessing the components of the atmospheric circulation. The relative vorticity distribution is an important index of synoptic motion in mid latitudes. Regions of positive relative vorticity are associated with cyclonic storms in the Northern Hemisphere. Thus the distribution of relative vorticity is an excellent diagnostic for weather analysis. Vorticity, the microscopic measure of rotation in a fluid, is a vector field defined as the curl of velocity. Relative vorticity is a measure of the intensity and direction of spin in a circular movement, which is performed by a unit volume of air around the vertical axis perpendicular to the plane over which this rotation occurs. Relative vorticity is a good quantity for studying atmospheric changes, because it presents the main order of magnitude of daily cyclonicity or anticyclonicity.
    Keywords: Relative Vorticity, Temperature, Canonical correlation analysis, Middle East, Iran
  • Fardin Saberi Louyeh *, Bohlol Alijabni, Shahriar Khaledi Pages 13-30
    Conclusion The present research, through applying precipitation and temperature extreme events, illustrates that percent of forecasted precipitation and temperature changes in comparison with the average base period of 1981-2010, in 2030, 2060, 2100 will increase procedurally. Spatial variability, and annual coefficient of variation in various regions are different. north, western north, eastern north and east will include the least temperature fluctuations, and the highest percent of precipitation with the highest coefficient of variation which conveys chronological period precipitation distribution with disordered accumulation and more local difference in this region in comparison with other regions. Then, Ghafghaz mountainous region has the highest percent of precipitation rise with suitable scattering in a year. The southern region of Caspian sea will experience the most rise of temperature and lowest percent of precipitation rise. High coefficient of variation in this area illustrates abnormal and disordered pattern on the threshold of precipitation. Sea level rise with three estimation regression, low average, high, on the basis of sea level ascending pattern equation For both scenarios, fluctuations in sea level based on subsidence Caspian pit seabed was calculated. In general, average annual sea level is increasing which is about 1.22 cm each year for scenario RCP8.5, and 0.93 cm yearly for scenario RCP4.5. Through this article, it can be found that changes in coastal region is unavoidable. However, inhabitants in this region have no system or not yet developed which can help them to adopt themselves with climate change issue. This study illustrated the significant effect of coastal climate which through climate change how society and economical activities are influenced.
    Keywords: Sea level rise, Caspian Basin, Emission scenarios, Extreme event, Coefficient of variation
  • Mohammad Moradi * Pages 31-42
    In this paper, data of the annual fastest wind speed in Bushehr station in south of Iran were used and graphical and numerical methods were applied to compute scale and local parameters of the Gumbel Distribution Function (GDF). Then, different return periods for the annual fastest wind speed were estimated. In the estimation process of local and scale parameters, Standard analytical procedures such as Method of Moments (MOM), Method of Order Statistics Approach (OSA), Least Squares Method (LSM) and Maximum Likelihood Method (MLM), were used.  
    Numerical computations show that the Method of Moments (MOM) provides better results compared to other methods and computed values for the scale and local parameters in estimation of annual fastest wind speed in Bushehr station are the best estimation.
    Computations of the annual fastest wind speed for return periods of 25, 50,100 and 1000 years, estimated to 29.7 m/s, 32.8 m/s, 35.8 m/s and 45.9 m/s, respectively. Moreover, we can say that, in the confidence level of 95%, every 207.2 and 82.9 years, annual fastest wind speed of 39 m/s and 35 m/s can happen, respectively.
    Keywords: least square, Maximum Likelihood, Order Statistics Approach, Method of Moments
  • Fahimeh Khadempour, Abbas Khashei *, Mahdi Amirabadizadeh Pages 43-56
    Introduction
    The solar energy is the source of energy on the planet, one of the most important sources of clean energy, especially in Iran. The correct estimation of solar radiation is considered one of the important parameters in simulating plant growth and the estimation of evapotranspiration is very important. Measuring the intensity of solar radiation, although it has a relatively long history in Iran, due to the high cost, most of the stations in the country are not equipped with a radiation meter, and this defect is even seen in developing countries. Since the prediction of solar radiation intensity in across Iran daily range with GeneXproTools software has not been done so far, due to the wide range of different climates in different regions of the country, this research has been carried out across Iran. Therefore, the ability of the simulation model of GeneXproTools software to predict solar radiation based on altitude, longitude, latitude, precipitation, relative humidity, maximum temperature, minimum temperature, average temperature and solar radiation on a daily basis in 31 stations over Iran. The period from 2007 to 2016 has been taken.
    Materials and Methods
    Iran with an average annual rainfall of 241 mm is located in dry and semi-arid latitudes of the planet, between the two meridians of 44⁰ and 64⁰ in the east, and two circuits 25⁰ and 40⁰ in the north. About 94.84% of its surface is located in arid and semi-arid regions with low atmospheric rainfall and high evapotranspiration. Meteorological stations were selected based on climatic variation in this research. In this way, the study stations were divided into six submergence classification systems, hyper-arid, arid, semi-arid, Mediterranean, humid and very humid (A). Gene expression planning (GEP) is one of the newest methods of artificial intelligence. This is the generalized method of genetic algorithm (GA), which was based on Darwin's theory and invented by Ferreira in 1999 (Roshgar and Mirhidaryan, 2014). In order to compare the results of the gene expression algorithm were defined in the prediction of solar radiation, 8 scenarios are considered based on the parameters affecting solar radiation such as altitude (m), longitude (degree), latitude (degree), precipitation (mm) relative humidity (%), maximum temperature (°C), minimum temperature (°C) and average temperature (°C).
    Results
    The main objective of this research is to select the best model for predicting daily average solar radiation in Iran using meteorological parameters. The number of training and testing data for the GEP model is presented in Table 1. Among the data collected, 80% of the data were used for training (total 8000 data, 6400 data from every parameter for every 31 stations) for the model. The experiment was performed for 20% of the data (a total of 1600 data from every parameter for each 31 stations) for the model.
     
    Table 1. Training and testing data for study stations for the GEP model

    Statistical period

     
    The performance of the GEP model by evaluating the best fitness fittings, namely RMSE, MAE, NSE and R2, which at best (Best Mode) of the values of best fitness RMSE, MAE, NSE and R2were 1000, 0, 0, 1, and 1, respectively. Comparing the results of model assessment statistics in different scenarios, it was found that in all scenarios, c scenario due to consideration of latitude (degree) parameters, precipitation (mm), average relative humidity (%), maximum temperature (°C), minimum temperature (°C), average temperature (°C), and average solar radiation (MJ.m-2.d-1), the desired prediction has presented. The best model for this study is R2 for training and data testing of 0.72 and 0.73, respectively. The GEP model shows the optimal output model as a tree, and also the equation derived from this structure. Since the four genes in this study are composed, each gene has a subtree and its own equation, which ultimately results from the final equation bond function. The final equation is a relation 1.
    Rs= SUB ET1+ SUB ET2 + SUB ET3 + SUB ET4 = sin((((((2RH)-cos(Tmax))-cos(Tmean))-p)/altitude)) + Tmean + sin(altitude) + ((((altitude + 7.30)-(2 Tmean))-p)/((sin(-8.41)+ 5.76)+atan(Tmin))))
    Conclusion
    From the results of the analysis, it was found that the predicted values ​​are well adapted to the values ​​measured in the Gene Expression Programming Model (GEP). In the GEP model, the e scenario is due to the fact that it does not consider altitude (degree), longitude (degree), latitude (degree), precipitation (mm), relative humidity (%), maximum temperature (°C) and minimum temperature (°C), with a lower correlation coefficient and more error and less efficiency. In this model, in the c scenario, with the addition of average temperature, and solar radiation parameters, a more favorable estimation of solar radiation is obtained. The results of the c and e scenarios are very close together, but the e scenario is weaker in the prediction of solar radiation, and the best model in this study is c scenario. More precise results will be obtained when latitude (degree), precipitation (mm), relative humidity (%), maximum temperature (°C), minimum temperature (°C), average temperature (°C) and solar radiation (MJ.m-2.d-1) are predicted for solar radiation. Also, the criteria for assessment under c scenario are higher with R2= 0.72, RMSE = 3.59, MAE = 2.82, and NSE = 0.72 in the education section of other scenarios. In general, the GEP method has the most accurate results in estimating daily solar radiation in Iran.
    Keywords: GeneXproTools software, Scenario, Solar Radiation
  • Javad Salarvand, Zahra Asodolahi, Farhad Ghasemi * Pages 57-72
    Introduction
    Rainfall erosivity is defined as the aggregative power of the rain. If other effective features on soil erosion be considered constant then soil loss could be directly connected to rainfall erosivity. Rain erosion term was proposed by Wichmeier and Smith in 1978 to consider the effect of climate on raw erosion. Measurements of meteorological parameters by the traditional methods require a dense rain gauge network. But, due to the topography and cost problems, it is not possible to create such a network in practice. Given the significant change in rainfall in time and space on the one hand and low rain-gauge stations to record rainfall on the other hand, the necessity of using geostatistical methods for rainfall erosivity mapping is inevitable. Geostatistical methods use the spatial correlation between observations in the estimation processes. In these cases the spatial distribution pattern of rainfall erosivity can be produced using different methods of interpolation.
    Materials and Methods
    Study area Lorestan province is located in southwest of Iran and covers an area of 28249 square kilometers. It is located between the latitudes 32º 37' and 34º 22' N and the longitudes 46˚51ʹ and 50˚30ʹE. The main objective of this research were: (1) analyze the spatial distribution of rainfall erosivity using two different interpolation methods namely ordinary Kriging and simple Kriging; (2) put forward the best interpolation method through cross-validation, construct the high resolution grid data of rainfall erosivity and provide the reliable information for relevant researches
    Monthly rainfall erosivity model

     In the first step, the precipitation data collected from 53 precipitation stations and Modified Fournier Index (MF) calculated based on Eq. (1)
                                                                                                    (1)
     Where MF is the modified Fournier index value (mm), pi  is averagemonthly precipitation (mm) and P is average annual precipitation (mm). Then Eq. (2) and Eq. (3) were used to estimate rainfall erosivity or R-factor values (MJ mm ha -1 h -1 year -1).
     
    It is suggested that Eq. (2) be applied for locations with MF values less than 55 mm and Eq. (3) be used for locations with MF values greater than 55 mm.
     
    Geostatistical Methods
    In this article two interpolation techniques namely simple and ordinary Kriging were compared in GS+5.1.1 and ArcGIS10.3 software’s in order to determine which one describe better the spatial distribution of rainfall erosivity. Kriging methods assume that the spatial variation of a continuous climatic variable is too irregular to be modeled by a continuous mathematical function, and its spatial variation could be better predicted by a probabilistic surface. The predictions of Kriging-based methods are currently a weighted average of the data available at neighboring sampling points (weather stations). The weighting is chosen so that the calculation is not biased and the variance is minimal. A function that relates the spatial variance of the variable is determined using a semi-variogram model which indicates the semi variance between the climatic values at different spatial distances.
     Validation and techniques comparison
    The resulting maps from interpolation were compared by using a set of validation statistics include Mean error (ME), Mean Standardized Error (MSE) and the root mean square error (RMSE) by Eq. (4), (5) and (6)
     
    Results and Discussion
    Based on the results, rainfall erosivity values varied from 11.1 to 749.5 MJmmha−1 h−1 y−1. Differences between the simple and ordinary Kriging models regarding the validation statistics were narrow, but allowed for a comparison. The obtained results showed that ordinary Kriging with higher R2 and lower ME, MSE and RMSE had better precision in mapping rainfall erosivity. The spatial distribution of rainfall erosivity showed the areas along north-south of Lorestan province and central regions had higher values while lower rainfall erosivity was seen in the western and eastern areas of the study area.
     Conclusion
    The availability of high-quality environmental maps is a key issue for agricultural and hydrological management in many regions of the world.  Produced rainfall erosivity map in this research can be used for estimation of soil loss by USLE model. Rainfall erosivity maps also can be suitable as guidance for soil conservation practices and identifying areas with the high potential of soil retention as an ecosystem service. Further research may be directed to find reliable erosivity indices which can be computed from daily precipitation data.
    Keywords: Erosion, Kriging, Variogram, Interpolation
  • Ahmad Seifi, Narges Salehnia, Fatemeh Golzadeh Khargh * Pages 73-89
    The traditional approachof economics is based largely on the wasteful use of natural resources and the lack of attention to future generations' rights to these resources. However, the environmental impacts of such an approach will be irreparable and in order to achieve the sustainable development, it has to be changed. Considering the importance of protecting the environment as well as the need to pay more attention to the sustained growth and development, this study, examines the factors affecting air pollution from the economic and institutional perspectives with a greater focus on the two variables: trade openness and the control of corruption index. Researches in the field of environment confirm that there is an inverted U-form relationship between the economic growth and the air pollution. In fact, the Environmental Kuznets Curve (EKC) claims that in the first stage of development, pollution grows rapidly; because the rapid growth results in untapped use of natural resources, and therefore the air pollutants generation increases. But at the high levels of development, people value the healthy environment, more effective environmental regulations will be implemented, regulatory institutions become more effective, the economic structure moves toward clean industries, and finally the pollution level declines. Regarding the role of trade patterns in the transmission of pollution, from high income economies which have stringent environmental standards to the middle and low income economies that place high levels of production and employment as their top priorities, new models have been proposed that can cover the role of trade in the EKC. Two conflicting hypothesis emerge from the debate: the pollution haven hypothesis (PHH) and the factor endowment hypothesis (FEH). According to the pollution haven hypothesis, pollution control policy differences across countries drive pollution-intensive industries to the countries with weaker regulations. Another hypothesis, the factor endowment hypothesis (FEH), predicts that rich and developed countries will specialize in polluting goods. Since the abundance of capital supports the production of more capital goods and polluting industries. On the other hand, countries with a lot of workforce and land tend to be more specialized in less polluting sectors. Studies show that with the entry of commercial variables into the discussion, increasing or decreasing of pollution is different depending on the country and the environmental index which is used. It seems that only the entry of economic variables, is not a solution to the problem, because a growing attention of countries will be given to the condition of pollutants and the need to move toward a sustainable development. As a result, they will face the challenge of complying with environmental regulations. Concurrent importance of economic growth, and reducing poverty and unemployment has made some managers to ignore the neglecting of rules by companies. So they can also obtain benefits for themselves through the rents. In this regard, some recent studies have come up with models that apply the effects of institutional factors. Most of these researches have shown that improving the governance, and in particular controlling the level of financial and administrative corruption, play a crucial role in reducing emissions. However, it seems that the quality of the regulatory institutions that can ensure the correct implementation of environmental regulations is less considered in the existing literature. So, using the Fully Modified Least Squares (FMOLS) approach in selected countries of three income groups from 1996 to 2011, the purpose of this study is to investigate the effects of trade openness and control of corruption on the particulate matter index. Empirical results indicate that the effect of trade openness on air pollution in all three income groups is negative and significant. It means that increasing in trade openness will improve the air quality and reduce PM10 emissions. Furthermore, the empirical findings show that with an increase in control of corruption, there is a significant decrease in PM10 emissions in Middle and Low income groups; while this effect is not statistically significant in high income group. Therefore, reducing the level of corruption and combating of corrupt activities will enhance the air quality in Middle and Low income groups. Since the clean air and healthy environment play an important role in the survival of human beings, the results of this study can lead to policies for the sustainable development.
    Keywords: Control of Corruption, Fully Modified Ordinary Least Squares (FMOLS), Particulate matter (PM10), Trade openness
  • Fakhreddin Iranpour *, Hassan Zohreh Wandi Pages 91-105
    Introduction
    In Iran, precipitation is one of the key variables for assessing the potential of water resources access. But because the spatial distribution of this variable is very uneven, the distribution of water resources of the country is not uniform. The maintenance and management of water resources, while also being a function of rainfall, depends on the variability of precipitation. The smaller the spatial variation of the rainfall, the greater the homogeneity and consistency of water resources. On the other hand, the less variability of rainfall is, the more water resources will be more stable and water supply will be possible. For this reason, the variability of rainfall time in the assessment of water resources, watersheds and the relative study of local and regional water resources is important.   
     materials and methods
      The data used in this study are available statistics for 60 years of daily precipitation of the synoptic station of Tabriz from 1951 to 2010, which was obtained from the Meteorological Office of East Azarbaijan Province. For accuracy raising the modeling stage, data is considered weekly. The most commonly used model used to show the time series of discrete random variables is known as the Markov chain. Examining these uncertain or random modes and selecting the model is the probability knowledge. In this research, it is attempted to use this knowledge and based on the Markov chain method, the probability of rainfall occurrence in Tabriz city is obtained.                                                                                              
    Results and discussion
    The distribution of rainfall during the year can have a large impact on water and agriculture planning. Considering the fact that Iran is located in the dry world belt and most of its regions have dry and semi-arid climates, and agriculture in these areas It is also based on this climate. Changes in rainfall can cause irreparable damage to the agriculture and water resources of these areas. Therefore, recognizing its system of changes can be a great help in this regard. The number of rainfall days is an appropriate criterion for assessing the distribution of rainfall. The average rainfall days in Tabriz average 80 days per year. An average of 30 days was observed in the spring, 5 days in the summer, 19 days in the fall and 26 days in winter.   Conclusion The study of state change matrix shows that during the statistical period (1951-2010), from the total of 21960 days of the studied statistics, 14649 days of change from dry day to the next dry day and 2519 days of change in the rainy day, which after dry day It happened. Also, the transformation matrix from rainy day to dry day is 2520 days and rainy day to rainy is 2272 days. The results of the X2 test indicate that the data are not independent and there is enough confidence to adhere to the daily precipitation data of Tabriz station from the Markov ranking model. The results of the trend test using Spearman correlation method at 95% confidence level indicate that the data are not trendy and validated by the approved chain. The results of the fitted final model showed that the probability of changing the order of the rainy day to the other rainy days has a higher percentage than the change from rainy day to dry, and this state is in the weeks leading to the spring and autumn seasons, More than other weeks.
    Keywords: Consecutive dry days, Markov chain, precipitation probabilities, Tabriz