Assessing the Influence of Nonpoint Source and Discharge Changes on Water Quality in a Tidal River Estuary Using a Three-dimensional Model

The change of nonpoint source and freshwater discharge as a result of land use change might affect the water quality in a river. In the current study, a coupled three-dimensional hydrodynamic and water quality model was created and applied to the Danshuei River estuarine system and its adjacent coastal ocean. The hydrodynamic and water quality models were validated using observations of water surface elevation, salinity distribution and water quality state variables. The predictions of hydrodynamics, salinity, dissolved oxygen and nutrients from the model simulation quantitatively agreed with the observational data. The validated model was then used to investigate the possible effects of the nonpoint source and freshwater discharge changes at the upstream reaches on water quality conditions in the Danshuei River estuarine system. Three scenarios were investigated to predict the dissolved oxygen using model simulations. The simulated results indicated that increasing nonpoint sources at the upstream reaches degraded the dissolved oxygen under low flow conditions. However, increasing freshwater discharges at the upstream reaches would overcome the loadings of the nonpoint source, which would result in increasing the dissolved oxygen in the tidal river estuary. The model can provide a useful tool for developing management practices for nonpoint sources to protect the water quality in the estuarine system.


Introduction
An increasing population, changes in urbanization and agricultural practices and structural landscape transformations have considerably modified the natural biogeochemical cycles of essential bio-elements on both a local and large scale [1][2][3].Surface water is controlled by natural processes, such as weathering and precipitation inputs, by anthropogenic activities, such as industrial effluents and wastewater treatment facilities, and by diffuse inputs, such as runoffs from urban areas and farms.Several studies demonstrate that surface water quality has severely deteriorated in numerous countries over the past few decades because of poor land use, which is indicated by a strong relationship between the declining water quality and the increasing development of the catchment scale [4,5].
Water quality is generally linked to land use/land cover (LULC) in catchments [6], and studies have focused on their relationship with water quality variables, such as dissolved salts, suspended solids and nutrients [7][8][9][10].These studies discovered that agricultural and urban land use strongly affect nutrients in the river water.Different methodologies have been applied to assess the effects of changes in the LULC on biogeochemical cycles, biodiversity and water quality [11][12][13][14][15].For example, Wang et al. [16] assessed the spatial-temporal water quality to identify the LULC sources of water pollution.Kibena et al. [17] combined in situ measurements and remote sensing techniques to assess the impacts of land use activities on the water quality of the Upper Manyame River in ZimbabweBu et al. [18] used land use types and landscape matrices, as well as statistical and spatial analyses to determine the relationships between land use patterns and river water quality in the Taizi River Basin in China during the dry and rainy seasons of 2009.However, measuring the LULC and its rates and patterns of change correlated with the water quality are traditional measures.Erol and Randhir [19] used a watershed ecosystem model to assess the impact of land use on water quality in the Lake Egirdir watershed, Turkey.Recently, different kinds of models have been widely adopted to evaluate the influences of LULC on water quality.
Estuaries are complex ecosystems that can include wetlands and mangrove swamps along their intertidal shores.They serve as a meeting point between land and sea and occupy specific ecological niches of considerable importance.However, because of human activities, adverse effects on the estuarine ecosystem have been widely reported [20,21].Numerical water quality models are useful in helping to understand the biological processes, nutrient loading and managing water quality conditions in aquatic systems.Deterministic models for water quality conditions are based on mass balance equations for dissolved and particulate substances in the water column, which consist of physical transport processes and biogeochemical processes.Information on the physical processes is typically obtained by applying hydrodynamic models.Once calibrated and verified against the observational data, the deterministic models then can be applied to explore a number of issues that can provide information for management practices [22,23].
Lin et al. [24] assessed the land use changes and their effects on land use patterns and hydrological processes at the Wu-Tu watershed in the upstream of the Keelung River Basin, but did not focus on the water quality issues.Because of urban development, the LULC in the Danshuei River watershed of northern Taiwan has been significantly changed, which resulted in an increase in the nonpoint source and discharge.A coupled three-dimensional hydrodynamic and water quality model (SELFE-WQ) was implemented to probe the influence of the increases in nonpoint source and freshwater discharge at the upstream reaches on water quality in the estuarine system.The model was validated by observations of the water surface elevation, salinity and water quality state variables.Then, the validated water quality model was used to predict the water quality for different scenarios by varying the inputs of the nonpoint source and freshwater discharge under low flow conditions.

Description of Study Area
The Danshuei River, along with its tributaries, is the largest river system in northern Taiwan.Its watershed encompasses 2726 km 2 , and it has a combined length of 158.7 km.This river consists of three major tributaries: the Tahan Stream, the Hsintien Stream and the Keelung River (Figure 1).The mean river discharges are 38.99 m 3 /s, 69.72 m 3 /s and 25.02 m 3 /s in the Tahan Stream, the Hsintien Stream and the Keelung River, respectively.In addition to the primary stream of the Danshuei River, the lower reaches of the three major tributaries are also affected by the tide.The principal tidal constituents of the estuary lean toward semi-diurnal tides, with a mean tidal range of 2.1 m and a spring tidal range of 3.5 m.Seawater intrusion reaches into all three tributaries, except during the period of extremely high river inflows.Generally, saltwater intrusion reaches 25-30 km from the mouth of the Danshuei River.The hydrodynamic characteristics in the system are primarily controlled by the tide, river inflow and density gradient induced by the mixing of seawater and freshwater [25,26].
The Danshuei River flows through the metropolitan area of Taipei, which has a population of approximately 6 million.A large amount of treated and untreated domestic sewage was discharged into the river system and resulted in low dissolved oxygen and high nutrient concentrations.Viable biological activities were observed only in the lowest reaches of the estuary, where the pollutant concentrations are reduced as a result of being diluted with sea water.

Hydrodynamic Model
The numerical modeling of ocean circulation at scales ranging from estuaries to ocean basins is becoming a developed field.Most modern oceanic and estuarine circulation codes solve for various forms of the three-dimensional Navier-Stokes equations and can be complemented with conservation equations for a given water volume and salt concentration.In this report, a three-dimensional, semi-implicit Eulerian-Lagrangian finite-element model (SELFE, [27]) was implemented to simulate the Danshuei River estuarine system and its adjacent coastal sea.The SELFE model solves the Reynolds stress-averaged Navier-Stokes equations, which use conservation laws for mass, momentum and salt with hydrostatic and Boussinesq approximations, to determine the free-surface elevation, three-dimensional water velocity and salinity.The SELFE model addresses the advection in transport equations using the Eulerian-Lagrangian method.The transport equations are solved at nodes and side centers along each vertical column using a finite element method, with the mass matrix lumped to minimize numerical dispersion.In the vertical direction, models based on terrain-following coordinates suffer from hydrostatic inconsistencies.This follows from the fact that the terrain-following coordinates do not conform to the geopotential, and the pressure gradient is evaluated as the difference between two large components that tend to cancel each other out, thus producing large discretization errors.Furthermore, these inconsistencies can be viewed as a result of efficiently evaluating the pressure at a grid point using extrapolation when steep bathymetric slopes are presented.Several remedies have been used to prevent this result, including evaluating the gradient along the Z-coordinate or applying high-order schemes [28].In the current model, the use of a hybrid coordinate system in the vertical direction effectively prevents the hydrostatic inconsistency, because the Z-coordinates used in the deeper part of the vertical axis serve to "stabilize", while the S-coordinates are used in the upper part of the water column.
Because turbulent mixing plays a critical role in determining the stratification in the tidal estuary, several reports have documented the modeling results of turbulence mixing parameterizations.The SELFE model uses the generic length scale (GLS) turbulence closure of Umlauf and Burchard [29], which has the advantage of encompassing most of the 2.5-equation closure model (K − ψ).A detailed description of the turbulence closure model, the vertical boundary conditions for the momentum equation, the numerical solution methods and the numerical stability can be found in Zhang and Baptista [27].

Water Quality Model
The water quality model used in this study was based on a three-dimensional conventional water quality analysis simulation program, called WASP5, which was originally developed by Ambrose et al. [30].It is comprised of a complicated system of four interacting parts: dissolved oxygen concentration, nitrogen cycle, phosphorus cycle and phytoplankton dynamics.Eight water quality components are included: dissolved oxygen (DO), phytoplankton as carbon (PHYT), carbonaceous biochemical oxygen demand (CBOD), ammonium nitrogen (NH 4 ), nitrate nitrogen (NO 3 ), organic nitrogen (ON), ortho-phosphorus or inorganic phosphorus (OPO 4 ) and organic phosphorus (OP).
A mathematical formulation of the conservation of mass can be written as follows: where C is the concentration of the water quality components, u, v and w are the water velocity components corresponding to a Cartesian coordinate system (x, y, and z), A h and K v are the coefficients of the horizontal and vertical eddy diffusion, respectively, and S c is a function that represents the internal source or sink of the water quality component.

Model Implementation
In the current study, the bathymetry and topography data of the Taiwan Strait and the Danshuei River estuarine system were obtained from the Ocean Data Bank and Water Resources Agency in Taiwan.The deepest point within the study area is 110 m (below mean sea level) near the northeast corner of the computational domain (Figure 2).The model mesh for the Danshuei River estuarine system and its adjacent coastal sea consists of 5119 elements (Figure 2).To meet the accuracy requirements, a fine-grid resolution was used locally, while a coarse resolution was implemented away from the region of interest.In this computational domain, the mesh size varied from 6000 m in the Taiwan Strait to 40 m in the upper reach of the Danshuei River estuarine system.In the vertical direction, ten Z-levels and ten evenly spaced S-levels were specified at each horizontal grid, i.e., the thickness of the cell depended on the bottom elevation of each grid.A 120-second time step was used in our simulations without any signs of numerical instability.

Tidal Elevation
The calibrated and verified model procedures were performed for daily freshwater discharges at the upstream boundaries in the Tahan Stream, Hsintien Stream and Keelung River in 2010.The model was run for a one-year simulation.Before the one-year simulated run, the model was spun up fifteen days to reach an equilibrium state.These conditions served to investigate the model response to the interactions of tidal forcing and varying river discharge.First, the coupled hydrodynamic and water quality model was calibrated with hydrodynamics (i.e., water surface and salinity) and then was calibrated and verified with water quality variables.The manual calibration procedure was conducted in the model simulation.The comparison between the model simulation and the measured water surface elevations is presented in Figure 3 at the following tidal gauge stations: Danshuei River mouth, Tu-Ti-Kung-Pi, Taipei Bridge, Hsin-Hai Bridge, Chung-Cheng Bridge, Ba-Ling Bridge and Ta-Chih Bridge (the locations are provided in Figure 2) for the period of 27 February to 8 March 2010.Generally, the modeling results accurately reproduced the water level variations.The mean absolute errors and the root mean square errors of the difference between the measured hourly surface elevations and the computed surface elevations during the period of 27 February to 8 March 2010, are provided in Table 1.Through the model validation with tidal elevation, a constant bottom roughness height (z o = 0.005 m) is used in the model simulation.

Salinity Distribution
The salinity distributions reflect the combined results of all processes, including density circulation and mixing processes, which, in turn, control the density circulation and modify the mixing processes.In the current study, the salinity distribution along the Danshuei River-Tahan Stream, which was collected by the Water Resources Agency in Taiwan, was used for model validation.The salinities of the open boundaries in the coastal sea were set to 35 ppt.Furthermore, the upstream boundary conditions at the three tributaries (Tahan Stream, Hsintien Stream and Keelung River) were specified with daily freshwater discharges; thus, the salinity at the upstream boundaries was set to 0 ppt.The simulated salinity distribution compared favorably with the salinity measurements along the Danshuei River-Tahan Stream during the flood and ebb tides on 26 November 2010, as indicated in Figure 4.It should be noted that the salinity field data were measured at 0.5 m below the water surface, and the simulated salinity was instead presented as the top layer and the bottom layer.Because ten Z-levels and ten evenly spaced S-levels were specified in the vertical direction for model simulation, the water depths at the top and bottom layer were different along the Danshuei River mouth.The absolute mean error and the root mean square error are 2.71 ppt and 3.72 ppt during the flood tide, respectively, and 0.49 ppt and 0.67 ppt during the ebb tide, respectively.The simulated results indicate that the model performs better for the ebb tide compared to that of the flood tide.

Water Quality State Variables
Adjacent to the metropolitan city of Taipei, it is believed that the spatial trend of the deteriorated water quality is primarily attributed to the wastewaters directly discharged into the river channels [31].Successive efforts were made to estimate the point source loadings that Montgomery Watson Harza (MWH) [32] used in the following water quality simulations.The daily freshwater discharges from 2010 were applied at the upstream boundaries at the Tahan Stream, Hsintien Stream and Keelung River.The five-constituent tide used to generate the time series tidal level was used at the ocean boundaries.The concentrations of the water quality state variables, including ammonium nitrogen, total nitrogen, total phosphorus, carbonaceous biochemical oxygen demand, dissolved oxygen and chlorophyll a, at the river and ocean boundaries were established based on the monthly measurement performed by the Taiwan Environmental Protection Administration (TEPA).
Two measured datasets were collected on 1 March and 1 June 2010, used for model calibration, while the other two datasets were collected on 3 September and 2 December 2010, adopted for model verification.The model parameters were initially estimated from the literature [33].These parameters were adjusted and tuned until a reasonable reproduction of the field data was obtained at the observation stations.Table 2 indicates the coefficients adopted for water quality simulations.Because of the constraint of space, thee longitudinal water quality distributions predicted by the water quality model on 1 March 2010, for the Danshuei River-Tahan Stream, Hsintien Stream and Keelung River are provided only in Figure 5 to Figure 7, respectively.The observational stations shown in Figures 5-7 refer to Figure 2. The water quality distributions of the dissolved oxygen, carbonaceous biochemical oxygen demand, ammonium nitrogen, total phosphorus and nitrate nitrogen concentrations at the top and bottom layers along the river channels are presented in the figures along with the observations at the monitoring stations.Because there were no measured nitrate nitrogen data, the simulated results at the top and bottom layers are only presented in Figure 5e. Figure 5a indicates that both the concentrations predicted by the model and the observed dissolved oxygen concentrations along the Danshuei River-Tahan Stream depict a decreasing trend from the Danshuei River mouth to the Hsin-Hai Bridge.In the lower estuary, the dissolved oxygen concentrations increase toward the river mouth as a result of seawater dilution.The lowest dissolved oxygen concentration occurs at the Hsin-Hai Bridge.The carbonaceous biochemical oxygen demand, ammonium nitrogen and total phosphorous all indicate the same spatial trends along the river channel from the Danshuei River to the Tahan Stream (Figures 5b-d).The concentrations increase from the mouth of the Danshuei River, reach a maximum at the Hsin-Hai Bridge and then gradually decrease towards the Fu-Chou Bridge.The maximum concentrations of these three water quality state variables all occur at the Hsin-Hai Bridge, which results in the lowest dissolved oxygen.Figures 5a-d indicates that the model generally captured the spatial trends of the observed longitudinal distributions, but could not mimic the observed dissolved oxygen concentration at the lower reaches.
Figure 6 illustrates the longitudinal water quality distributions along the Hsintien Stream.Both the observed data and the model predictions indicate that the water quality conditions degrade as one approaches the river mouth, where it joins the main stem of the Danshuei River, the dissolved oxygen decreases and the organic carbon, ammonium nitrogen and total phosphorus increase (Figures 6a-d).The model accurately captures the observed carbonaceous biochemical oxygen demand, ammonium nitrogen and total phosphorous along the Hsintien Stream.The longitudinal water quality distribution along the Keelung River (Figure 7) indicates that there were significant nutrient loadings discharged into the river section near the Bai-Ling Bridge, which is where the lowest dissolved oxygen in the Keelung River was observed.The model cannot actually mimic this lowest observed dissolved oxygen concentration, but does accurately match the observed carbonaceous biochemical oxygen demand, ammonium nitrogen and total phosphorous quite well along the Keelung River (Figures 7a-d).The comparison of simulated and measured nitrate nitrogen concentration is not possible, because there were no measured nitrate nitrogen data.The model performance with statistical errors, including the absolute mean error and the root mean square error, on 1 March, 1 June, 3 September and 2 December 2010, are illustrated in Tables 3-6, respectively.

Model Application: System Response to Nonpoint Source and Discharge Changes
The Danshuei River watershed can be classified into nine land uses, which include forest, agriculture, stream and reservoir, residential, transportation, grass land, bare soil, bush, and others.The percentages of 89% and 4.5% are occupied by forest and agriculture, respectively.To examine the system's response to a change of nonpoint source loadings and discharges from the upstream watershed, the validated model was used to predict the dissolved oxygen distribution under low flow conditions.To perform the model prediction, the five-constituent tide at the ocean boundaries was used to force the model simulation.The concentrations of the water quality state variables, i.e., ammonium nitrogen, total nitrogen, total phosphorus, carbonaceous biochemical oxygen demand, dissolved oxygen and chlorophyll a, at the river and ocean boundaries were established based on the mean values calculated from the measured water quality data, which were observed and collected from 2003 to 2013 by the TEPA (Table 7).Because the low flow could have significant influence on the water quality condition in the river, it was specified for model simulation.A Q 75 (a flow that equals or exceeded 75% of the time) flow condition is widely used in Taiwan to represent the designed low flow for water quality model simulation.The river discharges at the tidal limits of the three major tributaries, i.e., Tahan Stream, Hsintien Stream and Keelung River, were determined using a Q 75 low flow condition.The Q 75 river flows at the upstream reaches of the Tahan Stream, Hsintien Stream and Keelung River are 3.36, 14.23 and 3.33 m 3 /s, respectively.In order to comprehend the impact of the change in nonpoint sources loadings due to a change in land use, the point source loadings were kept the same with model calibration and verification (Section 3.3).Table 8 presents the estimated nonpoint source loadings for different water qualities under the Q 75 flow condition.
Three scenarios, as presented in Table 9, were adopted to respond to the change of nonpoint sources and freshwater discharges and to predict the dissolved oxygen distribution in the tidal estuarine system.To operate the model prediction for the different scenarios, we used 10%-100% nonpoint source increases and 10%-100% Q 75 flow increases to represent the increasing ratio of water quality concentrations at the upstream boundaries and the river flows at the upstream boundaries, respectively, as a result of the change in land use.The incremental ratio for the nonpoint source and Q 75 flow was set to 10% for each model run.

Scenario Simulation Condition
Scenario 1 10%~100% increasing ratio of nonpoint source and Q 75 flow at upstream boundaries (the increment ratio is 10%) Scenario 2 50% increasing ratio of nonpoint source and 10%~100% increasing ratio of the Q 75 flow at upstream boundaries (the increment ratio is 10%) Scenario 3 100% increasing ratio of nonpoint source and 10%~100% increasing ratio of the Q 75 flow at upstream boundaries (the increment ratio is 10%) Figure 8 presents the prediction of the dissolved oxygen distribution along the Danshuei River to the Tahan Stream, the Hsintien Stream and the Keelung River for the baseline case and different cases.In the figure, Line 1, Line 2 and Line 3 represent a 50% increasing ratio for the nonpoint source, a 50% increasing ratio for the nonpoint source, a 50% increasing ratio for the Q 75 flow, a 100% increasing ratio for nonpoint source and a 50% increasing ratio for the Q 75 flow, respectively.It can be determined that the dissolved oxygen in Line 2 is higher compared to that of the baseline in the Danshuei River to Tahan Stream, while the dissolved oxygen distribution in both Lines 2 and 3 is higher compared to that of the baseline.The dissolved oxygen distribution in all of the lines is lower than the baseline.This result indicates that the nonpoint source and the freshwater discharge in the upstream reaches would dominate the dissolved oxygen distribution in the primary stream and the two tributaries.Increasing the nonpoint source worsens the water quality, while increasing the freshwater discharge enhances the dissolved oxygen.(Note that Line 1 represents a 50% increasing ratio for the nonpoint source, Line 2 represents a 50% increasing ratio for the nonpoint source and a 50% increasing ratio for the Q 75 flow and Line 3 represents a 100% increasing ratio for the nonpoint source and a 50% increasing ratio for the Q 75 flow).
The prediction of dissolved oxygen concentration for Scenario 1 at different locations is provided in Figure 9.A 10% incremental ratio in the nonpoint source was conducted for each model simulation, which demonstrated that the dissolved oxygen decreased as the nonpoint source increased.A linear regression was performed at each location.We determined that a high determination coefficient was exhibited at each location.Figure 10 provides the predicted dissolved oxygen concentration for Scenario 2 at three locations.A 10% incremental ratio in the Q 75 flow was applied to each model simulation.The dashed line in Figure 10 represents the baseline condition.Above and below the dashed line represent the increasing and decreasing dissolved oxygen concentration, respectively.The modeling results indicate that when the Q 75 flow increases more than 20% under an increasing nonpoint source of 50% at the Hsin-Hai Bridge and the Chung-Cheng Bridge, the dissolved oxygen concentration will be higher than the baseline condition.The Q 75 flow will increase by at least 50% to overcome a 50% nonpoint source increase and maintain a high dissolved oxygen concentration that is not lower than the baseline condition at the Nan-Hu Bridge.The difference between Scenario 2 and Scenario 3 is the increasing ratio of the nonpoint source.For Scenario 3, the increasing ratio of the nonpoint source is set to 100% (i.e., two-times the nonpoint source).The simulated results indicate that when the Q 75 flow increases up to 50% and 40% at the Hsin-Hai Bridge and the Chung-Cheng Bridge, respectively, the dissolved oxygen concentration can be higher than the baseline condition.At the Nan-Hu Bridge, the Q 75 flow increases up to 100% (i.e., two-times the Q 75 flow) to overcome a 100% increase in the nonpoint source.
Garnier et al. [1] used a biogeochemical modeling approach to simulate historical changes in the nutrient delivery and water quality in the Zenne River in Belgium.It was determined that the change in land use could influence the nonpoint source pollution and flow rate and result in a change in water quality, especially the dissolved oxygen concentration in the river.Based on the simulations of different scenarios in this study, management strategies for water quality can be determined to at least match the dissolved oxygen to the baseline condition.

Conclusions
A coupled three-dimensional hydrodynamic and water quality model (SELFE-WQ) was applied to the Danshuei River estuary and its adjacent coastal region.The model was validated against the water surface elevation, salinity distribution and water quality variables.The simulated results obtained using the three-dimensional hydrodynamic model indicated that the computed water surface elevation and salinity reproduced the observational data well.The modeling results obtained using the water quality model indicated that the computed water quality parameters, including dissolved oxygen, carbonaceous biochemical oxygen demand, ammonium nitrogen and total phosphorus, mimicked the trends of the observational data.The overall performance of the model is in qualitative agreement with the available field data.
Then, the validated model was used to assess the system response to a change of nonpoint sources and freshwater discharges from the drainage basin.Three scenarios were investigated to predict the dissolved oxygen distribution in the tidal estuarine system.We determined that an increasing nonpoint source produced a worse water quality, while an increasing freshwater discharge enhanced the dissolved oxygen.The increasing percentage of the nonpoint source and freshwater discharge could produce different simulated results of dissolved oxygen concentration at different locations.However, the effects of combining the nonpoint source and freshwater discharge on dissolved oxygen would result in a nonlinear relationship.The numerical model can provide a useful tool for developing management practices and protecting the water quality in the estuarine system in the future.
In a future study, a watershed model (i.e., Soil and Water Assessment Tool, SWAT) can be used to predict flow and nonpoint source loading and then coupled with a hydrodynamic and water quality model to more accurately assess the influence of land use change on water quality in a tidal estuarine system.

Figure 2 .
Figure 2. Topography of the Danshuei River estuarine system and its adjacent coastal and unstructured grid for the computational domain.

Figure 4 .
Figure 4. Comparison between measured and simulated salinities along the Danshuei River-Tahan Stream on 26 November 2010, during the (a) flood tide and (b) ebb tide.

Figure 6 .
Figure 6.Comparison between measured and simulated water quality distribution along the Hsintien Stream on 1 March 2010: (a) DO; (b) CBOD; (c) NH 4 ; (d) TP; and (e) NO 3 (simulated result shown only as no NO 3 measurement).

Figure 7 .
Figure 7.Comparison between measured and simulated water quality distribution along the Keelung River on 1 March 2010: (a) DO; (b) CBOD; (c) NH 4 ; (d) TP; and (e) NO 3 (simulated result shown only as no NO 3 measurement).

Figure 8 .
Figure 8. Predicting dissolved oxygen distribution for baseline and different cases along the (a) Danshuei River to Tahan Stream, (b) Hsintien Stream and (c) Keelung River.(Note that Line 1 represents a 50% increasing ratio for the nonpoint source, Line 2 represents a 50% increasing ratio for the nonpoint source and a 50% increasing ratio for the Q 75 flow and Line 3 represents a 100% increasing ratio for the nonpoint source and a 50% increasing ratio for the Q 75 flow).

Figure 9 .
Figure 9. Predicting dissolved oxygen concentration for Scenario 1 at different locations: (a) Hsin-Hai Bridge; (b) Chung-Cheng Bridge; and (c) Nan-Hu Bridge.Note that the incremental ratio of the nonpoint source is 10% for each run.The black point and the solid line represent the model prediction and the regression line, respectively.

Figure 10 .
Figure 10.Predicting dissolved oxygen concentration for Scenario 2 at different locations: (a) Hsin-Hai Bridge; (b) Chung-Cheng Bridge; and (c) Nan-Hu Bridge.Note that the incremental ratio of the Q 75 flow is 10% for each run.The black point represents the model prediction.

Figure 11
Figure 11 illustrates the predicted dissolved oxygen concentration for Scenario 3 at three locations.The dashed line in Figure 11 also represents the baseline condition.The difference between Scenario 2 and Scenario 3 is the increasing ratio of the nonpoint source.For Scenario 3, the increasing ratio of the nonpoint source is set to 100% (i.e., two-times the nonpoint source).The simulated results indicate that when the Q 75 flow increases up to 50% and 40% at the Hsin-Hai Bridge and the Chung-Cheng Bridge, respectively, the dissolved oxygen concentration can be higher than the baseline condition.At the Nan-Hu Bridge, the Q 75 flow increases up to 100% (i.e., two-times the Q 75 flow) to overcome a 100% increase in the nonpoint source.Garnier et al.[1] used a biogeochemical modeling approach to simulate historical changes in the nutrient delivery and water quality in the Zenne River in Belgium.It was determined that the change in land use could influence the nonpoint source pollution and flow rate and result in a change in water quality, especially the dissolved oxygen concentration in the river.Based on the simulations of

Figure 11 .
Figure 11.Predicting dissolved oxygen concentration for Scenario 3 at different locations: (a) Hsin-Hai Bridge; (b) Chung-Cheng Bridge; and (c) Nan-Hu Bridge.Note that the incremental ratio of the Q 75 flow is 10% for each run.The black point represents the model prediction.

Table 1 .
Statistical error between simulated and measured water surface elevations at seven tidal stations during the period of 27 February to 8 March 2010.
Notes: AME, absolute mean error; RMSE, root mean square error.

Table 2 .
Coefficients adopted in the water quality model.

Table 3 .
Statistical error between simulated and measured water quality variables on 1 March 2010.
Notes: AME, absolute mean error; RMSE, root mean square error.

Table 4 .
Statistical error between simulated and measured water quality variables on 1 June 2010.
Notes: AME, absolute mean error; RMSE, root mean square error.

Table 5 .
Statistical error between simulated and measured water quality variables on 3 September 2010.
Notes: AME, absolute mean error; RMSE, root mean square error.

Table 6 .
Statistical error between simulated and measured water quality variables on 2 December 2010.

Table 7 .
The water quality specified at the upstream and ocean boundaries for model simulation.

Table 8 .
Estimated nonpoint source loadings under Q 75 flow condition.

Table 9 .
Description of different scenarios for model prediction.