Acessibilidade / Reportar erro

A coupled model of hydrodynamics circulation and water quality applied to the Rio Verde reservoir, Brazil

Acoplamento entre modelos de circulação de hidrodinâmica e de qualidade da água aplicado ao reservatório do Rio Verde, Brasil

Abstract

This work applied the coupled horizontal two-dimensional hydrodynamic circulation model (2DH) and the vertically integrated water quality model for non-conservative and passive scalars to the Rio Verde reservoir in the state of Paraná, Brazil, to simulate flow, temperature and water quality parameters. The water quality model uses the same spatial grid applied for the hydrodynamics model. Flow velocities and turbulence coefficients previously defined in the hydrodynamics model can be used directly in the water quality model. Modeling results were compared to field data for a period of 308 days, from February 27th 2010 and December 31th 2010. Field data included water and air temperature, wind, relative humidity, radiation, discharges and concentration values of some substances in the tributaries of the reservoir. The results indicated that hydrodynamic circulation and, consequently, horizontal transport, are strongly dependent on the wind. Inflows/outflows generated a localized circulation. The results for the temperature and dissolved oxygen concentration were compared with field measures and a satisfactory consistency was achieved. Despite the errors associated with boundary conditions, the models demonstrated their potential to adequately simulate the data set collected from the reservoir.

Keywords:
hydrodynamic circulation; Rio Verde reservoir; water quality models.

Resumo

No presente trabalho, os modelos de circulação hidrodinâmico bidimensional horizontal (2DH) e de qualidade da água integrados na direção vertical, para escalares não conservativos e passivos, são resolvidos de forma acoplada e aplicados ao reservatório de Rio Verde, localizado no estado do Paraná, Brasil, para simular a circulação e parâmetros qualidade da água no reservatório de Rio Verde. O modelo de qualidade da água utiliza a mesma grade espacial aplicada ao modelo hidrodinâmico. As velocidades e os coeficientes de relacionados à turbulência, previamente definidos no modelo hidrodinâmico, podem ser utilizados diretamente no modelo de qualidade da água. Os resultados de modelagem foram comparados com dados de campo para um período de 308 dias, de 27 de fevereiro de 2010 a 31 de dezembro de 2010. Dados de campo incluem temperatura da água e do ar, vento, umidade relativa, radiação, vazões e valores de concentração de algumas substâncias nos afluentes do reservatório de Rio Verde. Os resultados mostram que a circulação hidrodinâmica e, consequentemente, o transporte horizontal, são fortemente dependentes do vento. As vazões afluentes/efluentes geram uma circulação localizada, restrita as regiões próximas às embocaduras. Os resultados de temperatura e da concentração de oxigênio dissolvido obtidos pelo modelo foram comparados com dados medidos, mostrando uma concordância satisfatória. Apesar dos erros associados às condições de contorno, os modelos demonstraram seu potencial para simular adequadamente o conjunto de dados coletados do reservatório.

Palavras-chave:
circulação hidrodinâmica; modelos de qualidade da água; reservatório do Rio Verde.

1. INTRODUCTION

Reservoirs are complex environments which create interdependent relationships between their tributaries (rivers and streams that contribute for reservoir discharge) and the watershed environmental conditions, which are strongly influenced by water resources uses. By reducing rivers’ velocities with dams, lentic environments can be created which can in turn drastically alter the original conditions of the rivers with greater flow (lotic environments). In this context, the importance of limnological studies aiming at understanding the transport and mixing mechanisms of substances in lentic environments is emphasized.

Quantitative water quality modeling is a major challenge of these studies, due to the need for interdisciplinary approaches that involve advective-diffusive mechanisms and physical, chemical and biological processes in substance transportation. As a consequence, the development of mathematical models to quantify these processes has been the subject of intense research in the last decades. System complexity should be taken into account during the modeling process, ensuring that it has the capacity to model the environmental behavior and evaluate its several effects. Numerical models have been used as a scientific and managerial tool for analyzing temporal-spatial distributions of non-conservative water quality parameters. Among many works, the following deserve to be mentioned: Deus et al. (2013)DEUS, R.; BRITO, D.; MATEUS, M.; KENOV, I.; FORNARO, A.; NEVES, R. et al. Impact evaluation of a pisciculture in the Tucuruí reservoir (Pará, Brazil) using a two-dimensional water quality model. Journal of Hydrology, v. 487, p. 1-12, 2013. https://doi.org/10.1016/j.jhydrol.2013.01.022
https://doi.org/10.1016/j.jhydrol.2013.0...
show the application of the CE-QUAL-W2 model in the Tucurui reservoir in Brazil for different pisciculture scenarios. Using the same model, Kuo et al. (2006)KUO, J.; LUNG, W.; YANG, C.; LIU, W.; YANG, M.; TANG, T. Eutrophication modelling of reservoirs in Taiwan. Environmental Modelling & Software, v. 21, p 829-844, 2006. https://doi.org/10.1016/j.envsoft.2005.03.006
https://doi.org/10.1016/j.envsoft.2005.0...
simulated water quality changes at temperate and sub-tropical reservoirs in Taiwan. Liu et al. (2008)LIU, X.; PENG,W.; HE, G.; LIU, J.; WANG, Y. A Coupled Model of Hydrodynamics and Water Quality for Yuqiao Reservoir in Haihe River Basin. Journal of Hydrodynamics, v. 20, p. 574-582, 2008. https://doi.org/10.1016/S1001-6058(08)60097-9
https://doi.org/10.1016/S1001-6058(08)60...
published a study describing a 2-D coupled model of hydrodynamics and water quality applied to Yuqiao reservoir in China, developed to compute the hydrodynamic field and the variations of total phosphorus and total nitrogen.

The Rio Verde reservoir is located in the Metropolitan Region of Curitiba (RMC), in the state of Paraná, Brazil. Based on collected data included in the work by Carneiro et al. (2014)CARNEIRO, C.; ANDREOLI, C. V.; CUNHA, C. L. N.; GOBBI, E. F. Reservoir Eutrophication: Preventive Management An Introduction to Eutrophication and the Studies. London: IWA Publishing, 2014. 478p. , one can notice that the Rio Verde reservoir presented a moderately degraded environment. In 2005, an algae bloom of Cylindrospermopsis reciborskii was recorded, reaching 96,489 cells.mL-1. The Rio Verde reservoir has a dynamic system of stratification characterized by a reasonable stratification in the summer and a mixed water column in the winter. Due to the shallow depths, it is expected that the currents flowing can be well represented by depth-averaged variables.

Therefore, the purpose of this study was to use the hydrodynamic and water quality models of the SisBAHIA®, in coupled form, to simulate flow, temperature and water quality parameters in the Rio Verde reservoir. The advantage of coupling the two models, as developed in the current work, appears in the determination of the flow velocities and turbulence coefficients, which is done previously in the hydrodynamic model, and can be used directly in the water quality model. For this reason, the size of the system to be solved and, consequently, the computational effort, are reduced when compared to other schemes.

The SisBAHIA® (in Portuguese, Sistema Base de Hidrodinâmica Ambiental) is the Hydrodynamic Environmental System developed by the Coastal and Oceanographic Engineering Department, Oceanic Engineering Program, Federal University of Rio de Janeiro (COPPE/UFRJ). In the present work, only the water quality model of SisBAHIA® is shown. Further information about the hydrodynamic model of SisBaHiA® can be obtained through: www.sisbahia.coppe.ufrj.br.

2. MATERIALS AND METHODS

The water quality model of SisBaHiA® (MQA) is an Eulerian advective-diffusive transport model vertically integrated for non-conservative scales; in other words, the variability of substances concentrations is a result of physical, chemical and biological processes.

In the discretization of time and space, the SisBaHiA® adopts the finite difference method and the finite element method, respectively. The MQA uses the same spatial grid applied for the hydrodynamics model and different time-step lengths can be employed in analyses. Flow velocities and turbulence coefficients, previously defined in the hydrodynamics model, can be used directly in the water quality model. The advective mechanism is related to flow velocities and can be estimated through hydrodynamic models, which represent an important previous procedure for water quality models’ input data. The time step of the two-dimensional hydrodynamic model was 30s with an average Courant number of 4.0. The time step of the MQA was 60s.

The MQA is able to simulate the oxygen, nitrogen and phosphorus cycles, as well as temperature and salinity parameters due to their strong influence in kinetic processes (Sellers, 1965SELLERS, W. D. Physical Climatology. Chicago: The University of Chicago Press, 1965.). Overall, 11 water quality parameters can be evaluated, namely: Salinity, Temperature, Dissolved Oxygen (DO), Biochemical Oxygen Demand (BOD), Organic Nitrogen, Ammonia Nitrogen, Nitrate Nitrogen, Chlorophyll_a, Herbivore Zooplankton, Organic Phosphorus and Inorganic Phosphorus. The MQA uses an Eulerian approach by the finite element mesh, previously defined in the hydrodynamic model. With this condition it is possible to share the same flow velocity components, wind velocity and geometric characteristics. The spatial discretization employs, mainly, squared biquadratic finite elements (Rosman, 2016ROSMAN, P. C. C. Referência Técnica do SISBAHIA - SISTEMA BASE DE HIDRODINÂMICA AMBIENTAL. Programa COPPE: Engenharia Oceânica, Área de Engenharia Costeira e Oceanográfica. Rio de Janeiro: UFRJ, 2016.).

2.1. Mathematical model

In the MQA, the mass-balance equation for a non-conservative substance is applied considering depth-integrated Cartesian System (x, y) aligning in the east, north and vertical directions for three transport mechanisms: the advective term, the diffusive term, and the kinetic processes. The mass-balance equation is expressed by Equation 1, as defined by Sheng and Villaret (1989)SHENG, Y. P.; VILLARET, C. Modeling the effect of suspended sediment stratification on bottom exchange processes. Journal of Geophysical Research, v. 94, p.1429-1444, 1989. https://doi.org/10.1029/JC094iC10p14429
https://doi.org/10.1029/JC094iC10p14429...
:

C m t + U i C m x i = 1 H x j H D i j δ j k + Λ k 2 12 U j x k C m x k + R m (1)

In Equation 1, C m is the concentration of m substance (mg L-1); t is the time (s); U i represents the depth-averaged components of the horizontal velocity (U and V) (m s-1); H is the water depth (m); D ij is the turbulent viscosity coefficient of mass (m2 s-1); δ jk is the Kronecker delta; Λk represents the widths of the spatial and temporal Gaussian filters; and R m represents the kinetic processes of substance production and consumption. In Equation 1, i, j = 1,2 and k = 1, 2, 3, with k = 3 corresponding to time t. The following interpretation is valid for the index coefficient C m: C 1 = ammonia nitrogen (mg N L-1), C 2 = nitrate nitrogen (mg N L-1), C 3 = inorganic phosphorus (mg P L-1), C 4  = Herbivore Zooplankton (mg C L-1), C 5 = biochemical oxygen demand (mg O2 L-1), C 6 = dissolved oxygen (mg O2 L-1), C 7 = organic nitrogen (mg N L-1), C 8 = organic phosphorus (mg P L-1), C 9 = Chlorophyll_a (μg L-1), C T = temperature (0C) and C S = salinity (psu).

The MQA is coupled with the SisBaHiA® hydrodynamic model to provide the advective and diffusive terms in the mass-balance equation. The kinetic processes in MQA were obtained according to the Equations 2, 3, 4, 5, 6, 7, 8, 9 and 10 below (Sheng et al., 1996SHENG, Y. P.; YASSUDA, E. A.; YANG, C. Modeling the effect of reduced nitrogen loading on water quality. In: SPAULDING, M. L.; CHENG, R. T. Estuarine and Coastal modeling. New York: American Society of Civil Engineers, 1996.).

C1: Ammonia nitrogen

Σ R 1 = r n a k r a 1 - f o n C 9 + r n a r c a k r z 1 - f o n C 4 + r n a k e a 1 - f o n C 9 + r n a r c a k e z 1 - f o n C 4 - C 1 k a m + C 1 r n a f u N k g C 9 + k 71 C 7 - k 12 C 1 (2)

C2: Nitrate nitrogen

Σ R 2 = k 12 C 1 - k 2 D C 2 - 1 - C 1 k am + C 1 r na f u N k g C 9 (3)

C3: Inorganic phosphorus

Σ R 3 = r p a k r a ( 1 - f o p ) C 9 + r p a r c a k r z ( 1 - f o p ) C 4 + r p a k e a ( 1 - f o p ) C 9 + r p a r c a k e z ( 1 - f o p ) C 4 - r p a f u P k g C 9 + k 83 C 8 - v f r H C (4)

C4: Herbivore zooplankton

Σ R 4 = r ca E z k gz C 9 C 4 - k rz C 4 - k e z C 4 - k g z c C 4 (5)

C5: BOD - Biochemical oxygen demand

Σ R 5 = - k D C 5 - v s 3 ( 1 - f D 5 ) H C 5 + r o c r c a ( 1 - E Z ) k g z C 9 + r oc r ca k e a C 9 + r o c k e z C 4 (6)

C6: DO - Dissolved oxygen

Σ R 6 = k a ( O s - C 6 ) - k D C 5 - r o n k 12 C 1 + r oc r ca k g C 9 - r o c r c a k r a C 9 - r o c k r z C 4 - S O D H (7)

C7: Organic nitrogen

Σ R 7 = - k 71 C 7 - v s 3 1 - f 7 H C 7 + r n a 1 - E z k g z C 9 C 4 + r n a k r a f o n C 9 + r n a r c a k r z f o n C 4 + r n a k e a f o n C 9 + r n a r c a k e z f o n C 4 (8)

C8: Organic phosphorus

Σ R 8 = r p a K r a f o p C 9 + r p a r c a K r z f o p C 4 + r p a K e a f o p C 9 + r p a r c a K e z f o p C 4 + r p a 1 - E z K g z C 9 C 4 - K 83 C 8 - v s 3 ( 1 - f D 8 ) H C 8 (9)

C9: Chlorophyll_a

Σ R 9 = k g C 9 - k r a C 9 - k e a C 9 - k g z C 9 - v s H C 9 (10)

The coefficients used in the model are given in Table 1. In the MQA two different types of horizontal boundaries are considered: land boundaries and open boundaries. The land boundaries, in general, represent the water body margins and points with inflows or outflows, such as rivers, streams, surface water abstraction points, among others. The open boundaries represent the water domain limits of the study area. The prescription of normal fluxes in this model is associated with land boundaries. Along open boundaries, it is usual to neglect the diffusive fluxes; the model computes the mass balance equation with no diffusive terms.

The numerical implementation of the two models, hydrodynamic and water quality, are not discussed here; the interested reader is referred to Rosman (2016)ROSMAN, P. C. C. Referência Técnica do SISBAHIA - SISTEMA BASE DE HIDRODINÂMICA AMBIENTAL. Programa COPPE: Engenharia Oceânica, Área de Engenharia Costeira e Oceanográfica. Rio de Janeiro: UFRJ, 2016. for additional information concerning the hydrodynamic model. The numerical model developed for the advective and diffusive transport is described in detail by Cunha et al. (2002)CUNHA, C. L. N.; MONTEIRO, T. C.; ROSMAN, P. C. C. Two-Dimensional Modelling of Scalars Transport. Revista Brasileira de Recursos Hídricos, v. 7, p. 63-79, 2002..

Table 1.
Parameters employed in the water quality model and values adopted in the Rio Verde reservoir analysis.

2.2. Description of the Rio Verde Reservoir

The Rio Verde reservoir is located in the Metropolitan Region of Curitiba (RMC) and inserted in the Rio Verde watershed, which has a total area of 242 km2. The reservoir has a plan area of 5.971 Km2, an average volume of 25.6 hm³, an average water depth about 5.6 m and receives the water contribution of 14 sub-watersheds. The Rio Verde reservoir is oriented in a Northeast (NE)-Southwest (SW) direction, with a length of approximately 7,500 m and a maximum width of 1,300 m. The ratio between the depth and width of the reservoir is about 1:230 and between the depth and length is about 1:1330; these features suggest a predominance of horizontal velocities.

Figure 1 shows the bathymetry data and also the water quality and discharge points used for measurement. The main tributary of the reservoir is the Verde River (F4), accounting for about 70% of the total discharge flow; the other 13 smaller tributaries contribute with low flow levels. Figure 2 shows the daily discharge of each tributary from March 2010 to December 2010, estimated by the SWAT model (Soil and Water Assessment Tool). In relation to outflows, two points can be set: one close to the dam, which represents a REPAR (President Getúlio Vargas Refinery) license, with value of 0.83 m3/s; and another at the spillway, which represents difference between the sum of inflows and outflows (Carneiro et al., 2014CARNEIRO, C.; ANDREOLI, C. V.; CUNHA, C. L. N.; GOBBI, E. F. Reservoir Eutrophication: Preventive Management An Introduction to Eutrophication and the Studies. London: IWA Publishing, 2014. 478p. ).

Figure 1.
Map of Rio Verde reservoir bathymetry, water quality measure points and discharge measure points during 2010.

Meteorological data were obtained from a station installed near the dam (Latitude 25º31´36,83" S and Longitude 49º31´39,07" W). This station provided wind direction, wind velocity, air temperature, solar radiation, relative humidity and precipitation from March 2010 to December 2010. The daily maxima values of solar radiation showed a gradual decrease from March to July, reaching 457 W m-2 in May (Figure 2); from August to December the values increased, reaching 851 W m-2 in November. In Figure 2, the air temperatures data show a variability related to seasons, with minimum moving average of 12.6°C in August and maximum moving average of 20.7°C in April. With respect to the relative humidity, Figure 2 shows a small moving average variation over 2010, with values between 74.1% and 88.2%.

Figure 2.
Tributaries daily discharge (m3 s-1), daily average of solar radiation (W m-2), air temperature and relative humidity measured at the meteorological station of the Rio Verde reservoir from March 2010 to December 2010. The continuous line represents the moving average for 30 days.

Carneiro et al. (2014)CARNEIRO, C.; ANDREOLI, C. V.; CUNHA, C. L. N.; GOBBI, E. F. Reservoir Eutrophication: Preventive Management An Introduction to Eutrophication and the Studies. London: IWA Publishing, 2014. 478p. show that only the Verde River (F4) and two tributaries (TE10 and TD4) account for 90% of the pollutant load that arrives at the Rio Verde reservoir. Figure 3 shows the concentration values of eight parameters in three tributaries (F4, TE10 and TD4) from March 2010 to December 2010: water temperature, ammonia, nitrate, organic nitrogen, organic phosphorus, inorganic phosphorus, DO and BOD.

The Rio Verde reservoir presents a dynamic system of stratification characterized by a reasonable stratification in the summer and a mixed water column in the winter. Figure 4 shows the temperature profiles measured at R1 and R4 stations (see Figure 1) from March 2010 to December 2010. At the R4 station between April and September, the water column shows no stratification and can be considered well mixed; from October the water column begins to stratify with the increasing air temperature, reaching largest gradients in December. The station R1 is located in the shallow portion of the reservoir, and it is not possible to observe any stratification phenomena there.

With respect to the wind characteristics, the measured data between March and December 2010 demonstrate a predominance of NE direction, which coincides with the alignment of the reservoir axis. The wind flow represents the major influential variable for hydrodynamic circulation and mass transport in this reservoir. In two areas the hydrodynamic flow and mass transport are dominated by inflows and outflows: one near the dam, where the main water intake and the spillway are located; and the other near the Verde River, the main tributary. Outside these areas, circulation and transport are determined by the wind.

Figure 3.
Concentrations values of some substances in the tributaries of the Rio Verde reservoir.

Figure 4.
Temperature profiles at station R1 (left) and R4 (right) during 2010.

3. RESULTS AND DISCUSSION

Two-dimensional hydrodynamic models of Rio Verde reservoir were developed using SisBaHiA®. The spatial discretization in the horizontal x-y plane was carried out through sub-parametric Lagrangian quadrilateral finite elements with 9 nodes. In a sub-parametric element, the geometry of the element is linear and is defined only by the vertices. However, the variables are quadratic, and in addition to the values at the vertices, the middle values at each side are also needed and, in the case of quadrilaterals, a node at the center of the element is also necessary.

The mesh consists of 507 elements and 2402 nodes. The simulation was performed for a period of 308 days, from February 27th 2010 to December 31th, 2010, with a 60s time step. These dates were chosen due to available measurements of water quality parameters, made at four locations in the reservoir, which enabled the calibration process in MQA evaluation.

In the numerical model, wind conditions are considered unsteady and spatially homogeneous; the input data used in the model were the time series of wind speeds and directions measured at the meteorological station in the reservoir. At the open boundaries domain, the elevations were imposed; at the land boundary, except for the nodes corresponding to the rivers, all nodes were considered impermeable with null tangential and normal velocity components imposed.

The bottom friction coefficient can be written in terms of the Chezy coefficient, which depends on the amplitude of the equivalent bottom roughness. The amplitude of the equivalent bottom roughness, ε, was defined according to the bottom sediment characterization and distribution (Rosman, 2016ROSMAN, P. C. C. Referência Técnica do SISBAHIA - SISTEMA BASE DE HIDRODINÂMICA AMBIENTAL. Programa COPPE: Engenharia Oceânica, Área de Engenharia Costeira e Oceanográfica. Rio de Janeiro: UFRJ, 2016.). Silty-clay sediments are predominant in the reservoir (ε ≈ 0.015m) and in regions with large amount of submerged trees, bottom roughness amplitude values were increased (with values around 0.130 m).

Figure 5 shows the hydrodynamic model results for two situations: in June 26th, 2010, with an average wind speed of 1,81 m s-1, an average wind direction of 73°, and a main river discharge of 2,24 m3/s; and in October 18th, 2010, with an average wind speed of 4,84 m s-1, an average wind direction of 22°, and a main river discharge of 1,63 m3 s-1. The velocity fields show a strong relationship with the local wind pattern through the presence of a clockwise vortex formed in the region near the dam, which intensifies depending on the wind pattern. On June 26th, 2010, there was a northeasterly wind at an average of 1.81 m s-1. The currents follow the wind direction, forming a large clockwise vortex in the region near the dam. In the regions near the spillway and the Verde river, there is an intensification of the currents due to local effects of inflows/outflows. On October 18th, 2010, the wind speed is stronger, at an average of 4.84 m s-1. Under these conditions, two vortices were formed: a large vortex, clockwise in the central region of the reservoir, and a smaller one, also clockwise. There was an intensification of the currents in the region near the Verde river, caused by the flooding observed during this period. Taking into account these patterns, it can be stated that the hydrodynamic circulation, and consequently the horizontal transport, is strongly dependent on the wind and that inflows/outflows generate a localized circulation.

Figure 5.
Velocity fields obtained with SisBaHiA® from June 26th, 2010. Velocity fields obtained with SisBaHiA® from October 18 th, 2010.

From March 2010 to December 2010, data related to Temperature and DO were collected. The measurements were made at R4 station in reservoir, see Figure 1; the measurements were obtained at a several positions along the water depth. The analyses were performed for the same dates as that of the hydrodynamic simulation, between February 27th 2010 and December 31th, 2010. This article shows just the results for the Temperature and DO, which have measured data for comparison.

The turbulent diffusion coefficient used from the water quality model are: D xx (m2 s-1) = 2.0, D yy (m2 s-1) = 0.5, with a 60s time step. Table 1 lists the model coefficients and constants employed in the numerical simulation. The values of the parameters and constants are well-defined in the literature and were used in the application considering that there are no available studies for water quality modeling in Rio Verde reservoir. This is the main difficulty in performing the analysis proposed in this article. The field data are not sufficient for a complete calibration, and some adjustments were carried out manually to obtain numerical results close to the measured data, within a specified variation limit; some adjustments proved to be necessary, mainly in the reaeration and deoxygenation coefficients. In order to quantify the accuracy of fit, correlation coefficient R2 of the regression line between the model results and the field data at R4 station was computed. The best calibration required a R2 value as close to 1.0 as possible.

The initial concentrations are: Salinity = 0.0 psu, Temperature = 21.08ºC; Organic Nitrogen = 0.30 mg N L-1; Ammoniac Nitrogen = 0.157 mg N L-1; Nitrate Nitrogen = 0.038 mg N L-1; Biochemical Oxygen Demand = 1.50 mg O2 L-1; Dissolved Oxygen = 5.80 mg O2 L-1; Chlorophyll-a = 12.80 μg L-1; Zooplankton = 0.0 mg L-1; Organic Phosphorus = 0.0189 mg P L-1; Inorganic Phosphorus = 0.0039 mg P L-1. Such parameters represent the measurement made at the R4 monitoring point in the reservoir, see Figure 1, on March, 16th, 2010.

The temperature is used as the calibration parameter for variables related to advective and diffusive transport. The long-term temperature data simulated with SisBaHiA® were compared with average values of vertical profiles, measured at R4 station from March 2010 to December 2010. The model reproduces the average temperature variation and, according to Figure 6, it is possible to state that variable boundary conditions are adequately reproduced. The R4 is not under the direct influence of the rivers; consequently, there is a better agreement between the results. Figure 6 shows the DO concentrations measured at the R4 station and computed numerically by SisBAHIA®. The differences between these results showed that SisBAHIA® can predict adequately the variations of DO in the process of oxidation of organic matter present in the reservoir. The average relative errors between the model results and the field data at R4 station are: 0.095 for temperature and 0.406 for DO. Figure 6 also depicts the regression line; the correlation coefficient values (R2) were 0.8973 for the temperature, and 0.4575 for the DO. This indicated good agreement between simulated and field data for the temperature. However, DO concentrations do not indicate a good fit. For a better calibration of the model, a larger temporal series would be necessary.

Figure 6.
Temperature and DO concentration measured at the R4 station and computed numerically by SisBAHIA®.

4. CONCLUSIONS

The results presented in this paper show that the environmental modeling of the Rio Verde reservoir is appropriate and can be used quantitatively to study the phenomenon of eutrophication. Whereas the advective and diffusive transport is independent from the parameter modeled, the correct definition of the chemical, physical and biological processes involved in each water quality parameter, made it possible to characterize the distribution of these parameters in the Rio Verde reservoir. Based on the results of the hydrodynamic model, it was possible to note the significant influence of the wind in the hydrodynamic circulation. Thus, the importance of the proper use of wind data when modeling reservoirs becomes clear. Inflows/outflows generate impacts in circulation only in their immediate surroundings and do not contribute to the circulation of the reservoir as a whole. The results for the temperature and dissolved oxygen concentration were compared with field measures and a satisfactory consistency was achieved.

5. REFERENCES

  • CARNEIRO, C.; ANDREOLI, C. V.; CUNHA, C. L. N.; GOBBI, E. F. Reservoir Eutrophication: Preventive Management An Introduction to Eutrophication and the Studies. London: IWA Publishing, 2014. 478p.
  • CUNHA, C. L. N.; MONTEIRO, T. C.; ROSMAN, P. C. C. Two-Dimensional Modelling of Scalars Transport. Revista Brasileira de Recursos Hídricos, v. 7, p. 63-79, 2002.
  • DEUS, R.; BRITO, D.; MATEUS, M.; KENOV, I.; FORNARO, A.; NEVES, R. et al. Impact evaluation of a pisciculture in the Tucuruí reservoir (Pará, Brazil) using a two-dimensional water quality model. Journal of Hydrology, v. 487, p. 1-12, 2013. https://doi.org/10.1016/j.jhydrol.2013.01.022
    » https://doi.org/10.1016/j.jhydrol.2013.01.022
  • KUO, J.; LUNG, W.; YANG, C.; LIU, W.; YANG, M.; TANG, T. Eutrophication modelling of reservoirs in Taiwan. Environmental Modelling & Software, v. 21, p 829-844, 2006. https://doi.org/10.1016/j.envsoft.2005.03.006
    » https://doi.org/10.1016/j.envsoft.2005.03.006
  • LIU, X.; PENG,W.; HE, G.; LIU, J.; WANG, Y. A Coupled Model of Hydrodynamics and Water Quality for Yuqiao Reservoir in Haihe River Basin. Journal of Hydrodynamics, v. 20, p. 574-582, 2008. https://doi.org/10.1016/S1001-6058(08)60097-9
    » https://doi.org/10.1016/S1001-6058(08)60097-9
  • ROSMAN, P. C. C. Referência Técnica do SISBAHIA - SISTEMA BASE DE HIDRODINÂMICA AMBIENTAL. Programa COPPE: Engenharia Oceânica, Área de Engenharia Costeira e Oceanográfica. Rio de Janeiro: UFRJ, 2016.
  • SELLERS, W. D. Physical Climatology. Chicago: The University of Chicago Press, 1965.
  • SHENG, Y. P.; YASSUDA, E. A.; YANG, C. Modeling the effect of reduced nitrogen loading on water quality. In: SPAULDING, M. L.; CHENG, R. T. Estuarine and Coastal modeling. New York: American Society of Civil Engineers, 1996.
  • SHENG, Y. P.; VILLARET, C. Modeling the effect of suspended sediment stratification on bottom exchange processes. Journal of Geophysical Research, v. 94, p.1429-1444, 1989. https://doi.org/10.1029/JC094iC10p14429
    » https://doi.org/10.1029/JC094iC10p14429

Publication Dates

  • Publication in this collection
    14 Nov 2018
  • Date of issue
    2018

History

  • Received
    20 Feb 2018
  • Accepted
    05 Oct 2018
Instituto de Pesquisas Ambientais em Bacias Hidrográficas Instituto de Pesquisas Ambientais em Bacias Hidrográficas (IPABHi), Estrada Mun. Dr. José Luis Cembranelli, 5000, Taubaté, SP, Brasil, CEP 12081-010 - Taubaté - SP - Brazil
E-mail: ambi.agua@gmail.com