Hydrological Modeling Report
|
|
- Austen Terry
- 6 years ago
- Views:
Transcription
1 HYDRO EUROPE 2011 Hydrological Modeling Report Team Team members: Manuel GOMEZ (supervisor), Liza ASHTON, Thibault DESPLANQUES, Jens Harold DRASER, Panagiota GKAVAKOU, Kyung Tae LEE, Maxime NARDINI, Jieun PARK, Jiyoung PARK, Gabor SZABO & Anastassi STEFANOVA
2 Table of Contents Table of Contents... 1 List of Figures... 2 List of Tables Introduction Rainfall distribution Comparison of rainfall distribution methods Lumped or quasi distributed method Thiessen Polygon Method Distributed Method (Best Linears Unbiaised Estimators) Inverse Distance Weights Spline Kriging Conclusion Land Use Analysis Lag Time & Time of Concentration Hydrological modeling with HEC-HMS Why a model? Model Construction Steps Transform Method Routing Method Loss Method Baseflow Method HEC HMS limitations: Results Final Parameters Final Hydrograph Validation Comparison between Quasi-Distributed and Lumped Model Hydrological Modelling with MikeSHE Purpose of the model Model Construction Steps Validation / Results
3 8.4 Remarks / Problems Hydraulic Modelling input Hydrograph References List of Figures Figure 1 : Thiessen polygons method, apply to the case of study Figure 2 : Comparison of Thiessen method efficiency, depending of sampling density Figure 3 : Comparison between Thiessen and BLUE s method Figure 4 : Rainfall distribution with Spline method. Figure 5 : Rainfall distribution with IDW method Figure 6 : Rainfall distribution with Kriging method Figure 7 : Rainfall distribution with IDW and Kriging method, made in ArcGIS Figure 8 : ArcGIS print screen of zonal statistic, of IDW rainfall distribution in all the subbasins Figure 9 : Graph representing the rainfall distribution with Kriging and Thiessen method Figure 10 : Influence of the rainfall distribution method on the flood peak Figure 11 : Land use repartition, from ArcGIS Figure 12 : Flood hydrograph for the case Figure 13 : Flood hydrograph for the case Figure 14 : Influence of time of concentration methods Figure 15 : Graph showing the comparison of Socose and SCS lag time Figure 16: Typical representation of watershed runoff Figure 17: Models used in modeling the Var Runoff at the napoleon bridge Figure 18: Var Basin Model Figure 19: Meteorologic Model Figure 20: Gage weights input screen Figure 21: Time-Series Data Figure 22: Estimated hydrograph at Napoleon Bridge Figure 23: Input precipitations Figure 24: Final Hydrograph at Napoleon Bridge Figure 25: Final Hydrograph Figure 26: Daily discharges at the napoleon bridge in Figure 27: Daily discharges at Napoleon Bridge in January Figure 28: Comparison between computed and observed hydrograph in January Figure 29: Comparison of quasi-distributed and fully lumped model hydrograph Figure 30:Typical presentation of the processes in MikeSHE Figure 31: Methods used in MikeSHE model Figure 32: Model Domain and Grid (up) Figure 33: Topography of 300m (left) and 75m grid (right) Figure 34: Thiessen polygons and sample precipitation rate time series
4 Figure 35: River network and domain Figure 36: Landuse and assigned Strickler coefficients Figure 37: Output Hydrograph obtained by MikeSHE models Figure 38: Output Hydrographs obtained from all Hydrological models
5 List of Tables Table 1: Rain gauge location in Lambert II extend Table 2: Thiessen weights for each rain gauge depending on their location Table 3: Rainfall distribution for each sub-basin depending on the Thiessen weights Table 4: Comparison between Kriging, IDW and Thiessen polygons method Table 5: Comparison between kriging and Thiessen rainfall distribution, with hourly rainfall Table 6: Repartition of land use in percentage, for each sub-basin Table 7: Impervious percentage depending on the land use Table 8: Impervious percentage of all the sub-basins Table 9: Slope and maximum flow length for each sub-basin Table 10: Results of different time of concentration methods Table 11: Sensitivity analysis for 6 different cases of CN values Table 12: Comparison of the influence of time of concentration methods Table 13: Influence of the lag time on the hydrograph Table 14: Curve Number for each sub catchments Table 15: Monthly average flows for each subcatchments Table 16: Final calibration parameters Table 17: Sensitivity Analysis Application on the Initial Abstraction parameter Table 18: New parameters based on mean values Table 19: Comparison of peak outflow, time of peak flow and volume between quasi-distributed and fully lumped model
6 1. Introduction The hydrological modeling analysis consists in building the flood hydrograph located at the Napoleon Bridge. In that optic, it s necessary to conduct a rainfall distribution analysis all over the Var catchment. The choice of the method depends primarily on the information available on the raingauge. In that sense, according to the number of raingauges covering the Var catchment (6 in total), we can affirm that this number should not allow interpolation methods application. To confirm this statement, we ll compare the Thiessen polygons method with other interpolation methods like Best Linear Estimators Unbiaised (SPLINE, Inverse Distance Weight and Kriging). Once the precipitations are spatially distributed, it is important to focus on the land use in order to better understand the occupation of soils in terms of percentage of imperviousness, slope and longest flow path. The loss parameters which are the time to peak and the time of concentration are inferred from this analysis for each sub-basin. A sensitivity analysis is then conducted using the Excel spreadsheet of Mr. Laborde which allows us to rebuild a flood effect by the method of unit hydrograph. This approach will permit to understand the influence of the rainfall runoff method and parameters of the loss method (CN). The relevance of time of concentration methods are next evaluated (Passini, Ventura and Nasch) by comparing the simulation results with estimated hydrograph. to compare the lag time. Finally we will see the stages of construction of the flood hydrograph in HEC-HMS and MIKE SHE. In the first place, HEC-HMS software will help us to build a quasi-distributed numerical model in order to predict runoff volumes, peak flows, and timing of flows by simulating the behaviors of the watersheds. After validating the model, a fully lumped model will be build with the purpose of comparing it with the quasi-distributed one. In the second place, another hydrograph will be also compute with MIKE SHE software in order to evaluate the interest of distributed model. 5
7 2. Rainfall distribution We choose to use the Thiessen polygons to make the rainfall distribution of the 5 sub catchments. We built our polygons using the software ArcGIS. At first we define our sub-basins in ArcGIS by merging different sub-basins according to drainage points. We then collected the rain gauges locations on website of METEO-FRANCE ( ). METEO-France is the French meteorological agency, a public administrative institution in charge of forecasting and the study of meteorological phenomena. The website allows us to access to different data about the instruments of hydrological measurements. We were therefore able to reap the Lambert II extended coordinates for each rain gauge. We have 6 rain gauges in the Var catchement: Lambert II extend Rain gauge X Y Z Carros Levens St Martin Vésubie Guillaumes Roquesteron Puget Theniers Table 1: Rain gauge location in Lambert II extend. The Thiessen method permit to converts input rain gauges points to an output feature class of Thiessen proximal polygons. These polygons have the unique property that each polygon contains only one input point, and any location within a polygon is closer to its associated point than to the point of any other polygon. 6
8 Figure 1 : Thiessen polygons method, apply to the case of study. So we could put the geo-referenced rain gauges on the basins that had the same projection systems. Thanks to ArcGIS we select the intersect area between the sub-catchments and the Thiessen Polygons. That permits to compute the Thiessen Weights: Up Var Low Var Tinée Vésubie Esteron Guillaumes Puget Roquesteron Levens Carros St Martin Vesubie Table 2: Thiessen weights for each rain gauge depending on their location. Once we have these coefficients, we can assign them to the rain gauges rainfall data and determine the sub basins rainfall. Below we can observe the rainfall for each sub-basin: Date Time Up Var Low Var Tinée Vésubie Esteron 3/11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/
9 3/11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ Total (mm) : Table 3: Rainfall distribution for each sub-basin depending on the Thiessen weights. 3. Comparison of rainfall distribution methods We then took the time to compare the pattern of rainfall distribution in different ways. Here we will try to compare Thiessen polygon method with Best Linears Unbiaised Estimators (Spline, Inverse Distance Weights and Kriging). 3.1 Lumped or quasi distributed method Thiessen Polygon Method Thiessen proposed a method to evaluate spatial averages over a domain from anecdotal. One way underlying the "method of Thiessen" is an interpolation technique based on the "law " of the nearest neighbor. Clearly, the main advantage of this method is its simplicity. There is no objective information on the representativeness of interpolations. As illustrated in the diagrams below where the problem is presented in one dimension, if the density of sampling points is very strong, representativeness may be about right (left picture below), in other cases Instead the results can be very disappointing (right picture below). Figure 2 : Comparison of Thiessen method efficiency, depending of sampling density. 8
10 Correlation coefficient 1,0 Gandin Krigeage 0,9 0,8 0,7 0,6 0,5 0,4 Spline Ariytmétic mean Of 5neighbour Thiessen 0,3 0% y 50% 100% Pourcentage reconstructed rainfall periode Figure 3 : Comparison between Thiessen and BLUE s method. Thiessen polygons assign a weight to the gauges and then distributes in a homogeneous distribution on the sub-basins. However, we can consider this traditional method of interpolation as quasi Distributed. Indeed we allocate different weights to rain gauges and we do not consider an equal distribution of rainfall across our main basin. Full lumped model would consider a uniform distribution on all the Var basin. By going a little deeper, we can say that lumped model for each sub catchment basins would consider that with no rain sub-basins would not undergo rainfall. 3.2 Distributed Method (Best Linears Unbiaised Estimators) We have chosen here to compare three methods to give us a distribution of rainfall much more accurate and with a better distribution. We will detail three methods of interpolation's BLUE (Best Linear Estimators Unbiaised). We will begin with the Spline method (minimum curvature), then we will focus on the method of inverse distance weight (IDW) and we will then finish with the Kriging method Inverse Distance Weights This is a deterministic interpolation techniques that create surfaces from measured points, based on either the extent of similarit. The inverse distance weight technique provides an interpolation using a linearly we combination of the rain gauge locations. To predict a value for any unmeasured location, IDW will use the measured values surrounding the prediction location. Those measured values closest to the prediction location will have more influence on the predicted value than those farther away. Thus, IDW assumes that each measured point has a local influence that diminishes with distance Spline The basic form of the minimum curvature Spline interpolation imposes the following two conditions on the interpolant: The surface must pass exactly through the data points. The surface must have minimum curvature. 9
11 The basic minimum curvature technique is also referred to as thin plate interpolation. It ensures a smooth (continuous and differentiable) surface, together with continuous firstderivative surfaces. Rapid changes in gradient or slope (the first derivative) can occur in the vicinity of the data points; hence, this model is not suitable for estimating second derivative (curvature) Kriging This is a geostatistical interpolation techniques using the statistical properties of the measured points. Geostatistical techniques quantify the spatial autocorrelation among measured points and account for the spatial configuration of the sample points around the prediction location. Kriging generates an estimated surface from the points introduced (rain stations). Kriging is a spatial interpolation method, sometimes considered the most accurate a statistical point of view, which allows a linear estimate based on the expectation and also on the variance of the spatial data. As such, kriging is based on the calculation, interpretation and modeling of the variogram, which is a measure of the variance depending on the distance between data. This interpolation method differs from other methods (inverse distance, Thiessen polygons) because it has two advantages. First, kriging is best linear unbiased: the means are identical and the minimum variance. We built the three rainfall distribution in SURFER 8 software: Figure 4 : Rainfall distribution with Spline method. Figure 5 : Rainfall distribution with IDW method. 10
12 Figure 6 : Rainfall distribution with Kriging method. We have built spherical variogram method (scale = 375, length = 15100) and we also wanted to reflect uncertainty measures by imposing a Nugget effect (N = 10). But we met problems when we had to build our variograms. Indeed it is difficult to establish a regression with so few data, (whether in a linear exponential or spherical for the variograms components). Finally we decide to use ArcGIS software that allows us, to self calibrate the correlation and so variogram could be more accurate. So we rebuilt our distribution of rainfall in ArcGIS. But we rejected the Spline method. This method does not give us smooth. This method gives us a massive distribution of rainfall spread over Tinée while the sub basins do not have a rain gauge. Figure 7 : Rainfall distribution with IDW and Kriging method, made in ArcGIS. We use in ArcGIS the interpolation to raster from spatial analysis. We the use IDW and Kriging methods with a spherical variogram for Krigging and we choose a inverse distance squared weighted interpolation. We then apply a zonal statistic to know the rainfall distribution in all the sub-basins. 11
13 Figure 8 : ArcGIS print screen of zonal statistic, of IDW rainfall distribution in all the sub-basins. We initially wanted to compare the sums of hourly rainfall (from , 12:00 to 11/06/1994 1:00 am) on each sub-basins. Comparisons of the Thiessen polygon results with IDW method and Kriging method could then be carried out, which are shown in Table 2 below: Comparison: Kriging (Spherical) vs. IDW vs. Thiessen Polygons Kriging IDW (Power = 2) Thiessen Polygons Catchments mean (mm) Vesubie Esteron Lower Var Tinee Upper Var VAR Table 4: Comparison between Kriging, IDW and Thiessen polygons method. Then the hourly rainfall for the method of Thiessen polygons and krigging were compared. See below the summary table giving rainfall in mm for each sub-basin: Thiessen Kriging Thiessen Kriging Thiessen Kriging Thiessen Kriging Thiessen Kriging Date Time Up Var Low Var Tinée Vésubie Esteron 3/11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/ /11/
14 5/11/ /11/ /11/ /11/ /11/ /11/ /11/ Total (mm) : Table 5: Comparison between kriging and Thiessen rainfall distribution, with hourly rainfall. We can now compare the rainfall on each sub-basin by editing graphics. We note the similarity of rain on the various sub-basins. See below the rain by the two methods considered in the basin Up Var: Figure 9 : Graph representing the rainfall distribution with Kriging and Thiessen method. In order to confirm our observations, we compare the peaks of flood generated by these different methods. For this, we used our model calibrated for HEC-HMS. We have not changed any other settings, on the purpose to only observe the influence of the rainfall distribution method. 13
15 Figure 10 : Influence of the rainfall distribution method on the flood peak. 3.3 Conclusion The choice of interpolation method depends primarily on the information available. One can distinguish between Climate methods (where there are several observations of the field, or using purely spatial (this is our case). In our case the number of rain gauges is too small for the numerical values to be reliable. Indeed we note that there are not enough points to be able to assimilate the information for the variogram. Consequently whe the variogram was set up to perform a method of Best Linear Estimators Unbiaised (IDW, spline, Krigging) the gridding was not accurate and we will lose here all the interest of the correlation. It would be very interesting to have more data and to make a correlation between rainfall and the distance to the sea and altitude. It would require us to apply the exponential law for the distance to the sea and make a log normal regression to the altitude. This work could be done with SURFER software. By observing the three different methods we can notice that the results are almost similar. The impact of the distribution of rainfall over the flood peak is negligible. We can therefore conclude that in our case study the rainfall distribution method should be Thiessen polygons. This method is faster to implement and gives us more extreme rainfall, which may be interesting to take into account to model the event having the greatest impact. 14
16 4. Land Use Analysis We further discussed the land use. The land use will allow us to know the distribution of soil occupation characteristics. We therefore use the text file "land use" in ArcGIS. We then grouped the different land uses in 5 parts: Artificial surfaces, Agricultural Areas, Forest and semi natural areas, Wetlands, Water Bodies. Figure 11 : Land use repartition, from ArcGIS. Knowing the areas of all sub-basins and the occupied areas for each land use we have established a table: Esteron Low Var Up Var Tinée Vésubie Areas (km²) Artificial surfaces (%) Agricultures areas (%) Forest & Semi natural areas (%) Wetlands (%) Water bodies (%) Table 6: Repartition of land use in percentage, for each sub-basin. The table allowed us to identify our coefficients for the loss, SCS Curve Number, method. This table was also used to determine the coefficient of soil permeability. Here we have assigned a coefficient of impermeability to each land use (Source: U.S. Soil Conservation Service): 15
17 % Impervious = f (land use) Artificial surfaces 0.78 Agricultures areas 0.05 Forest & Semi natural areas 0.05 Wetlands 1 Water bodies 1 Table 7: Impervious percentage depending on the land use. Knowing what percentage of each land use for each sub basins, we have determined the permeability coefficients for each of our five sub-basins. % imperviousness Esteron Low Var Up Var Tinée Vésubie Table 8: Impervious percentage of all the sub-basins. To calculate the parameters of our loss and transform methods we need to determine in ArcGIS the slopes and the maximum flow length for each sub-basin. These parameters will allow us to calculate our time of concentration and rise times by different methods. 4.1 Lag Time & Time of Concentration The rise time depends on the concentration time of our basins. Time of concentration on a watershed is defined as the maximum required to a drop of water to traverse the path between hydrological points of the basin to the outfall. Area (km²) Slope (m/m) Maximum flow length (km) Esteron Low Var Up Var Tinée Vésubie Table 9: Slope and maximum flow length for each sub-basin. We after choose different methods for the time of concentration, and also for the lag time. Tc (hr) Kirpich Ventura Passini Nasch Esteron Low Var Up Var Tinée Vésubie Table 10: Results of different time of concentration methods. 16
18 We immediately dismiss the formula Kirpich. Indeed, this method is only applicable for basins with areas less than 0.8km ². The equation below concern Lag Time methods: Triangular Unit Hydrograph: We can also use the Nash lag time: We then conduct a sensitivity analysis using the Excel file of Mr. Laborde. This file allows us to rebuild a flood effect by the method of unit hydrograph.we decide to use the SCS Curve number, loss merhod. This will therefore enable us to understand the influence of the Curve Number on our flood peak but also understand the influence of our time of concentration and our lag time. We will then detail the flood hydrograph build the SCS CN method with a time of concentration corresponding to Ventura formula, and a lag time corresponding to the method of SCS (triangular hydrologic units),. Once we assess the CN we change our time of concentration by comparing with Ventura and Nasch Passini. Finally we use the lag time of the Nasch formula. Once all our parameters are clearly understood and defined by the senility analysis we will implement them in HEC HMS. Using the excel sheet we were able to simulate for each sub-basin, six cases, allowing us to understand the influence of CN on our hydrographs. Below we can see the results of different scenarios without changing the curve number. SCS Curve Number Difference with observed hydrograph Up Var Low Var Tinée Vesubie Esteron Time (hours) Discharge (%) Volume (%) Case Case Case Case Case Case Table 11: Sensitivity analysis for 6 different cases of CN values. 17
19 Figure 12 : Flood hydrograph for the case 1. Figure 13 : Flood hydrograph for the case 6. We chose to keep the parameters of case study 6. Indeed, this scenario is one that allows us to get as close as the volume observed. It should be noted that we are aware that the observed flood hydrograph is interpolated. During the flood of 1994 measuring devices was destroyed above 600 m 3 /s. So we will try to get closer to the curve on the peak time, which is the only parameter known with certainty. We will then in HMS HEC try to get closer to the volume and flood peak. Subsequently we will modify the time of concentration using different methods (Ventura, Passini, Nasch). Difference with the observed hydrograph Ventura Passini Nasch Time (hours) Discharge (%) Time (hours) Discharge (%) Time (hours) Discharge (%) Cas Table 12: Comparison of the influence of time of concentration methods. 18
20 Figure 14 : Influence of time of concentration methods. By comparing the concentration time, we note that this setting has no effect on volumes. However we note that the peak discharge and peak time varies. The formula for time of concentration of Ventura seems most appropriate here. Indeed it remains close to the curve observed at the peak time and at the maximum rate as well. Knowing our CN and therefore our method for determining the concentration time, they we have to understand the influence of lag time. To do this we will compare two methods. The first method is the rise time of the SCS triangular unit hydrograph, the latter being the method of SOCOSE. Difference between the observed hydrograph Time Discharge (%) (hours) Socose SCS Triangular UH Table 13: Influence of the lag time on the hydrograph. 19
21 Figure 15 : Graph showing the comparison of Socose and SCS lag time. We note that SOCOSE method allows us to obtain a flood peak very similar to the observed one, either in time or flow. To conclude this sensitivity analysis allowed us, to understand the various parameters involved in the construction of the hydrograph. We will use later in HEC-HMS, the CN determined, the concentration time of Ventura and SOCOSE Lag Time. By inputting these parameters at the beginning of our set up we will make a more accurate calibration and it will save time. We'll also be able to apply different methods in HEC-HMS that either loss or transform. The software will also allow us to integrate the routing settings and more. 20
22 5 Hydrological modeling with HEC-HMS Built by the US Army Corps of Engineers, HEC-HMS is a numerical model which can help to predict runoff volumes, peak flows, and timing of flows by simulating the behaviors of the watersheds, channels, and reservoirs. 5.1 Why a model? We propose in our project to use HEC-HMS in order to compute the runoff volume drained to the sea all over the Var catchment. This software use a simple view of the hydrologic process as illustrated in the figure 16 took from the technical reference manual. Each components of the runoff process are integrated in the software as models: Runoff-volume models that compute effective precipitations; Direct-runoff Models corresponding to water that moves over or beneath the soil; Baseflow Models that calculate the subsoil water volume drained into the channel; Routing Models that take into account storage and energy flux through the stream channels. The following illustration gives the final methods used in the project: Figure 16: Typical representation of watershed runoff. Rainfall distribution: THIESSEN Loss model : SCS Curve Number Direct runoff: Clark Unit Hydrograph Observed Hydrograph Channel flow routing : Muskingum Figure 17: Models used in modeling the Var Runoff at the napoleon bridge. 5.2 Model Construction Steps The model has been built following the steps below: 21
23 Enter basin, time-series and meteorological model data; Enter parameters values for loss, transform, baseflow and routing model; Enter control specifications; Run the model. The few following screenshots illustrate the precedent steps: Figure 18: Var Basin Model. Figure 19: Meteorologic Model. The Thiessen method involves choosing first for each subcatment the raingauges concerned and next to input depth weights (%) for each of them. Figure 20: Gage weights input screen. Precipitation gages correspond to the hourly rainfalls from November 3rd h to November 6 th 1994 at midnight. We also added the napoleon bridge estimated hydrograph. Figure 21: Time-Series Data. 22
24 Figure 23: Input precipitations. Figure 22: Estimated hydrograph at Napoleon Bridge. 5.3 Transform Method The transform parameters that have to be assumed in Clark Transform method are time of concentration (Tc) and the storage coefficient (R). Respectively, these 2 parameters represent two critical processes in the transformation of excess precipitation to runoff: Translation or movement of the excess precipitation from the watershed origin throughout the drainage to the outlet. Attenuation or reduction of the discharge peak as the excess is stored throughout the watershed [3]. The Ventura method has been chosen determing the time of concentration Tc regarding the peak shape with the Jean Pierre Laborde spreadsheet. The basin storage coefficient R, is an index of the temporary storage of excess precipitation in the watershed as it drains to the outlet point. Clark (1945) indicated that R can be computed as a flow at the inflection point on the falling limb of the hydrograph divided by the time derivative of flow (slope at that point). After several tests, this method of computation don t give relevant results. We can just assume that the storage coefficient varies from 0.1 to 0.9, depending upon the characteristics of the area of study. An average value of 0.05hours is for sandy geographies while a 0.95hours is for impervious areas. In our case, according to the land use analysis, the Var watershed is mostly made up of forest and semi natural areas, and overall charactherized by karst. Hence, the initial value is assumed to be 0.3. Besides, because of the lack of information relative to the volume of water drained by karst aquifers, the representation for this coefficient is not relevant and it could have an impact on the model. 23
25 For this reason, we decide to compare this method with another lumped model, the SCS Unit Hydrograph which depends only on the time to peak (SOCOSE Lag Time). This transform method revealed relevant results (see final calibration 6.2) conforting our final decision to use it for our hydrological modeling project. 5.4 Routing Method The Muskingum routing method uses a simple conservation of mass approach to route flow through the stream reach. By the continuity equation: time. Where I is inflow, O is the outflow rate, S is storage and t the Storage is modeled in the reach as the sum of prism storage (volume defined by a steady-flow water surface profile) and wedge storage (additional volume due to flood wave) multiplied by the travel time through the reach, K. This lead to the equation of the storage S t : Where the required parameters are K, the travel time of the flood wave through routing reach and X, a dimensionless weight ranging from 0 to 0.5. In order to apply a rational Muskingum model, we must also take into account the computational time step Δt respecting the ratio 0 < Δt / K < 2. Computing with an hourly time step, we can affirm that the travel time can t go below 0.5. Regarding the X parameter, experience has shown that for steeper streams, it will be closer to 0.5 whereas for channels with mild slopes and over-bank flow, it will decrease gradually to 0. Hence, the dimensionless weight X will be close to 0,5 in the upper part of the Var watershed and around 0 for the Lower Var reach. 5.5 Loss Method The Soil Conservation Service Curve Number model estimates precipitation excess as a function of cumulative precipitation, soil cover, land use and antecedent moisture by using equation: Where: P e = accumulated precipitation access at time t; P = accumulated rainfall depth at time t; S = potential maximum retention; I a = initial loss = 0,2*S (from SCS) The maximum retention, S is related with watershed characteristics through SCS Curve Number (CN) as: 24
26 CN for a specific watershed can be estimated using tables published in the Technical Reference Manual of HEC-HMS, appendix A. The soil type estimated range to moderate to high runoff corresponding to soil group C. For the case of the Var catchment which consists of several soil types and land uses, a composite CN was calculated as: in which CN composite is the CN used for runoff volume computations; i = an index of watersheds subdivisions of uniform land use and soil type; CN i = the CN for subdivision I, and A i = the drainage area of subdivision i. Esteron Low Var Var Up Var Tinée Vésubie CN Artificial Surfaces (km²) 1,80 15,81 2,40 1,80 2,36 90 Agricultures areas (km²) 34,28 33,44 36,30 5,83 13,00 82 Forest & Semi natural areas (km²) 414,83 99, ,67 740,00 378,56 77 Wetlands (km²) 0,09 0,00 1,64 0,00 0, Water Bodies (km²) 0,00 3,34 0,00 0,37 0, Sum Areas (km²) 451,00 152, ,00 748,00 394,00 CN composite computed Table 14: Curve Number for each sub catchments. Looking to the table 14, the CN composite will be most represented by forest and semi natural areas corresponding to CN = 77. Coupling these results with the intermediate calibration using the spreadsheet and the observed hydrograph, we finally chose the values shown in the final parameters table Baseflow Method Like the direct runoffs of precipitation, the baseflow is also a component of the streamflow hydrograph. Baseflow is the sustained or fair weather runoff of prior precipitation that was stored temporarily in the watershed, plus the delayed subsurface runoff from the current storm. Three alternative methods are proposed in the model: Constant, Monthly-Varying Baseflow; Exponential Recession Model; Linear-reservoir volume accounting model. We chose in our model the simplest one, constant monthly, because of the lack of data required for the other methods. However, it is important to notice that this baseflow method is intended normally to continuous simulation. Hence, it could be interesting to implement the linear-reservoir baseflow, which is more appropriate to model the recession after a storm event. We used constant daily flows over 36 years (1974 to 2000) available in the excel sheet var-discharge and compute the monthly average discharges for each subcatchments based on area weighting. The table below presents the final discharges: January February March April May June July August September October November December Area Monthly average flows 48,06 38,71 43,56 60,41 73,87 57,46 33,75 24,96 32,67 64,40 61,30 44,80 Upper Var 18,48 14,88 16,75 23,23 28,40 22,09 12,98 9,60 12,56 24,76 23,57 17, ,00 Vesubie 6,68 5,38 6,05 8,40 10,27 7,99 4,69 3,47 4,54 8,95 8,52 6,23 394,00 Tinee 12,68 10,21 11,49 15,94 19,49 15,16 8,91 6,59 8,62 16,99 16,17 11,82 748,00 Esteron 7,64 6,16 6,93 9,61 11,75 9,14 5,37 3,97 5,20 10,24 9,75 7,13 451,00 Lower Var 2,58 2,08 2,34 3,24 3,96 3,08 1,81 1,34 1,75 3,45 3,29 2,40 152,00 Table 15: Monthly average flows for each subcatchments. 25
27 5.7 HEC HMS limitations: The use of this software must be done keeping in mind the limitations. Indeed, the fluctuations of the tide, dams, bridges, and the existence of lateral flow phenomena can create backwaters that reduce and shift over time the peak flood. In that sense, we are aware that being based on an approximation of uniform flow, patterns of kinematic wave and Muskingum cannot model their effects on the flow. In our case of study, the flood overflowing the bed of the river, the wave will flood the surrounding plains. In some cases, these plains will thus constitute storage basins. The model used must take into account the possible passage of the flow of the riverbed floodplain. Often in these cases, we calculate the hydraulic properties of the riverbed and surrounding plains separately, and then combine them with a composite formula. We are conscious that this cannot be achieved by models Kinematic wave or Muskingum, because the respective parameters they use are constant and cannot take into account a change of hydraulic properties of flow. 6 Results 6.1 Final Parameters The table below presents the final parameters implemented in HEC-HMS after calibration. Hydrologic Elements Upper Var Tinee Vesubie Esteron Lower Var Reach1 Reach2 Reach3 Runoff volume Parameters Direct Runoff Parameters Channel Flow Parameters SCS Curve Number Clark Unit Hydrograph SCS Unit Hydrograph SCS CN Percent Initial Muskingum K Muskingum Storage Impervious Abstraction Tc (hr) Lag Time (min) (hr) X Coeff. (hr) (%) (mm) 77 5,3 15,17 28,934 0, , ,22 16,93 21,459 0, , ,46 21,77 13,323 0, , ,31 16,93 19,364 0, , ,6 12,70 19,885 0, , ,9 0, ,7 0, ,2 0,1 Table 16: Final calibration parameters. In order to better understanding the influence of each parameter on the final hydrograph, a sensitivity analysis has been done. The main purpose consists in computing the percentage difference between computing values with initial values in terms of peak discharge, time of peak and volume of water reaching the outlet of the Var basin. The table below gives the example for the initial abstraction influence: No. of Run Input Parameters Percent of change 1 initial value 2 5% 3 10% 4 30% 5 60% Initial Abstaction, Ia (mm) 6-5% 7-10% 8-30% 9-60% Calculated Values Percentage Difference Between Calculated Value to the initial Value Peak Discharge (m3/s) Time of Peak mm/dd/yr Volume (MM) Peak Discharge (%) Time of Peak (hr) Volume (%) 3215,7 11/05/1994, 16:00 51, ,1 11/05/1994, 16:00 50,67-0,89 0-1, ,7 11/05/1994, 16:00 49,85-1,80-0,33-3, ,6 11/05/1994, 16:00 46,62-5,63-0,33-9, ,3 11/05/1994, 17:00 41,93-11,95-0,33-18, ,8 11/05/1994, 16:00 52,32 0,87 0 1, ,3 11/05/1994, 16:00 53,15 1,73 0 3, ,2 11/05/1994, 16:00 56,52 4,99 0 9, /05/1994, 16:00 61,7 9,43 0,17 19,83 Table 17: Sensitivity Analysis Application on the Initial Abstraction parameter 26
28 6.2 Final Hydrograph 4000 Final Calibration Observed hydrograph (cm/s) Clark Unit Hydrograph (cm/s) SCS Unit Hydrograph (cm/s) Peak Discharge 05/11/ :00; 3680 cm/s 05/11/ :00; 3126 cm/s 05/11/ :00; 3586 cm/s Figure 24: Final Hydrograph at Napoleon Bridge The hydrograph built with the "SCS Unit Hydrograph" method will be finally the final one used for river modeling because of its better representation of reality. The main argument depends more on the required parameter than on the hydrograph shape as the observed one was estimated after 600cm/s. The final hydrograph, located at the junction 3, is presented below: Figure 25: Final Hydrograph. 27
29 6.3 Validation Model validation is the process of determining if a proposed model accurately represents a physical process and provides predictive capability [5]. In other words, the calibrated model built for the Var river project must be able to reproduce field observations from an independent data set. In that optic, a new model has been built using daily precipitation datas and observed hydrograph at the napoleon bridge from January 1 st to January 31 st The choice of January month is mostly explained by the fact that it constitutes the second most important peak flow in Concerning the other parameters values, they remain exactly the same. The table and graph below describe in detail the event [6] : Figure 26: Daily discharges at the napoleon bridge in 1994 Figure 27: Daily discharges at Napoleon Bridge in January 1994 The first observation that can be done concerns the 1994 flood in the figure 26 that doesn t represent the estimate peak flow as it stops at 1450cm/s. The figure below presents the final result after computation: Model Hydrograph (cm/s) Observed flow (cm/s) Figure 28: Comparison between computed and observed hydrograph in January
30 It appears clearly that the model built for the 1994 flood event doesn t match with this event. Indeed, it s difficult to have a relevant comparison knowing that we built the model of November 1994 based on extreme discharges and hourly whereas the observed flow corresponds to mean daily values. However, knowing that the maximum discharge for January 1994 was on 7 and reach 884cm/s, we can affirm that the trend of the curve is right and hence that the model should give relevant results for an episode closed to November 1994 one. In that sens, it could be interesting to validate the model on another event but based this time on extreme values. 7 Comparison between Quasi-Distributed and Lumped Model In this final part, we propose to compare the results of the final model (quasi-distributed) with a fully lumped model based on mean values for each subcatchment and parameter. The average values are presented in the table 17: Hydrologic Elements Upper Var Tinee Vesubie Esteron Lower Var Fully Lumped Model Reach1 Reach2 Reach3 Fully Lumped Model Runoff volume Parameters Direct Runoff Parameters Channel Flow SCS Curve Number Clark Unit Hydrograph Unit Hydrog SCS CN Percent Initial Muskingum Muskingum Storage Lag Time Impervious Abstraction Tc (hr) K (hr) X Coeff. (hr) (min) (%) (mm) 77 5,3 15,17 28,934 0, , ,22 16,93 21,459 0, , ,46 21,77 13,323 0, , ,31 16,93 19,364 0, , ,6 12,70 19,885 0, , ,40 7,18 16,70 20,59 0,50 411, ,9 0, ,7 0, ,2 0, ,9 0,3 Table 18: New parameters based on mean values. Regarding precipitation, mean hourly values have been computed gathering all the raingauge giving at the end only one rainfall event for the entire basin. After running the model, the results are the following: Observed hydrograph (cm/s) SCS Unit Hydrograph (cm/s) SCS fully Lumped (cm/s) Figure 29: Comparison of quasi-distributed and fully lumped model hydrograph. 29
31 SCS Unit Hydrograph (cm/s) SCS fully Lumped (cm/s) Observed hydrograph (cm/s) Total Residual SCS U H Total Residual SCS fully Lumped Peak Outflow (cm/s) 3585,8 3411, ,20 268,60 Time of peak flow (h) 5/11/94 21:00 5/11/94 21:00 5/11/94 18:00 3h 3h Volume (mm) 94,43 67,61 91,56-2,87 23,95 Table 19: Comparison of peak outflow, time of peak flow and volume between quasi-distributed and fully lumped model. Finally, the results show differences between the two models but not alarming. Indeed, the gap in peak outflow is cms below the estimated hydrograph peak against 94.2 cms for the quasi-distributed model. The volume decreases consequently to 24mm whereas the time of peak flow stands quite similar. To conclude, this study shows that using mean value give also relevant results. However, the choice depends more on the level of understanding of both hydraulics and hydrological processes, and also numerical modeling of physical processes by the user than the project purpose in itself. Indeed, even if we have switched for mean values, a detailed analysis of each parameter is essential in order to catch up with the reality. So for example, we cannot affirm after this study that a fully lumped model will not be adapted to modeling a retention basin as the volume is far from the reality. In our case, in view of the fact that the Var catchment brings complex processes, it is advisable to use a quasi-distributed model. Concerning the question of distributed model, it will be conducted in the next part related to Mike She. 30
32 8 Hydrological Modelling with MikeSHE MikeSHE is a complete software package by DHI, used as a modelling tool for a wide range of water resources and environmental problems related to surface water and groundwater and can be applied on scales ranging from local infiltration studies to regional watershed studies. It is a dynamic, and simulates all the major processes of the hydrological cycle. 8.1 Purpose of the model After computing the runoff volume of the flood event in the Var catchment with HEC-HMS, we chose to model and compare the output hydrograph with a fully distributed model that could be provided by MikeSHE. This software uses a more or less, according to the needs of the model, complex view of the hydrological processes, as illustrated in figure 30. Flow Model: computes the runoff by simulating all the processes involving water movement: Overland flow Channel flow Evapotranspiration Unsaturated Flow Saturated Flow Water Balance Calculation: computes the storage depths, changes and model errors for individual model components. Figure 30:Typical presentation of the processes in MikeSHE Overland Flow: Finite Difference Rivers and lakes (OC) Flow Model Rainfall distribution: THIESSEN Observed Hydrograph Figure 31: Methods used in MikeSHE model For our project we created a flow model using the Overland flow (Finite Difference method) and the Channel flow (Rivers and Lakes). The choice of these water movements was based on the assumption that after the 3 day event that proceeded the flood event, there is only saturated soil and there were no data provided that could indicate the opposite or assist in the calculation of saturated or unsaturated flow. Moreover, in this way we consider the full effect of rainfall to the flood event and model an extreme event as such of Especially for the overland flow, the finite difference was preferred over the Subcatchment Based since our interest was focused on the entire basin and particularly the lower Var. 31
33 An overview of the model and methods are illustrated in figure Model Construction Steps The model was built following the steps below and in this report is presented with the final parameters used: Fill simulation specifications; Choose water movement; Input grid, topography, precipitation rate, landuse; Input river network from Mike11 Determine the overland flow parameters and HD parameters; Run the model Simulation specification The simulation parameters are illustrated in the table below. Simulation Period 04/11/ :00 06/11/ :00 Time Step Control Initial time step 0.01 hours (6 minutes) Max allowed OL time step 0.1 hours Max precipitation depth per time step 0.5 mm Max infiltration amount per time step 1 mm Input precipitation rate requiring its own 0.5 mm time step OL Computational Control Parameters Maximum number of iterations 20 Maximum head change per iteration m Threshold gradient for applying lowgradient flow reduction Manning equation (using OL flow (checked) manning numbers) HD Parameters Upstream boundary condition 0 Initial conditions water depth 0.5m discharge 10m 3 /s Bed resistance 30 (uniform for the whole river) Table 19: Simulation specification Model Domain and Grid, Topography Both of the files used were derived in GIS from the ascii file (topography) and the hydro-analysis using ArcHydro Tools (domain). All files were re-projected in the same coordinate system (France III) and then were input to the model. The topography was converted to dfs file, in order to ensure speed of the model. 32
34 Figure 32: Model Domain and Grid (up).. Figure 33: Topography of 300m (left) and 75m grid (right) Climate For the Climate, the station based precipitation rate was used to ensure the spatial distribution in the basin, according to Thiessen polygons. The shape file with the polygons that was created in ArcGIS was imported and the time series of the six rain gages were edited in MikeSHE. The software calculates the Thiessen polygons without further interaction of the user. Figure 35: Thiessen polygons and sample precipitation rate time series 33
35 The net rainfall fraction was considered uniform and equal to 0.9, considering a very small part of losses since there was no infiltration parameter set. In fact, the sensitivity analysis provided a very small difference in the result using 0% or 10% infiltration, but with the 10% to reach the observed hydrograph better Rivers and Lakes The Mike11 files of the network and cross sections were used as input. The projection had to be changed to the system our data was projected in. A simulation file was created in Mike11, with the network, cross sections and the coupling with MikeSHE shifted to the right location. Upstream boundary condition (inflow) was considered to avoid numerical errors, since the river is considered to be fed only by the rainfall. Figure 35: 36: River network and domain Overland flow In the overland flow, no separated flow areas were used, as the only spatial distributed information we had was the Strickler coefficient through the landuse information. For the Strickler coefficient, the 34
36 values were obtained from the mean values of empirical tables for the manning coefficient and adjusted to the 4 types that occur in Var basin. Figure 37: Landuse and assigned Strickler coefficients For the Detention Storage and Initial Water Depth uniform values were used, of zero and 0.5m (as mentioned earlier), respectively. It has to be mentioned, that for the river bed, we used the maximum Strickler coefficient that can be used in the software (30 M). We would like to be able to use different coefficients for different parts of the river, using at least the only option given to divide the river in 3 parts, but during calibration, the results were not good. Therefore we used the uniform option of 30M, to approximate the mean value throughout the river Storing of results The desired results can be viewed in MikeSHE or MikeView (especially for the hydrograph). In order to indicate the outlet, points located near the mouth were put with an interval of 300m and 75m for the respective grids. 8.3 Validation / Results The sensitivity analysis, using logical ranges of the variables, did not provide much diversity. The best used variables were mentioned above. The validation of the 75m model originally would be done with the 300m one, but for reason explained in the Problem section, the 300m model was validated with the 300m model obtained by the data on HydroEurope archive (referred to as MikeSHE Tutorial). The main difference between the two MikeSHE models is the use of distributed Strickler according to the landuse, instead of uniform (Tutorial) and the initial conditions that were other than zero in our model. The simulation of flood event of 1994, compared to the observed hydrograph at Napoleon III Bridge and the Tutorial is illustrated in the following figure. 35
HYDRO EUROPE ISIS Report. Team
HYDRO EUROPE 2011 ISIS Report Team 6 23.02.2011 Team members: Liza ASHTON, Thibault DESPLANQUES, Jens Harold DRASER, Manuel GOMEZ (supervisor), Kyung Tae LEE, Maxime NARDINI, Jieun PARK, Jiyoung PARK &
More information1 n. Flow direction Raster DEM. Spatial analyst slope DEM (%) slope DEM / 100 (actual slope) Flow accumulation
1 v= R S n 2/3 1/2 DEM Flow direction Raster Spatial analyst slope DEM (%) Flow accumulation slope DEM / 100 (actual slope) 0 = no cell contributing 215 = 215 cell contributing towards that cell sqrt (actual
More informationWhat is runoff? Runoff. Runoff is often defined as the portion of rainfall, that runs over and under the soil surface toward the stream
What is runoff? Runoff Runoff is often defined as the portion of rainfall, that runs over and under the soil surface toward the stream 1 COMPONENTS OF Runoff or STREAM FLOW 2 Cont. The types of runoff
More informationURBAN FLOODING: HEC-HMS
1.0 Introduction URBAN FLOODING: HEC-HMS -Sunil Kumar, Director, NWA All major ancient civilisations were developed in the river valleys because river served as source of water, food, transportation and
More informationRainfall - runoff: Unit Hydrograph. Manuel Gómez Valentín E.T.S. Ing. Caminos, Canales y Puertos de Barcelona
Rainfall - runoff: Unit Hydrograph Manuel Gómez Valentín E.T.S. ng. Caminos, Canales y Puertos de Barcelona Options in many commercial codes,, HMS and others HMS Menu Transform method, User specified,
More information2
1 2 3 4 5 6 The program is designed for surface water hydrology simulation. It includes components for representing precipitation, evaporation, and snowmelt; the atmospheric conditions over a watershed.
More informationHydrologic Engineering Center Hydrologic Modeling System (HEC-HMS) Sunil KUMAR Director, National Water Academy
Hydrologic Engineering Center Hydrologic Modeling System (HEC-HMS) Sunil KUMAR Director, National Water Academy 22 April 2015 NWA, Pune Exercise Objective: To determine hydrological Response of the given
More informationRAINFALL - RUNOFF MODELING IN AN EXPERIMENTAL WATERSHED IN GREECE
Proceedings of the 14 th International Conference on Environmental Science and Technology Rhodes, Greece, 3-5 September 2015 RAINFALL - RUNOFF MODELING IN AN EXPERIMENTAL WATERSHED IN GREECE KOTSIFAKIS
More informationThe Texas A&M University and U.S. Bureau of Reclamation Hydrologic Modeling Inventory (HMI) Questionnaire
The Texas A&M University and U.S. Bureau of Reclamation Hydrologic Modeling Inventory (HMI) Questionnaire May 4, 2010 Name of Model, Date, Version Number Dynamic Watershed Simulation Model (DWSM) 2002
More informationRAINFALL-RUNOFF STUDY FOR SINGAPORE RIVER CATCHMENT
10 th International Conference on Hydroinformatics HIC 2012, Hamburg, GERMANY RAINFALL-RUNOFF STUDY FOR SINGAPORE RIVER CATCHMENT CHI DUNG DOAN (1)(3), JIANDONG LIU (1), SHIE-YUI LIONG (1), ADRI VERWEY
More informationDefinitions 3/16/2010. GG22A: GEOSPHERE & HYDROSPHERE Hydrology
GG22A: GEOSPHERE & HYDROSPHERE Hydrology Definitions Streamflow volume of water in a river passing a defined point over a specific time period = VxA discharge m 3 s -1 Runoff excess precipitation - precipitation
More informationMODULE 1 RUNOFF HYDROGRAPHS WORKSHEET 1. Precipitation
Watershed MODULE 1 RUNOFF HYDROGRAPHS WORKSHEET 1 A watershed is an area of land thaaptures rainfall and other precipitation and funnels it to a lake or stream or wetland. The area within the watershed
More informationHydrology and Water Management. Dr. Mujahid Khan, UET Peshawar
Hydrology and Water Management Dr. Mujahid Khan, UET Peshawar Course Outline Hydrologic Cycle and its Processes Water Balance Approach Estimation and Analysis of Precipitation Data Infiltration and Runoff
More informationTexas A & M University and U.S. Bureau of Reclamation Hydrologic Modeling Inventory Model Description Form
Texas A & M University and U.S. Bureau of Reclamation Hydrologic Modeling Inventory Model Description Form JUNE 18, 1999 Name of Model: MIKE 11 RR (Rainfall Runoff) Model Type: The MIKE 11 RR model is
More informationUncertainty in Hydrologic Modelling for PMF Estimation
Uncertainty in Hydrologic Modelling for PMF Estimation Introduction Estimation of the Probable Maximum Flood (PMF) has become a core component of the hydrotechnical design of dam structures 1. There is
More informationComputation of Hydrographs in Evros River Basin
European Water 31: 33-42, 2010. 2010 E.W. Publications Computation of Hydrographs in Evros River Basin P. Angelidis 1, G. Mystakidis 2, S. Lalikidou 2, V. Hrissanthou 2 and N. Kotsovinos 3 Department of
More informationWMS Tools For Computing Hydrologic Modeling Parameters
WMS Tools For Computing Hydrologic Modeling Parameters Lesson 9 9-1 Objectives Use the drainage coverage as a basis for geometric parameters as well as overlaying coverages to compute important hydrologic
More informationUNIT HYDROGRAPH AND EFFECTIVE RAINFALL S INFLUENCE OVER THE STORM RUNOFF HYDROGRAPH
UNIT HYDROGRAPH AND EFFECTIVE RAINFALL S INFLUENCE OVER THE STORM RUNOFF HYDROGRAPH INTRODUCTION Water is a common chemical substance essential for the existence of life and exhibits many notable and unique
More informationHydrologic Study Report for Single Lot Detention Basin Analysis
Hydrologic Study Report for Single Lot Detention Basin Analysis Prepared for: City of Vista, California August 18, 2006 Tory R. Walker, R.C.E. 45005 President W.O. 116-01 01/23/2007 Table of Contents Page
More informationCLARK COUNTY REGIONAL FLOOD CONTROL DISTRICT HYDROLOGIC CRITERIA AND DRAINAGE DESIGN MANUAL
CLARK COUNTY REGIONAL FLOOD CONTROL DISTRICT HYDROLOGIC CRITERIA AND DRAINAGE DESIGN MANUAL SECTION 600 STORM RUNOFF TABLE OF CONTENTS 601 INTRODUCTION 603 601.1 - Basin Characteristics 603 602 TIME OF
More informationWater Budget III: Stream Flow P = Q + ET + G + ΔS
Water Budget III: Stream Flow P = Q + ET + G + ΔS Why Measure Streamflow? Water supply planning How much water can we take out (without harming ecosystems we want to protect) Flood protection How much
More informationWater Budget III: Stream Flow P = Q + ET + G + ΔS
Water Budget III: Stream Flow P = Q + ET + G + ΔS Why Measure Streamflow? Water supply planning How much water can we take out (without harming ecosystems we want to protect) Flood protection How much
More informationHydrologic Modeling Overview
Hydrologic Modeling Overview Chuck Downer, PhD, PE Hydrologic Systems Branch Coastal and Hydraulics Laboratory Engineer Research and Development Center Vicksburg, Mississippi Hydrologic processes Hydrologic
More informationHydrologic Engineering Center. Training Course on. Hydrologic Modeling with HEC-HMS. Davis, CA. Course Description
Hydrologic Engineering Center Training Course on Hydrologic Modeling with HEC-HMS Davis, CA Course Description The course provides an introduction to HEC-HMS for new users, focusing both on using the program
More information1. Stream Network. The most common approach to quantitatively describing stream networks was postulated by Strahler (1952).
1. Stream Network The most common approach to quantitatively describing stream networks was postulated by Strahler (1952). First Order Streams streams with no tributaries. Second Order Streams begin at
More informationASSESSMENT OF DRAINAGE CAPACITY OF CHAKTAI AND RAJAKHALI KHAL IN CHITTAGONG CITY AND INUNDATION ADJACENT OF URBAN AREAS
Proceedings of the 4 th International Conference on Civil Engineering for Sustainable Development (ICCESD 2018), 9~11 February 2018, KUET, Khulna, Bangladesh (ISBN-978-984-34-3502-6) ASSESSMENT OF DRAINAGE
More informationHydrologic Engineering Center. Training Course on. Hydrologic Modeling with HEC-HMS April 2018 Davis, CA. Course Description
Hydrologic Engineering Center Training Course on Hydrologic Modeling with HEC-HMS 23-27 April 2018 Davis, CA Course Description The course provides an introduction to HEC-HMS for new users, focusing both
More informationAPPENDIX IV. APPROVED METHODS FOR QUANTIFYING HYDROLOGIC CONDITIONS OF CONCERN (NORTH ORANGE COUNTY)
APPENDIX IV. APPROVED METHODS FOR QUANTIFYING HYDROLOGIC CONDITIONS OF CONCERN (NORTH ORANGE COUNTY) Hydromodification design criteria for the North Orange County permit area are based on the 2- yr, 24-hr
More informationImpacts of Rainfall Event Pattern and Land-Use Change on River Basin Hydrological Response: a Case in Malaysia
International Proceedings of Chemical, Biological and Environmental Engineering, Vol. 93 (2016) DOI: 10.7763/IPCBEE. 2016. V93. 23 Impacts of Rainfall Event Pattern and Land-Use Change on River Basin Hydrological
More informationChapter H. Introduction to Surface Water Hydrology and Drainage for Engineering Purposes
Chapter H. Introduction to Surface Water Hydrology and Drainage for Engineering Purposes As seen in Figure H.1, hydrology is a complex science that deals with the movement of water between various stages
More informationNote that the Server provides ArcGIS9 applications with Spatial Analyst and 3D Analyst extensions and ArcHydro tools.
Remote Software This document briefly presents the hydrological and hydraulic modeling software available on the University of Nice Server with Remote Desktop Connection. Note that the Server provides
More informationHydrologic Calibration:
Hydrologic Calibration: UPDATE OF EFFECTIVE HYDROLOGY FOR MARYS CREEK October 2010 Agenda Background Hydrologic model Calibrated rainfall Hydrologic calibration 100 year discharges, Existing Conditions
More informationEvent and Continuous Hydrological Modeling with HEC- HMS: A Review Study
Event and Continuous Hydrological Modeling with HEC- HMS: A Review Study Sonu Duhan *, Mohit Kumar # * M.E (Water Resources Engineering) Civil Engineering Student, PEC University Of Technology, Chandigarh,
More informationCHAPTER 4 DISCUSSION. Total Flow. For the study watershed, the water budget equation is defined as inflow equals
CHAPTER 4 DISCUSSION Total Flow For the study watershed, the water budget equation is defined as inflow equals outflow minus any losses. Inflow consists of precipitation. Outflow consists of stream flow
More informationOverview of NRCS (SCS) TR-20 By Dr. R.M. Ragan
Overview of NRCS (SCS) TR-20 By Dr. R.M. Ragan TR-20 is a computer program for the simulation of runoff occurring from a single storm event. The program develops flood hydrographs from runoff and routes
More informationThe Islamic University of Gaza- Civil Engineering Department Sanitary Engineering- ECIV 4325 L5. Storm water Management
The Islamic University of Gaza- Civil Engineering Department Sanitary Engineering- ECIV 4325 L5. Storm water Management Husam Al-Najar Storm water management : Collection System Design principles The Objectives
More informationModeling catchment scale infiltration
Modeling catchment scale infiltration Serena Ceola, DICAM, University of Bologna Advanced Hydrology & Water Resources Management INTRODUCTION OBJECTIVE: Understanding the processes and the methods for
More informationAppendix 4 ADB. IN THE 6 CI's RIVER BASIN TERRITORY - PACKAGE B. Final Report B.2 - Strategic Spatial Planning
ADB TA 7189-INO: INSTITUTIONAL STRENGTHENING FOR INTEGRATED WATER RESOURCES MANAGEMENT (IWRM) IN THE 6 CI's RIVER BASIN TERRITORY - PACKAGE B Final Report B.2 - Strategic Spatial Planning Appendix 4 Modeling
More informationRunoff Processes. Daene C. McKinney
CE 374 K Hydrology Runoff Processes Daene C. McKinney Watershed Watershed Area draining to a stream Streamflow generated by water entering surface channels Affected by Physical, vegetative, and climatic
More informationEngineering Hydrology. Class 16: Direct Runoff (DRO) and Unit Hydrographs
Engineering Hydrology Class 16: and s Topics and Goals: 1. Calculate volume of DRO from a hydrograph; 2. Complete all steps to develop a. Class 14: s? HG? Develop Ocean Class 14: s? HG? Develop Timing
More information5/11/2007. WinTR-55 for Plan Reviews Small Watershed Hydrology Overview
WinTR-55 for Plan Reviews Small Watershed Hydrology Overview 1 Overview Course Outline Historical Background Program Description Model Capabilities and Limitations This is the Overview portion of the training.
More informationLecture 9A: Drainage Basins
GEOG415 Lecture 9A: Drainage Basins 9-1 Drainage basin (watershed, catchment) -Drains surfacewater to a common outlet Drainage divide - how is it defined? Scale effects? - Represents a hydrologic cycle
More informationA simple model for low flow forecasting in Mediterranean streams
European Water 57: 337-343, 2017. 2017 E.W. Publications A simple model for low flow forecasting in Mediterranean streams K. Risva 1, D. Nikolopoulos 2, A. Efstratiadis 2 and I. Nalbantis 1* 1 School of
More informationRunoff Hydrographs. The Unit Hydrograph Approach
Runoff Hydrographs The Unit Hydrograph Approach Announcements HW#6 assigned Storm Water Hydrographs Graphically represent runoff rates vs. time Peak runoff rates Volume of runoff Measured hydrographs are
More informationSuspended Sediment Discharges in Streams
US Army Corps of Engineers Hydrologic Engineering Center Suspended Sediment Discharges in Streams April 1969 Approved for Public Release. Distribution Unlimited. TP-19 REPORT DOCUMENTATION PAGE Form Approved
More informationModule 3. Lecture 4: Introduction to unit hydrograph
Lecture 4: Introduction to unit hydrograph (UH) The unit hydrograph is the unit pulse response function of a linear hydrologic system. First proposed by Sherman (1932), the unit hydrograph (originally
More informationFlood forecasting model based on geographical information system
doi:10.5194/piahs-368-192-2015 192 Remote Sensing and GIS for Hydrology and Water Resources (IAHS Publ. 368, 2015) (Proceedings RSHS14 and ICGRHWE14, Guangzhou, China, August 2014). Flood forecasting model
More informationNorman Maclean Snowmelt Flow rate Storm flows fs (c flow m a tre S
Eventually, all things merge into one, and a river runs through it. Norman Maclean Understanding Streamflow ADEQ SW Short Course June 13, 213 Phoenix, AZ Hydrographs Discharge (Q) USGS flow data & plots
More informationIndex. Page numbers followed by f indicate figures.
Index Aerodynamic method, 103, 110 111 Algae, 131, 173, 175 Alternate depth, 88 Alternating block method, 132, 140 141 Attenuation, 106, 107f, 118, 120 Page numbers followed by f indicate figures. Baseflow
More informationRainfall, Runoff and Peak Flows: Calibration of Hydrologic Design Methods for the Kansas City Area
Rainfall, Runoff and Peak Flows: Calibration of Hydrologic Design Methods for the Kansas City Area Bruce McEnroe, Bryan Young, Ricardo Gamarra and Ryan Pohl Department of Civil, Environmental, and Architectural
More informationWhat s so hard about Stormwater Modelling?
What s so hard about Stormwater Modelling? A Pugh 1 1 Wallingford Software Pty Ltd, ann.pugh@wallingfordsoftware.com Abstract A common misconception of stormwater modelling is that it is simple. While
More informationPhase 1 Part 2 CSO Control Plan Wellington Avenue CSO Facility. Hydraulic Modeling Software Selection
DRAFT Technical Memorandum Phase 1 Part 2 CSO Control Plan Wellington Avenue CSO Facility Hydraulic Modeling Software Selection Prepared for: City of Newport Public Works Department 70 Halsey Street Newport,
More informationChapter 4 "Hydrology"
Chapter 4 "Hydrology" Revised TxDOT Hydraulics Manual Jorge A. Millan, P.E., CFM TxDOT Design Division/Hydraulics 2012 Transportation Short Course October 17, 2012 Chapter 4 Sections 1 Hydrology s Role
More informationSection 600 Runoff Table of Contents
Section 600 Runoff Table of Contents 601 INTRODUCTION...600-1 602 RATIONAL METHOD...600-1 602.1 Rational Method Formula...600-2 602.2 Time of Concentration...600-2 602.3 Intensity...600-4 602.4 Runoff
More informationRadar-based flood forecasting: Quantifying hydrologic prediction uncertainty
Severe Storms Prediction and Global Climate Impact on the Gulf Coast, Rice University, October 29 31, 2008 Radar-based flood forecasting: Quantifying hydrologic prediction uncertainty Baxter E. Vieux,
More informationTHE RATIONAL METHOD FREQUENTLY USED, OFTEN MISUSED
THE RATIONAL METHOD FREQUENTLY USED, OFTEN MISUSED Mark Pennington, Engineer, Pattle Delamore Partners Ltd, Tauranga ABSTRACT The Rational Method has been in use in some form or another at least since
More informationThe use of soil moisture in hydrological forecasting
The use of soil moisture in hydrological forecasting V.A. Bell, E.M. Blyth and R.J. Moore CEH Wallingford, Oxfordshire, UK vib@ceh.ac.uk Abstract Soil moisture data generated by a land-surface scheme within
More informationBAEN 673 / February 18, 2016 Hydrologic Processes
BAEN 673 / February 18, 2016 Hydrologic Processes Assignment: HW#7 Next class lecture in AEPM 104 Today s topics SWAT exercise #2 The SWAT model review paper Hydrologic processes The Hydrologic Processes
More informationHypothetical Flood Computation for a Stream System
US Army Corps of Engineers Hydrologic Engineering Center Hypothetical Flood Computation for a Stream System December 1968 Approved for Public Release. Distribution Unlimited. TP-12 REPORT DOCUMENTATION
More informationAlternative Approaches to Water Resource System Simulation
US Army Corps of Engineers Hydrologic Engineering Center Alternative Approaches to Water Resource System Simulation May 1972 Approved for Public Release. Distribution Unlimited. TP-32 REPORT DOCUMENTATION
More informationEffect of Land Surface on Runoff Generation
Effect of Land Surface on Runoff Generation Context: Hydrologic Cycle Runoff vs Infiltration Infiltration: Process by which water on the ground surface enters the soil Runoff: Water (from rain, snowmelt,
More informationThe use of GIS to derive distributed unit hydrographs for streamflow prediction
8 th World IMACS / MODSIM Congress, Cairns, Australia 3-7 July 29 http://mssanz.org.au/modsim9 The use of GIS to derive distributed unit hydrographs for streamflow prediction Gibbs, M.S., Maier, H.R.,
More informationCHAPTER FIVE Runoff. Engineering Hydrology (ECIV 4323) Instructors: Dr. Yunes Mogheir Dr. Ramadan Al Khatib. Overland flow interflow
Engineering Hydrology (ECIV 4323) CHAPTER FIVE Runoff Instructors: Dr. Yunes Mogheir Dr. Ramadan Al Khatib Overland flow interflow Base flow Saturated overland flow ١ ٢ 5.1 Introduction To Runoff Runoff
More informationCONSISTENCY BETWEEN CUHP AND RATIONAL METHODS James C.Y. Guo, PhD, P.E., and Ben Urbonas, P.E., D.WRD 08/04/2008
CONSISTENCY BETWEEN CUHP AND RATIONAL METHODS James C.Y. Guo, PhD, P.E., and Ben Urbonas, P.E., D.WRD 08/04/2008 When estimating runoff for a storm event it is assumed that the storm runoff occurs from
More informationCIVE22 BASIC HYDROLOGY Jorge A. Ramírez Homework No 7
Hydrologic Science and Engineering Civil and Environmental Engineering Department Fort Collins, CO 80523-1372 (970) 491-7621 CIVE22 BASIC HYDROLOGY Jorge A. Ramírez Homework No 7 1. Obtain a Unit Hydrograph
More informationCEE6400 Physical Hydrology
CEE6400 Physical Hydrology Midterm Review Learning Objectives (what you should be able to do) Hydrologic data, the hydrologic cycle and water balance (HW 1) Work with hydrologic data, quantify uncertainty
More informationSan Antonio Water System Mitchell Lake Constructed Wetlands Below the Dam Preliminary Hydrologic Analysis
San Antonio Water System enhancement. This recommendation was based on limited water quality data provided by SAWS and using a free-water surface constructed wetland with approximately 112 acres of wetted
More informationTraining Course Brochure Building Capacity in Rural & Urban Water Management
Training Course Brochure 2015 Building Capacity in Rural & Urban Water Management Introduction The WastePro Academy seeks to encourage the use of desktop software applications in the water and environmental
More informationHydrologic Cycle. Water Availabilty. Surface Water. Groundwater
Hydrologic Cycle Hydrologic ydoogccyce cycle Surface Water Groundwater Water Availabilty 1 Hydrologic Cycle Constant movement of water above, on, and, below the earth s surface (Heath) Endless circulation
More informationAPPLICATION OF A HYDRODYNAMIC MIKE 11 MODEL FOR THE EUPHRATES RIVER IN IRAQ
2008/2 PAGES 1 7 RECEIVED 13.1.2008 ACCEPTED 26.4.2008 A. H. KAMEL APPLICATION OF A HYDRODYNAMIC MIKE 11 MODEL FOR THE EUPHRATES RIVER IN IRAQ Ing. Ammar H. Kamel Slovak University of Technology Faculty
More informationApplication of HEC-HMS for Flood Forecasting in Kabkian Basin and Delibajak Subbasin in Iran
IOSR Journal of Engineering (IOSRJE) e-iss: 2250-3021, p-iss: 2278-8719 Vol. 3, Issue 9 (September. 2013), V3 PP 10-16 Application of HEC-HMS for Flood Forecasting in Kabkian Basin and Delibajak Subbasin
More informationCreating a Flood Forecasting System for the San Diego River Watershed
Creating a Flood Forecasting System for the San Diego River Watershed Rand Allan County of San Diego Flood Control Martin J. Teal, P.E., P.H., D.WRE WEST Consultants inc. OUTLINE San Diego River Watershed
More informationModclark Model: Improvement and Application
IOSR Journal of Engineering (IOSRJEN) ISSN: 2250-3021 Volume 2, Issue 7(July 2012), PP 100-118 Modclark Model: Improvement and Application 1 Dr. Amartya Kumar Bhattacharya, 2 Dr.Bruce M. McEnroe, 3 Dr.Hongying
More informationCHOLLA POWER PLANT BOTTOM ASH POND INFLOW DESIGN FLOOD CONTROL SYSTEM PLAN CH_Inflowflood_003_
CHOLLA POWER PLANT BOTTOM ASH POND INFLOW DESIGN FLOOD CONTROL SYSTEM PLAN CH_Inflowflood_003_20161017 This Inflow Design Flood Control System Plan (Plan) document has been prepared specifically for the
More informationContext of Extreme Alberta Floods
Context of Extreme Alberta Floods Introduction Design of water management and stream crossing infrastructure requires determination of hydrotechnical design parameters. These parameters often consist of
More informationGwinnett County Stormwater System Assessment Program
Gwinnett County Stormwater System Assessment Program Jonathan Semerjian, PE Dept. of Water Resources Stormwater Management Sam Fleming, PE Dewberry Presentation Overview Project Background Drivers Enhanced
More informationContinuous Simulation Modeling of Stormwater Ponds, Lakes, & Wetlands: A BUILT-IN APPLICATION OF PONDS 3.2
Continuous Simulation Modeling of Stormwater Ponds, Lakes, & Wetlands: A BUILT-IN APPLICATION OF PONDS 3.2 PRESENTED AT THE SFWMD WORKSHOP PRE-DEVELOPMENT VERSUS POST DEVELOPMENT RUNOFF VOLUME ANALYSIS
More informationPeak discharge computation
Ia/P 4 Peak Dischage Method Graphical Peak Discharge Method This chapter presents the Graphical Peak Discharge method for computing peak discharge from rural and urban areas. The Graphical method was developed
More informationReview of existing simulation based flood frequency frameworks in Greece
European Procedures for Flood Frequency Estimation FLOODFREQ COST ACTION ES0901 3 rd Management Committee Meeting, Prague, 28 29 October 2010 WG3: Use of rainfall runoff models for flood frequency estimation
More informationAnalysis of the Sensitivity of the Variables Affecting Urban Flood through Two Rational and SWMM Methods by SSA
Bulletin of Environment, Pharmacology and Life Sciences Bull. Env.Pharmacol. Life Sci., Vol 4 [Spl issue ] 5: 36-37 4 Academy for Environment and Life Sciences, India Online ISSN 77-88 Journal s URL:http://www.bepls.com
More informationHydrologic cycle, runoff process
Hydrologic cycle, runoff process Motivation of hydrological modelling What happens at the catchment and in the stream when it rains? How does the increased/decreased runoff affect (not only) the landowners
More informationRepresenting the Integrated Water Cycle in Community Earth System Model
Representing the Integrated Water Cycle in Community Earth System Model Hong-Yi Li, L. Ruby Leung, Maoyi Huang, Nathalie Voisin, Teklu Tesfa, Mohamad Hejazi, and Lu Liu Pacific Northwest National Laboratory
More informationSoftware Applications for Runoff Hydrological Assessment
Bulletin UASVM Horticulture, 67(2)/2010 Print ISSN 1843-5254; Electronic ISSN 1843-5394 Software Applications for Runoff Hydrological Assessment Severin CAZANESCU 1), Sorin CIMPEANU 1), Oana GUI 2), Dana
More informationFlood Forecasting - What Can You Do With Your Data?
Flood Forecasting - What Can You Do With Your Data? C Druery 1, D McConnell,2 1 WorleyParsons, Sydney, NSW 2 WorleyParsons, Sydney, NSW Abstract Driven by the large scale flooding over the past several
More informationWatershed Management And Modeling. Gridded Surface Subsurface Hydrologic
Gridded Surface Subsurface Hydrologic Analysis What is GSSHA? GSSHA is a complete watershed simulation and management model used for hydrologic, hydraulic, sediment and quality simulation and management.
More informationSTORMWATER HYDROLOGY
..CHAPTER.. STORMWATER HYDROLOGY 3.1 Introduction to Hydrologic Methods Hydrology is the science dealing with the characteristics, distribution, and movement of water on and below the earth's surface and
More informationINFLOW DESIGN FLOOD CONTROL SYSTEM PLAN PLANT GREENE COUNTY ASH POND ALABMA POWER COMPANY
INFLOW DESIGN FLOOD CONTROL SYSTEM PLAN PLANT GREENE COUNTY ASH POND ALABMA POWER COMPANY Section 257.82 of EPA s regulations requires the owner or operator of an existing or new CCR surface impoundment
More informationChapter 3 Physical Factors Affecting Runoff
Chapter 3 Physical Factors Affecting Runoff Copyright 2003 David G Tarboton, Utah State University CHAPTER 3: PHYSICAL FACTORS AFFECTING RUNOFF The general climatic regime controls the total volume of
More informationComparison of Recharge Estimation Methods Used in Minnesota
Comparison of Recharge Estimation Methods Used in Minnesota by Geoffrey Delin, Richard Healy, David Lorenz, and John Nimmo Minnesota Ground Water Association Spring Conference Methods for Solving Complex
More informationFlood Modelling and Water Harvesting Plan for Paravanar Basin
International Journal of ChemTech Research CODEN (USA): IJCRGG, ISSN: 0974-4290, ISSN(Online):2455-9555 Vol.10 No.14, pp 01-08, 2017 Flood Modelling and Water Harvesting Plan for Paravanar Basin Dhinesh
More informationSTORMWATER MANAGEMENT: Emerging Planning Approaches and Control Technologies
STORMWATER MANAGEMENT: Emerging Planning Approaches and Control Technologies Chapter 5 " HYDROLOGIC SYSTEMS " CHAPTER 5 5.1 INTRODUCTION 5.1.1 Training Objectives The objectives of this module are:! to
More informationSimulation of Event Based Runoff Using HEC-HMS Model for an Experimental Watershed
International Journal of Hydraulic Engineering 213, 2(2): 28-33 DOI: 1.5923/j.ijhe.21322.2 Simulation of Event Based Runoff Using HEC-HMS Model for an Experimental Watershed Reshma T *, Venkata Reddy K,
More informationEngineering Hydrology Class 3
Engineering Hydrology Class 3 Topics and Goals: I.Develop s (estimate precipitation) II.Develop simple constant intensity design storm III.Develop SCS design storm Ocean s Why do we want to derive the?
More informationLinking hydrology to erosion modelling in a river basin decision support and management system
Integrated Water Resources Management (Proceedings of a symposium held at Davis. California. April 2000). I APIS Publ. no. 272. 2001. 243 Linking hydrology to erosion modelling in a river basin decision
More informationDRAINAGE & DESIGN OF DRAINAGE SYSTEM
Drainage on Highways DRAINAGE & DESIGN OF DRAINAGE SYSTEM P. R.D. Fernando Chartered Engineer B.Sc.(Hons), M.Eng. C.Eng., MIE(SL) Drainage Requirement of Highway Drainage System Introduction Drainage means
More informationFlorida s Deranged Hydrology and the National Water Model Challenges and Opportunities
Florida s Deranged Hydrology and the National Water Model Challenges and Opportunities Peter J. Singhofen, P.E. Streamline Technologies, Inc. psinghofen@icpr4.com Part 1 Florida is Different 2 2016, Streamline
More informationHYDROLOGICAL INVESTIGATIONS IN SEMI-ARID AREAS CASE STUDY: RAILWAY CULVERT DESIGN IN LIBYA
Ninth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 1027 YDROLOGICAL INVESTIGATIONS IN SEMI-ARID AREAS CASE STUDY: RAILWAY CULVERT DESIGN IN LIBYA S. A. Serrag 1, A. M.
More informationMeasuring discharge. Climatological and hydrological field work
Measuring discharge Climatological and hydrological field work 1. Background Discharge (or surface runoff Q s) refers to the horizontal water flow occurring at the surface in rivers and streams. It does
More informationHYDROLOGICAL IMPACTS OF URBANIZATION: WHITE ROCK CREEK, DALLAS TEXAS. Julie Anne Groening Vicars, B.A. Thesis Prepared for the Degree of
HYDROLOGICAL IMPACTS OF URBANIZATION: WHITE ROCK CREEK, DALLAS TEXAS Julie Anne Groening Vicars, B.A. Thesis Prepared for the Degree of MASTER OF SCIENCE UNIVERSITY OF NORTH TEXAS December 25 APPROVED:
More informationAvailable online at ScienceDirect. Procedia Economics and Finance 18 ( 2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Economics and Finance 18 ( 2014 ) 675 683 4th International Conference on Building Resilience, Building Resilience 2014, 8-10 September
More information