فهرست مطالب

پژوهش آب ایران - پیاپی 29 (تابستان 1397)

مجله پژوهش آب ایران
پیاپی 29 (تابستان 1397)

  • تاریخ انتشار: 1397/06/31
  • تعداد عناوین: 12
|
  • حسین صمدی بروجنی، روح الله فتاحی نافچی، داریوش ملکی، پری ملکی * صفحات 1-8
    رودخانه ها، در طی زمان برای رسیدن به پایداری، شکل های بستر متفاوتی را با رسوب گذاری یا فرسایش در بستر خود ایجاد می کنند. از راه های شناسایی رفتار رودخانه ها، مطالعه و بررسی ساختار جریان و شیوه ی شکل گیری شکل های بستر در آن ها است. در رودخانه های آبرفتی بخشی از تنش برشی مربوط به شکل بستر است. در این پژوهش به بررسی تنش برشی بستر، روی دو نوع شکل بسترشکنج (موازی و پولکی) پرداخته شده است. بدین منظور آزمایش هایی در چهار شیب و چهار دبی روی شکنج های مصنوعی ساخته شده در آزمایشگاه هیدرولیک دانشگاه شهرکرد انجام شد. برای اندازه گیری تنش و سرعت جریان از دستگاه سرعت سنج سه بعدی ADV استفاده شد. نتایج نشان داد میزان تنش در تمام آزمایش ها، روند ثابتی داشته است؛ به گونه ای که روی وجه بالادست شکنج، ماکزیمم مقدار را دارد و با نزدیک شدن به تاج از مقدار آن کاسته می شود و در تاج، کمترین مقدار را دارد. همچنین نتایج نشان داد شکنج های موازی و پولکی به ترتیب 26 و 23 درصد از میزان کل تنش برشی را تشکیل می دهند. میزان تنش برشی بستر، در روی شکنج های موازی، 27 درصد بیشتر از مقدار تنش برشی بستر روی شکنج های پولکی است.
    کلیدواژگان: تنش برشی بستر، رودخانه های آبرفتی، شکل بستر، شکنج
  • سعید سلطانی کوپایی، منیره فرامرزی، سعید پورمنافی، سارا آزادی * صفحات 19-28
    برای کمی سازی خشک سالی، شاخص های مختلفی توسعه داده شده اند که هر کدام دارای نقاط قوت و ضعفی هستند. ضریب های برخی از معادلات این شاخص ها، مربوط به شرایط مناطق استخراج آن ها است که برای دستیابی به نتایج قابل اعتماد در مناطق مختلف دیگر، لازم است مورد واسنجی قرار گیرند. هدف این پژوهش واسنجی ضرایب شاخص پالمر برای ایران است. به این منظور از داده های دما و بارش در مقیاس ماهانه در دوره آماری 2002-1990 و رطوبت ابتدایی خاک مربوط به 506 حوضه مطالعاتی کشور استفاده شد و میانگین وزنی اطلاعات برای هفده حوضه کشور به دست آمد. سپس شاخص پالمر برای کل ایران مورد واسنجی قرار گرفت و به صورت جداگانه برای شرایط خشک و تر به ترتیب ضرایب دوگانه (892/0، 002/0) و (946/0، 002/0) در مقایسه با معادله پالمر با ضرایب (897/0، 333/0) که تنها برای شرایط خشک بود، به دست آمد. نتایج این پژوهش انطباق بسیار خوبی را در اکثر حوضه ها برای سال های 2000 و 2001 به عنوان خشک ترین سال ها و سال های 1992 و 1993 به عنوان مرطوب ترین سال های مورد بررسی با توجه به نتایج حاصل از تحلیل شرایط بارشی سال های مورد نظر نشان داد. نتایج این پژوهش زمینه ای را برای واسنجی شاخص پالمر و ارزیابی خشک سالی در مناطق اقلیمی مختلف ایران فراهم می کند.
    کلیدواژگان: ایران، حوضه هیدرولوژیکی، خشک سالی، شاخص شدت خشک سالی پالمر، واسنجی
  • شبنم وکیلی * صفحات 41-49
    هدف اصلی از این مطالعه، بررسی عددی پروفیل سطح آزاد جریان روی سرریز اوجی و پلکانی با روش های حجم محدود و المان محدود و مقایسه این دو روش است. بدین منظور کد Fluent به عنوان نماینده ای از روش حجم محدود با کد ADINA به عنوان نماینده ای از روش المان محدود با نتایج آزمایشگاهی مقایسه شده است. شایان ذکر است برای مدل سازی آشفتگی از مدل بهره گرفته شده است. با بررسی نتایج، هماهنگی خوبی بین این دو روش و نتایج آزمایشگاهی مشاهده شده، اما با توجه به نتایج آزمایشگاهی موجود، در برخی نقاط روش حجم محدود در مقایسه با روش المان محدود، نتیجه قابل قبولی را ارائه کرده است. همچنین حداکثر خطا در پروفیل سطح آب در مقایسه با نتایج آزمایشگاهی به اندازه 2 میلی متر در نزدیکی دیواره سرریز مشاهده شده است.
    کلیدواژگان: سریز اوجی و پلکانی، Fluent و ADINA، المان محدود، احجام محدود
  • فاطمه یوسفوند، محمدرضا کاویانپور، محمد جواد منعم * صفحات 51-58
  • فرزین پرچمی عراقی *، فاطمه سمیع پور صفحات 59-69
    در پژوهش حاضر، امکان تولید نقشه های ابرآزاد روزانه و هشت روزه سطح برف (SC) و عمق آب معادل برف (SWE) با تفکیک مکانی 500 متر از طریق یک الگوریتم ادغام-ریزمقیاس سازی که مبتنی بر ترکیب مزیت تفکیک مکانی مطلوب تصاویر حاصل از سنجنده های نوری MODIS نصب شده بر روی ماهواره های Terra و Aqua (به ترتیب، MODIS Terra و MODIS Aqua) با مزیت پشت نمایی پوشش ابر سنجنده میکروویو AMSR E می باشد، مورد بررسی قرار گرفت. این الگوریتم برای حوضه کارون، سال آبی 1383-1382، به کار گرفته شد. ادغام نقشه های SC حاصل از سنجنده های MODIS Terra و MODIS Aqua (TAC) موجب شد میانگین سهم پیکسل های ابری در نقشه TAC طی دوره مطالعاتی در مقایسه با نقشه های MODIS Terra و MODIS Aqua به ترتیب، به میزان 29/17 و 43/14 درصد کاهش یابد. متوسط دقت کلی نقشه های SC در شرایط آسمان صاف (مقایسه به ازای پیکسل های غیرابری) برای تصاویر حاصل از سنجنده های MODIS Terra و MODIS Aqua به ترتیب، برابر با 20/97 و 86/95 درصد بود. متوسط دقت کلی نقشه های SC در شرایط جوی عام (مقایسه به ازای پیکسل های ابری و غیرابری) ، شامل نقشه های MODIS Terra، MODIS Aqua، TAC و MAC SC (حاصل از ترکیب نقشه TAC و AMSR E) به ترتیب، برابر با 38/49، 47/51، 48/65 و 15/90 درصد بود.
    کلیدواژگان: سنجش از دور، ماتریس اغتشاش، AMSR E، Aqua، MODIS، Terra
  • سیدناصرالدین بدیسار *، محمدصادق احمدی صفحات 71-80
    کمبود منابع آبی و خشک سالی های پی در پی، از مهم ترین مشکلات ناشی از تغییرات اقلیمی و فعالیت های دولت ها است که بدون شک در آینده موجب اختلافات بسیار زیادی خواهند شد. بحران آب در سطح جهانی و تلاش دولت ها برای مبارزه با خشک سالی و قحطی و آثار ناشی از آن، دولت ها را بر آن داشته است تا نهایت تلاش خود را برای کنترل و استفاده از منابع آب های زیرزمینی به کار گیرند. آب های زیرزمینی به عنوان یکی از اصلی ترین و مهم ترین منابع آبی در سطح کشور، تامین کننده عمده آب شرب، کشاورزی، صنعت و غیره است و نقش مهمی در شکوفایی اقتصاد کشور داشته و دارد. سیر تاریخی تغییرات قوانین آب از سال 1307 تاکنون، نه تنها متناسب با شرایط روز نبوده، بلکه موجب افزایش تخلفات و کاهش قدرت و اقتدار حاکمیت بر آب نیز شده است. روندی که با قانون توزیع عادلانه آب و قوانین بعدی شروع و تقویت شده است، موجب گردیده که علی رغم تعیین تکلیف تعداد زیادی از چاه های غیر مجاز طی سال های 1373 تا 1384، هنوز در شرایط کنونی بیش از 250 هزار حلقه چاه غیرمجاز موجود بوده و انگیزه برای تداوم این اقدام و افزایش تخلفات، وجود داشته باشد. علی رغم وجود ضعف هایی در قانون گذاری، از ظرفیت و اختیارات قوانین موجود نیز به خوبی استفاده و قسمت قابل توجهی از اهداف و مقاصد قانون گذار تامین نشده است. این مقاله با هدف بررسی مقررات حاکم منابع آب های زیرزمینی (چاه ها) و به روش توصیفی- تحلیلی به رشته تحریر در آمده است.
    کلیدواژگان: پروانه بهره برداری، چاه ها، مقررات، منابع آب زیرزمینی، مالکیت، مناطق ممنوعه
  • سوده قبادی خسرو، محمدرضا گودررزی، محمدرضا هاشمی، مهدی کماسی * صفحات 81-90
    خشک سالی، به عنوان یک حادثه ناگوار طبیعی، به طور مستقیم جوامع را از طریق تغییرات در دسترسی به منابع آب تحت تاثیر قرار می دهد. برای درجه بندی شدت خشک سالی، شاخص های متعددی وجود دارد که از میان آن ها، شاخص خشک سالی موثر (EDI) و شاخص بارش استاندارد (SPI) کاربرد فراوانی دارند. در این پژوهش، برای پیش بینی خشک سالی از داده های بارش دو حوضه واقع در استان های همدان و لرستان برای محاسبه شاخص SPIو شاخص EDI بهره گرفته و برای پیش بینی پدیده خشک سالی از دو مدل برنامه ریزی ژنتیکی و مدل درختی M5 استفاده شد. بررسی های انجام گرفته نشان داد که این مدل ها توانایی خوبی در پیش بینی پدیده خشک سالی داشتند و از دقت مناسبی برای مسائل پیش بینی برخوردار هستند. از دیگر مزیت این مدل ها، ارائه فرمول های ساده و صریح برای پیش بینی پدیده موردنظر است. ضریب تبیین در مدل M5 برای شاخص EDI در بهترین حالت 97/0 و در مدل برنامه ریزی ژنتیکی 95/0به دست آمد. همچنین ضریب تبیین برای شاخص SPI در بهترین حالت در مدل M5، 93/0 و در مدل برنامه ریزی ژنتیکی 83/0 حاصل شد. این مهم نشان دهنده این نکته است که مدل درختی M5 در مقایسه با مدل برنامه ریزی ژنتیکی از دقت بالاتری برخوردار بوده و به دلیل سادگی و قابل فهم بودن نسبت به مدل برنامه ریزی ژنتیکی دارای برتری نسبی است.
    کلیدواژگان: دشت سیلاخور، شاخص SPI، شاخص EDI، مدل سازی، داده کاوی
  • ابراهیم ابراهیمی، محمد ابراهیم بنی حبیب، عباس روزبهانی * صفحات 91-100
    پیش بینی سطح آب زیرزمینی، نقش مهمی در مدیریت منابع آب زیرزمینی دارد. در این پژوهش، با استفاده از شبکه های بیزین پویای احتمالاتی و در نظر گرفتن 10 رویکرد پیش بینی ماهانه، با استفاده از داده های واسنجی 12 سال آبی، سطح ایستابی آبخوان بیرجند، در یک دوره آبی 5 ساله و در دو رویکرد صریح (پیوسته) و خوشه بندی (دسته بندی) صحت سنجی و پیش بینی شد. در رویکرد پیش بینی صریح، از داده های صریح ماهانه استفاده و در رویکرد خوشه بندی با استفاده از شاخص عرض سیلهوت تعداد خوشه های بهینه تعیین و با خوشه بندی داده ها، از داده های خوشه بندی شده برای آموزش مدل بیزین استفاده و پیش بینی مدل با استفاده از داده های خوشه بندی انجام شد. با اعمال واسنجی در مدل بیزین و بررسی نتایج سناریوهای تحلیل حساسیت مختلف، کارایی بالای شبکه های بیزین پویا به ویژه در رویکرد صریح، مشاهده شد. با بررسی 10 رویکرد تحلیل حساسیت اعمال شده برای پیش بینی سطح ایستابی آبخوان، مشاهده شد که از بین عوامل پیش بینی کننده، پارامتر سطح آب زیرزمینی در ماه فعلی اثرگذاری بالایی را در افزایش دقت مدل برای پیش بینی سطح ایستابی آبخوان ماه بعدی دارد؛ به گونه ای که با حذف این پارامتر در مدل بیزین، ضریب تبیین از 9925/0 به 0004/0 کاهش یافت و میانگین مجذور مربعات خطا نیز از 2654/0 به 9709/0متر افزایش یافت. همچنین بالاترین ضریب تبیین در بین سناریوهای پیش بینی 995/0 است.
    کلیدواژگان: تحلیل حساسیت، تئوری بیز، خوشه بندی، مدیریت آب زیرزمینی
  • حسین ملکی نژاد، زهرا خسروانی، امیرحسین پارسامهر * صفحات 101-110
  • فریدون رادمنش، غلامعلی پرهام، رسول میرعباسی نجف آبادی، فرشاد احمدی * صفحات 111-120
    در مطالعه حاضر، تحلیل فراوانی سیلاب حوضه آبریز دز در دوره آماری 91-1335 با استفاده از رویکردهای قانون توانی (PL) و توزیع های متداول آماری انجام شد. ابتدا سری های جزیی سیلاب با استفاده از قانون توانی با فواصل زمانی 7، 14، 30 و 60 روزه استخراج و با سری حداکثر سالانه مقایسه شد. نتایج نشان داد که سری حداکثر سالانه، برای تحلیل فراوانی سیلاب حوضه آبریز دز مناسب نیست و سری جزیی 30 روزه حاصل از قانون توانی، مطابقت بهتری با وضعیت جریان این حوضه دارد. استقلال و ایستایی سری جزیی سیلاب 30 روزه با استفاده از آزمون والد- ولفوویتز (W-W) بررسی و استقلال سری موردنظر تایید شد. در مرحله بعد، توزیع توانی و توزیع های متداول آماری بر داده های سیلاب حوضه آبریز دز برازش داده و عملکرد هر توزیع با استفاده از آماره های جذر میانگین مربعات نرمال شده (NRMSE) و معیار نش- ساتکلیف (NS) بررسی شد. مقدار آماره های NS و NRMSE برای توزیع توانی در ایستگاه تنگ پنج بختیاری به ترتیب برابر با 99/0 و 75/2 محاسبه شد که نسبت به سایر توزیع ها از عملکرد بسیار بهتری برخوردار بوده و در ایستگاه های تنگ پنج سزار (96/0=NS و 12/5=NRMSE) و تله زنگ (98/0=NS و 40/3=NRMSE) نیز با اختلاف بسیار اندکی در رتبه دوم بهترین توزیع قرار گرفت. با توجه به اینکه عملکرد توزیع توانی در برآورد سیلاب حوضه آبریز دز بسیار مناسب و محاسبه ضرایب و کاربرد آن نسبت به توزیع های احتمالاتی آسان است، می تواند به عنوان توزیع برتر حوضه آبریز دز پیشنهاد شود.
    کلیدواژگان: ایستایی، سیل، سری جزیی، فرکتال
  • ابراهیم امیدوار، فرزانه فدایی، رضا قضاوی * صفحات 121-129
|
  • pari maleki *, hossein samadi boroujeni, Rohollah Fatahi nafchi, Darioush maleki Pages 1-8
    Interaction between water flow characteristics and bed erodibility plays an important role in sediment transport process. In order to reach stability, rivers with deposition or bottom erosion make different forms in the bed .The bed forms create extra resistance, which is called the bed resistance. The mutual interaction between the flow and the erodible bed through sediment transport phenomena in a sand-bed channel causes a variety of bed forms. Starting with ripples and gradually increasing in shear stress or water velocity, dunes, washed out dunes, flat bed, anti-dunes, and standing waves are formed. The most common boundary conditions in alluvial rivers are the mobile beds covered with Ripple and Dune. These forms, in many alluvial systems, play a critical role in contrast between the flow, discharge of sediment and morphology of bed. One way to identify the behavior of the rivers is to study the structure and the formation of bed forms within them. Ripples are among the smallest of the bed forms. The longitudinal cross-sections of ripples are usually asymmetrical. The upstream face of ripple is long and has a gentle slope, and the downstream face is short and steep. The height of ripples is usually between 0.5 cm and 2 cm, but not more than 5 cm. The wave lengths normally do not exceed to 30cm, and they are usually within the range of 1 cm to 15 cm. Some ripples that form in deep-water regions are symmetrical. Ripples are the smallest of the bed configurations. They are related to physical parameters near the river bed and have little correlation with the water depth. Their occurrence is the result of the unstable viscous layer near the boundary. They can form in both shallow and deep water. In plan, they either are parallel to each other or have a shape like fish scales. With increasing the flow velocity, the plan form of the ripples gradually develops form straight line to curves and then to a pattern like fish scales, symmetrical or unsymmetrical Resistance is a function of the geometrical dimensions of the bed forms and depth of water. Estimating of the flow resistance is one of the most important matters in planning, designing and operating of water resources projects, including water transfer and river system management. In this research, the effects of two different types of ripples (parallel and flake shape) on the hydraulic characteristics of flow were experimentally studied. The experiments flume located at the hydraulic laboratory of Shahrekord University, Iran. The flume used in this study was a Straight type that had the dimensions of 0.4 m wide and depth and 12 m long. This flume has vertical PVC sidewalls. Generally 48 tests in variety slopes of 0.0005 to 0.003 and variety discharges of 10 to 40 lit/s were conducted. Velocity and the shear stress were measured by using an Acoustic Doppler Velocimeter (ADV). Velocity measurements were performed with a frequency of 200 Hz, which provided accurate statistics on the mean flow and turbulence characteristics. Detailed velocity measurements were performed in 9 cross-sections in the Straight flume. The cross-sections were determined at: before crest, crest and after crest. Generally, in each section velocity was measured in 12 point (5, 20 and 35 cm from flume side and 4 points from bed). For All tests, flow depth was kept constant. In this study generally 48 test were tested under different hydraulic conditions. It was observed that the peak value of the bed shear stress appeared on the midpoint of upstream surface (before crest ripple) and the crest of the ripples had the lowest value of the bed shear stress. From the crest to trough (after crest) of the ripples, the general bed shear stress was in an increasing trend. It was generally found that with increasing Froude number and the bed shear stress increased. Also, in the case of parallel ripple bed form, the shear stress was about 26% more than that of plane bed and in the case of flake ripple bed form the increasing rate was 23%. It indicated that the shear stress was much affected by the parallel shape of ripple bed form compared with the flake shape, as it was 27% more than the flake shape.
    Keywords: Alluvial river , Bed form , Bed shear stress , Ripple
  • Sara Azadi *, Saeid Soltani kopaei, Monireh Faramarzi, Saeid Pourmanafi Pages 19-28
  • shabnam vakili * Pages 41-49
    The spillway is one of the most important elements of dam building. So studying and designing of the spillway hydraulic should be considered in dams and other similar projects. The ogee spillway, because of its superb hydraulic characteristics, has been widely studied. Its ability to pass flow efficiently and safely, when has been properly designed, with relatively good flow measuring capabilities, has enabled engineers to use it in a wide variety of situations. In the last two decades, there has been an increasing interest in the stepped spillways in various laboratory experiments around the world. This is partly because of technical advances in the construction of Roller Compacted Concrete (RCC) dams and considerable amount of energy dissipation along the chute, leading to reduction in the size of the stilling basin. Flow analysis using laboratory experiments usually involves considerable time and cost requirements .Nowadays, availability of computational fluid dynamic (CFD) programs and powerful computers has resulted in increasing usage of numerical methods of flow analysis. Free surface flows are encountered in hydraulic engineering problems including water jets, weirs and around gates. For the numerical modelling of a steady flow, the area of calculation and therefore the position of the free surface have to be known. The main object of this study is numerical analysis of free surface profiles on ogee and stepped spillways by finite volume and finite element methods and comparison of the two methods. The physical laws governing a fluid flow problem are represented by a system of partial differential equations regrouping the continuity equation, the Navier-Stokes equations and any additional conservation equations. The numerical analysis resolves these equations by accurate and complex numerical schemes. A program or code, where the numerical algorithm is implemented, is then solved on a computer. In recent study experimental results offered by Tabara and CHatila (2005) have been used to investigate stepped spillway physics models of Tabara and Chatara (2004) for studying ogee spillway. FLUENT and ADINA software have been used to simulate flow field of types of stepped and ogee spillways. It should be noted that k-£ model has been applied in order to modeling turbulence. Gambit preprocessor software is a tool for networking flow field. In recent study we have used PRESTO scheme for discretization of pressure, quick plan for discretization of terms of displacement of momentum equations, turbulence formulas and PISO algorithm for coupling velocity and pressure. In finite element method, Galerkin relation has the second order accuracy. SUPG method is used in momentum equations and turbulence for discretization of term advection which is effective in problem convergence. Meshing type has been selected as triangle form in ADINA code, but quadrilateral in FLUENT code regarding the spillway geometry. In some area structured mesh and in some other ones non-structured mesh have been chosen. Due to existence of rotatory flows over the steps, more fine meshes have been used in both codes near to body of stepped spillway as well as ogee spillway, because of the velocity gradient nearby wall, where flow enters spillway channel from reservoir. In recent study, water surface profile has been simulated in four types of stepped spillways and 3 types of ogee spillways on the basis of finite element method using FLUENT code and finite volume method using ADINA code. Relatively close consistency has been observed between water surface profiles in both codes by comparing experimental results. Note that free surface profile simulated over the ogee spillway was closer to water surface profile evaluated in laboratory than the profile by finite element method. There were observed better consistency in results from free surface profile over stepped spillway with finite volume method. Maximum differences between results of water surface profile in numerical model of FLUENT and ADINA were nearly 1.2 cm and 1.6 cm, comparing with laboratory results respectively. This value in numerical simulation is acceptable regarding to different meshing in model networking. On this basis, numerical methods of FLUENT and ADINA are proper in order to simulate similar structures for saving costs due to construction of physical model. Flow surface can be determined considering the lower error of FLUENT model in measuring free surface flow; thereby, designing lateral walls of spillway is realized with more precise.
    Keywords: Ogee spillway, Finite element, ADINA, FLUENT, Stepped spillway, Finite volume
  • Dr M. Monem *, Fatemeh Yousofvand, Mohammad Reza Kavianpour Pages 51-58
  • Farzin Parchami, Araghi *, Fatemeh Samipour Pages 59-69
    Snow cover, as an important component of land cover, is one of the most active natural and plays an important role in hydrological processes and climate. Variability in snow covered area has a significant influence on water and energy cycles, as well as socioeconomic and environmental repercussions. Frequent and long-term snow observation, accurate snow cover (SC) mapping and snow water equivalent (SWE) estimation are crucial for operational flood control, water delay planning, and resource management in snowmelt-dominated basins. Today, satellite-derived snow products obtained from visible and infrared imagery, as well as passive microwaves are available on the Internet, with few recently availability in near-real time . Optical sensors (e.g., Landsat, Advanced Very High Resolution Radiometer-AVHRR, Moderate Resolution Imaging Spectroradiometer-MODIS, Systeme Probatoire d’Observation de la Tarre-SPOT) have been well developed to provide snow information with good temporal and/or spatial resolution. But, cloud cover is a major factor in optical remote sensing that limits our capability to map the Earth’s surface. It is often not an easy task to collect a time-series cloud-free images for a particular area of interest, using optical remote sensing, which limits their wider applications for SC monitoring. Space-borne passive-microwave radiometers (e.g., Scanning Multichannel Microwave Radiometer-SMMR, Special Sensor Microwave/Imager-SSM/I, Advanced Microwave Scanning Radiometer-Earth Observing System-AMSR-E), can penetrate cloud to detect microwave energy emitted by snow and ice and provide information on SC and SWE. These passive microwave data are well suited to snow cover monitoring because of the characteristics such as all-weather imaging, large swath width with frequent overpass times. But, the coarse spatial resolution (e.g., 25 km of AMSR-E daily SWE product) hinders their applications in operational hydrological modeling. It seems that the combination of MODIS and AMSR-E can take advantage of both high spatial resolution of optical data and cloud transparency of passive microwave data. In this study, daily cloud-free SC and SWE maps at the 500-m resolution were produced for Karun watershed, Ahwaz, Iran (January 25–32, 2004). The daily MODIS-Terra, MODIS-Aqua, and AMSR-E snow data products were used via a fusion-disaggregation algorithm. The developed SC and SWE maps were evaluated, using total accuracy of snow mapping in clear-sky (Os) and all-sky (Oa) conditions, underestimation (UEc) and Overestimation (OEc) of snow covered area in clear-sky condition, snow accuracy in clear-sky (Sc) or in all-sky (Sa) conditions, and no snow accuracy in clear-sky (NSc) or in all-sky (NSa) conditions. The results of this study showed that the combination of MODIS-Terra and MODIS-Aqua considerably reduced the cloud coverage in such high resolution optical data. Although MAC-SC and MAC-SWE products have been developed to have 500 m spatial resolution, the massive and continuous cloud cover (larger than 25 km in size) in the MODIS-Terra and MODIS-Aqua and, hence, TAC products were simply replaced by the coarse AMSR-E pixels. In this case, although those cloud coverages were removed in the MAC-SC and MAC-SWE products, the actual resolution of the snow or no-snow pixels kept 25 km and such pixels had the false spatial resolution of 500 m. The SWE redistribution of AMSR-E based on MAC-SC products enhanced (to some extent) the spatial resolution of the AMSR-E SWE products. However, there was no measured data to evaluate the accuracy of the enhanced SWE products. It can be concluded that for pixels with scattered cloud cover (less than 25 km in size) in the TAC products, the MAC-SC and MAC-SWE products indeed improve the spatial resolution of those pixels to 500 m, while for massive cloud cover (larger than 25 km in size), the actual resolution of those pixels in the MAC-SC and MAC-SWE products are 25 km, even in 500 m pixel size. Despite of these limitations, the MAC-SC and MAC-SWE maps are suitable for hydrological, meteorological modeling on a daily basis in the study area.
    Keywords: AMSR E _Aqua_Confusion matrix _MODIS _Remote sensing_Terra
  • seyed naseraldin badisar *, mohammad sadegh ahmadi Pages 71-80
    Shortage of water resources and successive droughts are one of the most important types of problems caused by climate changes and governmental activities that will undoubtedly cause great disputes in the future. Global water crisis, drought and famine and their effects have made governments to make their utmost effort to control and use groundwater resources. Groundwater, as one of the main and most important water resources in Iran, is a major supplier of drinking, agricultural, and industrial water and has an important role in economic growth. Historical changes of water laws from 1928 so far, not only have not been in conformity with modern circumstances, but also have increased violations and reduced power and authority of sovereignty over the water. The process began and strengthened with the law of equitable distribution of water and subsequent laws. The laws has led to a situation that despite the legal decisions, a large number of illegal wells was made during the period of 1994 to 2005. There are still more than 250 thousand illegal wells in the current situation (Mazaheri and Abdulmanafi, 2015) and there is an incentive to continue these activities and increasing violations. In spite of the weaknesses in the legislation, the capacity and powers of the existing laws are not used properly and a remarkable part of the intentions of the legislator is not met. In the law of equitable distribution of water, no new decree has been issued regarding to the private ownership of water resources, except groundwater. Whether beneath a private land or public land, water is absolutely a joint asset that people can own by obtaining permission from the government. In fact, what has changed with the adoption of new rules is the issue of water resource allocationand private ownership of water resources still remains. However, the assumption of ownership of the water in the law of equitable distribution of water does not mean absolute ownership of holders. Since the fluid nature of groundwater prevents the physical possession of the underground source by the holder, it can be said that traditional rules of use of groundwater do not generally provide any right of ownership that the holder can appeal to the legal system. Based on the studies of ministry of energy, different regions of the country have been divided into two forbidden and free zones in terms of the possibility of exploiting underground water. In free zones, individuals can submit a written request to the regional water companies in each province to get a license for digging wells or constructing subterranean and exploiting underground water resources for any purposes. The use of groundwater resources up to 25 cubic meters per day for domestic, drinking, sanitation, and watering use does not need a permit in these areas and people can only use up to the abovementioned figure by informing the Ministry of Energy. In the forbidden zone, licensing well digging for industrial and civil projects of government with the approval of the Ministry of Energy and drinking water supply for cities, towns and villages with the offer of relevant ministries or their affiliated organizations. Industrial and drinking water, livestock and poultry health water supply with the certificate issued by the competent authorities will be possible based on the water facilities and the approval of the Ministry of Energy. Obvious point about the law is that tens of thousands of different beneficiaries didn’t see it necessary to make use of law for getting legal license. Because offenders, within three decades of experience, have found that regulated well digging leads to government’s supervision over their destiny and endless harassment upon them. This is the process that began and strengthened with the equitable distribution of water and subsequent laws. One of the opportunities created in forbidden areas for well digging licensing is related to livestock and greenhouse units. Experience has proven that these licenses are solely a certificate to permit legal well digging and investigations show an overdraft of more than 60 to 70 percent of these units. It can be asserted that more than 60% of these units have virtually no greenhouse or livestock activities and just tremendously extract water and continue traditional agriculture and high-consuming activities (in term of water). In many cases, the owners have used the value added of converting rain-fed land around wells to water-fed lands and try to sell with high prices.
    Keywords: Groundwater resources , Wells , ownership , Forbidden regions , Operation permit
  • MEHDI KOMASI *, SOUDEH GHOBADI KHOSRO, mohammad reza goodarzi, mohammadreza hashemi Pages 81-90
    Drought is a temporary and recurring meteorological event, originating from a lack of precipitation over an extended period of time. The success of drought preparedness and mitigation depends on timely information about drought onset and forecasting. This information may be obtained through continuous drought monitoring, which is normally performed using drought indices. Drought is an unpleasant, naturally occurring event caused by climate change that directly affects societies through changing their access to water resources. Among the numerous indices for drought intensity rating, the EDI and SPI have widespread applications. The SPI was computed by fitting a probability density function to the frequency distribution of the monthly precipitation records of each station. A drought event is considered to occur at a time when the value of the SPI is continuously negative and ends when the SPI becomes positive. The computation of the SPI drought index for any location is based on the long-term precipitation record (at least 30 years) cumulated over a selected time scale. This long-term precipitation time series is then fitted to a gamma distribution, which is then transformed through an equal probability transformation into a normal distribution. Positive and negative SPI values respectively indicate wet conditions (greater than median precipitation), and dry (lower than median precipitation). In most cases, the probability distribution that best models observational precipitation data is the Gamma distribution. Unlike most other drought indices, the EDI in its original form is calculated with the daily. The resulting EDI value represents standardized value for currently utilizable water resources, considering the continued dry period. If a negative DEP continues for more than 1day, the addition period of EDI will increase as long as the continued days. This variable addition period is limitless. The nature of genetic programming allows the user to gain additional information on how the system performs, i.e., gives insight into the relationship between input and output data. The GP is similar to genetic algorithm (GA) but unlike the latter, its solution is a computer program or an equation as against a set of numbers in the genetic algorithm. So, GP is more attractive than traditional GA for problems that require the construction of explicit models. The GP thus transforms one population of individuals into another one, in an iterative manner by applying operators. In evolutionary computation, it can distinguish between three different types of operators which are named crossover, reproduction, and mutation. M5 model tree approach is based on the principle of information theory that makes it possible to split the multi-dimensional parameter space and generate the models automatically according to the overall quality criterion. It allows for variation in the number of models created. The splitting in the M5 modal tree approach follows the idea of decision tree, but instead of the class labels, it has linear regression functions at the leaves, which can predict continuous numerical attributes. Model trees generalize the concepts of regression trees, which have constant values at their leaves. Therefore, they are analogous to piece-wise linear functions (and hence nonlinear). Computational requirements for model trees grow rapidly with increase in the dimensionality of the data set. Model trees learn efficiently and can tackle tasks with very high dimensionality. The major advantage of model trees over regression trees is that model trees are much smaller than regression trees and regression functions do not normally involve many variables. This research used precipitation data on two basins in Hamedan and Lorestan Provinces to calculate the SPI and EDI indices for monitoring drought. The genetic programming model and M5 model trees were used to predict the occurrence of drought in these two basins. It was found these models had good capability in predicting drought and enjoyed high accuracy in solving prediction problems. Another advantage of these models is that they use simple equations for predicting the phenomena under study. In the best-case scenario, the coefficients of determination for the EDI index in the M5 model trees and in the genetic programming model were 0.97 and 0.95, respectively. Moreover, the coefficients of determination for the SPI index in the M5 model trees and in the genetic programming model, in the best-case scenario, were 0.93 and 0.83, respectively. This suggests the M5 model trees are more accurate compared to the genetic programming model and enjoy relative superiority because they are simpler and more understandable than the genetic programming model.
    Keywords: Data driven , EDI Index , Modeling, SPI Index , Silakhor plain
  • Ebrahim Ebrahimi, Mohammad Ebrahim Banihabib, Abbas Roozbahani * Pages 91-100
    The aquifers are about four percent of the water on the earth, but they are considered as the best and most accessible source of fresh water. In recent years, they have been faced with severe water withdrawal, therefore some plains was considered as forbidden plains that it means the water withdrawal from these aquifers is unauthorized. At some point, plains have been faced with land subsidence that showed the severity of the disaster. Given such a critical situation in aquifers, management of groundwater resources in the form of tools such as monitoring the level of the aquifers is essential. One of the plains in Iran that has a critical groundwater resource is Birjand plain which requires management measures to be protected from future water resources crisis. Prediction of groundwater level in future periods is a useful tool to enforce management measures before a crisis occurs. Thus, in this study, groundwater level was predicted in Birjand aquifer taking 10 monthly forecasting scenarios in a period of 5 years and in both crisp (continues) and clustering approaches using probabilistic Dynamic Bayesian Networks (DBNs). Nowadays, various tools are used to predict the aquifer level including mathematical models, artificial neural network, neuro-fuzzy, Bayesian networks, and time series and so on. In recent years, due to the flexible and simple structure, Bayesian networks have been used for predictions of different parameters, especially in forecasting of hydrological parameters. Bayesian network as a modern forecasting probabilistic method shows probabilistic relationships between a set of variables by graphical model. It represents the dependence structure among several factors, that affecting on each other, and is based on Bayesian theory. Dynamic Bayesian Networks have been extended from Bayesian Networks which are created for two
    purposes
    first, as the cycle of dependency detector over the time, quite similar to Markov model; second, as the fixed process which is repeated in fixed-time interval. Another feature of Dynamic Bayesian network is their willingness to approximately structured changing. In this study, the input data (predictor parameters) of the model includes: temperature, rain, evaporation, monthly aquifer recharge in each Thiessen Polygon (Recharge), monthly withdrawals of groundwater in each aquifer Thiessen Polygon (discharge), groundwater levels in the current month and the groundwater level in the next month (predicted parameter). The first step in modeling by the dynamic Bayesian network is determining the dependent and independent data for calibration and validation. Model calibration data in both crisp and clustering approach include a 12-year period (1998 to 2009) and data validation include 5-year period (2010 to 2014), in monthly time step. Depending on the type of input data, both crisp and clustering approach was used. In the crisp approach, the probability functions were used and the predicted data were obtained by using the training data. In the clustering approach, by assigning each of the numbers to the right cluster, the modeling was done. A cluster is collection of objects that their relative distance to each other is low and relative to other members is high. In the first approach by using crisp data and NPC training structure in confidence level of 5%, the training was applied. In the second approach, the Silhouette index was calculated by using MATLAB software and by using the validation Silhouette index, and then, the number of clusters was determined. Finally, the clustering was performed by using k average method. Then, training was done in the second approach, using clustered data and NPC training structure in confidence level of 5%. Considering 10 different scenarios to predict groundwater level in the next month, the uncertainty of predicted parameters in both crisp and clustering approach were assessed. In fact, using these scenarios, sensitivity analysis was performed to check the accuracy of the model with respect to the existence or absence of different predictor parameters. In addition, the uncertainty of the model output is evaluated using dynamic Bayesian network probabilistic analysis. The results of the selected scenario in crisp approach showed the high prediction accuracy of Bayesian networks. For example, in piezometers 4 and 9, the coefficient of determination was estimated about 0.98. According to the results, crisp dynamic Bayesian networks approach predicted hydrograph aquifer more accurate than clustering approach. Due to the low efficiency of clustering approach in predicting groundwater hydrograph, to obtain the accuracy of 10 scenarios predictions in this approach, instead of coefficient of determination (R^2) and root-mean-square error (RMSE), the percentage of correct predicted clusters was determined. According to the results, the clustering approach predicts clusters with high accuracy for different piezometers. The scenario 6 had the best prediction which all predictor parameters except evaporation were used for the prediction. In this scenario, R^2 and RMSE were showed good accuracy as 0.9946 and 0.1275, respectively. The other scenarios had also the accuracy in their predictions very close to scenario 6, except scenario 9 which the groundwater levels were not used as input. Thus, in the crisp approach, the groundwater levels have a substantial impact on the accuracy of the prediction. Generally, in crisp approach, all predictor scenarios have acceptable accuracy of over than 90% except scenario 9. In clustering approach, by examining the accuracy of the scenarios in predicting clusters, most scenarios were accurate except scenario 9. Dynamic Bayesian Network model in clustering approach correctly predicted clusters but, unlike the crisp approach in predicting hydrograph, it was not able to present acceptable results. The results showed the ability of the proposed model in planning and management of groundwater resources, reducing the risk of aquifer level declining by applying short term management scenarios and predict its effects on rehabilitation. Moreover, this model can be used in the similar plains for aquifer management.
    Keywords: Bayes' Theorem, Clustering , Groundwater management, Sensitivity analysis
  • Amir Hossein Parsamehr *, Hossein Malekinezhad, Zahra Khosravani Pages 101-110
  • Farshd Ahmadi *, Feridon Radmaneh, Gholam Ali Parham, Rasoul Mirabbasi Pages 111-120
    Estimation of flood discharge and return period of extreme events is one of the most important factors in design of hydraulic structures. Therefore, conducting researches for improving the accuracy of estimation methods is very important. Nowadays, the fractal behavior of all earth science phenomena, including precipitation and runoff has been accepted. This method is used for modeling of natural phenomena because of its strong physical background. As power law is derived from fractal theory, it can be a suitable method for estimation of design flood. In the annual maxima method, a maximum value is extracted among all the recorded discharge occurred during a year. Therefore, in this method, out of a number of extreme events happened at the end of the year, only one value is chosen and other values are eliminated. The positive point of this method is that the extracted data are independent of each other. However, this method has its own limitations. For example, in years with high flood occurrences, the values that stand in the second and third places in terms of magnitude lose the chance to participate in the annual maxima series. On the other hand, in dry years, the occurred maximum values are far lower than the second or third ranked discharge of wet years. Thus, correspondence of the annual maximum time series with reality is not suitable, and such a selection may result in neglecting the real pattern in the recorded values. In another method, partial time series (the time series greater than a given threshold) is employed. In this sampling method, a certain threshold of flow rate is considered and all flow rates that are above the threshold, are taken into consideration. In this study, power law has been employed to select the flood partial series. Therefore, the main objective of this study is to use power law in a) selection of partial series of flood in Dez Basin, and b) flood frequency analysis. In this study, flood frequency analysis of Dez watershed during the period of 1956- 2012 was performed using two approaches of power law (PL) and ordinary statistical distributions. At first, the partial-duration flood series were extracted using power law with the intervals of 7, 14, 30, and 60 days, and then compared with the annual maxima series. The results indicated that the annual maxima series were not suitable for frequency analysis of the flood in Dez Basin, and the 30-day partial duration series obtained from the power law has a better correspondence with the flow and properties of the Dez Basin. Examining of the independence and stationarity of the 30-d partial duration series by Wald-Wolfowitz (W-W) test, confirmed the independence of the considered series. In the next step, the power distribution and ordinary statistical distributions were fitted onto the flood data of the Dez Basin, and performance of each distribution was investigated using Normalized Root Mean Square Error (NRMSE) and Nash-Sutcliffe (NS) criteria. The results showed that in the Tange Panj Bakhtiari station (TPB), power distribution had a better performance than other considered distributions, and in the Tange Panj Sezar (TPS) and Tale Zang (TZ) stations the power distribution stood in the second rank in terms of the best distribution with a slight difference. The results indicated that the annual maxima within shorter return periods estimate lower values for the flood of the studied stations. However, the partial series extracted using power law solved the problem in the annual maxima, and are in full correspondence with the physics of Dez Basin. Following selection of suitable statistical distributions, flood quintiles were calculated for return periods of 2, 5, 10, 20, 50, and 100 years in each station. Based on the calculated quintiles, it was observed that TPB station is more vulnerable to the flood occurrence, where the amount of flood increases in it with the prolongation of return period. In other words, in basins culminating in TPB station, development of warning systems and implementation of flood control programs are more needed. Since the performance of power distribution in estimating the flood in Dez Basin has been very satisfactory and calculation of its parameters and its application is easier than ordinary probability distributions, it can be suggested as the superior distribution for flood frequency analysis in Dez Basin.
    Keywords: Flood , Fractal , Partial Series, Stationarity
  • Reza Ghazavi *, Ebrahim omidvar, farzaneh Fadaii Pages 121-129