Numerical Simulation 2019/2020: Choosing the right model in applied hydraulics, 11 Jan 2020, Marlin Shlewet, Alaa Abusamra

CALIBRATED ROUGHNESS VALUES FOR DIFFERENT FLOW SCENARIOS

Marlin, Shlewet Mailing address [email protected]

Alaa, Abusamra Mailing address [email protected]

SUB TITLE

CASE STUDY: RIVER SECTION - , 1D MODELLING

KEY WORDS

Manning coefficient, Validation, Calibrition, Water Level, Unsteady Simulation, Cross Section, Boundary Conditions, Flood Period, Dry Period.

ABSTRACT

Roughness coefficient is the main important parameter in numerical surface water models, where flow discharge, water level and velocity are calculated from surface water eequations depending on it. This parameter´s value is non-measurable directely from the observation water field, but related to different other factors. The relevant concerned question in this paper is how can we define the exact value which can approximate our simulated model to the real case with enough accuracy?. In this present case study, 1D hydrodynamic model using HEC-ras software has been setted-up -to section of River Rhine approximately 34.1 km length, consists of 69 cross- sections from Ruhrort in the South to Wesel in the North in - using observation time series flow data and H/Q relationship of a flood return period 10 years. The model is used to calibrate the optimized roughness value using evaluation coefficients method , comparing with the observation water level values down stream the river, after running the unsteady flow simulation for three different flow discharge scenarios (high,mean and low flow discharge) and using different manning values between the range 25-40 for each. Also validation for the calibrated manning value of each scenario is implemented. As a result, there is a focus point to find out rather if the flow discharge conditions have an impact on the choosen roughness coefficient value to model the unsteady simulation and how could this impact be explained. This study could be as pre-step to prepare a well-flood model for management plaining or other flow condition.

1. INTRODUCTION

The one-dimensional (1D) hydrodynamic model has been used for river engineering and mapping flood risks which will help engineers to take precautionary measures to minimize flood damage. To obtained a realible and accureate results from the application of hydraulic models modeler must identify errors and inconsistency in the input data parameters1. To do so calibration is plays an important role to adjust a model’s parameters such as roughness to reproduce the suimulation results to an acceptable accuracy. The overall target of this task is to analyse potential flooding in the Lower Rhine from Ruhrort (the up-stream guage in South) to Wesel (the down-stream guage in North). Using 1-D hydraulic analysis model in Hec Ras to simulate and identify where and when the flooding may occure. This paper particulary investigate the dependency of the Manning’s roughness coefficient on the flow water condition by simulating defferent numerical unsteady model to check and calibrate Manning’s roughness coefficient for three diferent flow discharge scenarios, high flow during flooding, mean flow during normal conditions and low flow during driest conditions. Generally roughness decreas with increased stage and flow. However the Numerical Simulation 2019/2020: Choosing the right model in applied hydraulics, 11 Jan 2020, Marlin Shlewet, Alaa Abusamra findings of this paper showes different behavire for Rhine rever as the n values increase with increased stage due to effect of the vegitation which make the banks of the river are rougher than the channel. 2.

2. FIRST CALIBRATION SCENARIO FOR THE MANNING`S VALUE

2.1 Create and Run 1-D Model Simulation with Flood Flow Scenario

2.1.1 Geometry& Hotstart Model The geometry of the river is represented by defining topographic available data the coordinates (xyz) of the cross sections along the river and roughness coefficient parameter which is specified as a Manning´s value-n from a range values in 1-D HEC-ras model so it is associated with high degree of uncertainties. It represents the friction of water flow with the river bed and bank and depends on some direct factors like slope and bed type(high roughness associated with low n value)3. The concerned aim of this study is the roughness parameter to be calibrated to minimize errors in the simulation results. Thus , three simulations will be run for this scenario by estimating for each a uniform Manning roughness coefficient of 0.04, 0.029 and 0.025 respectively and that for all sections of the river and embankments, in order to do calibrartion. Hence, the available time series flow discharge of the studied flood event and the given rating curve are considered as boundary conditions up and down stream respectively. Result of HOTSTART simulation (unsteady simulation with steady state for 24h) is defined as initial conditions for the unsteady simulations.

2.1.2 Set-up and Run the Flood Model Simulation The choosen period to simulate the unsteady model for the flood event is continued for two months from 1/10/98 to 30/11/98, so boundary conditions of a flow hydrograph for upstream is used for that period and taken from the available time series flow discharge of Ruhrort station for that period. For the downstream the same rating curve is used. To run this unsteady simulation, a suitable time step for numerical computations of 5 min has been chosed according to the Courant Friedrichs Lewy (CFL) by taking into account the distance between cross sections (500m ) to avoid instability. This simulation run 3 times for each n value.

3.2 Calibration with Results and Discussions

In this study the calibration will be done for the roughness coefficient of Manning´s value for the river bed by using comparison period during the flood (high discharge values) from (30-Oct to 9-Nov)/98 and compare between the observation water levels downstream of Wesel station and the three different simulation data of water levels downstream resulting from changing Manning´s value (0.025, 0.029, 0.04). Down stream simulation data has chosen to be compared in order to involve all information from flood and upstream. This comparison to define the optimized manning´s value for this scenario will be evaluated by using different evaluation coefficients:  error index - Root Mean Square Error (RMSE) n obs sim 2 √∑ (Vi −Vi ) RMSE = i=1 , (1) n

- Percent bias (PBIAS) n obs sim 100∗∑ (Vi −Vi ) i=1 PBIAS = n , (2) obs ∑ Vi i=1  dimensionless evaluation - Nash-Sutcliffe efficiency (NSE)

Numerical Simulation 2019/2020: Choosing the right model in applied hydraulics, 11 Jan 2020, Marlin Shlewet, Alaa Abusamra

n obs sim 2 ∑ (Vi −Vi ) i=1 NSE = 1 − n , (3) obs obs 2 ∑ (Vi −Vmean) i=1

obs sim Where Vi is the observed water level at time i, Vi is the simuated water level at time i.

n=0.04 n=0.029 n=0.025 Evaluation Coefficient K=25 K=34 K=40 RMSE 0.085 0.093 0.100 PBIAS 0.080 0.077 0.083 NSE 9.995E-01 9.994E-01 9.993E-01 Table 1: Result of Coef evaluation for water level downstream from calibration

The best match for this scenario with the observed H values during flood discharge was with a K value of 25 (n=0.04) for the river Rhein, where the RMSE and PBIAS have shown the minimum value and the NSE has shown the best fit with observed plot as showen in figure 1. The interpretation for this result is because of that the high water level during flooding makes the friction between the water flow and the river bed on a higher elevation where normally the distribution of the vegetation and other rough things are existed intensively if we do a site visit, means big Manning´s value n (small Strickler´s value-K).

Figure 1: The graph shows Water Level at Wesel for different n value with high flow simulation

3.3 Validation

Validation is important to check the result from calibrartion of the chosen roughness-coe for the high flow scenario and determine whether the model provides reliable and accepted result for other similar events also. The chosen time of validation is another flooding period from 1/12/2002 to 31/01/2003, where an unsteady simulation is run with same conditions as before and calibrated for the period(01/01-11/01)2003, using the chosen Manning´s value n=0.04.

Evaluation Coefficient RMSE PBIAS NSE n=0.04 - K=25 0.043 0.006 9.998E-01 Table 2: Result of Coef evaluation for water level downstream from validation for n=0.04

The evaluation of the selected flood period gives even better result of RMSE and PBIAS, and a perfect fit with observation data from NSE value. Model is validated with n=0.04 during flood period.

Numerical Simulation 2019/2020: Choosing the right model in applied hydraulics, 11 Jan 2020, Marlin Shlewet, Alaa Abusamra

Figure 2: the graph represent validation of Water Level at Wesel for n=0.04 for with high flow simulation.

4. SECOND CALIBRATION SCENARIO FOR THE MANNING`S VALUE

4.1 Create and Run 1-D Model Simulation with Mean Flow Scenario

Aim of this scenario is to check wether the chosen Manning´s value n=0.04 from previous scenario is also could be an optimized value for mean flow discharge values. A period of mean discharge values is needed to be simulated which is already obviously existed on the same period of the previous two months have been simulated from 1/10/98 to 30/11/. Hence, results from previous scenario is used for this scenario.

4.2 Calibration with Results and Discussions

In this scenario the calibration will be done using same strategy as before, but using comparison period during the mean discharge values from (17-Oct to 24-Oct)/98,

n=0.04 n=0.029 n=0.025 Evaluation Coefficient K=25 K=34 K=40 RMSE 0.021 0.019 0.020 PBIAS -0.059 -0.009 0.001 NSE 9.999E-01 9.999E-01 9.999E-01 Table 3: Result of Coef evaluation for water level upstream from calibration mean flow scenario

The plot presented in the figure (3) and NSE values for the three different K show best fit, means for the comparison period for mean discharge there is no changing of the behavior itself from the simulated data with changing K. Beside, RMSE and PBIAS which express the error index between values have shown the minimum value for K=34, even if it is very close to K=40. The interpretation is that mean water level makes the friction between the water flow and the river bed on a normal elevation where normally the bed is smoother and distribution of vegetation is less. As a result, when this kind of scenario is simulated, K could be assumed 34 as a primary value. Wherase for K= 25 the error is higher due to less roughness with less water level.

4.3 Validation

For this scenario the chosen time of validation is another mean discharge period from 1/12/2002 to 31/01/2003, the unsteady simulation is calibrated for the period (15/12-24/12)2002, using n=0.029.

Evaluation Coefficient RMSE PBIAS NSE n=0.029 - K=34 0.033 -0.072 9.997E-01 Table 4: Result of Coef evaluation for water level downstream from validation for n=0.029 Numerical Simulation 2019/2020: Choosing the right model in applied hydraulics, 11 Jan 2020, Marlin Shlewet, Alaa Abusamra

The evaluation of the selected flood period gives very good result of RMSE and PBIAS and a perfect fit with observation data from NSE value. Model is validated with n=0.029 during mean discharge period.

Figure 3: The graph shows Water Level at Wesel for different n value with mean flow simulation.

5. THIRD CALIBRATION SCENARIO FOR THE MANNING`S VALUE

5.1 Create and Run 1-D Model Simulation with Low Flow Scenario

The third scenario testing which of the chosen Manning´s values n=0.04, 0.029 or 0.025 could be an optimized value for low flow discharge values. A period of low discharge values (Discharge less than 1300 m^3/sec) have been simulated for a two monthes from 1/08 to 30/09/2003.

5.2 Calibration with Results and Discussions

In this scenario the calibration will be done by comparing period during the low discharge values from (26- 28 Sep 2003).

n=0.04 n=0.029 n=0.025 Evaluation Coefficient K=25 K=34 K=40 RMSE 0.811 0.812 0.811 PBIAS -7.147 -7.156 -7.150 NSE 3.8E-02 3.8E-02 3.8E-02 Table 5: Result of Coef evaluation for water level downstream from calibration low discharge

The plot shown in the figure (4) and NSE values for the three different K show a similar fit, means for the comparison period for mean discharge there is no changing of the behavior itself from the simulated data with changing K. Beside, RMSE and PBIAS error index have shown the minimum value for K=25 and K=40, even if it is very close to K=34. As a result, when this kind of scenario is simulated, K could be assumed 25,34 or 40 as a primary range values.

5.3 Validation

Validation is done by taken similer low flow discharge from another period 1/09-30/10/1990 and calibrated with the observation values from the period (16/10-25/10)1990, using n=0.025.

Evaluation Coefficient RMSE PBIAS NSE Numerical Simulation 2019/2020: Choosing the right model in applied hydraulics, 11 Jan 2020, Marlin Shlewet, Alaa Abusamra

n=0.025 - K=40 0.039 -0.149 1.0E+00 Table 6: Result of Coef evaluation for water level downstream from validation for n=0.025, low flow scenario

The evaluation of the selected period gives very good result of RMSE and PBIAS and a perfect fit with observation data from NSE value. Model is validated with n=0.025 during low discharge period.

Figure 4: The graph shows Water Level at Wesel for different n value with low flow simulation

5. CONCLUSIONS

In conclusion, estimating manning’s n values by using hydrolic model and calibirate n values to observed profile data is considered as one of the best approach to estimate roughness values2. Results show increasing of the Rihne river roughnees with increasing stage and flow as an explanation could be the banks of the river are rougher than the channel bottom (vegetation). For the three flow conditions, appropriate n values was 0.04 for high flow and 0.029 for mean flow. However for the low flow scenario all tree values show reliable results.

ACKNOWLEDGEMENTS

The study was primarily undertaken during the third semester as a requirement to complet ( module 3.1: Numerical Simulation). The authers are grateful to Prof. Frank Molkenthin- lecturer of the cource and academic supervisor - for his continuous support to finish this paper. Tow authers have worked together on writing the structure of the first calibration scenario.the first auther: Marlin worked specifically on writing Abstract, Overview and the second scenario. The second author: Alaa worked on writing Introduction, Conclusion and the third scenario.

REFERENCES AND CITATIONS

1 Thamer Ahmed Mohammed et al 2011 IOP Conf. Ser.: Mater. Sci. Eng. 17 012040

2 Gary W.Brunner, P.E.: Calibriation of unsteady flow model

3 Wright, N., and Crosato, A. 2011. “2.07 - The Hydrodynamics and Morphodynamics of Rivers.” In Treatise on Water Science, 135–156. Oxford:Elsevier.