Frequency Analysis of Low-Flows in the Large Karoun River Basin in Iran

Similar documents
Prediction of the Weekly Mean Discharge into Reservoir of Doroudzan Dam Using HEC-HMS Model and its Comparison with Observed Data ( )

A peaks over threshold approach for hydrological drought severity assessment.

SEES 503 SUSTAINABLE WATER RESOURCES. Floods. Instructor. Assist. Prof. Dr. Bertuğ Akıntuğ

THE QUANTIFICATION OF THE THRESHOLD LEVEL METHOD ON LOW FLOWS STUDIES

Conference Proceedings Paper Prediction of Annual Inflow to Karkheh Dam Reservoir using Time Series Models

A simple model for low flow forecasting in Mediterranean streams

Comparison of Probability Distributions for Frequency Analysis of Annual Maximum Rainfall

Examining the Effectiveness of Check Structures in reducing Flood Discharge of Kotok Watershed in Khuzestan Province in Iran

PART IV WATER QUANTITY MONITORING, TECHNOLOGICAL ADVANCES AND CONCLUSIONS

Relationships between low-flow characteristics of South African streams

Determination of Flood Plain Zoning in Zarigol River Using the Hydraulic Model of HEC-RAS

CHAPTER ONE : INTRODUCTION

Hydrological Flood Routing in Rivers

Assessing ecological flow conditions for wetlands fed from ungauged stream reaches

ANALYSIS OF TIME OSCILLATIONS OF WATER ON LAKE MODRIC AS A MULTI-PURPOSE RESERVOIR

Europe - IMBO th European Conference on the Implementation of the WFD

Monitoring of water levels and flood discharge

Impacts of Rainfall Event Pattern and Land-Use Change on River Basin Hydrological Response: a Case in Malaysia

Water Resources in Korea

LARGE DAM EFFECTS ON FLOW REGIME (CASE STUDY: KARKHEH RIVER, DOWNSTREAM OF RESERVOIR DAM)

Study of Hydrology based on Climate Changes Simulation Using SWAT Model At Jatiluhur Reservoir Catchment Area

RIVER DISCHARGE PROJECTION IN INDOCHINA PENINSULA UNDER A CHANGING CLIMATE USING THE MRI-AGCM3.2S DATASET

Low Flow Frequency Analysis on Selected Rivers in Malaysia

PARAMETERS OF LOW FLOW AND DATA ON LOW FLOWS IN SELECTED IRISH RIVERS

Are we water rich or water poor. Levente Ronczyk: The Challenge of Water Resources Management in Hungary

HYDROLOGICAL DROUGHT INDEX AT THE IRRIGATION AREA IN PEMALI-COMAL RIVER BASIN

STORAGE-YIELD CURVES WITH INFLOWS FROM A DIVERSION CHANNEL

Long-term change of stream water quality as a consequence of watershed development and management

CE 2031 WATER RESOURCES ENGINEERING L T P C

Impact of Future Climate Change on the Water Resources System of Chungju Multi-purpose Dam in South Korea

Flood Forecasting using Adaptive Neuro-Fuzzy Inference System (ANFIS)

Journal of Water Resource and Hydraulic Engineering Jul. 2015, Vol. 4, Iss. 3, PP

CHAPTER FIVE Runoff. Engineering Hydrology (ECIV 4323) Instructors: Dr. Yunes Mogheir Dr. Ramadan Al Khatib. Overland flow interflow

Ontario Flow Assessment Tool

Simulation of Climate Change Impact on Runoff Using Rainfall Scenarios that Consider Daily Patterns of Change from GCMs

ACTIVITIES FOR COMBATTING AGRICULTURAL DROUGHT IN TURKEY. I. Introduction

Review of Hydrological Potential in Combined Gin and Nilwala River Basins of Sri Lanka

Simulation of Rainfall-Runoff Using WEAP Model (Case Study: Qaraso Basin)

Low flow estimation in Austria

Assessment of Climate Change Impacts on Water Resources and How to Adapt

International Journal of Advance Research in Engineering, Science & Technology

COMPARISON OF ESTIMATED FLOOD HYDROGRAPHS USING LUMPED AND DISTRIBUTED MODELS

Comparison of Streamflow and Precipitation in the Upper Provo River Watershed

Latest tools and methodologies for flood modeling

Study of drought with SPI index (case study: Ghareh Chai and Karkheh basins)

Karun River Iran. May 2016 Sepideh Nayemi International cooperation coordinator Geneva Water Bureau- (Direction générale de l'eau )- DETA

BIABAN Journal Vol 11, No 1 (2006) PP , Salajegheh and Dastorani

What we are looking for through utilizing the global earth observation. 14 January 2012 Japan International Cooperation Agency (JICA) Yusuke AMANO

Including the human factor in adaptive actions to face climate change

Flood control and safety evaluation of Banqiao reservoir

Analysis of the Sensitivity of the Variables Affecting Urban Flood through Two Rational and SWMM Methods by SSA

of Surface Water Data (Environment Canada, 1994). The primary database consisted of

CASE STUDY OKUYOSHINO, JAPAN

Distribution Restriction Statement Approved for public release; distribution is unlimited.

Comparison of Artificial Neural Networks ANN and Statistica in Daily Flow Forecasting

6-5 Analysis of water storage state by the subsurface dam

HYDRODYNAMIC SIMULATION OF SURFACE WATER CONTROL SLUICE GATES BY HEC-RAS MODEL

The Fourth Assessment of the Intergovernmental

Analysis of Future Discharge in Nakdong River Basin According to Effect on Climate Change

Chapter 1 Introduction

Background and History of Water Management of Okanagan Lake and River

Water balance and observed flows in the Anllóns river basin (NW Spain).

PRELIMINARY DESIGN OF THE HYDRAULIC STRUCTURES DAM IN THE PISÃO RIVER

A comparative study of the methods for estimating streamflow at ungauged sites

Risk based decision making in Kaunas hydropower system during flood period

UNIT HYDROGRAPH AND EFFECTIVE RAINFALL S INFLUENCE OVER THE STORM RUNOFF HYDROGRAPH

Module 3. Lecture 6: Synthetic unit hydrograph

A hydro-ecological assessment method for temporary rivers. The Candelaro river case study (SE, Italy)

Hydrological Modelling of Narmada basin in Central India using Soil and Water Assessment Tool (SWAT)

4. Present Activities and Roles

Oregon Water Conditions Report January 11, 2017

Brannen Lake Storage Feasibility Potential Effects on Water Levels

Water Balance Assessment of the Roseires Reservoir

ANALYSIS OF RAINFALL DATA FOR DROUGHT INVESTIGATION AT BRAHMAPURI (MS)

The Development of Hydrological Services

CENTRAL ASSINIBOINE INTEGRATED WATERSHED MANAGEMENT PLAN SURFACE WATER HYDROLOGY REPORT

Federal Democratic Republic of Ethiopia STATE PARTY REPORT IN RESPONSE TO THE WORLD HERITAGE COMMITTEE DECISION WHC 34 COM 7B.44

Probabilistic Analysis of Peak Daily Rainfall for Prediction purposes in Selected Areas of Northern Nigeria.

Hydrologic frequency analysis: History and contemporary research. Cameron Bracken University of Colorado at Boulder

Module 2 Measurement and Processing of Hydrologic Data

August 6, Min-Ji Park, Hyung-Jin Shin, Jong-Yoon Park Graduate Student. Geun-Ae Park. Seong-Joon Kim

MODULE 1 LECTURE NOTES 2 MODELING OF WATER RESOURCES SYSTEMS

Climate Change in the Columbia Basin. Stephanie Smith Manager of Hydrology, BC Hydro

Hydrology and Water Management. Dr. Mujahid Khan, UET Peshawar

Flood Inundation Mapping of Tadi River

Second meeting of the Task Force on the Water-Food-Energy-Ecosystem Nexus, 8-9 September 2014,

ESTIMATING THE EFFECT OF CLIMATE CHANGE ON THE HYDROLOGY OF THE NILE RIVER IN THE 21 ST CENTURY

Simulation of a reservoir with standard operating rule

Optimization of long-term Han River basin hydropower scheduling using linear programming

Comparison of hedging rule curves depending on available water and current storage in reservoir operation under droughts

NAP-GSP REGIONAL TRAINING WORKSHOP FOR ASIA Mainstreaming Climate Change Adaptation into Water Resources Seoul, September 2017

WATER QUALITY MODEL FOR EL CAJON AND AGUAMILPA RESERVOIRS. Edith G. Lugo. Rushit Hila. Oliver Obregon. Final Report

Regional Association VI, Forum Hydrology, Koblenz May 8 10, Statements of the national hydrological service of the Czech Republic

Analysis of Flood Routing

Country proposal - SRI LANKA

Water Supply Development in Iran during Recent Drought

RAINFALL-RUNOFF RELATIONSHIPS AND FLOW FORECASTING, OGUN RIVER, NIGERIA

Closure Water Balance Model to Support Closure Designs for a Mine in Laos

1 THE USGS MODULAR MODELING SYSTEM MODEL OF THE UPPER COSUMNES RIVER

Hydrological And Water Quality Modeling For Alternative Scenarios In A Semi-arid Catchment

Transcription:

International Research Journal of Applied and Basic Sciences 2013 Available online at www.irjabs.com ISSN 2251-8X / Vol, 4 (1): 146-151 Science Explorer Publications Frequency Analysis of Low-Flows in the Large Karoun River Basin in Iran kazem Hamadi 1, Samaneh Abdovis 2, Narges Zohrabi 3 * 1. Department of Basic Sciences of Water Resources, Khouzestan Water and Power Authority, Iran. 2. Department of Basic Sciences of Water Resources, Khouzestan Water and Power Authority, Iran. 3. Department of irrigation, Science and Research Branch, Islamic Azad University, Khouzestan Iran. * Corresponding author email: n.zohrabi@khouzestan.srbiau.ac.ir ABSTRACT: Frequency analysis is one of the suitable statistical methods in assessing the amounts of low-flow in various return periods which could be obtained using statistical distribution theories and choosing proper probability law or distribution function. The purpose of this study is a frequency analysis for low-flows of Karoun, Dez, and the Large Karoun river basin in Iran, find the dominant probability distribution functions and the low-flows' values until their behaviors in 1-day to 30-day periods could be determined. The observed statistics and information from three hydrometric stations were used in this investigation, namely, Polshaloo as the inflow to the Karoun-3 storage dam, Talezang as the inflow to the Dez storage dam, and Ahvaz as the performance location of the Large Karoun system. The results from this study showed the low-flow generation in Karoun basin (Polshaloo) and the Dez basin (Talezang) are closed together in low probabilities. Analysis of daily to monthly low-flows of the Large Karoun River in Ahvaz station varies between 190 to 250 m 3 /s. These values for 80% frequencies are between 143 to 183 m 3 /s. The time series analysis of 1-day to 30- day low-flows of this station shows that flow during the first 15-year period prior 1969 had been normal and during the following years it was affected by the basin s soil and water resources development plans. Keywords: Frequency Analysis, Low-Flows, Water Shortage, Large Karoun River. INTRODUCTION International glossary of hydrology defines low flow as flow of water in a stream during prolonged dry weather. This definition does not make a clear distinction between low flows and droughts. Low flows is a seasonal phenomenon, and an integral component of a flow regime of any river. Drought, on the other hand, is a natural event resulting from a less than normal precipitation for an extended period of time. McMahon and Arenas (1982) presented a compilation of methods used in different countries to compute low-flow characteristics, illustrated some of these methods with case studies and addressed some theoretical aspects of natural and man-induced factors affecting low flows. Smakhtin (2001) presented a paper with matter "Review Low flow hydrology" that further focuses on the techniques for low-flow estimation in ungauged river catchments. This paper presents a summary of recent international low-flow related research initiatives. The review is largely based on the research results reported during the last twenty years. Analyzing rivers discharge and flood water is one of the most important statistical topics in the field of the applied hydrology science. The observed data is a primary and substantial requirement for every frequency analysis. whenever there are sufficient statistical samples and the possibility to fit a proper distribution function, the frequency analysis is a suitable statistical method in assessing discharges in different return periods. Rao and Hamed (2000) for calculating the statistical parameters of density or distribution functions suggested several methods could be used including: method of moments, maximum likelihood, probability weighted moments, least squares, etc. Usually, the first three mentioned methods are common. The lowflow is the lowest mean flow in one day or 3, 5, 7, 30 or 60 consecutive days during one year. Modarres (2008) Stated the reason for calculating the low-flow for some consecutive days is to minimize the flow s small oscillations and variations. Tusker (1987) with using statistical distributions studied the hydrometric statistics of about 20 rivers in Virginia State of US. He suggested the 3-parameter Weibull and Log-Pearson distributions for analyzing the 7-day low-flow series. Yureklh and Kurunc (2005) properly used the low-flows

analysis for drought analyses. For study drought they applied the daily low-flow information and statistics of each month observed in three hydrometric stations located at the Cockrick basin in Turkey. The objective of this study is the frequency analysis of Karoun and Dez rivers low-flow at inflows to Karoun-3 and Dez dams reservoirs to while clarify probability distribution functions also to determine the low-flow amounts and their behavior during 1-day and 30-day periods. Also, the low-flow time series for Ahvaz hydrometric station at the large Karoun system is discussed. MATERIALS AND METHODS In this study, the information and statistics from three hydrometric stations located at large Karoun basins including Karoun-Polshaloo, Dez-Talezang, and large Karoun-Ahvaz were used. The observed flow in Polshaloo and Talezang hydrometric stations are as inflows to the reservoirs of Karoun-3 and Dez dams. The common statistical period for stations was 1957-2008. The aforementioned data were acquired from Ministry of Energy [Iran Water Resource Management Company and Khouzestan Water and Power Authority (KWPA)].The geographic coordinates for the selected hydrometric stations show in table (1). The equipments and number of stations in during sampling were sufficient and ranked as grade 1 in terms of hydrometrics. Regularly, discharge measurement, sectioning, and water sampling for qualitative tests and sediments are performed in these stations. In 1955, the first sampling from Khouzestan surface water resources began at the Gotvand section, and during the following years, Polshaloo and Talezang hydrometric stations became operational. By oprating of Karoun-3 storage dam, the Polshaloo station was placed at the dam s reservoir (close to the dam s axis) where that the inflow is considered as the recorded flow in this station. In fact, recently, due to opration of the Karoun-4 storage dam, the flow normalization requires a series of hydrological modifications. The general characteristics of flows in studied hydrometric stations are appended to table and figure(1).

Figure1. Location of selected hydrometric stations. Table 1. Coordinates and general characteristics of flows in selected hydrometric stations River s name Karoun Dez Large Karoun Station s name Polshaloo Talezang Ahvaz Longitude (UTM) 417051.32 291061.3 280701.07 Latitude (UTM) 3515497.8 3633751 3469350 Height (m) 690 444 9.78 Statistical duration (Year) 50 53 54 Basin s area (km 2 ) 24210 16213 60737 Annual flow volume (million m 3 ) 9429 7852 19962 Annual flow variations coefficient 0.37 0. 0.40 Purpose of frequency analysis is obtaining the probable law of density function or distribution function dominant on the studied random variable according to the observed sample and calculating the return period or x T occurrence value (Rao and Hamed, 2000). Generally, the frequency analysis is the probable investigation of the flow s random variable including: sample collection (minimum daily or more consecutive days discharge), sample data mending and completion, required statistical test such as random, independency, integrity, outliers and stationary, choosing proper statistical distribution (density function f(x) or distribution function F(x)), assessing the density function parameters and the required tests for density function fitting and also choosing the best function. In order to achieve the research s goal, including the low-flows in 1-day or more consecutive days of a year, hydrometric statistics from the above-mentioned stations were considered in daily scale for the statistical period mentioned in table (1). Provided the daily flow statistics of the mentioned sections for more than 50 years, the following steps were taken: a) Using the daily discharge, the minimum daily discharge for each station-water year was determined. In this method the minimum discharge from among discharges of 365 days of each year was determined and considered for same year. Thus, for each hydrometric station a 50-years time series of low-flow was developed. b) Since 3, 5, 7, 15, and 30 consecutive days low-flows are aimed in this study, using the daily discharge and the moving average technique, the discharge of 3, 5, 7, 15, and 30 consecutive days low-flows was determined for each station. c) For the series frequency analysis, the Hydrological Frequency Analysis software(hyfa), and such statistical distribution as normal, log-normal (II), log-normal (III), gamma (II), Pearson (III), log Pearson (III) and Gumball were used. for calculate the distribution parameters, two methods were applied, namely, the moment and the maximum likelihood method. The best statistical distribution corresponding to data was selected based on the minimum mean relative deviations of the flow observed values and their calculated values in distribution, relative deviations mean square, K-square (K 2 ), and comparing the distribution fitting on observed data separately for each station. The low-flow with different return periods was determined based on this distribution. RESULTS AND DISCUSSIONS In this study, the minimum discharges of Karoun, Dez, and large Karoun rivers were considered. The statistical period is 52 year. First the preliminary analysis on time series of low-flows in one day or more consecutive days of each year including mean, maximum and minimum parameters, standard deviation, range and coefficient variations for each hydrometric station was conducted. The results are given in table (2). As shown in table(2), the time series ranges of one-day, three-day and to monthly time series of Polshaloo and Talezang hydrometric stations have slightly changed. This is due to flow consistency in dry seasons. Moreover, the annual flow coefficient of variations for these stations are about 37%, while the low-flow time series coefficient of variations for Karoun River is 27% and for Dez River is 21%. The Ahvaz hydrometric station does not follow the above-mentioned trend due to regulate upstream flow and there is a significant difference between one-day and monthly time series flow. The HYFA hydrology software was used for frequency analysis. The moments and maximum likelihood methods assessed the statistical distribution parameters. The best statistical distribution was determined based on the mean relative deviations of observed and calculated flows in the distribution, relative deviations mean squares, k-square and comparing the distribution fitting on observed data separately for each station and time series. The analysis results are presented in table (3). 148

Table 2. the preliminary analysis of low-flow time series of hydrometric stations (m3/s)- Pol shalo Karoun river-pol shalo station Parameter/time series 1-day 3-day 5-day 7-day 15-day Monthly Mean 96 97 98 99 100 103 Maximum 164 166 168 168 169 179 Minimum 29 36 36 39 40 Standard deviation 26 27 26 29 27 27 Range 135 130 132 130 131 139 Coefficient of variations 0.27 0.27 0.27 0.29 0.27 0.27 Dez river-talezang station Parameter/time series 1-day 3-day 5-day 7-day 15-day Monthly Mean 56 56 56 57 58 61 Maximum 101 101 101 101 102 104 Minimum 34 35 35 36 37 39 Standard deviation 12 12 12 12 12 13 Range 67 66 66 65 65 65 Coefficient of variations 0.22 0.21 0.21 0.21 0.21 0.22 Large Karoun river-ahvaz station Parameter/time series 1-day 3-day 5-day 7-day 15-day Monthly Mean 189 199 209 215 232 250 Maximum 328 333 344 353 410 434 Minimum 86 90 91 92 100 107 Standard deviation 54 57 60 64 72 79 Range 242 244 253 261 310 327 Coefficient of variations 0.28 0.28 0.29 0.30 0.31 0.31 Table 3. The analysis results for frequency distribution functions fitting of low-flows in the studied hydrometric stations Karoun river-pol shalo station Time series Type of function Method of assessment K2 Sq.Rel.Dev Rel.Dev 1-day log Pearson (III) Momentum 1.1 44.44 3.56 3-day Normal Momentum 0.9 13.70 2.12 5-day Normal Momentum 0.90 13.63 2.25 7-day Normal Momentum 1.91 18.57 2.43 15-day Normal Momentum 0.08 12.04 2.18 Monthly Normal Momentum 1.30 12. 2.05 Dez river-talezang station 1-day gamma(ii) maximum likelihood 4.45 12.09 2.14 3-day Gumball(I) maximum likelihood 2.00 8.76 1.96 5-day log Normal(II) Momentum 2.18 9.72 2.09 7-day log Normal(II) Momentum 2.75 8.69 2.75 15-day Gumball(I) maximum likelihood 0.87 6.28 1.70 Monthly log pearson (III) Momentum 2.00 4.27 1.65 Large Karoun river-ahvaz station 1-day log Pearson (III) Momentum 3.59 11.44 2.66 3-day gamma(ii) maximum likelihood 1.74 16.17 2.39 5-day log Normal(III) Momentum 0.81 14.51 2.85 7-day gamma(ii) maximum likelihood 1.37 18.96 2.67 15-day log Normal(III) Momentum 1.37 9.00 1.91 Monthly log Pearson (III) Momentum 0.63 9.70 2.00 According to the information provided in table (3) and also using the best distribution function, the low-flow discharge with various return periods was determined. Table (4) indicates the frequency analysis results for low-flows in the studied hydrometric stations. The area of Karoun at Polshaloo and Dez at Talezang catchment-basin is 24210 and 16213 km 2, respectively. Hence, the area ratio of generating basin of Dez River flow to Karoun in the mentioned stations is equal to 0.67. But the mean ratio of Dez branch to Karoun annul flow is 80%. This ratio could be 58% for lowflows average discharges, which as the occurrence probability decreases, this ratio increases up to 96%. That is, the minimum flow generation of the Karoun basin (Polshaloo) gets closer to the Dez basin (Talezang). Analyzing the daily to monthly low-flows of Karoun-e-Bozorg River at the Ahvaz station varies between 190 and 250m 3 /s. These numbers for 80% frequencies are between 143 to 183m 3 /s. It should be noted that the analyzed figures for Ahvaz station are related to regular and normal data for the last 50 years and the upstream consumptions as well as returned water are effective on them. Figures (2) to (4) compare the one-day and thirty-day time series of low-flows for studied hydrometric stations. As shown in these figures, these series are close to each other in Karoun (Polshaloo) and Dez (Talezang) rivers. This means that the low-flows for at least thirty days are consistent in each day. Indeed, this is slightly different for dry seasons. As shown in table (1), the applied statistical period in analyzing the large Karoun river at Ahvaz hydrometric station is 54 years. The 149

time series analysis of one-day and thirty-day low-flows of this stations shows that during the first 15 years of the statistical period before 1969 the flow was relatively normal and for the following years it was affected by basin s soil and water resources development plans figure (4). Table 4. The results from frequency analysis of low-flows at the studied hydrometric stations (m /s) Karoun river-pol shalo station Time sequence Mean 50% 60% 70% 75% 80% 90% 95% 97.5% 1-day 96 1000 93 85 80 75 61 50 3-day 5-day 97 98 97 98 91 96 83 92 79 81 75 76 63 65 54 55 45 47 7-day 99 999 92 85 81 77 65 56 15-day 100 100 93 86 82 78 66 56 48 Monthly Dez river-talezang station 1-day 103 56 103 555 96 52 89 49 85 47 80 46 68 41 58 3-day 5-day 56 56 54 55 51 52 49 49 47 48 46 46 43 42 40 39 37 7-day 57 56 55 53 48 47 43 39 15-day 58 56 54 51 50 48 45 42 40 Monthly 61 59 56 53 51 50 46 43 Large Karoun river-ahvaz station 1-day 189 185 172 158 151 143 123 108 3-day 5-day 199 209 194 205 180 190 166 175 159 167 151 158 131 135 117 117 105 102 7-day 215 209 194 178 170 161 139 123 15-day 232 226 209 191 182 172 146 126 110 Monthly 250 2444 225 205 195 183 156 135 99% 42 32 35 37 47 49 39 35 32 36 34 37 34 41 95 81 92 85 110 96 92 119 101 Figure 2. Time series comparison of low-flows Karoun River Polshaloo Figure 3. Time series comparison of low-flows Dez River Talezang 150

Figure 4. Time series comparison of low-flows large Karoun Ahvaz ACKNOWLEDGEMENT We thank the Ministry of Energy, Water Researches and Engineering Standards of Khouzestan Water and Power Authority in Iran for supporting this research. REFERENCES McMahon TA, Arenas AD. 1982. Methods of computation of low streamflow. Paris, UNESCO Studies and reports in hydrology., 36:107. Modarres R. 2008. Regional Frequency Distribution Type of Low Flow in North of Iran by moment. Journal Water Resource Manage.,22:823-841 Rao AR, Hamed KH. 2000. Flood Frequency Analysis, CRC Press, Boca Raton, 350 p. Smakhtin VU. 2001. Review Low flow hydrology: a review. Journal of Hydrology., 240(3) : 186. Smith J, Jones MJ, Houghton LD. 1999. Future of health insurance. N Engl J Med. 965: 325 329. Tasker GD. 1987. A comparison of metods for estimating low flow characteristics of streams", JAWRA Journal of the American Water Resources Association., 23(6) : 1077 1083. WMO: World Meteorological Organization, 1974. International Glossary of Hydrology, WMO, Geneva. Yureklh K, Kurunc A. 2005. Performances of Stochastic Approaches in Generating Low Stream flow Data for Drought Analysis. Journal of Spatial Hydrology Spring., 5 (1). 151 14 47