Pollutant Flux Releases During Summer Monsoon Period based on Hydrological Modeling in Two Forested Watersheds, Soyang Lake
Pollutant Flux Releases During Summer Monsoon Period based on Hydrological Modeling in Two Forested Watersheds, Soyang Lake
Environmental Engineering Research. 2009. Mar, 14(1): 13-18
Copyright ©2009, Korean Society of Environmental Engineering
  • Received :
  • Published : March 31, 2009
Export by style
Cited by
About the Authors
S. H. Kang

In this study, specific pollutant releases during the Asian monsoon season were estimated and the information was applied to the non-point pollutant sources management from two forested watersheds of the Soyang Lake. The two watersheds are part of the 2,703 km 2 Soyang Lake watershed in the northern region of the Han River. The outlets of the two watersheds were respectively analyzed for continuous water quality concentration and for discharge during various single rainfall events. Statistical power function methods are utilized to compare stream discharge and pollutant flux release during the study period. Based on the monitoring data during the study period, the specific load flux method using simulated discharge was conducted and validated in the two watersheds. The model predictions corresponded well with the measured and calculated pollutant releases. The modeling approach taken in this study was found to be applicable for the two forested watersheds.
1. Introduction
Since the pioneering research that involved estimating the sediment transport rate for predicting the life span of a dam in a river, various rating curves have been used to estimate pollutant release from a watershed. 1 , 2) The reason we need an empirical relationship is because most of the rivers have not been gauged or stream flow and suspended sediment flux data for these rivers is not readily available. Vogel et al. 3) reported that the amount of pollutant constituents including sediments from river flows depends on the form of the relationship between the effective river discharge and the pollutant flux as well as on the form of the frequency distribution of the river discharges based on over 300 journal articles. Cooper et al. 4) examined how opportunistic sampling in addition to weekly sampling of dissolved organic matter affected regression and ratio based on estimates of loading. Furthermore, Christiansen et al. 5) used a real time water quality monitoring and regression model to estimate the constituent concentration, loads, and yields in the Little Arkansas River in South Central Kansas. The study revealed that these indicators yielded errors greater than 25 percent. However, many of these studies have been carried out either in watersheds in terms of sustainable urbanized effects on surface water quality 6- 9) or on land use effects with an emphasis on the observation of stream flow. 10 12) Furthermore, the study of a large scale (>100 km 2 ) and non-impacted watershed for controlling a Non Point Source (NPS) has seldom been analyzed, 13) despite the longterm studies carried out in these fields. The process of determining the daily pollutants load that runs into a lake is both challenging and complex. Many lakes receive a high portion of their pollutant loading from non-point sources, which are highly variable and are difficult to quantify. Lakes manifest pollutant loading on a seasonal scale, not on a daily basis. The loading during the winter months may have little effect on summer algal densities. While non-point pollutant loads tend to flow out with extreme rainfall, forming a watershed, little information is available on how to control the flow effectively.
In this study, field monitoring in undisturbed watersheds are examined and an evaluation is performed of the characteristics of specific pollutant release generating flow discharge by using the hydrological modeling approach. The study consists of two stages. In the first stage, the hydrological processes for stream discharge were estimated by comparing the observed data. The second stage of this work was presented with a model validation to estimate the specific pollutant release based on the calculated flow discharge during the rainy season. Originally, this study was designed on the basis of the long-term rainfall phenomenon over 30 years in the two non-impacted watersheds of SoyangLake.
2. Method and Tools
- 2.1. Watershed Description
The Soyang reservoir is an artificial dam located about 100 km to the east of Seoul. The reservoir is the deepest and largest reservoir in Korea; it was constructed on the North Han River in 1973 for flood control of down stream areas and water supply for Chunchon city including metropolitan cities. The hydraulic retention time is the longest (0.7 year) in Korea because the outflow rate is small, considering the reservoir’s large volume (Table 1) . The water depth varies from about 110 m on the outer part of the Dam to 30 m on the inner part and about 42 m on average. The Watershed of Soyang Lake has two main subwatersheds: the Inbuk River (578 km 2 watershed area), of which part of the watershed is included in the north Kangwon province, Democratic People’s Republic of Korea (North Korea), and the Naelin River (1,075 km 2 watershed area). The elevations above sea level ranged from 215 m to 1,348 m for the Inbuk River watershed and from 201 m to 1,567 m for the Naelin River watershed, which rapidly changed to steep slopes.
Hydrological characteristics and land use in watershed of Soyang reservoir
PPT Slide
Lager Image
Hydrological characteristics and land use in watershed of Soyang reservoir
The two watersheds consist of over 60% of the total watershed area. The water quality conditions in the lake maintain an oligotrophic to mesotrophic state. The water blooms with cyanobacteria, mainly Anabaena spp., which has also been annually reported in some tributaries of the lake. 14 , 15)
- 2.2. Methods
Data obtained during the rainy season over the past 2 years, from 2004 to 2005, for a rating curve with a long term water quality were analyzed when available. The annual mean precipitations of the Inbuk River watershed and the Naelin River watershed were 1,083 and 1,118 mm, respectively, and over 80% of the annual mean rainfalls in the two main watersheds of the Soyang River were concentrated in heavy showers that occurred several times during the rainy season from May to September. This rainfall phenomenon results in a deposit of pollutive components during dry seasons.
The higher rainfall that occurs from May to September will mobilize the build-up of contaminants, creating a larger discharge. The basic water quality parameters such as Biochemical Oxygen Demand (BOD), Chemical Oxygen Demand (COD), SS, Total Nitrogen (T-N), and Total Phosphorus (T-P) were chosen for the derivation of the flux rating curve, because they are major indicators for the estimation of organic, turbid, and nutritive conditions, respectively. The samplings in the two watersheds are concentrated in single rainfall events, because the samples in the period tend to be sufficient to characterize the seasonal pollutant load. A limited number of samples are also monitored during low discharge periods to allow the pollutant load with a baseflow. This method assumes that river flows vary significantly and that concentrations increase with flow during storm runoff events. Based on the observed data, experimental relations such as the pollutant regression curve are applied. The most commonly used basic pollutant rating curve is the power function: 16)
PPT Slide
Lager Image
where C is the concentration (mg/L), Q is the discharge (m 3 /s), and a and b are the regression coefficients. Large scale watershed models for NPS lack a variable source area hydrology routine for stream flow, and are unable to identify spatially dynamic runoff. In this regard, hydrological efforts to evaluate stream discharge will be described in the next session.
An important issue for the water resource manager is the determination of a specific pollution flux on a daily basis. If continuous concentration and flow measurements could be made for the entire study period, the mass flux could be estimated by multiplying by the appropriate unit conversion. For these comparisons, load (kg/day) L , which is integrated over the study period and given as a day of observations, is defined by the following Eq. (2);
PPT Slide
Lager Image
where C is the concentration (mg/L), Q is the discharge (m 3 /s), and F is the conversion factor. In order to summarize the specific load, Eq. (3) uses a byproduct of load estimates, normalizing them to the monitored watershed area, A (km 2 ). The specific pollutant flux which is a modified power function equation was used as follows:
PPT Slide
Lager Image
In Eq. (3), Q is not provided as a continuous measurement due to various reasons including cost restrictions. In order to capture the seasonal variation of a specific pollutant release during the rainy period, continuous samples of single rainfall events in two watersheds of the Soyang Lake were corrected and validated considering the antecedent rainfall, the rainfall intensity, and the term of the sampling periods. In order to assess the difference in efficiency between the observed and the predicted rating curve, the Root Mean Square Error (RMSE) method was used according to the following formula:
PPT Slide
Lager Image
where n: the number of observations.
In an ungauged watershed, the first step for estimating NPS is to generate stream discharge. The TOPMODEL selected for this study is a rainfall-runoff model in which distributed predictions of watershed responses are made based on the hydrological similarity of any point in a watershed. TOPMODEL is a topographically based hydrological model which reproduces the hydrological behavior of watersheds in a semi-distributed way. 17-19) Since a detailed explanation of TOPMODEL has been reported elsewhere, only a brief description will be provided as follows: The simplicity of TOPMODEL is in the topographic index, the ln(α/ tanβ) index, as an index of hydrological similarity derived from the digital elevation models as shown in Fig. 1 , where α is the cumulative drainage though a point, and the hydraulic gradient of the saturated zone can be approximated by the local surface topographic slope, tanβ,
PPT Slide
Lager Image
where T 0 represents the local value of the saturated transmissi-
PPT Slide
Lager Image
Spatial distribution of topographic index in two catchments.
vity of the soil. β i is the slope angle, S i is the local storage deficit, and m is the change in transmissivity with depth. The local deficit is defined as:
PPT Slide
Lager Image
where S is the average storage deficit and A is the local average area of ln(α / tan β).
The vertical drainage q v from the unsaturated store at any point i is controlled by the local saturated zone deficit S i , which depends on the depth of the local water table:
PPT Slide
Lager Image
Where S UZ is the storage in the unsaturated zone and t d is a time delay constant. The evaporation for the TOPMODEL description is allowed at the full potential rate for water draining freely in the unsaturated zone and for predicted areas of surface saturation. When the gravity drainage zone is exhausted, actual evapotranspiration E a as a function of potential evaporation E p may continue to deplete the root zone store at the rate E a given by:
PPT Slide
Lager Image
where S rz and S r max are the root zone storage deficit and the maximum allowable storage deficit, respectively. The watershed flux of water entering the water table, Q v is calculated by summing the q v of each topographic index class:
PPT Slide
Lager Image
The output from the saturated store is represented by the baseflow term, Q b which can be calculated using a subsurface storage deficit-discharge function of the form:
PPT Slide
Lager Image
where Q 0 = Ae is the discharge in the case when S is zero.
3. Results and Discussion
The single event rainfall amount for two watersheds ranged from 24 mm to 198 mm. The fitted relationships between specific discharge and pollutant flux releases using Eq. 3 are summarized in Table 2 . Most of the regressions for the two watersheds during the rainy season that are given in Table 2 have positive values. This may be due to the discharge through the wash off process in which rainfall results in pollutant deposition during the light rainy season from January to May or October to December. In certain cases, the regression values were high during the consistently short sampling period term in the cases of the 2 nd and 3 rd rainfall events or the same monthly data in the cases of the 3 rd and 4 th rainfall events in the Inbuk River watershed. 20) It has been reported in a number of studies that runoff quality during the initial storm of the season has higher pollutant concentrations than in later rainfall events. Lee et al. 21) defined this phenomenon as seasonal first flush. Also, their study reported that peak pollutant concentrations tend to precede peak runoff flow rates and thus result in the runoff of a greater pollutant release during early storm stages. However, in this study, a seasonal first flush phenomenon of specific pollutant fluxes in two watersheds did not occur for selected indicators. This phenomenon of pollutant release is important in the management of the pollutant load from the diffused watershed. As a higher pollutant mass load is delivered during the rainy season, controlling the entire pollutant load during the rainy season will be more effective than trying to treat the pollutant mass of all storms throughout the season. Because the months of July and August usually have the greatest rainfall, the bulk of the discharge in the two river watersheds is expected to occur during the rainy season. Conversely, the long dry period from April or May to October allows contaminants to build up in these areas.
The watershed-scale hydrological model for pollutant flux release could assist environmental managers to focus their limited discharge data set. In order to estimate and validate pollutant load flux release, two single rainfall events were selected for the Inbuk River and the Naelin River respectively. Fig. 2 illustrates the observed pollutant flux using Eq. (3) versus the calculated pollutant flux from the predicted discharge using distributed TOPMODEL in two rivers and demonstrates the close relationship of over 80% of all the selected parameters. This modeling approach allows seasonal pollutant fluxes to be produced effectively, when the discharges in the insufficient rating curve areas are generated well. Fentie et al. 22) reported that the relative differences between total suspended sediment loads using rating curves varied depending on the scale of the watershed areas. However, further data sets are required to confirm the significant relationship between them. The magnitude of RMSE values depends on the variances between a true discharge and pollutant fluxes distribution. The values of five indicators including discharge were 0.14 for the Inbuk River and 1.58 for the Naelin River, which implies lower RMSEs. The calculated results indicate that the pollutant releases that control discharge at the large scale non-urbanized forest area are well captured by the hydrological-link model. This approach can also simulate the trend as well as the variability in the daily or monthly pollutant release at any point of outlet of the two watersheds during the rainy season.
4. Conclusions
The rainfall phenomenon is an important factor in the control of pollutant load from diffuse areas. The procedure demonstrated in this study is a frame approach for simulating specific pollutant flux releases in two forested watersheds where detailed river flow data is not available during the rainy season. Extensive water sampling data in two forested watersheds was collected to derive specific pollutant flux regression curves during various rainfall events. For simulating a specific pollutant release, a conceptual model was applied in the large forested watersheds. The results of this study indicate that the relationship between discharge and specific pollutant load have a high regression coefficient. The flows during the rainy season provide most of the pollutant fluxes in this study, providing evidence of the important link between discharge and pollutant load. The five pollutant load rating curves for this study were developed using continuous sampling data during the heavy rainfall season. The simulation results indicate that the proposed approach can be used to predict both the discharge and pollutant release from the forested watershed, while the values of RMSE’s are 0.14 for the Inbuk River watershed and 1.58 for the Naelin River watershed, respectively. This approach allows the estimation of other pollutive indicators following the same procedure. Concerning the treatment of pollutant loads from diffused watersheds, there is no agreement about at which level the contaminated waters should be removed. This result, considering the long term rainfall phenomenon, suggests that treating pollutant fluxes from a
Correlation coefficients with pollutant rating coefficient a and exponent b for indicators in Inbuk River, 2004 and Naelin River, 2005
PPT Slide
Lager Image
Correlation coefficients with pollutant rating coefficient a and exponent b for indicators in Inbuk River, 2004 and Naelin River, 2005
PPT Slide
Lager Image
The predicted specific pollutant fluxes versus observed pollutant fluxes based on calculated discharge (o: observed,? : calculated in Inbuk-R watershed; x: observed, ... : calculated in Naelin-R watershed, respectively).
watershed during the rainy season may be more effective than treating a randomly released load during the dry season. Though further modeling testing based on other field data is needed, the results shown above will be helpful for controlling the pollutant flux release from forested watersheds during the rainy season.
The author would like to thank Prof. Keith Beven of the University of Lancaster (UK) for providing various reports and the TOPMODEL computer code. In addition, part of this research was supported by a grant from the Eco-STAR Project (07-Ⅰ-3) of the Korean Ministry of Environment.
Komai Y , Umemoto S , Inoue T , Imai A 2003 Evaluation of annual loading of major ionic species and nutrients in forested watershed Japan Diffuse Pollution Conference Dublin part 5 13 - 18
McKee L. J , Ganju N. K , Schoellhamer D. H 2006 Estimates of suspended sediment entering San Francisco Bay from the Sacramento and San Joaquin Delta San Francisco Bay California Journal of Hydrology 323 (1) 335 - 352    DOI : 10.1016/j.jhydrol.2005.09.006
Vogel R. M , Stedinger J. R , Hooper R. P 2003 Discharge indices for water quality loads Water Resour Res SWC 39 (10) 1 - 9
Cooper D. M , Watts C. D 2002 A comparison of river load estimation techniques: application to dissolved organic carbon Environmetrics 13 (7) 733 - 750    DOI : 10.1002/env.525
Christiansen V. G , Rasmussen P. P , Ziegler A. C 2000 Regression analysis and real time water quality monitoringto estimate constituent concentration loads and yields in the Little Arkansas River South-central Kansas 1995-99:US Geological Survey Water Resources Investigations Report 00-4126
Endreny T. A , Hassett J. M 2005 Robustness of pollutant loading estimators for sample size reduction in a suburban watershed Intl J River Basin Management 3 (1) 53 - 66
Ide J , Nagafuchi O , Kume A , Otsuki K , Ogawa S 2003 Runoff nutrients from an afforested watershed Chamaecyparis Obtusa during rain events Diffuse Pollution Conference Dublin part 5 30 - 34
Solo-Gabriele Helena M 1998 Generation of long-term record of contaminant transport Journal of Environmental Engineering 124 (7) 619 - 627    DOI : 10.1061/(ASCE)0733-9372(1998)124:7(619)
Tagg A. F , Millington R. J 2003 Development of an integrated catchment planning system for the urban environment Diffuse Pollution Conference Dublin part 4 18 - 23
Horowitz A. J , Elrick K. A , Smith J. J 2001 “Estimatingsuspended sediment and trace element fluxes in large riverbasins: methodological considerations as applied to theNASQAN program” Hydrological Processes 15 1107 - 1132    DOI : 10.1002/hyp.206
Inoue T , Ebise S , Mastui Y , Matsushita T 2003 Evaluation of organic pollutant and nutrient loading in a rural river Diffuse Pollution Conference Dublin part 5 19 - 22
Novotny V 1999 Integrating diffuse/nonpoint pollutant control and water body restoration into watershed management JAWRA Journal of the American Water Resources Association 35 (4) 717 - 727    DOI : 10.1111/j.1752-1688.1999.tb04169.x
Inamdar S 2006 Challenges in Modeling Hydrologic and Water Quality Processes in Riparian Zones JAWRA Journal of the American Water Resources Association 42 (1) 5 - 14    DOI : 10.1111/j.1752-1688.2006.tb03818.x
Park S. J , Lee E. J , Lee D. H , Lee S. H , Kim S. J 1995 Spectrofluorometric assay for rapid detection of total and fecal coliforms from surface water Appl Environ Microbiol 61 (5) 2027 - 2029
Kim B. C , Choi K. S , Kim C. G , Lee U. H , Kim Y. H 2000 Effects of the summer monsoon on the distribution and loading of organic carbon in a deep reservoir Lake Soyang Korea Water Research 34 (14) 3495 - 3504    DOI : 10.1016/S0043-1354(00)00104-4
Asselman N. E. M 2000 Fitting and interpretation of sediment rating curves Journal of Hydrology 234 (3-4) 228 - 248    DOI : 10.1016/S0022-1694(00)00253-5
Beven K. J 1997 TOPMODEL: A critique Hydrological Processes 11 (9) 1069 - 1087    DOI : 10.1002/(SICI)1099-1085(199707)11:9<1069::AID-HYP545>3.0.CO;2-O
Kirkby M. J 1997 TOPMODEL: A personal view Hydrological Processes 11 (9) 1087 - 1097    DOI : 10.1002/(SICI)1099-1085(199707)11:9<1087::AID-HYP546>3.0.CO;2-P
Campling Paul , Gobin Anne , Beven Keith , Feyen Jan 2002 Rainfall-runoff modelling of a humid tropical catchment: the TOPMODEL approach Hydrological Processes 16 (2) 231 - 253    DOI : 10.1002/hyp.341
Kim L. H , Kayanian M , Lau S. L , Stenstrom M. K 2003 A new modeling approach for estimating first flush metal mass loading Diffuse Pollution Conference Dublin part 4 38 - 46
Lee H. J , Lau S. L , Kayhanian M , Stenstrom M. K 2004 Seasonal first flush phenomenon of urban stormwater discharges Water Research 38 (19) 4153 - 4163    DOI : 10.1016/j.watres.2004.07.012
Fentie B , Joo M , Yu B , Hunter H , Marsh H , Carroll C , Dougall C 2000 Comparison of Mean Annual SuspendedLoads Estimated by the Sednet Model and Rating Curve in the Fitzroy Catchment Australia International Congress on Modelling and Simulation Advanced and Applicationsfor Management and Decision Making Melbourne December 1133 - 1139