Microsoft Word - 40-3703l_s_ETASR_V10_N4_pp6126-6130 Engineering, Technology & Applied Science Research Vol. 10, No. 4, 2020, 6126-6130 6126 www.etasr.com Khalifa: Simulation of Dissolved Oxygen and Dissolved Substrate for Hasel River Simulation of Dissolved Oxygen and Dissolved Substrate for Hasel River Walid M. A. Khalifa Civil Engineering Department Hail University, Saudi Arabia and Fayoum University, Egypt khalifawalid@yahoo.com Abstract—Hasel is considered a moderately polluted river in Germany. This study investigated its water quality, examining Dissolved Oxygen (DO) and dissolved substrate (COD) with the use of AQUASIM. The calibration procedure used observed data from various locations along the river. The model’s calibration was used to study the response of Hasel River to the effluents of wastewater treatment plants and sewer overflow emissions. Results revealed that high emissions from sewerage systems may reduce the oxygen concentration to low levels. Furthermore, joined sewer overflows may disrupt the oxygen levels for a long period. In addition, oxygen was over saturation in some periods of the calibration period. The proposed model can be utilized in future analyses, improving the functional understanding of ecological processes in rivers and the identification of ecological effective management strategies. Keywords-water quality models; rivers; dissolved oxygen; eutrophication; impact assessment I. INTRODUCTION Most water quality models are flexible to several environments and subjected to proper definitions of boundary conditions, dimensional variation, and factor characterization. The basic goal of River Water Quality Models (RWQMs) is to simulate the spatial-temporal effect of organic pollution on oxygen level, because of its significance for aquatic life. The Streeter-Phelps model, which describes the balance between deoxygenation and reaeration, was presented in [1]. Another important example of RWQM is the simulation of toxic substances, such as organic chemicals or heavy metals, via the WASP model package [2]. The Aquatic Simulation (AQUASIM) program [3] was designed for creating standardized and consistent river water quality models. Besides the fundamentals of biochemical conversion processes, the structure of compartments of flowing water ecosystems was examined in [4], including the longitudinal, vertical, and lateral zonation patterns. This aspect has not gained much attention in river water quality modeling so far [5], while it should be examined in order to attain an ecologically convenient choice of model compartments and state variables [6, 7]. Streaming waters are linked elements within the hydrological continuum [8-10]. As a result, the hydrological interactions between water compartments are important, as they influence the transport and storage of water, chemical compounds, and nutrients [11-13]. Quantifying the prominence of the exchange processes is complex, due to high spatial and temporal variations of the hydrological system based on a set of factors such as river morphology and hydraulic gradients [14-17]. Furthermore, rivers’ dynamics are shaped by a complicated temporal style, due to daily and seasonal changes in trophic and respiration activities [18-22]. This complexity may be narrowed by analyzing the relevance of temporal dynamics of water constituents in a eutrophic shallow river [5]. This case study considered the Hasel River in Germany. Its main goal was to analyze the oxygen concentrations as one of the water quality parameters using an AQUASIM model [23], and discuss the future river water quality management. II. MATERIALS AND METHODS A. Study Site Hasel River is a right-sided tributary in the middle reach of Werra River, located in Thuringia, Germany between 10°28ʹ and 10°47ʹE longitude and 50°32ʹ and 50°38ʹN latitude. Its total length is 25Km, and it is divided into four reaches coming from four catchment areas of 321Km 2 , having an average gradient of 2‰. The first reach has a length of 6.0Km and an average discharge of 1.54m³/s. The second reach has a length of 9.25Km, and the catchment area discharges 0.02m³/s in average. The third reach has 2.75Km length, and its catchment area brings an average discharge of 3.01m³/s. The fourth reach has 7Km length, and its catchment area has an average discharge of 0.01m³/s. This study selected five trapezoidal cross-sections in the boundaries of the four reaches, with surface width ranging from 0.70m at the spring to 11.00m at the mouth. Hasel River is affected by the influx of the wastewater treatment plants at Suhl and Rohr, located at 19.3Km and 6.1Km upstream (US) its mouth, respectively. The Suhl plant, serving 48,000 inhabitants, discharges directly US the first reach terminal, while the Rohr plant, serving 300 inhabitants, discharges downstream (DS) the third reach terminal [24]. The temporal and spatial dynamics of flow and physico-chemical parameters were measured during 1990-2001 by Thüringen Landesamstadt für Umwelt und Geologie, TLUG (National Thuringia Town of Environment and Geology, Germany). The surface water body was measured nearly every month in monitoring stations at Ellinghausen, Diezhausen, and Simon (2.99, 14.6, 17.9Km upstream the mouth respectively) for pH, Oxygen, COD, BOD, ON, NH4, NO2, NO3, OP, OPO4. Temperature and light intensity were measured daily in 2001. Corresponding author: Walid M. A. Khalifa Engineering, Technology & Applied Science Research Vol. 10, No. 4, 2020, 6126-6130 6127 www.etasr.com Khalifa: Simulation of Dissolved Oxygen and Dissolved Substrate for Hasel River B. Modeling Approach One-dimensional river hydraulics can be described by a set of two partial differential equations representing mass and momentum balance. The two most important approximations to these equations, the kinematic and diffusive wave [25], are implemented in AQUASIM to describe river hydraulics. The equations for river hydraulics are coupled with advection- diffusion equations to describe the transport of dissolved or suspended substances in the water, leading to the following set of differential equations. The first equation describes the water flow through the compartment as: �� �� � �� �� � � (1) The temporal change in the immersed cross-sectional area A, is determined by the spatial gradient of the discharge Q and the lateral inflow q. The behavior of the substances transported with the water flow is described by: �� �� � � �� �� �� � �� ���� �� ������� (2) The concentration is affected by four terms: advection with the water flow, longitudinal dispersion, transformation processes, and lateral inflow or outflow. In order to combine the above system of differential equations, one boundary condition is necessary for (1), and two boundary conditions are required for (2). The boundary condition for (1) that describes the discharge through river’s sections at the start point xs is: ����� � ��� (3) The boundary conditions for (2) are given by the continuity of the substance mass flows entering the river section and a transmission boundary condition [26] at the compartment end: ����� �� � �� � ���, (4) �� ��� � 0 (5) where Iin,C is the total mass input of substances, described by the concentration C per unit of time. Equation (5) is omitted for dispersion-free transport. The longitudinal dispersion coefficient in (2) is estimated according to [27], as: � � �� �!� "∗$ (6) where cf is a non-dimensional dispersion coefficient, w is the surface width of the river, d is the mean river depth, and %∗ � &'(�� is the friction velocity. The empirical expression (7) is used to calculate the friction slope: �� � ) *+,� ) -./0 12 (7) where KSt is the friction coefficient according to Strickler, R is the hydraulic radius of the river, and v is the average cross- sectional flow velocity. C. Modeling of Oxygen Series The oxygen balance of Hasel River is characterized by high daily and seasonal temporal dynamics triggered from radiation and intense production-respiration processes. According to (2), it can be modeled with a simple version of AQUASIM using: � � 32� �4� 52�� 67 $ - $ 3$89 :� *:�; :� <= (8) where r is the net oxygen production rate, K2 is the reaeration rate constant, Csat is the oxygen saturation concentration, CO2 is the oxygen concentration, I is the light intensity, d is the mean river depth, P and R are production and respiration parameters respectively, Kdeg is the degradation rate constant, KO2 is the half saturation concentration with respect to oxygen, and COD is the concentration of substrates. III. RESULTS AND DISCUSSION A. Verification of Hydraulic Model The study of Hasel River was modeled using cross-section profiles throughout river’s length, a constant slope, and an effective coefficient through every reach. All simulations were performed using an extended version of AQUASIM [3, 28-30]. The calibration of the effective Strickler friction coefficient in the four reaches of Hasel River (KSt,I, KSt,II, KSt,III, and KSt,IV) was performed according to the TLUG estimated values. The non-dimensional dispersion coefficient cf had a constant value based on the tracer transients in a water column [5]. The estimated parameters are given in Table I. Figure 1 shows the water level profile of Hasel River. Hydraulic height gradients through the four reaches and the calculated longitudinal water surface profile were in acceptable agreement. TABLE I. HYDRAULIC PARAMETERS ESTIMATES Parameter Unit Value KSt,I m 1/3 s -1 40.0 KSt,II m 1/3 s -1 35.0 KSt,III m 1/3 s -1 30.0 KSt,IV m 1/3 s -1 32.5 cf 0.006 Fig. 1. Profile of water levels of Hasel River. B. Verification of Dissolved Oxygen Model Estimates for physical reaeration were based on empirical assessments [31, 32], considering the given flow velocities and water depths for the investigated time periods. Empirical formulas resulted in K2-values ranging between 17.0 and 29.0d -1 . A value of K2=20d -1 was selected for model calculations. The oxygen saturation concentration formula is given in [33]. Based on these boundary conditions, production and respiration rate parameters were estimated from continuously measured oxygen time series (Table II). In [34], an approximate value of Kdeg was reported ranging between 0.5 and 3.0d -1 . Engineering, Technology & Applied Science Research Vol. 10, No. 4, 2020, 6126-6130 6128 www.etasr.com Khalifa: Simulation of Dissolved Oxygen and Dissolved Substrate for Hasel River TABLE II. PARAMETER ESTIMATES FROM OXYGEN TIME SERIES Parameter Unit Value Std. dev. P g/(Wd) 0.10097 0.0006 R g/(m 2 d) 11.27 0.06 Fig. 2. Oxygen time series for Hasel River. The effect of the effluents on oxygen concentration in the river is based on a degradation rate coefficient of 1.0d -1 . Figure 2 shows the concentration of oxygen at five sections of Hasel River at the boundaries of the river reaches. This method achieved a meaningful agreement between measured and calculated oxygen concentrations. Despite the enhanced gas exchange with the atmosphere, due to low depth and turbulent flow, daily oxygen amplitudes highly exceeded saturation up to 6mg/l in May-August 2001, while oxygen deficits were almost equal pronounced. These patterns were still existent in the remaining period of 2001, but in a more narrow range of oxygen concentrations. In addition, the concentrations of oxygen were validated for the mouth with data from [35]. When considering the oxygen balance in shallow eutrophic rivers, one important concern is the over-saturation which leads to gas-blisters in fishes. According to the European Classification of Water Quality in water bodies, Table III shows the trophic state of Hasel River [36-38]. As it is evident, Hasel River is moderately polluted except some locations such as near the waste water treatment plant at Suhl. TABLE III. TROPHIC STATUS FOR HASEL RIVER Cross-sectional locations* O2-saturation % Trophic state 19.30Km (Suhl) 132.53 Eutrophic 19.00Km 120.37 Mesotrophic 9.75Km 116.74 Mesotrophic 7.00Km 119.08 Mesotrophic 0.00Km (mouth) 107.92 Mesotrophic * The locations are measured US of the mouth. C. System Response to Inputs of Organic Matter The sewage treatment plant at Suhl is at x=19300m, having a mean dry weather discharge of 0.07m 3 /s and a mean COD of 77.3mg/l. The sewage treatment plant at Rohr is at x=6100m, having a mean dry weather discharge of 0.004m 3 /s and a mean COD of 50mg/l. In addition, there is a sewer overflow at x=20600m with a mean COD of 125mg/l. In order to examine the pollution load impact, oxygen concentrations were calculated at x=19200m, 100m DS the sewage treatment plant effluent at Suhl, under the assumption of constant COD concentrations of 90mg/l (upper legislation limit) combined with a sewer overflow duration of 2.4h, a discharge of 2m 3 /s, and COD concentration of 125 mg/l (about the yearly average). Figure 3 shows the actual and hypothetical load cases. The comparison shows that during a dry period COD effluent concentration leads to a slight decrease in DO concentrations, with no significant differences between the actual and hypothetical loads of COD. The decrease is between 0.3 and 0.5mg/l while daily fluctuations do not change significantly. The effect of a sewer overflow in August 26, 2001 is shown in Figure 3(b). The simulation is in agreement with the documented effects of immediate and delayed oxygen depletion in the surface flow of running waters [39]. Figure 4 shows the time series for COD and oxygen at 100m DS Suhl. The joined sewer overflow increased COD concentrations for more than 4 hours. Meanwhile, the oxygen concentrations decreased below 4mg/l for several hours, followed by a slow recovery. Therefore it is concluded that joined sewer overflows Engineering, Technology & Applied Science Research Vol. 10, No. 4, 2020, 6126-6130 6129 www.etasr.com Khalifa: Simulation of Dissolved Oxygen and Dissolved Substrate for Hasel River have the potential to disrupt oxygen balances for extended time periods. Those single events suggest a real endangerment for populations of macro-organisms and fishes [5]. (a) (b) Fig. 3. Oxygen time series at 100m DS Suhl of Hasel River: (a) actual, b) hypothetical. Fig. 4. Hypothetical time series for COD and oxygen at 100m DS Suhl. IV. CONCLUSION The proposed AQUASIM model was used to quantify the oxygen balance of a highly eutrophic shallow river, using a systematic procedure of river compartmentalization including identification of flows and a simplified description of biochemical conversion processes in Hasel River. The modeling approach quantified the temporal dynamic of oxygen in surface flow for spatially averaged scales, with consideration of distinct patterns at the boundaries of the river reaches, using physico-chemical data of 2001. The proposed model opens new perspectives in the study of ecological system responses to anthropogenic impact. The model’s potential is demonstrated by simulations showing the effect of emissions from wastewater treatment plants operating according to emissions standards, which may disrupt oxygen concentrations. Moreover, combined sewer overflows may cause severe oxygen depletion for extended time periods, endangering macro-invertebrate and fish populations. Further research is needed in ecologically meaningful river water quality modeling. In these efforts, this AQUASIM model could be utilized during system analysis at several levels of intricacy and human impact, improving the functional understanding of ecological processes in rivers and the identification of ecological effective management strategies. ACKNOWLEDGMENTS Measurements and data from the National Thuringia Town of Environment and Geology in Germany were a key factor for this study. Most of the work was performed during a research visit in 2004 at Kassel University financed through the German Research Foundation (DFG). Finally, the staff members at Kassel University, especially Prof. Dr. Borchardt, Prof. Dr. Tönsmann, Mr. Funke, and Mr. Weiss offered valuable advice, support, and great hospitality. REFERENCES [1] S. C. Chapra, Surface water quality modeling. New York, NY, USA: McGraw-Hill, 1997. [2] R. B. J. Ambrose, T. A. Wool, and J. L. Martin, “WASP-Water Quality Analysis Simulation Program, Version 5.2-MDEP, Model documentation,” Environmental Research Laboratory - AScI Corporation, Athens, Georgia, USA 2001. [3] P. Reichert, AQUASIM 2.0: computer program for the identification and simulation of aquatic systems. Dübendorf, Switzerland: Swiss Federal Institute for Environmental Science and Technology (EAWAG), 1998. [4] U. von Gunten, M. Elovitz, and H. P. Kaiser, “Characterization of ozonation processes with conservative and reactive tracers: prediction of the degradation of micropollutants,” Analusis, vol. 7, no. 25, pp. 29–31, 1997. [5] D. Borchardt and P. Reichert, “River Water Quality Model no. 1 (RWQM1): Case study I. Compartmentalisation approach applied to oxygen balances in the River Lahn (Germany),” Water Science and Technology, vol. 43, no. 5, pp. 41–49, Mar. 2001, doi: 10.2166/wst.2001.0247. [6] J. J. Kuiper, “Making eco logic and models work: an integrative approach to lake ecosystem modelling,” Wageningen University & Research, Wageningen, Netherlands, 2016. [7] Kenneth Irvine et al., “Water Framework Directive - An assessment of Mathematical Modelling in its Implementation in Ireland (2002-W-DS- 11),” Environmental Protection Agency, Wexford, Ireland, 2005. [8] Kelsey G. Jencso, Brian L. McGlynn, Michael N. Gooseff, Steven M. Wondzell, Kenneth E. Bencala, and Lucy A. Marshall, “Hydrologic connectivity between landscapes and streams: Transferring reach- and plot-scale understanding to the catchment scale,” Water Resources Research, vol. 45, Art. no. W04428, Apr. 2009. [9] M. Hrachowitz et al., “Transit times—the link between hydrology and water quality at the catchment scale,” WIREs Water, vol. 3, no. 5, pp. 629–657, May 2016, doi: 10.1002/wat2.1155. [10] S. G. Leibowitz, P. J. Wigington, K. A. Schofield, L. C. Alexander, M. K. Vanderhoof, and H. E. Golden, “Connectivity of Streams and Wetlands to Downstream Waters: An Integrated Systems Framework,” JAWRA Journal of the American Water Resources Association, vol. 54, no. 2, pp. 298–322, Mar. 2018, doi: 10.1111/1752-1688.12631. [11] K. M. Fritz et al., “Physical and Chemical Connectivity of Streams and Riparian Wetlands to Downstream Waters: A Synthesis,” JAWRA Journal of the American Water Resources Association, vol. 54, no. 2, pp. 323–345, Mar. 2018, doi: 10.1111/1752-1688.12632. [12] J. Harvey and M. Gooseff, “River corridor science: Hydrologic exchange and ecological consequences from bedforms to basins,” Water Resources Research, vol. 51, no. 9, pp. 6893–6922, Jul. 2015, doi: 10.1002/2015WR017617. [13] G. Weigelhofer, T. Hein, and E. Bondar-Kunze, “Phosphorus and Nitrogen Dynamics in Riverine Systems: Human Impacts and Management Options,” in Riverine Ecosystem Management: Science for Governing Towards a Sustainable Future, S. Schmutz and J. Sendzimir, Engineering, Technology & Applied Science Research Vol. 10, No. 4, 2020, 6126-6130 6130 www.etasr.com Khalifa: Simulation of Dissolved Oxygen and Dissolved Substrate for Hasel River Eds. Cham, Switzerland: Springer International Publishing, 2018, pp. 187–202. [14] Y. Schindler Wildhaber et al., “Effects of river morphology, hydraulic gradients, and sediment deposition on water exchange and oxygen dynamics in salmonid redds,” Science of The Total Environment, vol. 470–471, pp. 488–500, Feb. 2014, doi: 10.1016/j.scitotenv.2013.09.100. [15] C. Anibas et al., “A hierarchical approach on groundwater-surface water interaction in wetlands along the upper Biebrza River, Poland,” Hydrology and Earth System Sciences, vol. 16, no. 7, pp. 2329–2346, 2012, doi: 10.5194/hess-16-2329-2012. [16] M. V. Japitana and M. E. C. Burce, “A Satellite-based Remote Sensing Technique for Surface Water Quality Estimation,” Engineering, Technology & Applied Science Research, vol. 9, no. 2, pp. 3965–3970, Apr. 2019. [17] M. V. Japitana, M. E. C. Burce, and C. Ye, “A Geoinformatics-based Framework for Surface Water Quality Mapping and Monitoring,” Engineering, Technology & Applied Science Research, vol. 9, no. 3, pp. 4120–4124, Jun. 2019. [18] M. Iwanyshyn, M. C. Ryan, and A. Chu, “Separation of physical loading from photosynthesis/respiration processes in rivers by mass balance,” Science of The Total Environment, vol. 390, no. 1, pp. 205–214, Feb. 2008, doi: 10.1016/j.scitotenv.2007.09.038. [19] W. K. Dodds, A. M. Veach, C. M. Ruffing, D. M. Larson, J. L. Fischer, and K. H. Costigan, “Abiotic controls and temporal variability of river metabolism: multiyear analyses of Mississippi and Chattahoochee River data,” Freshwater Science, vol. 32, no. 4, pp. 1073–1087, Dec. 2013, doi: 10.1899/13-018.1. [20] R. Adrian et al., “Environmental Impacts—Lake Ecosystems,” in North Sea Region Climate Change Assessment, M. Quante and F. Colijn, Eds. Cham, Switzerland: Springer International Publishing, 2016, pp. 315– 340. [21] W. M. A. Khalifa, “Simulation of water quality for the El-Salam canal in Egypt,” presented at the Water Polution 2014, The Algarve, Portugal, May 2014, pp. 27–37, doi: 10.2495/WP140031. [22] W. M. A. Khalifa, “Evaluation of water quality parameters using numerical modeling approach for the El-Salam Canal in Egypt,” International Journal of Advanced and Applied Sciences, vol. 7, no. 2, pp. 99–112, Feb. 2020, doi: 10.21833/ijaas.2020.02.014. [23] Md. J. B. Alam, M. R. Islam, Z. Muyen, M. Mamun, and S. Islam, “Water quality parameters along rivers,” International Journal of Environmental Science & Technology, vol. 4, no. 1, pp. 159–167, Dec. 2007, doi: 10.1007/BF03325974. [24] “Thüringer Landesamt für Umwelt, Bergbau und Naturschutz.” http://tlubn.thueringen.de/. [25] B. C. Yen, “Unsteady Flow Mathematical Modelling Techniques,” in Modelling of Rivers, H. W. Shen, Ed. New York, NY, USA: John Wiley, 1979, p. 13.1-13.33. [26] U. Y. Shamir and D. R. F. Harleman, “Numerical solutions for dispersion in porous mediums,” Water Resources Research, vol. 3, no. 2, pp. 557–581, 1967, doi: 10.1029/WR003i002p00557. [27] H. B. Fischer, J. E. List, C. R. Koh, J. Imberger, and N. H. Brooks, Mixing in Inland and Coastal Waters. San Diego, CA, USA: Academic Press, 1979. [28] P. Reichert, “AQUASIM – A Tool for Simulation and Data Analysis of Aquatic Systems,” Water Science and Technology, vol. 30, no. 2, pp. 21–30, Jul. 1994, doi: 10.2166/wst.1994.0025. [29] P. Reichert, “Design techniques of a computer program for the identification of processes and the simulation of water quality in aquatic systems,” Environmental Software, vol. 10, no. 3, pp. 199–210, Jan. 1995, doi: 10.1016/0266-9838(95)00010-I. [30] “Software, Eawag” https://www.eawag.ch/de/abteilung/siam/software/ (accessed Jul. 17 2020) [31] M. Owens, R. W. Edwards, and J. W. Gibbs, “Some reaeration studies in streams,” Air and water pollution, vol. 8, pp. 469–486, Sep. 1964. [32] P. Wolf, Simulation des Sauerstoffhaushaltes in Fliessgewässern. Vol. 53. Forschungs- und Entwicklungsinstitut für Industrie- und Siedlungswasserwirtschaft sowie Abfallwirtschaft, 1974. [33] A. E. Greenberg, R. R. Trussell, and L. S. Clesceri, Standard methods for the examination of water and wastewater, 16th ed. Washington, DC, USA: APHA, 1985. [34] R. V. Thomann and J. A. Mueller, Principles of Surface Water Quality Modeling and Control. Harper & Row, 1987. [35] J. Arle and F. Wagner, “Effects of anthropogenic salinisation on the ecological status of macroinvertebrate assemblages in the Werra River (Thuringia, Germany),” Hydrobiologia, vol. 701, no. 1, pp. 129–148, Jan. 2013, doi: 10.1007/s10750-012-1265-z. [36] A. Nędzarek, A. Tórz, and J. Kubiak, “Oxygen conditions and trophic state of Lake Głębokie (Szczecin, Poland) in the years 2008-2010,” Limnological Review, vol. 10, no. 3–4, pp. 163–172, Jan. 2010, doi: 10.2478/v10194-011-0019-z. [37] Z. Xu and Y. J. Xu, “Determination of Trophic State Changes with Diel Dissolved Oxygen: A Case Study in a Shallow Lake,” Water Environment Research, vol. 87, no. 11, pp. 1970–1979, 2015, doi: 10.2175/106143015X14362865226716. [38] H. Siwek, M. Włodarczyk, and R. Czerniawski, “Trophic State and Oxygen Conditions of Waters Aerated with Pulverising Aerator: The Results from Seven Lakes in Poland,” Water, vol. 10, no. 2, Art. no. 219, Feb. 2018, doi: 10.3390/w10020219. [39] P. Harremoës, “Immediate and delayed oxygen depletion in rivers,” Water Research, vol. 16, no. 7, pp. 1093–1098, Jan. 1982, doi: 10.1016/0043-1354(82)90124-5. AUTHORS PROFILE Walid M. A. Khalifa is an Assistant Professor at Hail University, Saudi Arabia, and Fayoum University, Egypt. He received his B.Sc. in Civil Engineering and his Ph.D. in Water Resources and Environmental Hydrology at Cairo University. His research interests include Water Quality and Hydrodynamics Modeling, Pressurized Irrigation Modeling, and Design of RC Water Structures.