SIMULATION OF PHYSICAL AND CHEMIC AL PROCESSES IN RESERVOIRS: TWO CASE STUDIES
Selma L. Garcia Iturbe, B.A., M.S.
Dissertation Prepared for the Degree of
DOCTOR OF PHILOSOPHY
UNIVERSITY OF NORTH TEXAS
December 2005
APPROVED:
Miguel Acevedo, Major Professor Samuel F. Atkinson, Committee Member Kenneth Dickson, Committee Member William T. Waller, Committee Member Thomas La Point, Committee Member and Graduate Program Coordinator for Environmental Sciences Art Goven, Chair of the Department of Biological Sciences Sandra L. Terrell, Dean of the Robert B. Toulouse School of Graduate Studies García Iturbe, Selma L., Simulation of physical and chemical processes in reservoirs:
Two case studies. Doctor of Philosophy (Environmental Science), December 2005, 206 pp., 64
tables, 65 illustrations, references, 61 titles.
Managing water quality aspects requires the use of integrative tools that allow a holistic
approach to this problem. Water quality models coupled to hydrodynamic models are these tools.
This study presents the application of the water quality model WASP coupled to the
hydrodynamic model DYNHYD for two distinct reservoirs: Lake Texoma and Tocoma
Reservoir. Modeling the former included simulations of water velocities, water level, and four
chemical and physical compounds: chlorides, dissolved oxygen (DO), biochemical oxygen
demand (BOD), and total suspended solids (TSS); and validation of the results by comparing
with observed values during March – May, 1997. The latter is still under project status and the
simulation was performed in a prospective way. The analysis included simulations of water
velocities under current and for expected conditions, DO and BOD. Both models, DYNHYD and
WASP, fitted pretty well to observed conditions for Lake Texoma and for where Tocoma
Reservoir has been planned.
Considering management and decision support purposes, the role of boundary and
loading conditions also was tested. For Lake Texoma, controlling boundary conditions for
chlorides is a determinant factor for water quality of the system. However, DO and TSS in the reservoir are governed by additional process besides the condition of the boundary. Estimated
loadings for this system did not provided significant effects, even though the allocation of a load
for chlorides resulted in significant changes in the trend for expected chloride concentrations at
the Washita River Arm of Lake Texoma. For Tocoma Reservoir, the expected concentration of
DO all over the reservoir is going to driven by boundary conditions, as well as by the management of autochthonous BOD loadings provided by vegetation decomposition. These two factors will be determinant for the resulting water quality of the future reservoir.
Copyright 2005
by
Selma L. Garcia Iturbe
ii
ACKNOWLEDGEMENTS
I would like to thank her dissertation adviser Dr. Miguel Acevedo for his guidance and support on this study. Many thanks are due Dr. William Thomas Waller for providing invaluable advice during all the stages of this work. I also wish to acknowledge Dr. Samuel Atkinson, Dr.
Kenneth Dickson and Dr. Thomas La Point, members of the committee, for their advice. I would also like to thank my fellow graduate student Alexandra Christina Upton for her friendship and support during most of the stages of this study.
Finally, I am grateful for the financial support of this research provided by CVG
EDELCA in Venezuela, the U.S. Environmental Protection Agency and the U.S. Army Corp of
Engineers.
iii
TABLE OF CONTENTS
Page
ACKNOWLEDGMENTS ...... iii
LIST OF TABLES...... vi
LIST OF ILLUSTRATIONS...... x
Chapters
1. INTRODUCTION ...... 1
2. HYDRODYNAMIC SIMULATION...... 7 Model Description ...... 9 Lake Texoma: Description of the Study Area ...... 10 Model Parameterization for Lake Texoma ...... 15 Hydrodynamic Modeling Results for Lake Texoma ...... 26 Tocoma Reservoir: Description of the Study Area...... 33 Model Parameterization for Tocoma Reservoir...... 37 Hydrodynamic Model Results for Tocoma Reservoir...... 44
3. RESERVOIR WATER QUALITY MODELING...... 50 Model Description ...... 51 Previous Water Quality Studies of Lake Texoma...... 54 Water Quality Model Parameterization for Lake Texoma ...... 56 Water Quality Modeling Results for Lake Texoma...... 78 Water Quality Studies for Tocoma Reservoir...... 95 Water Quality Model Parameterization for Tocoma Reservoir...... 98 Water Quality Modeling Results for Tocoma Reservoir ...... 108
4. ROLE OF BOUNDARY CONDITIONS AND LOADINGS IN WATER QUALITY SIMULATIONS...... 115 Sensitivity Analysis to Upstream Boundary Conditions ...... 115 Sensitivity Analysis for Loading Conditions...... 122
5. CONCLUSIONS...... 130
iv
Appendices
A. GEOMETRY DATA FOR LAKE TEXOMA MODEL ...... 136
B. ESTIMATIONS OF TRAVEL TIME FOR LAKE TEXOMA MODEL ...... 139
C. GEOMETRY DATA FOR TOCOMA RESERVOIR MODEL...... 141
D. ESTIMATIONS OF HYDRAULIC TRAVEL TIME FOR TOCOMA RESERVOIR ...... 145
E. INPUT DATA FOR WATER QUALITY SIMULATION OF LAKE TEXOMA ...... 148
F. STATISTICAL ANALYSIS FOR DEVELOPING INITIAL CONDITIONS FOR LAKE TEXOMA WATER QUALITY MODEL...... 168
G. STATISTICAL ANALYSIS FOR DEVELOPING DAILY BOUNDARY CONDITIONS FOR LAKE TEXOMA WATER QUALITY MODEL ...... 174
H. VALUES FOR TESTING SENSITIVITY TO BOUNDARY CONDITIONS ...... 182
I. VALUES FOR TESTING SENSITIVITY TO LOADING CONDITIONS.... 193
REFERENCES ...... 201
v
LIST OF TABLES
Page
1. Some morphometric characteristics of Lake Texoma, Tocoma Reservoir, and their watersheds...... 6
2. Distribution of annual volumes in Lake Texoma for different usages (Compiled from TWDB 2003) ...... 10
3. Some morphometric characteristics of Lake Texoma and its basin ...... 11
4. Main creeks that drain directly into Lake Texoma ...... 15
5. Estimation of maximum time step (in seconds) for different conditions of water velocity and depth of the channels in the developed hydrodynamic network for Lake Texoma ...... 22
6. Time step values for different junction depths of the hydrodynamic network for Tocoma Reservoir...... 42
7. Predicted and measured velocities for Tocoma under river condition ...... 45
8. Relative change of water velocities for Caroni River when comparing river and reservoir condition at Tocoma site...... 46
9. General characteristics of segments used in WASP for Lake Texoma ...... 58
10. Relationship between water quality segments identification and station identification used by Atkinson et al. (1999) ...... 59
11. Some constants and rates used by WASP for water quality simulation ...... 60
12. Different loading conditions evaluated for Lake Texoma ...... 65
13. Identification of routine stations monitored by Atkinson et al. (1999) in Lake Texoma and their estimated location...... 66
14. Estimated initial concentrations of chlorides for Washita and Red River arms for WASP calibration ...... 69
15. Equations used to calculate daily concentrations for TSS at the boundary of the Washita River arm for Lake Texoma...... 73
16. Percentage distribution of solids per segment used in WASP...... 75
17. Equations for estimations of loadings of solids using estimations made by Upton (2004) ...... 77
vi
18. EMCs selected for estimating loadings for Lake Texoma. Nitrogen values are from Upton (2004)...... 78
19. Estimated EMCs for each sub-watershed around Lake Texoma...... 79
20. Statistical analysis of predicted values for TSS concentrations in Lake Texoma using WASP ...... 87
21. Geographic coordinates for water quality sampling stations nearby Tocoma Dam site ...... 97
22. Annual average value for the main limnological characteristics of some water quality sampling stations at the lower Caroni River...... 99
23. Geometric data for water quality segments developed for the WASP network of Tocoma Reservoir...... 101
24. Time functions used for DO and BOD simulation of Tocoma Reservoir ...... 105
25. Constants and rates used by WASP for water quality simulation of Tocoma Reservoir ...... 105
26. BOD loading conditions for adjacent watershed of the Tocoma Reservoir ...... 106
27. Phytomass and degradation factors used for BOD loading estimates during the creation of Tocoma Reservoir...... 108
28. DO values for Guri Dam discharge channels and Tocoma Dam site measured during January and August, 2005...... 110
29. Estimated values used for DO balance at water quality segments at Necuima Canyon under river condition...... 112
30. Results of the sensitivity test for chlorides and DO showing percent change in daily, volume-averaged concentrations for some WASP segments of the Lake Texoma model compare to the base case in March-May, 1997...... 120
31. Results of the sensitivity test for solids 1, 2, and 3 showing percent change in daily, volume-averaged concentrations of TSS for some WASP segments of the Lake Texoma model compare to the base case in March-May, 1997...... 121
32. Output percent change for TSS concentrations at several segments of the water quality network of Lake Texoma during the sensitivity test to Solid 2 loading conditions ..... 124
33. Output percent change for chloride concentrations at several segments of the water quality network for Lake Texoma during the sensitivity test of the model to chloride loading conditions...... 125
34. Junction geometry data for Lake Texoma model ...... 137
vii
35. Channel geometry data for Lake Texoma model...... 138
36. Estimations of travel time per junction for Lake Texoma model using DYNHYD ..... 140
37. Junction geometry data for Tocoma Reservoir model...... 142
38. Channel geometry data for Tocoma Reservoir model ...... 143
39. Estimated travel time per segment in Tocoma Reservoir model ...... 146
40. Water quality segment definition for Lake Texoma model...... 149
41. System Data for Lake Texoma model ...... 150
42. Initial concentrations for chlorides and fractions of suspended solids for Lake Texoma model...... 151
43. Exchange functions defined for Lake Texoma model ...... 152
44. Settling and resuspension functions defined for solids type 1 (clay) in Lake Texoma model...... 152
45. Settling and resuspension functions defined for solids type 2 (silt) in Lake Texoma model ...... 155
46. Settling and resuspension functions defined for solids type 3 (sand) in Lake Texoma model...... 157
47. Red River boundary condition for calibration of WASP for Lake Texoma ...... 159
48. Washita River boundary condition for calibration of WASP for Lake Texoma ...... 160
49. Big Mineral Creek boundary condition for calibration of WASP for Lake Texoma ... 161
50. Water temperature function defined for calibration of WASP for Lake Texoma ...... 162
51. Print interval used in WASP for Lake Texoma water quality simulation ...... 162
52. Time step used in WASP for Lake Texoma water quality simulation ...... 162
53. Chloride loadings (Kg/day) for calibration of WASP for Lake Texoma during March to May, 1997 ...... 162
54. Solid 2 loadings for water quality simulation of Lake Texoma during March to May, 1997...... 163
55. BOD loadings for water quality simulation of Lake Texoma during March to May, 1997 ...... 166
viii
56. Values used to test change in 10% of boundary conditions of chlorides for Lake Texoma ...... 183
57. Values used to test change in 10% of boundary conditions of Solids 1 (Clay) for Lake Texoma ...... 185
58. Values used to test change in 10% of boundary conditions of Solids 2 (Silt) for Lake Texoma ...... 187
59. Values used to test change in 10% of boundary conditions of Solids 3 (Sand) for Lake Texoma ...... 189
60. Values used to test change in 30% of boundary conditions of DO for Lake Texoma ...... 192
61. Values of chloride loads at WASP segment WRT22 to test sensitivity of Lake Texoma model to changes in loading conditions...... 194
62. Values of TSS loads at several WASP segments to test sensitivity of Lake Texoma model to changes in plus 50% of loading conditions ...... 195
63. Values of TSS loads at several WASP segments to test sensitivity of Lake Texoma model to changes in minus 50 % of loading conditions ...... 197
64. BOD loadings for testing sensitivity to 50 % change of loading conditions in Tocoma Reservoir model...... 199
ix
LIST OF ILLUSTRATIONS
Page
1. Lake Texoma indicating different zones in the lake based on water salinity and turbidity (Adapted from Atkinson et al. 1999) ...... 4
2. Tocoma Reservoir site in the Caroni River...... 5
3. Lake Texoma watershed (Adapted from Upton 2004) ...... 12
4. Lake Texoma and surrounding counties (Adapted from Upton 2004)...... 13
5. Image of the DEM file for Lake Texoma showing different depths. Lighter colors correspond to shallower areas (Adapted from: USGS 1999)...... 18
6. Junction network of the hydrodynamic model of Lake Texoma ...... 19
7. Channel segmentation of Lake Texoma for the hydrodynamic model...... 20
8. Variability of error in pool elevation estimations of Lake Texoma using different percentages of accumulated volume in a 24-hour period as inflow in Big Mineral Creek and applying DYNHYD ...... 21
9. Precipitation/Evaporation balance between March,1997-May,1997 using daily reported values for Denison Dam rain gage station ...... 23
10. Boxplot of mean daily total inflow in Lake Texoma during 1997 represented by month. The highest dispersion of values is observed in May and the lowest occurs in August ...... 24
11. Graphical comparison between observed and predicted elevations in Lake Texoma during March, 1997 and May,1997, using DYNHYD and estimated runoff for Big Mineral Creek provided by Upton (2004) ...... 26
12. Graphical comparison between observed and predicted elevations in Lake Texoma during March, 1997 and May,1997, using DYNHYD and a 50% of thedaily accumulated volume in the reservoir not explained by the Red and Washita rivers inflows based on mass balance estimations as a daily inflow of Big Mineral Creek ...... 27
13. Variation of water flow in different zones of Lake Texoma. 0 to 1 corresponds to the Red River Zone, 19 to 28 corresponds to the Red River Transition Zone, 0 to 20 corresponds to the entrance of Washita River to Lake Texoma, 27 to 28 represents Washita River Transition Zone, 28 to 29 and 31 to 0 represents the Main Lake Zone, 0 to 32 corresponds to the entrance of Big Mineral Creek to Lake Texoma, and 10 to 36 is the junction between the Big Mineral Creek and the Red River Transition Zone...... 29
x
14. Spatial distribution of predicted water velocities for Lake Texoma under the steady-state condition ...... 30
15. Graphical representation of travel time in Lake Texoma when a pulse of tracer is applied in the upper portion of the Red River arm ...... 31
16. Graphical representation of travel time in Lake Texoma when a pulse of a tracer is applied in the upper portion of the Washita River arm...... 32
17. Estimation of travel time for Lake Texoma using a combined loading at the Red and Washita River...... 33
18. Scheme of hydroelectrical development in the Lower Caroni River (Adapted from CVG EDELCA 2005b) ...... 34
19. Geographic location for the future Tocoma Reservoir, in the Caroni River ...... 35
20. Mean monthly turbine-used and spilled flow at Guri Dam discharge channels under average conditions ...... 36
21. Average monthly precipitation at Las Babas weather station, close to Tocoma Reservoir site (Period: 1986-2004) ...... 37
22. Schematic location for bathymetric cross-sections in the Caroni River for the Tocoma Reservoir area (Adapted from CVG EDELCA 1987) ...... 39
23. Bathymetric map for Tocoma Reservoir...... 40
24. Hydrodynamic network for Tocoma Reservoir ...... 41
25. Comparison between hydrodynamic network for river and reservoir condition for Caroni River at Tocoma Reservoir site...... 43
26. Spatial variation of simulated velocities for Tocoma Reservoir...... 47
27. Peak concentration of the tracer at segments close to the future dam, in Tocoma reservoir, when a pulse is applied at Guri Dam discharge channels...... 48
28. Temporal variation of a tracer concentration when applied at a) Cunaguaro River, and peak concentration of the tracer at b)the segment next to the future dam, in Tocoma reservoir ...... 49
29. Graphical representation of the configuration of water quality segments for the application of WASP for Lake Texoma...... 53
30. Graphical representation of the configuration of water quality segments for the application of WASP for Lake Texoma...... 54
31. Segmentation of Lake Texoma used in WASP6...... 57
xi
32. Estimated water velocities in every segment of the water quality network developed for Lake Texoma using DYNHYD ...... 62
33. Settling velocities used for the water quality network developed for Lake Texoma using WASP ...... 63
34. Spatial variation of the average concentration of chlorides in the Red River arm, based on measurements made by Atkinson et al. (1999)...... 67
35. Spatial variation of the average concentration of chlorides in the Washita River arm, based on measurements made by Atkinson et al. (1999)...... 68
36. Local watershed defined by Upton (2004) for Lake Texoma...... 76
37. Temporal variation of chlorides concentration in some stations located in the Red River arm of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997...... 81
38. Temporal variation of chlorides concentration in some stations located in the Washita River arm and Main Lake Zone of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997...... 82
39. Temporal variation of total suspended solids concentration in some stations located in the Red River arm of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997...... 85
40. Temporal variation of total suspended solids concentration in some stations located in the Washita River arm of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997...... 86
41. Temporal variation of total suspended solids concentration in some stations located in the Main Lake Zone of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997...... 88
42. Temporal variability of particle resuspension velocity applied to solid type 2 in segment 27 of the water quality network developed for Lake Texoma using WASP, obtained by model calibration ...... 89
43. RMS of errors of TSS concentration for different water quality segments in a) the Red, and b) the Washita rivers arms of Lake Texoma using WASP ...... 91
44. Temporal variation of DO concentration in some stations located in the Red River arm of Lake Texoma in comparison with simulated values for the period between March, 1997 and May, 1997 ...... 92
45. Temporal variation of DO concentration in some stations located in the Washita River arm and Main Lake Zone of Lake Texoma in comparison with simulated values for the period between March, 1997 and May, 1997 ...... 93
xii
46. RMS of errors of DO concentrations for different water quality segments in the Red and Washita River arms of Lake Texoma using WASP...... 95
47. RMS of errors of DO concentrations for some water quality segments at the: a) Red River arm, and b) Washita River arm of Lake Texoma under different conditions of salinity and sediment oxygen demand...... 96
48. Spatial location of water quality measurement sites at the lower Caroni River considered for this study ...... 98
49. Correspondence between water quality and hydrodynamic networks developed for Tocoma Reservoir...... 102
50. Temporal variation of BOD loadings for two segments of the water quality network, while simulating the filling up process for Tocoma Reservoir...... 107
51. Spatial distribution of water quality segments with autochthonous BOD loadings ..... 110
52. Spatial distribution of DO concentrations for Caroni River between Guri Dam discharge channels and Tocoma Dam site ...... 111
53. Expected spatial distribution of DO for Tocoma Reservoir when simulating the reservoir a) with, and b) without BOD loadings...... 114
54. Spatial distribution of percent change of chloride concentrations in Lake Texoma when an increment of +10% of boundary condition for chlorides was applied...... 119
55. Spatial distribution of percent change obtained during the sensitivity test to changes in the boundary condition of DO for Tocoma Reservoir...... 122
56. Chloride concentrations simulated for WASP segment WRT22 while testing sensitivity of Lake Texoma model to changes in chloride loadings ...... 126
57. Comparative spatial variation of potential DO change during the flooding phase of Tocoma Reservoir by reducing autochthonous BOD loading in 50% in comparison with the effect of the total BOD load...... 128
58. Comparative expected spatial distribution of DO for Tocoma reservoir after flooding with 50% and 100% of autochthonous BOD loadings ...... 129
59. Linear regression curve between chloride concentrations and distance at the Washita River arm in Lake Texoma ...... 170
60. Linear regression curve between chloride concentrations and distance at the Red River arm in Lake Texoma ...... 171
61. Linear regression curve between chlorides and conductivity values measured at the gaging station at the Red River near Gainesville, TX ...... 173
xiii
62. Analysis of data distribution for conductivity values measured in the Washita River at the gage station near Dickson, OK ...... 176
63. Analysis of data distribution for conductivity values measured at Station 25 at Lake Texoma ...... 177
64. Linear regression curve between conductivity values measured at Washita River gaging station near Dickson, OK and chloride concentrations measured at Station 25 in Lake Texoma ...... 179
65. Linear regression curve between conductivity values and chloride concentrations measured at Station 25 in Lake Texoma...... 181
xiv CHAPTER 1
INTRODUCTION
Traditionally, monitoring water quality in rivers and lakes has been used as the primary source of information for making decisions and management of water bodies.
Empirical analysis of water quality data has resulted in useful descriptions of aquatic ecosystems (Hakanson and Peters 1995). This approach, however, has been considered insufficient for prediction and quantification of secondary effects, especially those associated to human activities and their effects on water bodies, and the sustainability of this important resource (Goodwing and Hardy 1999).
A new paradigm in water quality management is evolving, and a broad range of variables and relationships among them are considered before defining a management objective for water bodies. This approach requires extensive data and simulation tools to assist in the decision-making processes that managers and communities must make in order to select strategies that contribute to conservation of the considered water system
(Goodwing and Hardy 1999). Consequently, developing predictive models for rivers and lakes is a powerful tool in applied ecology and environmental management.
Water quality simulation is a representation of the natural ecosystem integrating the temporal and spatial variability of the physical, chemical, and biological characteristics of aquatic systems (Ford 1999). Aquatic ecosystem models, and particularly those applied to reservoirs, provide an integrated analysis of the morphometry, hydrology, ecology, and internal and external forcing functions acting on a water body, based on monitoring physical, chemical, and biological features of the water body (Tufford and McKellar 1999). Water quality measurements, normally made at a
1 finite number of specific localities during specified periods of time, provide the information needed to infer process dynamics and their effects over much broader spatial and temporal scales using water quality models. Thus, the decision-making process and water quality management are supported on more integrated and complete information, taking into account not only observed but also expected water quality conditions.
Using simulation tools to study reservoirs is particularly important because these impoundments may be more complex than natural lakes (Hakanson and Peters 1995,
Tufford and McKellar 1999). In a reservoir, river inflows determine a circulation pattern based on lake morphometry, defining advective transport and density gradients (Ford
1999). In addition, hydrology, land use, and weather establish temporal patterns and variations on water quality. As a consequence, the biotic component of the ecosystems respond, developing many different functional habitats (Kalff 2002). These processes can be better evaluated in their whole dimension (temporal and spatial) with a tool that allows their integration. Water quality simulation models represent this tool. When studying a reservoir with multiple uses, this approach can be particularly useful and desirable.
Modeling can also be used for prediction of impacts resulting from the creation of reservoirs. As a predictive tool, models provide a prospective vision of the processes that can result in the environment as a consequence of human developments (Canter 2002). In this study, the hydrodynamic model DYNHYD and the water quality simulation program
WASP have been used for providing a better understanding of physical, chemical, and biological processes developed in a created reservoir, Lake Texoma in the United States, and those expected for a new projected reservoir, the Tocoma Reservoir, in Venezuela.
2 Lake Texoma is a reservoir located on the border between States of Oklahoma and Texas (USA) at the confluence of the Red and Washita rivers. The U.S. Army Corps of Engineers created it in 1945 initially to generate hydroelectrical energy (Breeding
1997). However, at the present time it is also used for flood control, recreation, water supply, and fisheries, which is a very important source of revenue for the local economy
(USACE 2003a).
The Red and Washita rivers provide the main inflows to the lake. Variations of the mean monthly flow in the Red River, at the Gainesville gage station, in Texas, are between 36 and 226 m3/s (period: 1936-2004). In the Washita River, mean monthly flows
at Dickson, in Oklahoma, fluctuate between 18 and 117 m3/s (period: 1928-2004) (USGS
2005). These fluctuations are associated with variations in water quality of the reservoir,
especially during spring time (Work and Gophen 1999).
Analysis of data resulting from intensive monitoring of physical and chemical
characteristics of the lake reflects a spatial variation of conductivity and turbidity
(Atkinson et al. 1999, Waller et al. 2002). Four zones have been defined in those studies:
the Red River Arm and the Washita River Arm, considered as lotic zones, and the Red
River Transition Zone, the Washita River Transition Zone, and the Main Body,
considered as lentic zones (Figure 1). Other studies also consider the Big Mineral Creek
as an additional zone because of its differentiated characteristics (Mabe 2002).
Additional studies on water quality of the reservoir and its relationship with land
use of the shoreline and the watershed also show that there are biological gradients,
affecting biological productivity of the reservoir (Work and Gophen 1999, An et al. 2002,
3 An and Kampbell 2003), and the effects of additional human influence in the water body are still unknown.
Figure 1. Lake Texoma indicating different zones in the lake based on water salinity and turbidity (Adapted from Atkinson et al. 1999) Tocoma Reservoir will be created in Bolivar State (Venezuela) on the lower basin of the Caroni River. CVG EDELCA, a Venezuelan government-owned energy company that produces nearly 70% of the total energy consumed in the country, will create and operate this reservoir mainly for hydroelectrical generation. Located between two impoundments, Guri Reservoir and Caruachi Reservoir, the main source of water will come from the turbinated and spilled water of Guri Reservoir, on the Caroni River. Five tributaries to the Caroni River above the future dam site will contribute an additional, but small amount of flow. These tributaries are small, intermittent creeks, being the most important Cunaguaro River, flowing into the Caroni River close to the future dam site
(Figure 2).
4
C C
a a
r r
u u
a a
c c
h h
i i
R
es
er
v
o
i
r
Tocoma Dam Site
C a r o C unagua n ro Ri í ver Dam R iv Town k. e e re r C Principal river
a is Island u L a í Reservoir r a . M k erey Cree e El M k Direct Influence Area Limit e r C Rivers o p a r Permanent River a u G Guri Dam Intermittent River Unknown Flow Regime
GURI RESERVOIR
Figure 2.Tocoma Reservoir site in the Caroni River The reservoir will cover an area of 85 Km2, and it will impound a volume of water of 1.758 Km3. Considering an annual average inflow of 4,250 m3/s, theoretical residence time has been estimated in 4.8 days
Comparative values of some characteristics related with morphometry of both reservoirs and their watershed are presented in Table 1.
This document is organized into five chapters. The first chapter corresponds to this introduction. Chapters 2 and 3 provide a brief explanation of the characteristics of each reservoir, the calibration and simulation of hydrodynamic and water quality processes for both reservoirs. Then, in the fourth chapter, the role of loadings and boundary conditions for each of these water bodies is studied by means of sensitivity analysis. Lastly, the fifth chapter is dedicated to conclusions, as well as opportunities and needs of additional research.
5 Table 1. Some morphometric characteristics of Lake Texoma, Tocoma Reservoir, and their watersheds. Characteristics Texoma Lake Tocoma Reservoir Area (at normal pool elevation) 264,901,462 m2 (1) 85,000,000 m2 (3) Volume (at normal pool elevation) 3,064,601,914 m3 (1) 1,758,000,000 m3 (3) Mean depth 11.6 m 14.7 m Maximum depth 29.3 m 77 m Drainage basin area 87,498 km2 (2) 88,067 km2 Theoretical Residence Time 255 days 4.8 days Mean annual inflow 139 m3/s (1) 4,261 m3/s Mean annual outflow 168 m3/s (2) 4,261 m3/s (1) USACE 2003b (2) Atkinson et al. 1999 (3) CVG DELCA 2005
6 CHAPTER 2
HYDRODYNAMIC SIMULATION
Modeling and simulation of flow and water quality in rivers and lakes is a representation of the natural ecosystem integrating the temporal and spatial variability of the physical, chemical, and biological characteristics of aquatic systems (Ford 1999).
Aquatic ecosystem models, and particularly those applied to reservoirs, provide an integrated analysis of the morphometry, hydrology, ecology, and internal and external forcing functions acting on a water body, based on monitoring physical, chemical, and biological features of such a water body (Tufford and McKellar 1999). Water quality measurements, normally made at a finite number of specific localities during specified periods of time, provide the information needed to infer process dynamics and their effects over much broader spatial and temporal scales using water quality models. Thus, the decision-making process and water quality management are supported on more integrated and complete information, taking into account not only observed but also expected water quality conditions.
Using simulation tools for studying a reservoir is particularly important because these impoundments may be more complex than natural lakes (Hakanson and Peters 1995,
Tufford and McKellar 1999). In a reservoir, river inflows determine a circulation pattern based on lake morphometry, defining advective transport and density gradients (Ford
1999). In addition, hydrology, land use, and weather establish temporal patterns and variations on water quality. As a consequence, the biotic component of this ecosystem responds selecting many different functional habitats (Kalff 2002). Water quality models
7 allow the integration of these processes and their evaluation in their whole dimension.
Moreover, a preliminary step in this process is the understanding of hydrodynamic
processes through simulation. In the case of studying a reservoir with multiple uses, this
approach is particularly useful and desirable.
The U.S. Environmental Protection Agency (USEPA), in conjunction with the
Tulsa District of the U.S. Army Corps of Engineers (USACE), has been studying Lake
Texoma and its watershed. This dissertation is a contribution to the study and represents
one step in developing a decision support system, based in environmental modeling, to
manage the water quality of that reservoir. Hydrodynamic studies of Lake Texoma are
not available; consequently, in this study the hydrodynamics of Lake Texoma were
evaluated using a simple hydrodynamic model. The DYNHYD model was selected for its
flexibility and simplicity, and for the fact that it can be coupled with the Water Analysis
Simulation Program (WASP) (Ambrose et al. 1993).
On the other hand, creating a reservoir for hydroelectrical generation can require
exhaustive evaluations to establish the magnitude and intensity of environmental effects
of the project. In this sense CVG EDELCA, a state-owned company in charge of
approximately 70% of the electrical generation in Venezuela, has promoted the
evaluation of Tocoma Reservoir and its environmental impact. This study is a part of the
environmental impact evaluation of the project.
The general objective of this chapter is to describe the hydrodynamics of both
Lake Texoma and the future Tocoma Reservoir using the simulation model DYNHYD.
The specific objectives include calibrating and validating DYNHYD predictions of these
water bodies, and estimations of the hydraulic travel time and residence time in those
8 reservoirs. This information will also be used as an input for the analysis of physical, chemical and biological processes that determine the water quality of these reservoirs, using the water quality modeling components of WASP, presented in the following chapter.
Model Description
The model selected for this research is DYNHYD5, the hydrodynamic model of the WASP system. WASP (Water Analysis Simulation Program) is a system integrated by two independent, computational programs: DYNHYD5 and WASP6. Both of these are dynamic compartment models, very flexible, and are used to simulate water transport and fate and transport of water quality constituents and toxic compounds. Particularly,
DYNHYD5 simulates hydrodynamics, based on the momentum and mass conservation equations. The model system is developed and supported by the US Environmental
Protection Agency (USEPA) and it has been used to study a variety of issues in many different water bodies, such as rivers, lakes, reservoirs, and estuaries, including, but not limited to, physical and ecological characterization, evaluation of the effects of human activities, and development of watershed management strategies (Lang and Fontaine
1990, Lung and Larson 1995, Branch et al. 1996, De Smedt et al. 1998, Tufford and
McKellar 1999, Lee 2000).
To simulate the hydrodynamic of a water body, DYNHYD solves the equations of motion and continuity in a computational network integrated by junctions and channels, which are defined for the system to be simulated. The model provides volume and velocity in every junction and channel of the hydrodynamic network, respectively.
Channels define the movement of water between junctions, and the junctions hold the
9 volumes where chemical and biological components of the water systems are going to be distributed. These values will be also used as input for the water quality model.
Lake Texoma: Description of the Study Area
Lake Texoma is a man-made lake located on the border between Oklahoma and
Texas. It was created in 1945 at the confluence of the Red and Washita rivers by the U.S.
Army Corps of Engineers in order to generate hydroelectrical energy (Breeding 1997) and for flood control. However, other beneficial uses have been considered and at the present time, it is also used for recreation, water supply, and fisheries. These uses provide a very important source of revenue for the local economy (USACE 2003a). Only about
5.5% of the storage capacity of the reservoir is available for consumptive usages, and the rest is used for energy generation. The present distribution of water volumes by different usage is presented in Table 2. The rest of usages are a consequence of consumptive usage of the water of the reservoir.
Table 2. Distribution of annual volumes in Lake Texoma for different usages (Compiled from TWDB 2003).
USAGES ADJUDICATED TO VOL/YEAR (m3) Red River Authority of Texas 2,342,700 City of Denison 30,085,200 Greater Texoma Utility Authority 18,495,000 Municipal Purpose North Texas Municipal Water District 95,310,900 State of Oklahoma 16,029 SUB-TOTAL 146,249,829 State of Oklahoma 5,395,608 Irrigation Red River Authority of Texas 308,250 SUB-TOTAL 5,703,858 Greater Texoma Utility Authority 12,330,000 Texas Utilities 12,330,000 Industrial - Mining Red River Authority of Texas 123,300 SUB-TOTAL 24,783,300 TOTAL 176,736,987
10 Denison Dam was built to create Lake Texoma, and it is located on the Red River,
at river mile 725.9, five miles northwest of Denison in Grayson County, Texas (USACE
1996b). Two large arms depict the physiography of this reservoir: Red River and Washita
River arms. Lake Texoma has a very dendritic shape with a shoreline estimated to be
more than 970 Km. Table 3 summarizes some characteristics related with the
morphometry of the reservoir and its watershed.
Table 3. Some morphometric characteristics of Lake Texoma and its basin. Area (at normal pool elevation)(1) 264,901,462.4 m2 Volume (at normal pool elevation)(1) 3,064,601,913.8 m3 Mean depth 11.6 m Maximum depth(1) 28.8 m Length of the shoreline (at 620 ft a.s.l.) (4) 971.2 km Drainage basin area(2) 87,498 km2 Mean annual inflow (3) • Red River at Gainesville (1936-2000) 93 m3/s • Washita River near Dickson (1928-2000) 49 m3/s 3 • Total(1) 185 m /s Mean annual outflow (1) 168 m3/s (1) USGS 1999 (2) Atkinson et al. 1999 (3) USGS 2005 (4) TWDB 2003
Watershed Characteristics
The Lake Texoma watershed extends to the Red River headwaters in New
Mexico, and the Washita River headwaters, located in NW Oklahoma (Figure 3). About
70 million years ago, much of the western reaches of the Red River, draining into Lake
Texoma, were covered with a shallow saline sea. Consequently, soils and surface geology
of the watershed provide significant amounts of calcium carbonate, calcium sulfate, and
sodium chloride, which determine most of the ionic composition of the water in the
reservoir today (RRA 2003b).
11 OKLAHOMA
O
C
I
X
E
M
EW
N TEXAS
Figure 3. Lake Texoma watershed (Adapted from Upton 2004) The local watershed of the reservoir encompasses Cooke and Grayson counties, in the Texas reach; and Carter, Love, Marshall, and Bryan counties in the Oklahoma reach
(Figure 4). Several small communities characterize much of the area of the local watershed of the reservoir. However, the area around the City of Denison has been reported as one of the fastest growing areas in Texas due to the overflow from the
Dallas/Fort-Worth area (RRA 2003a). Cities within the local watershed of Lake Texoma with a population over 10,000 include Denison and Gainesville. A significant number of farms that produce mainly hay, wheat, corn, peanuts and fruits are also located in the local watershed of Lake Texoma. Some farms in this area also raise beef cattle, poultry, and horses. Mining activity is also identified in this area and it includes oil and gas, sand, and gravel (RRA 2003a).
12 Figure 4. Lake Texoma and surrounding counties (Adapted from Upton 2004). Considering the geology of the region, adjacent watersheds of Lake Texoma are mainly located in the ecoregion of the Central Oklahoma-Texas Plain, integrated by sandstone and shale beds of the Triassic Age. Some Permian and Pennsylvania Age sandstone, shale, limestone, dolomite and gypsum are found in that ecoregion (Patrick
1994).
The topography surrounding Lake Texoma varies from gently sloping flats to rocky cliffs and steep, wooded hillsides. Many areas along the shoreline are forested, predominantly covered by different species of oak, hickory, American elm, and juniper.
Non-forested vegetation consists largely of many annual grasses, forbs, and wild legumes
(USACE 1996b).
13 Hydrology
The total drainage area of Lake Texoma is estimated to be 87,498 km2 (Atkinson et al. 1999). The surface area of the reservoir is 264.9 km2 at normal pool elevation (195 m a.s.l), and it holds a volume of approximately 3.03 km3 (USACE 2003a). The maximum depth has been estimated at 34 m (Atkinson et al. 1999), and the mean depth is calculated to be 11.6 m. The reservoir has a highly developed dendritic shoreline, with a length estimated to be more than 970 km.
The hydrodynamics of the reservoir are driven by two main inflows, from the Red and Washita rivers, and one main outflow at the Denison Dam. Daily inflow to Lake
Texoma coming from the Red River is measured near Gainesville, TX (USGS 07316000
- latitude 33o43’40”N, longitude 97o09’35”W) and the inflows of the Washita River are measured near Dickson, OK (USGS 11130303 – latitude 34o14’00”N, longitude
96o58’32”W). Variations in the mean monthly flow in the Red River, measured at
Gainesville gage station, are between 36 and 226 m3/s (period: 1936-2004). In the
Washita River, mean monthly flows at Dickson fluctuate between 18 and 117 m3/s
(period: 1928-2004) (USGS 2005). The highest variability of inflow to the reservoir occurs between May and June, and the lowest during August. These fluctuations are associated with variations in the water elevation and water quality of the reservoir, especially during spring time (Work and Gophen 1999). Other inflows come from numerous small tributaries draining to Lake Texoma that are not measured.
In addition, several small, intermittent creeks drain directly into the reservoir, providing additional inputs of water estimated at less than 10% of the total annual accumulated volume of the reservoir. The most important of the smaller inputs are
14 Hickory Creek and Big Mineral Creek. A list of some of the creeks draining directly to
Lake Texoma is presented in Table 4.
Table 4. Main creeks that drain directly into Lake Texoma Location Margin(1) Identified Creeks Hickory Creek Boggy Creek Wilson Creek Little Hauani Creek Left House Creek Red River arm Blair Creek Buncombe Creek Caney Creek Big Mineral Creek Right Mill Creek Little Mineral Arm Glasses Creek Washita River arm Right Little Glasses Creek (1): In reference to the downstream direction
Outflows of the reservoir are mainly at the dam. There, there are eight 20-ft diameter reinforced concrete conduits that are 800 ft long with an invert elevation of
523.0 ft. Five of them are dedicated for the delivery of hydroelectric power, and the other three are used for the required downstream releases and floodwater releases (USACE
1996b). Channel capacity below Denison Dam is approximately 1,700 m3/s. However, the mean annual outflow through the dam is estimated in 168 m3/s.
Model Parameterization for Lake Texoma
The simulation process involves parameterization or calibration of the selected model and validation process, including adjustment of initial conditions in the system to be simulated (Ford 1999). DYNHYD was set to describe water movement in the lake and periodic variations of pool elevations. Input parameters adjusted for DYNHYD to simulate the hydrodynamic of Lake Texoma included values associated with the
15 geometry of the reservoir, hydrology of the main inflows and outflows of the reservoir,
and meteorological information from weather stations located in the reservoir area.
The validation process was performed by estimation of the root mean square of the error (RMSE). The error is the difference between predicted and observed daily elevations of the reservoir, using reported values by USACE (2003b). The RMSE is a
measure of the differences between observed and predicted (Davis 2002):
2 ⎡ ⎛ ∧ ⎞ ⎤ ⎢∑⎜ yi − yi ⎟ ⎥ ⎢ ⎝ ⎠ ⎥ RMSE = ⎣ ⎦ , where: n
yi = Observed value
∧ yi = Predicted value
n = Number of paired values
The relative RMSE was calculated relating the RMSE with the average of the
observed values.
Calibration of hydrodynamic models is normally performed by adjusting
simulated and real values of water velocity (Ambrose et al. 1993). However, there were
no available measurements of water velocities for the reservoir. Consequently, calibration
of DYNHYD5 was performed by adjusting the estimated daily pool level of the reservoir to real pool level registers, generated under selected conditions of inflow and outflow.
The objective of this calibration process was to generate the best approximation between observed and simulated pool levels of the lake. A range of +/- 2.0 m. of the reported daily pool level value of a selected period of time was considered acceptable for this objective, using a dynamic condition for simulation. This number represents less than 10% of the
16 mean depth of the Main Lake Zone in Lake Texoma where water elevation of the
reservoir is measured. Similar studies have reported an RMSE less than 20 % of daily
pool elevation estimations (Crain et al. 2000, Tufford and McKellar 1999).
Parameterization of DYNHYD required considering parameters associated with
channel geometry, hydrology and meteorology (precipitation/evaporation balance).
Geometry of the reservoir was estimated using Geographic Information Systems (GIS).
Inflows were set to the measured values of daily mean discharge, using the information
from two USGS gage stations located at Gainesville, on the Red River, and at Dickson,
on the Washita River. Additional inflow was adjusted using runoff estimations for the
Big Mineral Creek sub-watershed (Upton 2004) or by a mass balance procedure.
Model Geometry
DYNHYD requires the definition of the network of channels and junctions for the
water body to be analyzed. The main objective of this segmentation process was to improve the understanding of the hydrodynamics in the reservoir and its relationship to
reported lake zonation. Atkinson et al. (1999), Clyde (2002), and Mabe (2003) have
reported the existence of five zones in Lake Texoma, determined by the chemistry of the
water. These zones are: the Red River arm, the Washita River arm (considered as lotic
zones), the Red River Transition Zone, the Washita River Transition Zone, and the Main
Body (lacustrine zone). The network of channels and junctions for the hydrodynamic modeling of Lake Texoma was developed taking into account this zonation. In addition, this segmentation was developed considering smaller sub-watersheds around the lake.
The bathymetry of the reservoir was used to develop the model segmentation needed for the hydrodynamic and water quality simulation models. The number of
17 segments that the system to be simulated can have is constrained, in part, by the scale at
which reasonable bathymetric estimates can be made (Tufford and McKellar 1999).
Some bathymetric surveys have been made in the reservoir recently (USGS 1999,
TWDB 2003). In the study performed by the USGS (1999), a Digital Elevation Model
(DEM) file of the lake and its surrounding area with a grid of six by six meters each cell
was produced. This file was used to define the network of channels and junctions for the
hydrodynamic simulations as well as determine the required segments parameters (Figure
5).
Figure 5. Image of the DEM file for Lake Texoma showing different depths. Lighter colors correspond to shallower areas (Adapted from: USGS 1999) Surface areas, volumes, and maximum depth for each segment designed for the
hydrodynamic network of the lake were adjusted to this source of information and were
calculated using GIS software ArcGIS 8.0TM (Environmental Systems Research Institute
Inc., Redlands, CA). In addition, because the model assumes every junction and channel to be a rectangular receptacle or conveyor, and because the geometry of the network
18 determines the time step for the simulation, consideration on water depths, length and surface area for each junction and channel were also made as criteria for segmentation of the reservoir.
After combining all mentioned constraints, the resulting network was set up in a one-dimensional framework and it resulted in 40 junctions and 39 channels (See Figures
6 and 7). Data on the geometry of each channel and junction are presented in Appendix A.
21 22 24 23
25 26 27 3 2 4 28 1 5 29 6 18 19 20 30 7 17 31 8 10 9 16 32 11 15 33 40 12 14 39 13 34 38 37 36 35
# = Junction Number
Figure 6. Junction network of the hydrodynamic model of Lake Texoma
Meteorological and Hydrologic Data
Three tributaries were defined in this simulation: Red River, Washita River, and
Big Mineral Creek. A review of the available literature showed that hydrological
19 measurements of the Washita and Red rivers have been performed for more than 60 years
continuously at Dickson, OK and Gainesville, TX, respectively (USGS 2005). These
locations represent the main inlets of water to the lake.
21
22 23 24 25 26 2 3 1 4 27 28 5 19 18 29 6 17 20 7 30 9 16 31 8 10 15 11 14 32 33 39 12 38 13 37 36 35
34
# = Channel Number
Channel
Figure 7. Channel segmentation of Lake Texoma for the hydrodynamic model Additional inputs for the reservoir were adjusted using estimated runoff for the
Big Mineral Creek sub-watershed (Upton 2004). Another estimation of additional inputs was made using the results of a mass balance of volumes in the reservoir (Tortorelli
personal communication), considering accumulated volume during a day period, reported
by USACE (2003b), and precipitation and evaporation data registered at Denison Dam weather station. During the calibration process it was determined that approximately one half of the accumulated volume during a 24-hour period, estimated from the reported
20 pool level by USACE and the lake stage-volume curve, and the consideration of
evaporation and precipitation daily data from Denison Dam station provided the lowest
error in pool level estimations of the reservoir and a good representation of the variability
of daily pool elevation of the reservoir (See Figure 8). The boundary downstream of the
system can be set both stipulating outflows and comparing resulting water elevations, or
defining surface elevation and using predicted velocities for validation. For this
simulation, the downstream boundary was defined specifying outflows in order to get
estimations of daily pool elevations and these numbers were used in the validation
process.
1.40 3.00
1.20 2.50
1.00 2.00
0.80 1.50 0.60 Absolute Error (m) RMS of Errors (m) 1.00 0.40
0.50 0.20
0.00 0.00 0 253040455060100 % additional inflow at Big Mineral Creek
RMS Max.Abs.Err.
Figure 8. Variability of error in pool elevation estimations of Lake Texoma using different percentages of accumulated volume in a 24-hour period as inflow in Big Mineral Creek and applying DYNHYD Meteorological data were also available for several years, in a daily frequency, for
some rain gage stations around the lake (NOAA 2003), as well as at Denison Dam
(USACE 2003b). A daily precipitation/evaporation balance was performed and was
21 included as an input of the model (Figure 9). Wind directions and wind velocities were not considered in this simulation.
The time step selected for this simulation was 60 seconds, based on the Courant condition (Ambrose et al. 1993). Considering this condition, the time step (Dt) should be less than the shortest channel length divided by the wave celerity, obtained by the equation:
L Δ=t i , where: ()gyii± U
Li = Length of channel i, m
yi = mean depth of the channel i, m
Ui = velocity in channel i, m/sec
g = acceleration of gravity (9.8 m/sec2)
The minimum length of the channels in the developed network was approximately
1,200 m, and the mean depth of the channels ranked between approximately 2 and 22 m.
Using different values of water velocity in the channels, various values for the time step were calculated, as shown in Table 5. The lowest calculated value for the time step was
85 sec, higher than the selected time step for this simulation.
Table 5. Estimation of maximum time step (in seconds) for different conditions of water velocity and depth of the channels in the developed hydrodynamic network for Lake Texoma
Water Velocity (m/s) Depth 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 2 270 270 269 269 268 267 267 266 5 171 171 171 170 170 170 170 169 10 121 121 121 121 121 120 120 120 15 99 99 99 99 99 98 98 98 20 86 86 86 85 85 85 85 85
22 March/1997
35
30
25
20
15
10
mm/day 5
0 1 2 3 4 5 6 7 8 9 10111213141516171819202122232425262728293031 -5
-10
-15 April/1997
35
30
25
20
15
10
mm/day 5
0 1 2 3 4 5 6 7 8 9 101112131415161718192021222324252627282930 -5
-10
-15 May/1997
35
30
25
20
15
10
mm/day 5
0 1 2 3 4 5 6 7 8 9 10111213141516171819202122232425262728293031 -5
-10
-15
Figure 9. Precipitation/Evaporation balance between March,1997-May,1997 using daily reported values for Denison Dam rain gage station
23
Two 92-days periods were simulated. One of them corresponded to March, 1997 to May, 1997, considered a period of high inflow variability (Figure 10); and the second one was between July, 1997 and September, 1997, when the main inflows were in the lowest hydrologic condition. These two ranges were selected because they provided the most complete, longest, and continuous set of physical, chemical and biological observations. Validation was performed for the interval between March, 1997 and May,
1997.
Figure 10. Boxplot of mean daily total inflow in Lake Texoma during 1997 represented by month. The highest dispersion of values is observed in May and the lowest occurs in August
24 Hydraulic Travel Time and Residence Time
To complete the analysis of the hydrodynamics of Lake Texoma an estimation of
the hydraulic travel time and the residence time was performed. The hydraulic travel time
is considered as the amount of time that it takes for some volume of water to move
between the segments of the hydrodynamic network, considering the flow. The residence
time is understood as the amount of time that it takes for a conservative component to
travel through the water body during a specific hydrologic condition (Tufford and
McKellar 1999).
Estimations of the theoretical hydraulic travel time can be performed using the
equation t = V/Qin, where V represents the volume of the reservoir at normal pool
elevation and Qin represent the annual average of the total inflow. However, it also can be
evaluated using hydrodynamic simulation performed in a steady-state condition and also
using the average flow for each river. Using this approach it is possible to identify portions of the considered water body where the highest hydraulic travel time is expected.
Residence and travel time were calculated using a hydrodynamic simulation under
steady-state setting at average hydrologic conditions for each river. Mean flows in each
channel were simulated and these values were related to volumes. Estimations of travel
time were performed using the water quality model and introducing a simulated pulse of
a tracer (i.e., a conservative constituent) as a point source load at the first segment of the
water quality model for each inflow on the first day of the simulation period. The
difference between the day of occurrence of the peak of tracer concentration in water
quality segment 31 (closest to the dam) and the first day of application of the pulse was
considered the travel time.
25 Hydrodynamic Modeling Results for Lake Texoma
Calibration of the hydrodynamic model using estimated runoff for Big Mineral
Creek resulted in a maximum difference between observed and predicted elevation of
2.03 m, with an overestimation of the reservoir level during the period of highest flow variability (Figure 11). The estimated RMSE in this prediction was 1.06 m.
Predicted Observed
Figure 11. Graphical comparison between observed and predicted elevations in Lake Texoma during March, 1997 and May,1997, using DYNHYD and estimated runoff for Big Mineral Creek provided by Upton (2004)
Setting Big Mineral Creek flows as a half of the daily accumulated volume in the reservoir which is not explained by the Red and Washita rivers inflows, in a mass balance estimation provided a better estimation of daily elevations of the reservoir, with a maximum difference between observed and predicted elevation of 1.93 m (Figure 12). In this case, the estimated RMSE was 0.94 m. During low flow conditions, both simulations provided similar results. This suggests that some additional losses should be included in the runoff model to adjust flows for the Big Mineral Creek sub-watershed during high precipitation condition. However, both methods provided results considered to be in the accepted range of lake elevation values for this simulation.
26 One aspect to take into account is the fact that these estimations assume negligible
seepage to groundwater and therefore the difference is attributed entirely to surface
inflows. This assumption was made to simplify the analysis and it does not affect the
resulting prediction of bulk movement of the water column.
Predicted Observed
Figure 12. Graphical comparison between observed and predicted elevations in Lake Texoma during March, 1997 and May,1997, using DYNHYD and a 50% of the daily accumulated volume in the reservoir not explained by the Red and Washita rivers inflows based on mass balance estimations as a daily inflow of Big Mineral Creek. In general, the model tracks daily variations of the elevation of the reservoir closely. Elevation of the reservoir directly responds to a combination of periods of
highest discharge and regulation of the outflows. This corresponds to the function of
flood control assigned to this dam. Precipitation/evaporation balance is also an important
factor in estimates of the reservoir elevation. If this factor is not considered,
overestimations of water elevation are expected. The effect of this parameter is similar
during both conditions, high variability of flow or lowest inflow with low variability.
The model also suggests the existence of a longitudinal gradient of velocity.
Water velocity suddenly declines once the Red River gets into the lake. This effect is not
observed in the Washita River, probably because of the effect of the Cumberland Pool,
located upstream of the boundary set for this simulation. Results of the steady-state
27 simulation seem to support the idea proposed by Atkinson et al. (1999), Clyde (2002) and
Mabe (2002). However, the identified zones in this hydrodynamic simulation are slightly different from those reported by Atkinson et al. (1999). A clear lotic zone is identified in
the proximity of the entrance of the Red River, and four lentic zones can be defined as the
Red River transition zone, the Washita River transition zone, the Main Lake zone, and
the Big Mineral Creek zone, with the lowest simulated velocities and even with some
backflow during peak flows of the Red River (Figure 13). A lotic zone for the Washita
River Arm is not clearly identified with this simulation.
Hydraulic Travel Time
Using the theoretical method to estimate the hydraulic travel time for Lake
Texoma, it resulted in approximated 255 days. Using the proposed hydrodynamic
segmentation, the largest travel time per segment is about 56 days, and it occurs in the
Big Mineral Creek area, where some back flows and lower velocities occur due to the
damming effect generated by Red River flows. Figure 14 presents a spatial zonation of
predicted water velocities for the reservoir. The average travel time is estimated around
15 days per segment, which is similar to the average hydraulic travel time in the
transition zones. Estimations of travel time in the Red River zone are the lowest and the
average is around 3 days. In transition zones the average travel time varies between 5 and
35 days, and in the Main Lake area the average travel time per segment is around 19 days.
Results of these estimations for each segment are presented in Appendix B.
28 Residence Time
A simulated pulse of load of 500 Kg/d of a tracer was applied in segment 1 during the first day of the simulation. It took 78 days to reach the maximum concentration in segment 31, using a dynamic condition with high flow variability. Dilution process and the effects of transport processes are also examined in this simulation (See Figure 15).
Figure 13. Variation of water flow in different zones of Lake Texoma. 0 to 1 corresponds to the Red River Zone, 19 to 28 corresponds to the Red River Transition Zone, 0 to 20 corresponds to the entrance of Washita River to Lake Texoma, 27 to 28 represents Washita River Transition Zone, 28 to 29 and 31 to 0 represents the Main Lake Zone, 0 to 32 corresponds to the entrance of Big Mineral Creek to Lake Texoma, and 10 to 36 is the junction between the Big Mineral Creek and the Red River Transition Zone.
29
Figure 14. Spatial distribution of predicted water velocities for Lake Texoma under the steady-state condition The same procedure was applied using a simulated pulse loading of 150 Kg/d of a tracer in segment 20, at the boundary in the Washita River Arm (See Figure 16).
Residence time in this case resulted shorter, considering also the longitude of this portion of the lake. The length of the Washita River Arm was estimated to be 30.8 Km, while the length of the Red River Arm is 54.9 Km, so the path of the tracer in this portion of the reservoir was almost a half of the path through the Red River Arm. It took about 48 days for the peak of concentration of the tracer to reach segment 31.
30 Figure 15. Graphical representation of travel time in Lake Texoma when a pulse of tracer is applied in the upper portion of the Red River arm A combined loading of a tracer next to both boundaries Washita and Red River was also tested and there is an overlap of the effect of both loadings. However the maximum peak is reached approximately 77 days after the occurrence of loading in segment 1. This suggests that residence time during high-variable-flow conditions is mainly driven by hydrologic conditions in the Red River (See Figure 17).
31 Figure 16. Graphical representation of travel time in Lake Texoma when a pulse of a tracer is applied in the upper portion of the Washita River arm
32 Figure 17. Estimation of travel time for Lake Texoma using a combined loading at the Red and Washita River
Tocoma Reservoir: Description of the Study Area
Tocoma Reservoir has been projected between miles 368 and 379 in the Caroni
River, Venezuela. Caroni River is a blackwater river, right-sided tributary in the lower reach of the Orinoco River, in Venezuela, with a total length of 431 miles approximately.
It flows to the Orinoco River in a North-South direction and its watershed is about 92,170
Km2. Among the very last 62 miles of the Caroni River, the most important hydroelectrical company in Venezuela, CVG EDELCA, has built three reservoirs for hydroelectrical generation: Guri, Macagua and Caruachi. Tocoma Reservoir is going to be the fourth and it will be located between two of these three reservoirs: Guri and
Caruachi, as presented in Figure 18. Consequently, the very last portion of river conditions in the lower Caroni will be flooded once the reservoir is created.
33 GURI DAM
250
200 TOCOMA CARUACHI 150 DAM MACAGUA 100 DAM Orinoco Elevation (masl) Elevation (masl) 50 River
0 Caroní River
100 90 80 70 60 50 40 30 20 10 0
Distance from the Orinoco River (km)
Figure 18. Scheme of hydroelectrical development in the Lower Caroni River (Adapted from CVG EDELCA 2005b)
Hydrograph and Hydrology
The total length to be flooded for the construction of the reservoir extends from
Guri Dam discharge channels to the La Pollera Island, a place located 17 Km downstream of Guri Dam, on the Caroni River (Figure 19). Under current conditions, water flows from Guri Dam for approximately 6 Km in a very straight pattern through a narrow canyon named Nekuima Canyon, approximately 200 to 500 m wide. Then, the river opens widely, reaching approximately 7 Km wide, as a result of a geologic fault.
Consequently, numerous islands and rocks emerge at this portion of the river, developing
rapids and small falls. Once the fault area ends, the river gets narrower again, with
approximately 3.5 Km wide. At this point a 6,250 m length concrete dam is going to be
erected to create a reservoir with a total area of 85 Km2.
The total drainage area of the future reservoir is 88,064 Km2, and the immediate
drainage area of the new reservoir has been estimated on 400 Km2. Other tributaries of
34 Tocoma Reservoir will be Cunaguaro River, and other three small creeks: El Merey,
Maria Luisa, and Guarapo. RESERVOIR CARUACHI
Tocoma Dam Site
C Cun aguaro ar River o n í R iv e r
. k e e r C Dam
a s i Principal river u L erey Creek Future flooded area a El M í r a . Island M k e e r Reservoir C Rivers Permanent River Intermittent River
Guarapo Guri Dam Unknown Flow Regime
GURI RESERVOIR
Figure 19. Geographic location for the future Tocoma Reservoir, in the Caroni River The main source of water for the new reservoir will come from the discharge of
Guri Reservoir. Since 1986, this reservoir controls the discharges of the lower portion of
the Caroni River. Therefore, most of the time there is no seasonal flow fluctuation, except
for those years with high precipitation in the watershed. The mean annual flow
downstream of Guri Dam is 4,630 m3/s (Period: 1986-2004) (CVG EDELCA 2005a), and
it is almost constant during the year. Only for very rainy years, the mean daily flow can ascend to 12,500 m3/s. A comparison of mean monthly flows discharged through the
turbines and the spillway during a typical year is presented in Figure 20.
35 4000
3500
3000
2500
2000 Flow (cms)
1500
1000
500
0 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
Turbines Spillway
Figure 20. Mean monthly turbine-used and spilled flow at Guri Dam discharge channels under average conditions However, the flow of those small tributaries draining directly to the Tocoma
Reservoir area fluctuates considerably during the rainy season, with monthly daily mean
values ranging between 1.2 and 18 m3/s (CVG EDELCA 2005a). The annual inflow of these small creeks directly draining to Tocoma Reservoir has been estimated on 11 m3/s
(Acevedo and Sherman 2002).
Mean annual rainfall has been estimated at 1,270 mm. The highest precipitation
occurs between June and August, and the lowest between March and April, as presented
in Figure 21. Differently from what happens to those small rivers draining directly to the
future reservoir, the main inflows would not vary significantly during the rainy season
because of the flow regulation that Guri Dam exerts at the lower portion of the watershed.
36 250
200
150
Rainfall (mm)
100
50
0 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
Figure 21. Average monthly precipitation at Las Babas weather station, close to Tocoma Reservoir site (Period: 1986-2004)
Model Parameterization for Tocoma Reservoir
Because the Tocoma reservoir has not been built yet, the hydrodynamic model
was examined as an expected quantitative prospective of potential water flow behavior.
For Tocoma Reservoir, model parameterization was based on validations of DYNHYD results describing river conditions at this portion of the Caroni River, as described by
Acevedo (2002) and Acevedo et al. (2002). Simulated velocities were compared with
measured values. Then, using topographic information for the area to be flooded by the
reservoir, and some bathymetric cross-sections of the Caroni River, a digital elevation model (DEM) was built for the new reservoir and a new hydrodynamic network was created. These results provided the parameters associated with channel geometry for describing the hydrodynamic network of the new reservoir. Inflows were set to the
37 monthly mean flow discharged by Guri Dam, using the information from the
hydrological data base of CVG EDELCA (CVG EDELCA 2005a). Additional parameters for the hydrodynamic simulation, such as Manning roughness coefficient, came from the literature review (Maidment 1992), considering topography and vegetation information.
Model results were evaluated for continuity and they should be compared with real values once the reservoir is completed and the area flooded.
Model Geometry
Based on the model network developed by Acevedo (2002) for river conditions of
the Caroni, a new water network was created, but this time considering the total area to
be flooded by the reservoir. According to project specifications, water level for the
simulation was set constant at 127 m.a.s.l. The approach used for developing the
hydrodynamic network described the flow in a quasi two-dimensional manner. This is to
say, the flow was considered to have a predominant one-dimension condition. However,
because of the topography of the zone, and the existence of multiple islands, some lateral flow components were also considered, resulting in a two-dimension network. Smaller sub-watersheds around the reservoir were also considered for defining location of junctions of the network.
The bathymetry of the new reservoir was built based on topographic maps, available at 1:10,000 scale, with 5-m contours. A digitalization of these data was
performed and a Digital Elevation Model (DEM) file for the new reservoir was created
using GIS software ArcGIS 9.0TM (Environmental Systems Research Institute Inc.,
Redlands, CA). This information was complemented by cross-section bathymetric data of the river acquired in previous studies (CVG EDELCA, 1987). The location of the cross
38 sections are indicated in Figure 22, and the resulting bathymetric model of the reservoir is presented in Figure 23.
In similar manner to the Lake Texoma parameterization previously described, additional criteria for developing the hydrodynamic network for Tocoma Reservoir were: homogeneous water depths, length, and surface area for each junction and channel of the resulting network. However, by combining these criteria, junction areas ranged between
267,973 and 667,290 m2.
1501 1501 16
872000
17
16
17
1701 1701
1702 1702 19 18
1902 1901 1904 1903 1905
BATHYMETRIC CROSS-SECTION
Figure 22. Schematic location for bathymetric cross-sections in the Caroni River for the Tocoma Reservoir area (Adapted from CVG EDELCA 1987)
39 The hydrodynamic network for Tocoma Reservoir resulted in 52 junctions and 66
channels, approximately twice than those described for Acevedo (2002) for this portion of the Caroni River functioning as a river. The resulting network is presented in Figure 24,
and a visual comparison of both, river and reservoir networks, is presented in Figure 25.
Data on geometry of each channel and junction are included in Appendix C. Values for
channel length and width also came from the analysis using ArcGIS9.0TM (Environmental
Systems Research Institute Inc., Redlands, CA).
Figure 23. Bathymetric map for Tocoma Reservoir
Hydrologic Data
The main inflow for the reservoir will come from the discharge of Guri Dam. At
this dam, there are two power houses. These two function differently from each other. In
40 addition, they release their discharges in two channels, one natural and one excavated, separated from each other by a rocky promontory. Average monthly flows estimated as a function of the energy produced for each power house were used for setting up the inflow for the new reservoir, considering data from 1986, when Guri started to work under these conditions. Two additional inputs, representing El Merey Creek and Cunaguaro River, were also considered. A constant level of 127 m was set as reservoir water level, considering that the reservoir will function as a constant-level reservoir. Meteorological data were not included under this analysis.
Island
Figure 24. Hydrodynamic network for Tocoma Reservoir The time step was set up considering the Courant condition (Ambrose et al. 1993), as explained for the hydrodynamic model calibration of Texoma Reservoir. A 20-second
41 time step was used, taking into account that the minimum length of the channels in the hydrodynamic model network developed for Tocoma Reservoir was approximately 700 m, and the mean depth of the channels ranged between 4.5 and 49.5 m. The lowest calculated value for the time step was 25, as presented in Table 6.
Table 6. Time step values for different junction depths of the hydrodynamic network for Tocoma Reservoir Depth Water Velocities (m/s) (m) 0.1 0.5 1 2 3 4 5 6 4.5 104 98 91 81 73 60 60 55 5 98 93 87 78 70 64 58 54 10 70 67 64 59 54 50 47 44 15 57 55 53 49 46 43 41 39 20 50 48 47 44 41 39 37 35 25 44 43 42 40 37 36 34 32 30 41 40 39 37 35 33 32 30 35 38 37 36 34 32 31 30 29 40 35 34 34 32 31 29 28 27 45 33 33 32 30 29 28 27 26 50 31 31 30 29 28 27 26 25
The model was set up as a steady-state condition and two flow conditions were simulated: 4,250 and 12,000 m3/s. These two conditions were selected based on the average discharge for Guri Dam and the maximum discharged allowed for this dam while hydroelectrical projects downstream of Guri are under construction. Validations for river conditions were performed using the flow reported for those days when velocity measurements were made at Tocoma dam site.
42 Tocoma Dam 32 site
31 30
29 28 27 25 26 24 23 21 22 19 20 18 17 16 1214 13 15
11 9 8 10 7 6
5
4
3
1 2 Guri Dam
Tocoma Dam site
55 40 52 45 39 53 38 44 49 54 46 25 23 37 50 22 21 36 51 43 47 20 24 19 35 48 42 18 34 41 33 17 27 26 16 31 15 32 13 30 29 28 14 10 12 11 9 8 7 6 5
4
3 2 1 Guri Dam
Figure 25. Comparison between hydrodynamic network for river and reservoir condition for Caroni River at Tocoma Reservoir site
43 Hydraulic and Travel Time
Residence and travel time were calculated under the steady-state condition for
both, as a river and as a reservoir, considering two hydrologic conditions: average and high flow. Estimates for travel time were performed aided by the water quality model, introducing a conservative substance as a tracer at the first water quality segment of the network. The difference between the day of occurrence of the peak of tracer’s concentration at the lower boundary of the network and the day of application of the load was considered the travel time. Estimates for these conditions, both river and reservoir conditions, were compared to establish the net change resulting by the construction of the reservoir on the hydrodynamics of the system.
Hydrodynamic Model Results for Tocoma Reservoir
Water velocity data collected during 2002 and 2005 were used to assess the “fit” of the DYNHYD model under river condition to the portion of the Caroni River where
Tocoma Reservoir will be created. During this period of time, water discharges from Guri
Dam did not substantially change. Therefore, it was not possible to evaluate the performance of the model under high flow conditions.
Model estimates for water velocity for channel 38 of hydrodynamic network under river condition were compared with observed values, and the RMSE was used as a measure of model fit. An RMSE of 9.63% (less than 20%) was achieved when simulating river conditions, and it was considered acceptable for these estimations. These data are presented in Table 7.
44 Table 7. Predicted and measured velocities for Tocoma under river condition Date Flow Observed Velocity Predicted Velocity (m3/s) (m/s) (m/s) 09/05/2002 5,242 0.63 0.77 09/10/2002 4,794 0.77 0.71 09/24/2002 4,453 0.65 0.66 10/01/2002 4,226 0.66 0.63 10/05/2002 3,967 0.60 0.60 10/22/2003 4,165 0.62 0.62
As a result of the simulation for Tocoma Reservoir, estimates of flow velocity and volumes for each channel and junction of the hydrodynamic network were obtained. The model suggests the reduction of water velocities in at least one order of magnitude at
Nekuima Canyon as a reservoir, when compared with simulation results for river conditions with average inflow (See Table 8).
Highest water velocity values are expected in the central channel of the new reservoir. For those littoral areas associated inflows of creeks and small rivers, like
Cunaguaro and El Merey, the model indicates that water velocities will have very low values. A graphical representation of water velocities results is presented in Figure 26.
Hydraulic Travel Time
A theoretical estimation of the hydraulic travel time for Tocoma Reservoir resulted in 4.8 days, considering the average inflow and the total volume for the reservoir.
Using the proposed segmentation, the largest travel time is about 32 days per segment, and it occurs in the Cunaguaro River area, where some back flows are expected, considering the amount of flow in the main channel of the reservoir in comparison with the flow of this small river. The average travel time is estimated in 1.7 days per segment, approximately, with the lowest values in segments located in the Nekuima Canyon and
45 followed by those in the central channel of the reservoir. Results of these estimations are presented in Appendix D.
Table 8. Relative change of water velocities for Caroni River when comparing river and reservoir condition at Tocoma site
River Condition Reservoir Condition Simulated Simulated Percent of Segment ID Water Velocity Segment ID Water Velocity Change (%) (m/s) (m/s) 1 3.45 4 0.82 23.81 2 5.07 5 0.50 9.88 3 5.04 6 0.20 4.04 4 3.66 9 0.11 2.89 5 1.75 11 -12 0.07 3.77 6 0.96 28 0.04 4.63 7 1.02 13 0.03 3.27 8 2.21 9 -14 0.08 3.65 9 0.93 29 - 30 0.01 1.58 10 0.81 31 0.03 3.55 11 1.32 15 - 26 0.03 2.51 12 1.09 27 0.02 1.86 13 1.11 15 -16 0.03 3.12 14 1.01 16 0.03 2.51 15 1.07 33 0.01 0.70 16 0.82 41 0.04 4.55 17 0.53 19 0.02 3.82 18 0.88 18 0.03 2.96 19 0.52 24 0.02 3.60 20 0.65 42 0.05 7.02 21 0.49 47 0.01 2.21 22 0.6 28 - 51 0.03 4.97 23 0.53 50 0.01 2.34 24 0.78 43 0.05 6.08 25 0.38 46 0.01 3.51 26 0.44 49 0.01 3.09 27 0.73 44 0.05 6.81 28 0.6 45 0.04 6.69
46 Figure 26. Spatial variation of simulated velocities for Tocoma Reservoir
Residence Time
A simulated pulse of load was applied at the discharge channels for Guri Dam which represents the first segment of the water quality network. It took 2 to 5 days to register the concentration peak at segments next to the future dam when the treatment was applied to segment 1. A graphical representation of temporal variation of the loading at these segments is presented in Figure 27. These values represent the residence time for the future reservoir. When testing the application of load at Cunaguaro River, it took 23 days for registering a peak concentration of the tracer at the end of the water quality network, indicating a high residence time for this portion of the new reservoir. Figure 28 shows temporal variation of tracer concentration when applying the treatment at
47 Cunaguaro River area. The difference in magnitude of these two values can be understood as a measure of the environmental sensitivity that Cunaguaro River arm at the future reservoir will have.
Tracer (ppt)
Segment 1 Segment 41 Segment 48
Figure 27. Peak concentration of the tracer at segments close to the future dam, in Tocoma reservoir, when a pulse is applied at Guri Dam discharge channels.
48 a)
b)
Figure 28. Temporal variation of a tracer concentration when applied at a) Cunaguaro River, and peak concentration of the tracer at b)the segment next to the future dam, in Tocoma reservoir.
49 CHAPTER 3
RESERVOIR WATER QUALITY MODELING
Water quality has become one of the topics of major environmental and health concern in recent years (De Villiers 2002). Traditionally, monitoring water quality in rivers and lakes has been used as the primary source of information for making decisions relating to the management of water bodies. Empirical analysis of water quality data has resulted in useful descriptions of aquatic ecosystems (Hakanson and Peters 1995). This approach, however, has been considered insufficient for prediction and quantification of secondary effects, especially those associated with human activities and their effects on water bodies, and the sustainability of these important resources (Goodwing and Hardy
1999). A broad range of variables and relationships among water resources and their usages are considered before defining a management objective for water bodies. This approach requires extensive data and simulation tools to assist in the decision-making process for managers and communities in order to select management strategies that support the conservation of the water system under consideration (Goodwing and Hardy
1999). Consequently, developing predictive models for rivers, reservoirs, and lakes is a powerful tool in practical ecology and environmental management.
Natural processes and anthropogenic effects in constructed impoundments have been better understood using simulation tools considering that these impoundments may be more complex than natural lakes (Hakanson and Peters 1995, Tufford and McKellar
1999). This is important to be considered when evaluating the environmental impact associated to a reservoir creation, like the Tocoma Reservoir. In the case of studying a
50 reservoir with multiple uses, such as Lake Texoma, this approach is particularly useful
and desirable.
The specific objectives of this chapter are to describe: the procedure used in the
calibration and validation of WASP6 for both reservoirs, Lake Texoma and Tocoma
Reservoir; and the evaluation of predictions for some chemical components for these two study cases. For Lake Texoma, predictions for chlorides, total suspended solids and dissolved oxygen were performed using the selected water quality simulation model; and dissolved oxygen-biochemical oxygen demand estimates were analyzed for Tocoma
Reservoir with the same model.
Model Description
The model selected for this research was the Water Analysis Simulation Program,
WASP (Ambrose et al. 1993). WASP is a dynamic compartment model, which is very flexible, and is used to simulate fate and transport of water quality constituents and toxic compounds. WASP6 simulates transformation of water quality constituents based on the mass conservation equation, using two kinetic sub-models: EUTRO and TOXI, the former is related to the problem of reduced dissolved oxygen content and eutrophication processes, and the latter is related to pollutants and toxic substances in the water
(Ambrose et al. 1993). The model system was developed and is supported by the US
Environmental Protection Agency (USEPA) and has been used to study a variety of issues in many different water bodies, such as rivers, lakes, reservoirs, and estuaries.
These studies typically provided physical and ecological characterization, evaluation of the effects of human activities, and development of watershed management strategies
51 (Lang and Fontaine 1990, Lung and Larson 1995, Branch et al. 1996, De Smedt et al.
1998, Tufford and McKellar 1999, Lee 2000).
Within TOXI, organic and inorganic chemicals, metals, and sediments can be modeled. A chemical tracer can also be analyzed with this sub-model of WASP by specifying no decay. In this study, the TOXI sub-model of WASP was used to simulate the processes associated with chlorides and total suspended solids concentrations in Lake
Texoma. The model was set to predict concentrations of these variables in the reservoir by simulating: the inputs of flow and loads of suspended solids and chlorides, the advective and dispersive transport of these components, and the processes of sediment deposition. Therefore, the sub-model TOXI was set to describe kinetic transport of one chemical without degradation (chlorides) and three types of solids (clay, silt, and sand) within multiple water quality segments and their corresponding benthic layers. A sub- benthic layer covering the entire bottom of the reservoir was also defined in this simulation (Figure 29).
In EUTRO, several processes determine the transport and fate of nutrients, phytoplankton, carbonaceous material, and dissolved oxygen in a water body. These processes can be physical, chemical and biological. EUTRO was used to simulate dissolved oxygen (DO) and biochemical oxygen demand (BOD) for Lake Texoma and
Tocoma Reservoir. Simulations may incorporate phosphorous and nitrogen cycle, and phytoplankton kinetics in a subsequent effort.
Four complexity levels are included in EUTRO to simulate carbonaceous biochemical oxygen demand (CBOD) and DO. In this dissertation the simplest complexity level of EUTRO was used considering that this was the first approach of
52 evaluation of water quality for these water bodies; that is to say, a system for dissolved
oxygen balance that solves the Streeter-Phelps BOD-DO equation in a slightly modified
form, allowing time-variable temperatures to be specified. Under this approach, sources
of oxygen are given by reaeration. The sinks are mainly driven by total biochemical demand. A salinity function for Lake Texoma was included considering the natural brackish condition of the system and its effects on oxygen solubility.
Configuration of the water quality network for Lake Texoma in EUTRO was the same used for TOXI. Water quality network for Tocoma Reservoir consisted on multiple water quality segments with a single benthic layer covering the entire bottom surface of the future reservoir (Figure 30).
Water Quality Segment B 1 2 3 Flow
Benthic Layer
Sub-benthic Layer
B = Boundary # = Water Quality Segment ID
Figure 29. Graphical representation of the configuration of water quality segments for the application of WASP for Lake Texoma
53 Inflow
Water quality B 123 Flow segments (Surface Water)
Sedimentation and Resuspension
Benthic Segment
B = Boundary # = Water quality segment ID Figure 30. Graphical representation of the configuration of water quality segments for the application of WASP for Lake Texoma
Previous Water Quality Studies of Lake Texoma
Studies on the water quality of Lake Texoma include a variety of issues.
Chemistry of the water has been one of the topics of major concern in this reservoir, considering the value of the reservoir as a source of water supply and the revenue associated with the production of striped bass (Morone saxatilis) (USACE 2003a).
Analysis of data resulting from intensive monitoring of physical and chemical characteristics of the reservoir reflects a spatial zonation based on its conductivity and turbidity (Atkinson et al. 1999, Waller et al. 2002). Four zones were defined in those studies: the Red River Arm and the Washita River Arm, considered as lotic zones, and the Red River Transition Zone, the Washita River Transition Zone, and the Main Body, considered as lentic zones. Other studies also consider the Big Mineral Creek as an additional zone because of its differentiated characteristics (Mabe 2002).
54 Additional studies on water quality of the reservoir and its relationship to land use
of the shoreline and the watershed also show that there are biological gradients, affecting
biological productivity of the reservoir (Work and Gophen 1999, An et al. 2002, An and
Kampbell 2003), and the effects of additional human influence in this water body are still
unknown.
Biological aspects that have been analyzed in the reservoir include effects of
pollution, biological productivity, ecology of fish population, zooplankton, algae, and
bacteria. Atkinson et al. (1999) define the reservoir as slightly eutrophic, based on
chlorophyll a values measured in the main body area. A more recent study (Mabe 2002) considers the possibility of mesotrophic status, based on lower chlorophyll a concentrations measured in the lake. Hunter and Carroll (1982) analyzed the concentration of herbicide residues in fish and sediment samples collected in Lake
Texoma, showing elevated concentrations of PCB in bottom-feeding fishes but less significant amounts of herbicide residues in the same population. Work and Gophen
(1999) studied the environmental factors in Lake Texoma and other U.S. reservoirs affecting the distribution of a cladoceran. They report the gradient of water conductivity of Lake Texoma as an important factor in determining the distribution of Daphnia
lumholtzi in the reservoir. An et al. (2002) identified a seasonal occurrence of E. coli and total coliforms at the marinas of Lake Texoma, and stated that it was possibly due to variations of anthropogenic stress or environmental conditions in those areas. A seasonal variation in chlorophyll a and nitrates is also reported by An and Kampbell (2003) at the
marinas of Lake Texoma.
55 Water Quality Model Parameterization for Lake Texoma
The simulation process involved parameterization, calibration of the selected
model, and validation of WASP for Lake Texoma, including adjustment of initial system
conditions and simulation time step (Ford 1999). The model was set in a one-dimension
condition in order to describe variations in the average concentration of total suspended
solids, chlorides, BOD, and DO in different sectors of Lake Texoma, during a period of
92 days, between March, 1997 to May, 1997. Selection of this period for the simulation
was done because of the availability of water quality data and the variability of the total
inflow to the reservoir and its influence on pool elevation. The WASP model was linked
to the hydrodynamic model DYNHYD, previously calibrated for the reservoir with a time
step adjusted to 30 seconds, as described in the previous chapter, covering the same
period of time. Consequently, time step for WASP was adjusted to 300 seconds which
corresponds to tenfold of the hydrodynamic model time step.
WASP Network Characteristics
In order to model the reservoir using WASP, it was necessary to divide the water
body into segments. In this simulation, every segment was assumed to be a completely
mixed reactor with uniform modeling parameters, such as surface area, volume, cross
sectional area and chemical or physical constants or rates. The WASP network utilized
the same segmentation as the hydrodynamic study. Every junction of the hydrodynamic
model was defined as a water quality segment, except the junctions at the very beginning
or end of the hydrodynamic network because these junctions were used to set boundary
conditions of the water quality model, as explained later. In addition, a benthic layer of
20-cm depth was defined for every water quality segment. This layer acted as a
56 convenient sink for settling BOD and its depth was selected to adequately reflect the thickness of the active layer that interacts with overlying water. Also, it was defined to reflect a reasonable time history in the sediment layer. Below this benthic layer, a sub- benthic layer covering the entire bottom of the reservoir was also defined in order to define a boundary for simulations of the dynamic of solids at the bottom of the reservoir.
The resulting network has 36 water quality segments, 36 benthic segments and one 50-cm- depth sub-benthic layer. Figure 31 shows the final segmentation used, while
Table 9 presents some general characteristics of each water quality segment.
B 20 22 21
23 24 25 2 1 3 26 B 4 27 5 17 18 19 28 6 16 29 7 9 8 15 30 10 14 31 36 11 13 35 12 B 34 33 32 B
# = Segment ID
Figure 31. Segmentation of Lake Texoma used in WASP6
57 Table 9. General characteristics of segments used in WASP for Lake Texoma Description Volume Mean Depth (m) Description Volume Mean Depth (m) (m3) (m3) RR1 10900000 4.17 Benthic1 67849.36 0.02 RR2 57000000 6.38 Benthic2 227809 0.02 RR3 62100000 7.11 Benthic3 172382 0.02 RR4 46500000 9.59 Benthic4 105700.4 0.02 RR5 38300000 9.73 Benthic5 85160.28 0.02 RR6 38600000 10.68 Benthic6 75636.02 0.02 RR7 55100000 11.63 Benthic7 100194.2 0.02 RR8 50100000 12.1 Benthic8 82174.78 0.02 RT9 53600000 11.51 Benthic9 96575.22 0.02 RT10 65200000 10.46 Benthic10 137206.4 0.02 RT11 62500000 21.22 Benthic11 145330.2 0.02 RT12 17500000 6.29 Benthic12 145747.9 0.02 RT13 1.03E+08 14.06 Benthic13 207399.2 0.02 RT14 1.5E+08 14.25 Benthic14 193493.2 0.02 RT15 1.76E+08 15.82 Benthic15 242573.1 0.02 RT16 1.63E+08 16.09 Benthic16 191682.2 0.02 ML17 2.41E+08 16.17 Benthic17 290756.5 0.02 ML18 1.2E+08 15.65 Benthic18 159279.8 0.02 ML19 1.3E+08 17.02 Benthic19 202149.3 0.02 WR20 59100000 9.5 Benthic20 128412.7 0.02 WR21 96500000 11.12 Benthic21 196998.6 0.02 WRT22 1.58E+08 11.27 Benthic22 279131.8 0.02 WRT23 1.33E+08 16.72 Benthic23 173223.5 0.02 WRT24 72000000 16.63 Benthic24 82775.06 0.02 WRT25 83000000 15.76 Benthic25 116072.2 0.02 WRT26 1.03E+08 14.5 Benthic26 138434 0.02 WRT27 1.37E+08 14.32 Benthic27 230113.4 0.02 ML28 3E+08 14.56 Benthic28 368154.9 0.02 ML29 3.11E+08 19.88 Benthic29 372119.1 0.02 ML30 1.37E+08 21.39 Benthic30 121963.5 0.02 ML31 1.58E+08 21.43 Benthic31 142306.1 0.02 BMC32 5610000 3.6 Benthic32 43163.36 0.02 BMC33 12300000 2.91 Benthic33 64745.04 0.02 BMC34 16200000 5.49 Benthic34 64745.04 0.02 BMC35 26800000 6.64 Benthic35 86326.74 0.02 BMC36 59400000 7.77 Benthic36 160437.3 0.02 Sub-benthic 1.42E+08 0.5
A correspondence between the segment identification in the model and the station identification number used in Atkinson et al. (1999) is presented in Table 10.
58 Boundaries
Those segments of the water quality model that import, export, or exchange water with the locations outside of the defined network for the reservoir were considered
“boundaries”. Therefore, four segments were defined as boundaries in this simulation:
o At the very beginning of the Red River arm, and represented by the gage
station located on the Red River near Gainesville, TX.
o At the very beginning of the Big Mineral Creek arm.
o At the very beginning of the Washita River arm, represented by the area of the
reservoir where Station 25 of Atkinson et al. study (1999) is located, and
o At the very end of the reservoir, in the Main Body zone next to the dam.
Table 10. Relationship between water quality segments identification and station identification used by Atkinson et al. (1999) Water Quality Segment Identification Station ID (1) Number of WASP WASP Segment ID Station No. Zone Segment 4 RR4 1 RRZ 8 RR8 3 RRZ 34 BMC34 7 RRTZ 11 RRT11 8 RRTZ 14 RRT14 9 RRTZ 31 ML31 17 MLZ 29 ML29 19 MLZ 27 WRT27 20 WRTZ 22 WRT22 22 WRTZ 20 WR20 24 WRZ Boundary Boundary WR 25 WRZ (1) According to Atkinson et al. (1999) RR = Red River ML = Main Lake BMC = Big Mineral Creek WRT = Washita River Transition RRT = Red River Transition WR = Washita River
59 Model Constants
Several constants, as well as physical and chemical rates were required for WASP calibration to Lake Texoma. A list of the required rates and constants included in the model for the variables simulated in this project is presented in Table 11.
Table 11. Some constants and rates used by WASP for water quality simulation. Constants & Rates Units Value Chlorides
• Dispersion coefficient m2/sec Segments 28 to 22 69 Segments 21 to 20 46 Total Suspended Solids • Settling velocity: m/day Fine-grained sediment (clay) 0.032 to 0.32 Medium-grained sediments (silt) 19 Coarse-grained sediments (sand) 94 • Resuspension velocity: m/day Fine-grained sediment (clay) 0 to 1.27 x 10-07 Medium-grained sediments (silt) 0 to 1.27 x 10-06 Coarse-grained sediments (sand) 0 3 • Solid density: g/cm Fine-grained sediment (clay) 2.7 Medium-grained sediments (silt) 2.0 Coarse-grained sediments (sand) 1.8
Dissolved Oxygen and Biochemical Oxygen Demand
• Dissolved Oxygen Oxygen to Carbon ratio mg O2/mg C 2.57 Reaeration rate day-1 Model calculated Reaeration temperature coefficient 1.04 • Biochemical Oxygen Demand BOD decay day-1 0.1 BOD decay temperature coefficient 1.04 Half saturation constant for oxygen limitation mg O2/L 1
A review of the literature and a calibration process were performed to obtain the corresponding values for dispersion coefficients, settling and resuspension velocities, and density of solids used for this simulation, as well as for constants and kinetic rates associated with DO and BOD (Tufford and McKellar 1999, Benaman et al. 1996, Schultz
2001, Wool et al. 2001). For level one complexity in EUTRO, only two constants had to
60 be set in the model: the reaeration rate and the deoxygenation rate. The reaeration rate helps to determine the rate of gas transfer of oxygen from the overlying atmosphere into the surface. Within EUTRO there were three options to define this rate, including its calculation by the model from water velocity, wind velocity (default set to 0.6 m/sec), water temperature, and air temperature (default set to 15oC). This option was the selected
for this simulation. The deoxygenation rate is set with the sediment oxygen demand and
they account for the rate at which BOD consumes oxygen to stabilize the organic matter
and for the effects of settling and oxidation on BOD in sediments. After a calibration
process, a value of sediment oxygen demand of 0.4 g/m2day was set for the entire
reservoir.
Initially, transport of the components in the water quality simulation was driven
by advection, described by the calibrated hydrodynamic model for Lake Texoma.
However, predictions of chloride concentrations in the Washita River arm did not agree with observed values. As a result, dispersive transport was also included for the transport
of chemical components in this portion of the reservoir. A dispersion function was
specified in the water quality model between segments 28 and 22 (69 m2/sec). Another
dispersion function was defined between segments 21 and 20 (46 m2/sec) after
considering gradients of conductivity reported for that portion of the reservoir (Upton
2004). These values also account for the action of the wind and exceed those reported for
the same parameter in a simulation made by Tufford and McKellar (1999) of Lake
Marion, a large reservoir in the South Carolina (USA) coastal plain, but are lower than
the values used by Benaman et al. (1996) in a water quality simulation in the Houston
ship channel (Texas, USA).
61 In order to set the density of each type of solid considered for this simulation, and
its corresponding settling and resuspension velocity, results of the hydrodynamic study
were considered. Velocity gradients can influence the sedimentation pattern in a water
body. Heavier particles are expected to settle faster than lighter material (Allen et al.
1999). Consequently, settling and resuspension velocities were adjusted taking into account different water velocity gradients in each arm of Lake Texoma. Higher settling velocity values were assigned to larger particles in the Washita River arm, where water
velocities were lower, while the lowest settling velocities were assigned to suspended
particles in the Red River arm (Figure 32 and Figure 33). In addition, the diameter and
density of each type of solid were considered in selecting settling velocity values, as
recommended by Wool et al. (2001).
0.040
0.035
0.030
0.025
0.020
Velocity (m/s) 0.015
0.010
0.005
0.000
2 5 1 R4 R6 10 16 23 3 RR2 R R RR8 T T1 T L18 R T29 L M ML20 W RT2 R M ML33 RR RR RRT14 RR W WRT27 W Water Quality Segment ID
Red River Zone Red River Transition Zone Main Lake Zone Washita River Transition Zone
Figure 32. Estimated water velocities in every segment of the water quality network developed for Lake Texoma using DYNHYD
62 35.00
30.00
25.00 y
20.00
15.00 Settling Velocity (m/da Velocity Settling
10.00
5.00
0.00
1 3 7 1 7 R1 T9 1 2 23 35 R RR3 RR5 RR7 L R T L29 L31 RR RT1 RT15 M ML19 W M M R RRT1 R WR WRT25 WRT2 BMC33 BMC Water Quality Segment ID Red River Zone Red River Transition Zone Main Lake Zone Washita River Transition Zone Big Mineral Creek Zone
Figure 33. Settling velocities used for the water quality network developed for Lake Texoma using WASP Resuspension velocities also were site- and time-dependent. However, in most of
the cases the same value was used for both river arms and no resuspension was defined
for the lacustrine area. The ranges of the selected values are presented in Table 11.
Model Calibration
The calibration process for the water quality model required determining constants and rates for several physical and chemical processes involved in the
sedimentation and resuspension rate and transport processes, as well as those processes
that affect dissolved oxygen concentration. Degradation rates were considered for the
simulation of DO and BOD. Since chloride is a very conservative substance, it was not
necessary to estimate degradation rates for the chloride simulation. All rates and constants utilized were obtained from the literature and adjusted during the calibration.
63 A calculation of the RMSE of the predicted values for chlorides and total
suspended solids concentrations was performed in order to select the constant or rate value that provided the best estimation for these concentrations. Chlorides were evaluated using specific conductance values, which are highly correlated with the concentration of this ion in Lake Texoma (R2 = 95.51%, from Atkinson et al. 1999).
Initial conditions for each cell were set to the measured values of the constituent
concentrations in the lake. A selected period from previous studies made by Atkinson et al. (1999) was used during this process. Additional constituent loadings were set using
the results of the simulation of the local watersheds, based on results provided by Upton
(2004).
The validation process for chlorides, DO, and total suspended solids was
performed using the relative RMSE, a two-sample Kolmogorov-Smirnov test, and a two-
sample t-test for those sets of values with homogeneous variances. One sample corresponded to simulated values of the analyzed variable, and the second sample
corresponded to field measured values reported in previous studies (Atkinson et al. 1999,
Waller et al. 2002).
Model Runs
Data used to adjust the model for the simulated conditions between March 1997
and May 1997 are presented in Appendix E. In this study, the average loading condition
was considered for chlorides and BOD, while two different loading conditions for total
suspended solids were evaluated: the predicted actual condition and the maximum
expected condition, considering estimates made by Upton (2004) using a calibration
based on the Event Mean Concentration method for the local watersheds of Lake
64 Texoma. The highest estimation of total suspended solids reported by Upton (2004) was considered in this simulation as the maximum expected condition. No variations of loadings for chlorides were analyzed in this report since it was assumed that most of the source of this ion is natural. Table 12 summarizes conditions evaluated with WASP in this report.
Table 12. Different loading conditions evaluated for Lake Texoma
Case Boundary Condition Loading Chlorides TSS 1 Measured or simulated Simulated to average Simulated to average to average 2 Measured or simulated Simulated to average Maximum simulated to average
Data Sets
Water quality data were needed to calibrate the model and validate the outputs of the simulation. A review of the available literature shows that hydrological measurements of the Washita and Red rivers have been done continuously for more than 60 years at
Dickson, OK and Gainesville, TX, respectively (USGS 2005). However, water quality in the lake and in the rivers has not been measured as continuously as flow. As a consequence, the availability of water quality data was used to establish the period of time to be simulated.
Based on a review of different studies reported for Lake Texoma, the period of time considered for water quality simulation corresponded to the hydrological year 1996-
1997. During this time, an intensive survey was conducted for the reservoir, and results of this study were included in a database available on the web (Atkinson et al. 1999).
Atkinson et al. (1999) conducted physical, chemical, and biological sampling at
11 routine stations in order to develop a baseline of water quality data for Lake Texoma
65 (August 1996 through September 1997). Also, five random stations were sampled to get
additional information on the spatial distribution of water quality in the reservoir. A list
with routine stations and their locations is presented in Table 13.
Chlorides
Selection of values for chlorides included:
o The analysis of reported data for setting boundary conditions and the initial
concentration for every water quality segment in the model, and
o Identification of concentrations to perform validation of predicted values.
Table 13. Identification of routine stations monitored by Atkinson et al. (1999) in Lake Texoma and their estimated location. Lake Latitude (N) Longitude (W) Station Zone Degrees Minutes Seconds Degrees Minutes Seconds 1 33 54 45 96 54 9 Red River 3 33 52 26 96 50 2 Big Mineral Creek 7 33 47 54 96 48 15 Red River 8 33 49 38 96 46 15 Transition 9 33 50 7 96 42 18 17 33 49 28 96 35 26 Main Lake 19 33 53 22 96 35 10 Washita River 20 33 55 49 96 33 21 Transition 22 34 0 3 96 37 36 Washita 24 34 1 43 96 34 10 River 25 34 2 38 96 34 16
Reported values for chlorides show a general pattern consistent with historic data.
The highest concentration of chlorides is found in the Red River Zone, with a decreasing pattern toward the Main Lake. In the Washita River arm, chloride concentrations decrease from the Main Lake Zone toward the Washita River (Atkinson et al. 1999); the
highest numbers are in the Red River Zone, with decreasing values in the direction of the
66 flow. Instantaneous measurements did not reflect this pattern when the flow of the Red
River, measured at Gainesville gaging station, increases. The spatial variation of the concentration of chlorides in each arm of Lake Texoma during March through May, 1997 is shown in Figure 34 (Red River arm) and Figure 35 (Washita River arm).
Considering this spatial pattern, a linear regression between distance from the boundary of the water quality network and chlorides concentration was performed to estimate the initial concentration for each water quality segment of the model network.
The R2 for the Red River arm was 0.82 (p-value < 0.05), and 0.90 for the Washita River arm regression (p-value < 0.01). Statistical results are included in Appendix F and the resulting concentrations using the regression model are presented in Table 14.
500
450
400
350 Station 17 Main Lake Zone Chlorides (mg/l) 300 Station 4 Red River Zone 250
200 0 10000 20000 30000 40000 50000 60000 Distance (m)
03/09-12/97 04/13-16/97 05/04-07/97 05/18-21/97
Figure 34. Spatial variation of the average concentration of chlorides in the Red River arm, based on measurements made by Atkinson et al. (1999)
67 300
250
200
150 Station 17 Main Lake Zone Chlorides (mg/l) 100
Station 25 50 Washita River Zone
0 0 5000 10000 15000 20000 25000 30000 35000 Distance (m)
03/09-12/97 04/13-16/97 05/04-07/97 05/18-21/97
Figure 35. Spatial variation of the average concentration of chlorides in the Washita River arm, based on measurements made by Atkinson et al. (1999) In order to set boundary conditions for chlorides, a regression between the historic
values of chlorides and conductivity reported by the USGS (2005) for the Gainesville
gaging station, on the Red River during the 1996-1997 hydrologic year, was performed.
This analysis showed that approximately 98% of the variation in chloride concentration
could be explained by its relationship with conductivity during the period considered,
with a p-value of 6.19 x 10-9. Consequently, this relationship was used to transform daily conductivity measurements made at Gainesville gaging station in the Red River, into chloride concentrations. A daily function for chlorides was defined for the Red River boundary using this method. Results of the regression analysis for chlorides and conductivity at Gainesville gage station are included in Appendix G.
68 Table 14. Estimated initial concentrations of chlorides for Washita and Red River arms for WASP calibration WASP Segment ID Corresponding Water Distance from the Estimated Chloride Quality Station ID(1) boundary (m) Concentration (mg/l) RR1 1020.5 360.91 RR2 2286.5 358.21 RR3 4353.5 353.80 RR4 1 6105.5 350.06 RR5 8021.5 345.97 RR6 9514.5 342.79 RR7 11033.5 339.54 RR8 3 12651.5 336.09 RRT9 14600.5 331.93 RRT10 17241.5 326.30 RRT11 8 18929.5 322.69 RRT12 21527.5 317.15 RRT13 24474.5 310.86 RRT14 9 28603.5 302.05 RRT15 31582.5 295.69 RRT16 34438.5 289.60 ML17 36397.5 285.42 ML18 38464.5 281.01 ML19 40625.5 276.40 WR20 24 1148 96.60 WR21 3065 109.18 WRT22 22 5826 127.29 WRT23 8019 141.67 WRT24 11033 161.45 WRT25 12668 172.17 WRT26 14322 183.02 WRT27 20 16920 200.07 ML28 21947 233.04 ML29 19 22951 239.63 ML30 25823 258.47 ML31 17 27560 269.86 (1) as indicated by Atkinson et al. (1999)
In the case of the Washita River boundary, values obtained during water quality sampling at Station 25 in the Atkinson et al. (1999) survey were used to define the boundary condition. A limitation of these measurements was there were only four values for the period considered for the simulation. In order to improve the definition of this boundary condition, a Welch two-sample t-test was performed between chloride concentrations measured by Atkinson et al. (1999) at Station 25 and those measured by
69 the USGS at Dickson gaging station in the Washita River. No significant differences were found between these two stations. Consequently, a regression between chloride concentrations measured at Station 25 in Lake Texoma, and conductivity values measured at the gage station near Dickson, was performed. This analysis showed that about 85% of the variation of chlorides at Station 25 in Lake Texoma could be explained by its relationship with conductivity measurements at Dickson gage station in the
Washita River. Consequently, using this regression, a daily function for chlorides was defined for the boundary in the Washita River, using daily conductivity measurements obtained at the gaging station at Dickson. Statistics are included in Appendix G.
Dissolved Oxygen and Carbonaceous Biochemical Oxygen Demand
The analysis of available data for DO and BOD was performed in order to determine concentrations for these two variables in order to:
o Set the boundary condition at Red River, Washita River, and Big Mineral
Creek
o Set initial conditions for each segment of the water quality network developed
for the simulation, and
o Perform validation of predicted values
Values of monthly measurements of DO corresponding to the simulated period were used to define a monthly function for the boundary condition at Red River, Washita
River, and Big Mineral Creek. Values reported by the USGS for stations 07316000, at the
Red River near Gainesville, and 07331000, at the Washita River near Dickson, were used to set the boundary condition in those two rivers, respectively. Values reported by the
70 Red River Authority for station 15320, located at Big Mineral Creek at FM 901 North of
Sadler, were used to set the boundary condition for Big Mineral Creek (RRA 2003a).
The model requires the definition of the boundary condition for BOD as the ultimate BOD. Consequently, BOD5 values reported by the USGS for the gaging station
USGS 07316000 at the Red River near Gainesville, TX, during the period of time
between 1978 and 1986 were used to define monthly average concentrations for that
river. These values were transformed using an empirical relationship proposed by Weijers
-1 (1999) for waters with low BOD decay rate (KBOD § 0.16 day ). There were no data for
BOD at the gaging station USGS 07331000 at the Washita River near Dickson, OK. So,
in order to set the boundary condition for BOD at the Washita River, reported values by
the USGS of COD at the mentioned gaging stations were analyzed to establish average
monthly concentrations of COD and then transformed into BOD using a factor of 0.7,
considering that between 60 and 80% of biodegradable matter measured as COD in
natural waters is considered biologically degradable (Grady et al. 1999).
The selected level of complexity for the simulation of DO and BOD also
considered the effects of salinity on DO concentrations. Consequently, a boundary
condition for salinity was also defined. Considering the relationship between chlorides
and conductivity, and conductivity and salinity, chloride concentrations defined for the simulation of chlorides in TOXI were also used to set the boundary condition for salinity in EUTRO, but expressed in parts per thousand (ppt).
Only initial conditions for DO and salinity were set for the simulation. The initial concentration of DO for each segment represented the average value of the measurements reported by Atkinson et al. (1999) for each sampling station and the same number was
71 used for all the segments located in the same zone of the corresponding sampling station.
Estimated initial chlorides concentrations developed for TOXI were used as initial
conditions for salinity at each segment of the network, expressed in ppt.
Total suspended solids.
An analysis of the available data on total suspended solids (TSS) was performed
in order to establish the boundary conditions and initial concentrations for every water
quality segment in WASP. Reported values of TSS were less numerous than those of
chlorides. For the Red River, at the gaging station near Gainesville, no determination of
TSS was available for the period of time between March, 1997 to May, 1997. In the case
of the Washita River, only four values were reported by Atkinson et al. (1999) for Station
25 during the same period of time selected for the calibration of the water quality model.
Consequently, it was necessary to use historical data in order to develop a daily function
for the boundary condition in these two points using empirical methods.
The sediment rating curve method (James et al. 2001, Jain 2001) was selected to
develop daily concentrations for the boundaries of the system. For the two gaging stations
(Gainesville in the Red River and Dickson in the Washita River) time series of river discharge and sediment concentration were downloaded from the web server of the
USGS. Data for September, 1930 to August, 1995 were chosen for calibration of Dickson gaging station, and data for October, 1966 to September, 1986 were chosen for the calibration of the Gainesville gaging station. Using the relationship between the recorded
sediment concentration and discharge of the river given by the equation:
LaQ= b , where:
72 L = load; Q = flow; a = Event Mean Concentration (EMC); and b = washout
coefficient. It was assumed that the EMC was equal to the average of the sediment
concentration (C) recorded for each month considered in the calibration process while the
washout coefficient (b), defined as the first-order rate coefficient that relates the rate of
removal of a particle from the soil as a consequence of the runoff to its concentration,
was estimated performing a linear regression of the logC on logQ. Once the coefficient b was obtained, the load of sediments for the period between March, 1997 through May,
1997 was calculated. Finally, a daily concentration of TSS was estimated using the calculated load and the observed mean daily discharge at each gaging station. The washout coefficient for the Red River was estimated to be 0.0267 and for the Washita
River it was 0.9912.
An additional procedure of interpolation of the concentrations of TSS reported by
Atkinson et al. (1999) at Station 25 for the period of time between March and May, 1997
was used to estimate daily concentrations of TSS for the boundary condition at the
Washita River. Equations developed to calculate daily concentrations for each period of
time using this method are included in Table 15. This second procedure provided better
results for predictions in the Washita River Arm, according to the RMS of errors for
selected water quality segments of this portion of the network implemented for WASP6.
Table 15. Equations used to calculate daily concentrations for TSS at the boundary of the Washita River arm for Lake Texoma Simulation Time Period Equation (d = days of simulation) 1 to 32 TSS = -1.3548d + 67.355 32 to 62 TSS = 0.2d + 17.6 62 to 92 TSS = -.0667d + 34.133
73 Another adjustment required for this calibration was the discrimination of three
different types of solids for the estimated TSS concentrations at each boundary.
Historical values were downloaded from the USGS web site, but they were scarce for the
purpose of this simulation. Consequently, an approximate distribution was set based on estimations found in the literature (Schultz 2001) and historic values were only used to compare to the selected ranges for these parameters.
The initial concentration defined for every segment in the water quality network was established using a liner regression of the distance of each arm of the reservoir and the average of observed concentrations, using values reported by Atkinson et al. (1999).
Because of the lack of sediment distribution studies in the reservoir, a calibration process
was used to define the distribution of different types of solids in every water quality
segment of the network. Distribution of solids was also based on estimates found in the
literature (Schultz 2001) and by taking into consideration that sedimentation rates of
suspended particles in Lake Texoma are influenced to some degree by the amount of total
dissolved solids, which in turn are directly related to chloride concentrations (Schroeder
and Toro 1996). As a consequence, different proportions of each fraction of solids
considered in the model were applied to different segments of the water quality network
in the model. The distribution of solids in every segment of the water quality network in
the model is presented in Table 16.
74 Table 16. Percentage distribution of solids per segment used in WASP
%SolidsType 3 %SolidsType 1 %SolidsType 2 WASP Segment (Sand + Fine (Clay) (Silt) Sand) 1 29.22 46.76 24.02 2 29.90 47.85 22.25 3 31.09 49.75 19.15 4 32.20 51.51 16.29 5 33.47 53.55 12.99 6 34.55 55.29 10.16 7 35.71 57.14 7.14 8 37.04 59.26 3.70 9 38.79 61.21 0.00 10 41.43 58.57 0.00 11 43.29 56.71 0.00 12 46.55 53.45 0.00 13 50.92 49.08 0.00 14 58.55 41.45 0.00 15 65.70 34.30 0.00 16 74.29 25.71 0.00 17 81.70 18.30 0.00 18 91.24 8.76 0.00 19 100.00 0.00 0.00 20 5.00 95.00 0.00 21 5.00 95.00 0.00 22 5.00 95.00 0.00 23 5.00 95.00 0.00 24 5.00 95.00 0.00 25 5.00 95.00 0.00 26 5.00 95.00 0.00 27 5.00 95.00 0.00 28 100.00 0.00 0.00 28 100.00 0.00 0.00 29 100.00 0.00 0.00 29 100.00 0.00 0.00 30 100.00 0.00 0.00 30 100.00 0.00 0.00 31 100.00 0.00 0.00 31 100.00 0.00 0.00 Loading
Initially, no loads were considered for the calibration process performed to establish the distribution of different fractions of solids. Once calibrated, estimated loads defined by Upton (2004) for each local watershed around Lake Texoma, as presented in
Figure 36, were included in some segments of the water quality network.
75 Figure 36. Local watershed defined by Upton (2004) for Lake Texoma Estimates of daily runoff and mean event concentrations in fifteen sub-watersheds defined by Upton (2004) were used to define daily loads for some water quality segments of the reservoir. Because these estimates were performed for larger local sub-watersheds than those considered in the water quality segmentation, an equal amount of loading was included in those segments corresponding to a larger sub-watershed. In other cases it was necessary to add estimates made for local sub-watersheds located at the left and right banks of the water quality segment. The algorithms used for calculations of loadings for each water quality segment, using estimations made by Upton (2004), are included in
Table 17. Considering that the largest amount of solids found in the literature for the
76 runoff corresponded to silt, all estimated loads in the water quality model were assumed to be type 2 solids. Daily values of loading for TSS in the model are included in
Appendix E.
Table 17. Equations for estimations of loadings of solids using estimations made by Upton (2004) WASP Segment Algorithm 1 to 8 Load = (NS-RR)/8 + (SS-RR)/8 9 Load = NS-RRT 17 Load = NS-ML 20 Load = (WS-WR2)+ (ES-WR2) 21 Load = (WS-WR)/2 + (ES-WR)/3 22 Load = WS-WRT + (ES-WR)/3 27 Load = ES-WRT 30 Load = SS-ML 34 Load = SS-RRT NS-RR = North side Red River sub-watershed SS-RR = South side Red River sub-watershed NS-RRT = North side Red River transition sub-watershed NS-ML = North side Main Lake sub-watershed SS-ML = South side Main Lake WS-WR = West side Washita River sub-watershed WS-WR2= West side Washita River sub-watershed 2 ES-WR= East side Washita river sub-watershed ES-WR2= East side Washita river sub-watershed 2 SS-RRT = South side Red River transition sub-watershed ES-WRT = East side Washita River transition sub-watershed
The loading estimate approach developed by Upton (2004) for TSS was used to obtain loading estimates for BOD and COD. In other words, COD and BOD Event Mean
Concentrations (EMC) were multiplied by the runoff estimate of each sub-watershed around Lake Texoma.
The Event Mean Concentration (EMC) data for COD and BOD were estimated by compiling literature and water quality data related to the 17 land use categories described by Upton (2004). For each land use, the lowest value of the range reported in the literature was selected. The results are summarized in Table 18. A weighted average of the EMC based on proportions of land use for each sub-watershed was calculated to estimate the corresponding value for each sub-watershed. A summary of the estimated values of EMC for COD and BOD for each sub-watershed is included in Table 19.
Loadings used for each sub-watershed are included in Appendix E.
77 Table 18. EMCs selected for estimating loadings for Lake Texoma. Nitrogen values are from Upton (2004). Event Mean Concentrations (mg/l) Land Use Category BOD COD Min Max Min Max Open water 3.5b 18.0b Low Intensity Residential 7.7b 25.5c 47.4b 54.3c High Intensity Residential 8.8c 9.1b 54.3c 54.3c Commercial 6.6c 23.0c 38.0c 116.0e Barren Rock 3.8d 13.0f 18.0b Quarry 3.5b 28.9d 18.0b Deciduous Forest 3.8b 6.0f 20.0b Evergreen Forest 3.8b 6.0f 20.0b Mxed Forest 3.8b 6.0f 20.0b Shrubland 0.5e 6.0f 20.0b Grass Herbaceous 0.5e 6.0f 20.0b Pasture, Hay 0.5e 6.0f 20.0b Row Crop 4.0e 18.3d 23.9b Small Grains 4.0e 18.3d 23.9b Urban/Recreational Grasses 3.0d 4.1b 21.9b Woody Wetland 3.5b 6.0f 18.0b Herbaceous Wetland 3.5b 6.0f 18.0b
a. Upton 2004 b. Calculated according to Osborne 2002 c. NCTCOG 2001 d. McConnel 1999 e. Baird et al. 1996 f. Naranjo 1997
Water Quality Modeling Results for Lake Texoma
Simulation of chlorides, DO, and total suspended solids of Lake Texoma using
WASP provided important insights of the spatial and temporal variations of these variables in a complex system. The integration of hydrodynamics with other transport processes and sedimentation was evaluated on a daily basis and results of these simulations follow.
78 Table 19. Estimated EMCs for each sub-watershed around Lake Texoma Estimated EMCs by Sub-watershed (mg/l) Sub-watershed Nitrogen BOD COD ESWR 0.99 2.9 21.1 ESWR2 0.97 2.4 21.1 ESWRT 1.38 3.2 22.1 NSML 0.90 2.6 21.5 NSRR 1.01 2.8 21.5 NSRRT 1.15 3.0 21.6 SSML 1.39 3.4 22.9 SSRR 0.98 2.6 20.9 SSRRT 1.95 3.4 23.4 WSWR 0.80 2.7 20.8 WSWR2 0.87 2.4 20.7 WSWRT 1.08 2.3 21.7
Chlorides
As a result of the geology of the watershed, the Red River carries a high
concentration of chlorides (Patrick 1994). Chloride is the most abundant anion in Lake
Texoma (Atkinson et al. 1999), and many of the chemical, physical, and biological
processes are influenced by this ion (Lerfeste et al. 1971, cited by Atkinson et al. 1999).
Considering the conservative properties of chlorides, it was used as a tracer in the water
quality simulation, providing additional information for validation of hydrodynamic
simulation. Initially, only those loads at the boundary of the system provided by the rivers
were included. No additional loads from the local watersheds of Lake Texoma were considered. Results of the simulation showed that the variation of chlorides in the Red
River arm is driven mainly by advection. As shown in , there was a good agreement between the trend followed by measured values of chloride concentrations and simulated values for the period comprised during March, 1997 and May, 1997.
However, in the Washita River arm, predicted values did not agree with observed values during the simulated period when transport was defined mainly by advection
79 forces. It was necessary to define a dispersion function, considering the gradient of
concentration between the Red and Washita rivers. In addition, given the information on
land use, a load was added at segment 22 of the water quality network using estimated
runoff for this sub-watershed (Upton 2004) and the average concentration reported for
chlorides for smaller creeks located in this area. Once these adjustments were applied,
correspondence between predicted and observed values improved. A representation of the
results is presented in Figure 38.
Some temporal variations of chloride concentrations that occur in the reservoir as a consequence of variations in the inflows of the rivers, as well as in the concentration of chlorides in those water bodies, are not accounted for by monthly measurements.
However, the model was able to reproduce those changes, as can be seen in the graphical
representation of predicted values for the Red River Zone (WASP Segments RR4 and
RR8) (Figure 37). These changes do not affect the trend followed by the concentration of
chlorides in the Main Lake Zone of the reservoir, as shown in Figure 38 for WASP
Segments ML29 and ML31. The presence of a large transition zone in the Washita River arm was also sustained by the pattern shown by simulated values for chlorides in WASP
Segments 22 through 27, even though graphical comparison was only performed with the
results obtained for the WASP Segment WRT27, as presented in the same Figure.
The relative RMSE for predictions for the Washita River arm was below 25%,
with a lowest value at WASP segment WRT27 (9%). In the Main Lake zone, the relative
RMSE of the predictions was around 7%. For the Red River arm, the relative RMSE of the predictions was around 10%.
80 WASPWASP Se gmSegment e nt RR4 RR4 WASPWASP Segment Segment RR8 RR8
800 700
700 600 600 ))) /// 500 500 400 400 300 300 Chlorides (mg (mg (mgChlorides Chlorides Chlorides Chlorides (mg/l) (mg/l) (mg/l) Chlorides Chlorides Chlorides ChloridesChloridesChlorides (mg/l (mg/l (mg/l 200 Chlorides (mg/l) (mg/l) (mg/l) Chlorides Chlorides Chlorides 200 100 100 0 0 7 7 7 7 7 7 7 7 7 7 7 99 997 /199 199 199 199 199 199 199 1 /1997 199 2/ 9/ 6/199 7/1 1/ /6/1997 12 19/ 26/ 4/ 4/ 1 23/ 5/ 14/ 11/1997 4/ 4/8 15/1997 29/1997 5 /13/1997 3/ 3/ 3/ 4/ 4/ 4/30/199 5/ 3/ 3/18/1997 3/25/1997 4/ 4/22/ 4/ 5
Predicted (mg/l) Observed(Top) Observed (Bottom) Predicted (mg/l) Observed(Top) Observed (Bottom)
WASP Segment Segment RRT11 RRT11 WASPWASP Se gme Segment nt RRT14 RRT14
600 600 500 500
))) 400 400
300 300 200 ChloridesChloridesChlorides (mg/l (mg/l (mg/l
Chlorides (mg/l) (mg/l) (mg/l) Chlorides Chlorides Chlorides 200 ChloridesChloridesChlorides (mg/l) (mg/l) (mg/l) Chlorides (mg/l) (mg/l) (mg/l) Chlorides Chlorides Chlorides 100 100 0 0 7 7 7 7 7 7 7 7 9 9 9 7 7 7 97 97 9 97 99 9 99 9 99 99 99 99 9 9 1 /1 /199 1 /1 /1 /1 19 /19 1/ 1 5/19 2 9 0 8/1 2/1 /6 3/1 1 4/1/ 4/8/1997 1 2 2 5/6/ 2 4/1/1997 4/8/ /15/1997 /2 /29/1997 5 /1 / / / / / 3/11 3/1 3/25 4 4 4 5 3 3/18/1997 3/25/ 1997 4 4 4 5/13/ 1997 5
Predicted (mg/l) Observed(Top) Observed (Bottom) Predicted (mg/l) Observed(Top) Observed (Bottom)
Figure 37. Temporal variation of chlorides concentration in some stations located in the Red River arm of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997
81 WASPWASP Segment Segment WR20 WR20 WASPWASP Segment Segment WRT2 WRT4277
180 300 160 250 140 ))) 120 200
100 150 80 100 Chlorides (mg/l) (mg/l) (mg/l) Chlorides Chlorides Chlorides
ChloridesChloridesChlorides (mg/l (mg/l (mg/l 60 ChloridesChloridesChlorides (mg/l) (mg/l) (mg/l) Chlorides (mg/l) (mg/l) (mg/l) Chlorides Chlorides Chlorides 40 50
20 0 0 97 7 97 97 7 7 7 7 7 7 7 7 7 7 9 - -9 97 9 9 9 9 9 9 9 9 9 r- y y - r- r- r- r- r- r- - - - ar-97 pr- ar a a a p p p ay ay ay -Mar- M -Ma -A -Ma -Ma M M M M A A A M M M 1 1 0 0 0 - 1- 1- 1- 0- 0- 0- - - - 11- 21-Mar-97 3 1 20-Apr-97 30-Apr-97 1 2 30-May-97 1 1 2 3 1 2 3 10 20 30
Predicted (mg/l) Observed(Top) Observed (Bottom) Predicted (mg/l) Observed(Top) Obs e r ve d ( Bot t om)
WASPWASP SegmSegment ent ML29 ML29 WASPWASP SeSegment gme nt ML31 ML31
300 350
250 300
lll 250 200 200 150 150 ChloridesChloridesChlorides (mg/ (mg/ (mg/ 100 ChloridesChloridesChlorides (mg/l) (mg/l) (mg/l)
Chlorides (mg/l) (mg/l) (mg/l) Chlorides Chlorides Chlorides 100 Chlorides (mg/l) (mg/l) (mg/l) Chlorides Chlorides Chlorides
50 50
0 0 7 7 7 7 7 7 7 7 7 7 9 97 9 9 9 9 9 -9 - -9 r- r- - y r- r- y a a a ar-97 a ar-97 pr ay-97 ay-97 a -M M -M M -A -Apr-97 M M -M -M -Mar-97 -Mar- -Ap -Apr M -May-9 1 1- 1 0 0 1 1 1 0 0 1 2 31- 1 20-Apr-9 30 10- 20- 3 11-Mar-97 2 3 1 20-Apr-97 30 10- 2 30-May-97
Predicted (mg/l) Observed(Top) Obs e r ve d (Bottom ) Predicted (mg/l) Observed(Top) Obs e r ve d ( Bot t om)
Figure 38. Temporal variation of chlorides concentration in some stations located in the Washita River arm and Main Lake Zone of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997
82 The largest errors were found in those portions of the reservoir where a clear difference between concentration at the top and at the bottom of the water column was evident. This can be explained by the configuration used in this simulation. With the
arrangement used during this project, the model does not account for vertical differences
in the concentration of the chemicals. Vertical segmentation of the water column could be
a future improvement to develop in the implementation of WASP as a management tool
for Lake Texoma.
Total Suspended Solids
The water quality model calibration for total suspended solids demonstrated the
capability of WASP to simulate the spatial heterogeneity of this variable in Lake Texoma.
Predictions in the Red River arm depict seasonal and spatial variation of the
concentration of total suspended solids quite well (see Figure 39). This cannot be said for
estimations made for the Washita River arm and the Main Lake zone where the model
does not fully represent the average concentration of suspended solids (Figure 39 and
Figure 40). Information about the distribution of particles in these portions of the water
system and multiple simulations where boundary conditions and initial concentrations for
every segment of the water quality network were adjusted could improve the accuracy of
these predictions.
Statistical tests were used to examine differences between observed TSS values
for every water quality monitoring station in Atkinson et al. (1999) and predicted values
using WASP. A Kolmogorov-Smirnov test was performed for all the stations listed in
Table 10, and a Welch two-sample t-test was performed in those sets of values with
83 homogenous variances (as evaluated with a F test). Results of these analyses are
summarized in Table 20.
The Kolmogorov-Smirnov test shows that there is not enough evidence for
significant differences between predicted and observed values for all the stations.
Variances of predicted and observed values for TSS concentrations were homogeneous for most of the sets of values tested. WASP segments 4, 22, 29 and 31 were found to be non-homogeneous and therefore the t-test was not applied. The lowest p-value was obtained for WASP segment 20. However, when performed, results of the t-test also supported that there is no evidence of statistically significant differences between observed and predicted values for TSS concentrations during the simulated period for most of the segments in the transition zones of the reservoir.
Results for WASP segment 4 and 20, close to the boundaries were affected by estimates that defined boundary conditions on a daily basis. These estimates were made because of the lack of information for the simulated period. In the transition zone, the relationship between changes in TDS and the amount of suspended solids reported for
Lake Texoma (Schroeder and Toro 1996) and its effect in settling velocities can also represent a source of uncertainty in this simulation.
The model was able to simulate the temporal variability of TSS in Lake Texoma,
and many predictions are closer to values reported for the surface layer of the water column than bottom values. This is particularly noticeable in those conditions when there is a large difference between observed values for TSS concentrations at the top and at the bottom of the water column.
84 WASP Segment RR4 WASP Segment RR8
200 35 30 150 25 20 100 15 TSS (mg/l)TSS (mg/l) 50 TSS (mg/l)TSS (mg/l) 10 5 0 0 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97
Predicted Observed(Top) Observed (Bottom) Predicted Observed(Top) Observed (Bottom)
WASP Segment RT11 WASP Segment RT14
20 25
15 20 15 10 10 TSS (mg/l) TSS (mg/l) TSS (mg/l) TSS (mg/l) 5 5 0 0 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97
Predicted Observed(Top) Observed (Bottom) Predicted Observed(Top) Observed (Bottom)
Figure 39. Temporal variation of total suspended solids concentration in some stations located in the Red River arm of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997
85 WASP Segment WR20 WASP Segment WRT22
50 20
40 15 30 10 20 TSS (mg/l)TSS (mg/l) TSS (mg/l)TSS (mg/l) 10 5 0 0 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97
Predicted Observed(Top) Observed (Bottom) Predicted Observed(Top) Observed (Bottom)
WASP Segment WRT27
10 8 6 4 TSS (mg/l)TSS (mg/l) 2 0 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97
Predicted Observed(Top) Observed (Bottom)
Figure 40. Temporal variation of total suspended solids concentration in some stations located in the Washita River arm of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997.
86 The distribution of particles and the settling and resuspension velocities applied for different segments seem to explain these results. However, it can be said that the general temporal and spatial pattern provided by simulations fits closely to the observed values for the two large arms of Lake Texoma.
It was necessary to set different settling and resuspension rates for every type of solid simulated in the model in order to improve predictions. This can be explained by varying temporal and spatial hydrodynamic conditions reported for the reservoir using simulation tools described earlier in this report. In a reservoir with longitudinal velocity gradients it can be expected to find differences in net settling rates among locations in the reservoir. This was considered in the simulation by adjusting different settling and resuspension rates not only by location, but also temporally.
Table 20. Statistical analysis of predicted values for TSS concentrations in Lake Texoma using WASP
WASP K-S test F-test t-test Location Segment p-value F p-value t-value p-value 4 Red River Zone 0.2286 25.33 0.0259 * * 8 Red River Zone 1 0.43 0.5099 0.2003 0.8488 11 Red River Transition Zone 0.7714 0.22 0.2454 0.1124 0.9156 14 Red River Transition Zone 0.7714 1.31 0.8307 1.3559 0.2248 20 Washita River Zone 0.2286 1.26 0.852 2.7192 0.0352 22 Washita River Transition Zone 0.2286 0.0905 0.0791 * * 27 Washita River Transition Zone 0.7714 0.5727 0.6584 0.5205 0.6227 29 Main Lake Zone 0.2286 30.6479 0.0189 * * 31 Main Lake Zone 0.6994 14.3698 0.0552 * *
* = Variances are not homogeneous
In general, settling velocities obtained by calibration increase as velocity declines across the reservoir. The same trend is found on a temporal basis, where reductions in the flow and in the velocity of the water promote higher settling velocities. Consequently, higher temporal variations in the concentration of total suspended solids in the reservoir were more remarkable in the river zones and in the transitions zones rather than in the
87 Main Lake Zone. In the Main Lake Zone, variations in water velocity were less
noticeable and the temporal pattern of total suspended solids in this portion of the
reservoir was more stable (See Figure 41).
WASP Segment ML29
15
10
TSS (mg/l) TSS 5
0 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97
Predicted Observed(Top) Observed (Bottom)
WASP Segment ML31
8 7 6 5 4 3 TSS (mg/l) TSS 2 1 0 1-Mar-97 31-Mar-97 30-Apr-97 30-May-97
Predicted Observed(Top) Observed (Bottom)
Figure 41. Temporal variation of total suspended solids concentration in some stations located in the Main Lake Zone of Lake Texoma in comparison with simulated values for the period between March, 1997 through May, 1997 Resuspension velocities were also dependant on predictions of water velocity in
different portions of the reservoir. In those areas where water velocities were found to be
very low, no resuspension velocities were applied. In Figure 42 temporal variation of
resuspension velocities applied in WASP6 for the calibration of solids type 2 in segments
20 and 27 is presented as an example. This means that there was a threshold water velocity affecting resuspension and settling velocities, and, consequently, the spatial distribution of suspended solids in the reservoir.
88 1.40E-06
1.20E-06
1.00E-06
8.00E-07
6.00E-07
4.00E-07 Resusupension VelocitiesResusupension (m/day)
2.00E-07
0.00E+00
1 5 9 13 17 21 25 29 33 37 41 45 49 53 57 61 65 69 73 77 81 85 89 Simulation Day
Figure 42. Temporal variability of particle resuspension velocity applied to solid type 2 in segment 27 of the water quality network developed for Lake Texoma using WASP, obtained by model calibration
89 Distribution of the size of the particles in different portions of the reservoir is
another important factor considered in the calibration process. Using only one type of solid does not account for the temporal variation of TSS concentrations in different locations in the reservoir.
During events of rapid washout, the distribution of particles changes and the
amount of total suspended solids in the reservoir is heavily affected, especially close to
the boundaries of the system. Rapid washout also affects predictions among different
locations in the reservoir, but by a smaller magnitude in absolute terms.
The RMSE was estimated for every water quality segment, and it was found that this value decreases in the direction of the flow (Figure 43). Close to the boundaries,
there is a clear differentiation between top and bottom concentration for TSS in the water
column. Vertical variation is not accounted for in the arrangement used for the model in
this simulation. However, lack of field data impeded the development of a more detailed
simulation.
The model also does not account for the effect of turbulence generated by wind, especially in shallow areas. This can be improved in future additional simulations of this variable for Lake Texoma. Another important factor to be improved in subsequent
simulations of this variable is the effect of the TDS gradient reported for Lake Texoma
and its relationship with the sedimentation of particles, mainly silt and clay.
DO and BOD
The model was successfully executed with loadings, flows, and parameters explained in previous sections of this report, using the Streeter-Phelps equation, but
representing a dynamic condition, as presented in Figure 44 and Figure 45. This approach
90 allowed a more realistic calibration to observed data. In addition, the sensitivity of the
model to boundary conditions, salinity, and sediment oxygen demand was also
investigated. Further research including nutrients and phytoplankton in the dynamic of oxygen for Lake Texoma could provide additional insights of this variable in the system.
45.00 a) 40.00 35.00 30.00 25.00 20.00 15.00 10.00 5.00 0.00 RMS of Errors for(mg/l) TSS Seg.4 Seg.8 Seg.11 Seg.14 WASP Segment
b) 12.00 10.00 8.00 6.00 4.00 2.00 0.00 RMS of Errors for TSS (mg/l) TSS for Errors of RMS Seg.20 Seg.22 Seg.27 Seg.29 Seg.31 WASP Segments
Figure 43. RMS of errors of TSS concentration for different water quality segments in a) the Red, and b) the Washita rivers arms of Lake Texoma using WASP. The model accounts for temporal and spatial variations of the average
concentration of dissolved oxygen in the reservoir. An overestimation of DO concentrations occurred during the first 15 days of the simulation, especially in some segments at the Red River Arm. It seems that the lack of simulation of other processes associated with oxygen consumption is affecting the results. Another possibility is related with the estimations made in order to set the boundary condition.
91 WASP Segment RR4 WASP Segment RR8
12 12 10 10 8 8 6 6 4 4 DO (mg/l) DO DO (mg/l) DO DO (mg/l) (mg/l) DO DO 2 2 0 0 3/1/19973/1/1997 3/8/19973/8/1997 4/5/19974/5/1997 5/3/19975/3/1997 3/1/19973/1/1997 3/8/19973/8/1997 4/5/19974/5/1997 5/3/19975/3/1997 3/15/1997 3/22/1997 3/29/1997 4/12/1997 4/19/1997 4/26/1997 5/10/1997 5/17/1997 5/24/1997 5/31/1997 3/15/1997 3/22/1997 3/29/1997 4/12/1997 4/19/1997 4/26/1997 5/10/1997 5/17/1997 5/24/1997 5/31/1997 3/15/19973/15/1997 3/22/19973/22/1997 3/29/19973/29/1997 4/12/19974/12/1997 4/19/19974/19/1997 4/26/19974/26/1997 5/10/19975/10/1997 5/17/19975/17/1997 5/24/19975/24/1997 5/31/19975/31/1997
Predicted Observed (Top) Observed (Bottom) Average Predicted Observed (Top) Observed (Bottom) Average
WASP Segment RRT11 WASP Segment RRT14
12 12 10 10 8 8 6 6 4 4 DO (mg/l) (mg/l) DO DO 2 (mg/l) (mg/l) DO DO 2 0 0 3/1/19973/1/1997 3/8/19973/8/1997 4/5/19974/5/1997 5/3/19975/3/1997 3/1/19973/1/1997 3/8/19973/8/1997 4/5/19974/5/1997 5/3/19975/3/1997 3/15/19973/15/1997 3/22/19973/22/1997 3/29/19973/29/1997 4/12/19974/12/1997 4/19/19974/19/1997 4/26/19974/26/1997 5/10/19975/10/1997 5/17/19975/17/1997 5/24/19975/24/1997 5/31/19975/31/1997 3/15/19973/15/1997 3/22/19973/22/1997 3/29/19973/29/1997 4/12/19974/12/1997 4/19/19974/19/1997 4/26/19974/26/1997 5/10/19975/10/1997 5/17/19975/17/1997 5/24/19975/24/1997 5/31/19975/31/1997
Predicted Observed (Top) Observed (Bottom) Average Predicted Observed (Top) Observed (Bottom) Average
Figure 44. Temporal variation of DO concentration in some stations located in the Red River arm of Lake Texoma in comparison with simulated values for the period between March, 1997 and May, 1997
92 WASP Segment WR20 WASP Se gme nt WRT 27
12 14 10 12 8 10 8 6 6 4 4 DO (mg/l) (mg/l) DO DO (mg/l) (mg/l) DO DO 2 2 0 0 3/1/19973/1/1997 3/8/19973/8/1997 4/5/19974/5/1997 5/3/19975/3/1997 3/1/19973/1/1997 3/8/19973/8/1997 4/5/19974/5/1997 5/3/19975/3/1997 3/15/19973/15/1997 3/22/19973/22/1997 3/29/19973/29/1997 4/12/19974/12/1997 4/19/19974/19/1997 4/26/19974/26/1997 5/10/19975/10/1997 5/17/19975/17/1997 5/24/19975/24/1997 5/31/19975/31/1997 3/15/19973/15/1997 3/22/19973/22/1997 3/29/19973/29/1997 4/12/19974/12/1997 4/19/19974/19/1997 4/26/19974/26/1997 5/10/19975/10/1997 5/17/19975/17/1997 5/24/19975/24/1997 5/31/19975/31/1997
Predicted Observed (Top) Observed (Bottom) Average Predicted Observed (Top) Observed (Bottom) Average
WASP Segment ML29 WASP Segment ML31
12 14 10 12 8 10 8 6 6 4 4 DO (mg/l) (mg/l) DO DO (mg/l) (mg/l) DO DO 2 2 0 0 3/1/19973/1/1997 3/8/19973/8/1997 4/5/19974/5/1997 5/3/19975/3/1997 3/1/19973/1/1997 3/8/19973/8/1997 4/5/19974/5/1997 5/3/19975/3/1997 3/15/19973/15/1997 3/22/19973/22/1997 3/29/19973/29/1997 4/12/19974/12/1997 4/19/19974/19/1997 4/26/19974/26/1997 5/10/19975/10/1997 5/17/19975/17/1997 5/24/19975/24/1997 5/31/19975/31/1997 3/15/19973/15/1997 3/22/19973/22/1997 3/29/19973/29/1997 4/12/19974/12/1997 4/19/19974/19/1997 4/26/19974/26/1997 5/10/19975/10/1997 5/17/19975/17/1997 5/24/19975/24/1997 5/31/19975/31/1997
Predicted Observed (Top) Observed (Bottom) Average Predicted Observed (Top) Observed (Bottom) Average
Figure 45. Temporal variation of DO concentration in some stations located in the Washita River arm and Main Lake Zone of Lake Texoma in comparison with simulated values for the period between March, 1997 and May, 1997
93 However, the average of the relative RMSE did not exceed 20% (approximately 1.6 mg
O2/L) of the average concentration of DO in the reservoir. A representation of the absolute
values for the RSME for each water quality segment is presented in Figure 46. The highest
discrepancy between observed and predicted values was detected at WASP segment 11,
corresponding to the water quality station No. 8 at the Red River Zone. This is attributed to
divergences in magnitude, though the pattern followed by the simulation fitted pretty well the
observed pattern of temporal variation. The lowest RMS of errors was obtained for the
simulations of DO concentrations at the Main Lake Zone. Considering these results, it seems that the average condition for oxygen in the reservoir can be described taking into account physical processes only. In order to make a distinction between different strata in the reservoir it should
also be necessary to consider biochemical processes.
The effects of salinity and sediment oxygen demand for DO simulations in Lake Texoma
were compared by activating and deactivating the salinity function using two different conditions
for sediment oxygen demand in the model. It was found that, under conditions of high sediment
oxygen demand, simulated values are closer to the average of observed values when the salinity
function is activated. The magnitude of changes is more remarkable in segments at the Washita
River Arm, where salinity values are lower, and it is expected to be more sensitive to changes in
this variable. Using a lower sediment oxygen demand, like the one applied for simulations of DO
in Lake Texoma (0.4 g/m2-day), the salinity function does not exert any detectable effect in the
results (Figure 47). In all cases, a reduction in sediment oxygen demand remarkably improved
simulation results
94 Red River Arm
1.80 1.60 1.40 1.20 1.00 0.80 0.60 0.40 0.20
RMS of Errors for DO (mg/L) DO for Errors of RMS 0.00 Seg.4 Seg.8 Seg.11 Seg.14
Washita River Arm
1.40 1.20 1.00 0.80 0.60 0.40 0.20
RMS of Errors for DO (mg/L) DO for Errors of RMS 0.00 Seg.20 Seg.22 Seg.27 Seg.29 Seg.31
Figure 46. RMS of errors of DO concentrations for different water quality segments in the Red and Washita River arms of Lake Texoma using WASP
Water Quality Studies for Tocoma Reservoir
Physical, chemical, and biological characteristics of the Caroni River in the portion where
Tocoma Reservoir has been planned are driven by the characteristics of the water drained from
Guri Reservoir. These are, in general sense, very similar to those reported for the lower Caroni
River.
95 a) 25.00
20.00
) S-P+High SOD 15.00 SP+High SOD+Salinity SP+Low SOD 10.00 RMSE (% SP+Low SOD+Salinity
5.00
0.00 Seg.4 Seg.8 Seg.11 Seg.14 Aver.
b) 35.00 30.00
25.00
) S-P+High SOD 20.00 SP+High SOD+Salinity
15.00 SP+Low SOD RMSE (% SP+Low SOD+Salinity 10.00
5.00
0.00 Seg.20 Seg.22 Seg.27 Seg.29 Seg.31 Aver.
Figure 47. RMS of errors of DO concentrations for some water quality segments at the: a) Red River arm, and b) Washita River arm of Lake Texoma under different conditions of salinity and sediment oxygen demand The Caroni River, a blackwater river as a result of the high content of humic substances, is poor in electrolytes (Paolini 1986). Published studies for the lower Caroni River and Guri
Reservoir report a very low conductivity, poor nutrient contents, low concentrations for total suspended solids, and high content of dissolved organic carbon for this portion of the river
(Paolini 1986, Weibezahn 1992, Latorraca 1999, CVG EDELCA 1999).
However, considering the proximity that Tocoma Reservoir will have with respect of
Guri Reservoir, it was necessary to evaluate not only water quality characteristics for this
96 reservoir in a water quality measurement site close to the dam, but also some characteristics of the dam itself, to gain information about specific water quality conditions developed for this portion of the river.
In this sense, five water quality stations were considered during the study, as represented in Figure 48. A list of the geographic location for each water quality station is included in Table
21, and the annual average value for the main limnological characteristics of these places is included in Table 22.
Table 21. Geographic coordinates for water quality sampling stations nearby Tocoma Dam site Water Quality Station ID (1) UTM Coordinates Sector 20 – La Canoa N E Guri Reservoir at Boya 800 805676 500130 Guri Dam discharge channel 1 859673 500153 Guri Dam discharge channel 2 859642 500398 El Merey Creek 862018 503136 Caroni River at Tocoma site 872069 494140 (1) = as CVG EDELCA 2004
Dissolved oxygen content in the water column of the sampling station located at Guri
Reservoir close to the dam, is around 5 mg/l as an annual average. However, during May to
November, the reservoir is stratified and dissolved oxygen concentration close to the bottom is around 3 mg/l. The water drained through the power house No 1 at Guri Dam comes from this stratum close to the bottom of the reservoir, approximately 80 to 100 m deep. On the other hand, the power house No 2 at Guri Dam draws water from the stratum close to the surface, with higher dissolved oxygen content. This condition explains differences found in the dissolved oxygen content for each discharge channel of Guri Dam. Downstream, at Tocoma Dam site the water recovers the amount of dissolved oxygen required for supporting aquatic life, with values closer to the saturation condition (CVG EDELCA 2005b).
97 Tocoma Site C A Cuanaguaro R O N I R Cunaguaro I River V E R
El Merey Qda. El Merey
Discharge Discharge Channel 2 Channel 1
Boya 800
GURI RESERVOIR
Limnological sampling stations
Figure 48. Spatial location of water quality measurement sites at the lower Caroni River considered for this study Other limnological variables show that this portion of the river has soft waters with a
very low buffer capacity, considering total hardness and alkalinity values. Redox potential describes an oxidative environment. Total suspended solid contents are low along the system, as well as turbidity values. Biologically, the system is oligotrophic, with very low phytoplanton biomass, low nutrient contents, and low chlorophyll a values.
Water Quality Model Parameterization for Tocoma Reservoir
Due to the low value of DO found at the inflow of Tocoma, the major question for the model was the impact of impoundment on DO and whether it recovers at the future dam site.
Therefore, in contrast to Lake Texoma, water quality simulation for Tocoma Reservoir only
98 described the DO-BOD processes. Then, parameterization of the water quality model for the future Tocoma Reservoir mainly consisted on adjustments of boundary conditions, loadings, and several physical time functions describing the expected real system for these two processes.
Table 22. Annual average value for the main limnological characteristics of some water quality sampling stations at the lower Caroni River LIMNOLOGICAL Water Quality Sampling Stations VARIABLES Boya 800 Boya 800 Guri Dam Guri Dam Tocoma (surface) (bottom) Channel 1 Channel 2 Temperature (oC) 27.4 27.1 26.4 27.3 27.5 Dissolve Oxygen 5.2 1.9 2.4 4.0 7.4 (mg/l) Transparency (m) 2.35 - - - 2.19 Turbidity (NTU) 2.21 22.2 2.10 2.05 2.08 Real Color (Co-Pt 46.18 38.02 41.30 32.87 43.07 Units) Redox Potential (mv) 198 205 200 197 210 TSS (mg/l) 1.86 45.96 1.85 1.31 1.85 pH 5.89 5.66 5.60 5.89 5.89 Total Alkalinity (mg 3.03 3.03 3.04 2.99 3.02 CaCO3/l) Total Hardness (mg 2.65 2.60 2.53 2.64 2.67 CaCO3/l) Specific Conductivity 8.93 9.35 8.44 8.54 8.21 at 25oC (mS/cm) Ca (mg/l) 0.47 0.41 0.46 0.50 0.55 Mg (mg/l) 0.22 0.19 0.20 0.23 0.24 Na (mg/l) 0.73 0.96 0.79 0.78 0.80 K (mg/l) 0.40 0.41 0.41 0.41 0.42 TSP (mg/l) 0.009 0.020 0.027 0.010 0.012 Nitrites (mg/l) ND ND ND ND ND Nitrates (mg/l) 0.18 0.51 0.48 0.50 0.44 Chlorides (mg/l) 0.62 1.53 0.56 2.64 0.64 Chlorophyll a (mg/l) 0.269 - 0.185 0.294 0.44 TSP = Total soluble phosphorus ND = Non detectable
While evaluating water quality for Tocoma Reservoir, a water quality network for the reservoir was adapted from the hydrodynamic network previously described in this document. In addition, water quality model was set up for describing the reservoir in a one-dimension condition. This means that model estimates represent the average condition for every water
99 quality segment and do not reflect differences in the water column as a consequence of
stratification, if any. A time period of 40 days was established as simulation time, considering
that the system stabilizes in a long period of time for the Cunaguaro River area, based on
estimations of hydraulic and travel time for this portion of the future reservoir. For the rest of the
system, hydrodynamic simulation demonstrated that stability is reached in a short period of time.
The WASP model was linked to hydrodynamic results of the DYNHYD model. Consequently,
time step for water quality simulations was set at 198 seconds, which is approximately 10 times
the hydrodynamic time step.
WASP network for Tocoma Reservoir
Water quality network for the reservoir used the hydrodynamic network prepared for the
simulation of the reservoir with DYNHYD, described in the previous chapter of this document.
Every hydrodynamic junction corresponded with one water quality segment, and inflow and
seaward boundaries defined for DYNHYD became boundaries for the water quality network.
Therefore, the WASP network for Tocoma Reservoir has five boundary segments: one for every discharge channel of Guri Dam, and the other three for El Merey Creek, Cunaguaro River and at the end of the network, next to the future dam. It also has 50 water quality segments. A 20-cm depth benthic segment was defined for describing bottom conditions. This value responds to the need of developing enough “memory” for describing processes at the bottom of the reservoir,
such as oxygen consumption, as well as bottom-water interactions. As presented previously in
this chapter, Figure 30 shows the schematic representation of the water quality network
developed for Tocoma Reservoir. The correspondence between junctions of the hydrodynamic
network developed for DYNHYD and water quality segments for WASP is presented in Figure
100 49. Geometric data for water quality segments used for simulation of the future Tocoma
Reservoir are included in Table 23.
Table 23. Geometric data for water quality segments developed for the WASP network of Tocoma Reservoir WASP Volume Mean Depth WASP Volume Mean Depth Segment ID (m3) (m) Segment ID (m3) (m) 1 3,500,000 14.0 27 57,000,000 34.5 2 8,260,000 18.1 28 64,700,000 23.9 3 13,500,000 22.5 29 57,000,000 19.2 4 8,190,000 28.1 30 4,640,000 11.1 5 6,350,000 17.6 31 22,600,000 12.2 6 32,400,000 30.0 32 19,600,000 9.7 7 26,700,000 27.1 33 11,900,000 8.7 8 18,600,000 34.9 34 5,690,000 9.0 9 60,400,000 31.9 35 24,800,000 12.1 10 41,000,000 37.0 36 22,800,000 22.0 11 39,200,000 25.9 37 81,200,000 37.7 12 57,700,000 24.7 38 86,900,000 38.9 13 71,300,000 24.2 39 59,900,000 40.9 14 18,000,000 19.4 40 67,000,000 43.2 15 51,000,000 24.5 41 85,400,000 44.5 16 69,300,000 26.2 42 74,100,000 42.0 17 12,700,000 12.9 43 57,100,000 38.9 18 27,500,000 14.7 44 73,000,000 37.8 19 43,600,000 15.3 45 48,600,000 45.6 20 95,500,000 34.4 46 80,500,000 32.4 21 14,800,000 12.7 47 50,900,000 38.8 22 40,000,000 34.7 48 42,500,000 35.3 23 92,000,000 37.2 49 24,000,000 21.0 24 25,500,000 33.5 50 21,200,000 15.6 25 29,200,000 25.3 Benthic 16,954,281 0.2 26 61,600,000 20.4 Segment
101 Tocoma Dam site
DYNHYD Junctions Network Guri Dam
Tocoma Dam site
WASP Segments Network
Guri Dam
Figure 49. Correspondence between water quality and hydrodynamic networks developed for Tocoma Reservoir
102 Boundary and Initial Conditions
Water quality measurements were used as the source of data for defining boundary
conditions for the system. Water quality sampling sites at discharge channels of Guri Dam were considered for establishing the boundary condition at this portion of the network. A simple mass balance for these two channels was performed for calculating DO concentration at this boundary.
Therefore,
()[]DO × Q + ()[]DO × Q []DO = 1 1 2 2 , where Q1 + Q2
[DO1] = average DO concentration at discharge channel 1
Q1 = annual average water flow at discharge channel 1
[DO2] = average DO concentration at discharge channel 2
Q2 = annual average water flow at discharge channel 2
Using 2 mg/l for dissolved oxygen concentration at discharge channel 1 and 4 mg/l for
dissolved oxygen concentration at discharge channel 2, the resulting concentration for the boundary condition was 3.3 mg DO/l, considering that the average flow for these channels is
1,500 and 3,000 m3/s, respectively.
Boundary concentrations for El Merey Creek and Cunaguaro River came from the annual average of dissolved oxygen at the sampling stations located in these two water bodies. A value of 6 mg/l was used for both Cunaguaro River and El Merey Creek.
There were no measured data for BOD concentrations. Consequently, considering land use of the immediate watershed, BOD concentrations were set at zero value for the boundary condition at water quality segments 1 and 2, representing Guri Dam discharges.
103 Model Constants
The simplest complexity level for describing DO-BOD dynamic for Tocoma Reservoir was used for this simulation. The Streeter-Phelps equation described this process as indicated below:
o For BOD:
∂C ⎛ C ⎞ v (1 − f ) BOD (T −20) ⎜ DO ⎟ oms DS = −k Dθ D ⎜ ⎟ − C BOD ∂t ⎝ K1/ 2 + C BOD ⎠ D
o For DO:
∂C ⎛ C ⎞ SOD DO (T −20) ⎜ DO ⎟ (T −20) = k R (Cs − C DO ) − k Dθ D ⎜ ⎟C BOD − θ s , where: ∂t ⎝ K BOD + C BOD ⎠ D
CBOD = Biochemical oxygen demand concentration (mg/l)
CDO = Dissolved oxygen concentration (mg/l)
o -1 KD = Deoxygenation rate at 20 C (day )
qD = Deoxygenation temperature coefficient
T = Temperature (oC)
KBOD = Half-saturation constant for limited oxygen (mg O2/l)
voms = Organic matter settling velocity (m/day)
fDS = Dissolved BOD fraction
D = Surface water column depth (m)
-1 kR = Reaeration rate (day )
Cs = Dissolved oxygen saturation (mg/l)
SOD = Sediment oxygen demand at 20oC (g/m2day)
qs = Temperature coefficient
104 This simulation only considered reaeration processes and sediment oxygen demand for
describing variations on dissolved oxygen concentrations in the reservoir. Using this approach,
oxygen gets into the system by reaeration, driven by hydraulic processes and wind action forces;
and it is consumed by BOD decomposition and sediment oxygen demand. Based on mean wind
velocity registered at Las Babas weather station (CVG EDELCA 2005a), close to Guri Dam,
wind velocity was set at 3 m/s for the entire system.
In addition, WASP was allowed to make spatial and temporal adjustment of the
reaeration rate, considering differences between water quality segment depths in the network.
Then, two time functions were defined, as presented in Table 24.
Table 24. Time functions used for DO and BOD simulation of Tocoma Reservoir Time Function Value Water velocity 1 m/s Water temperature 27oC
Oxygen demand considered the oxidation process. Equations that govern oxidation in the
DO-BOD balance for WASP include phytoplankton anaerobic decomposition and anaerobic decomposition of benthic organic carbon. The latter fraction was the only one considered for this simulation of Tocoma Reservoir. As a result of the calibration process, constants and rates selected for Tocoma Reservoir are presented in Table 25, considering values recommended by
Ambrose et al. (1993b) and Acevedo (2002).
Table 25. Constants and rates used by WASP for water quality simulation of Tocoma Reservoir Constants and Rates Units Value Sediment oxygen demand g/m2day 1.4 • Temperature correction factor 1.0477 BOD degradation rate day-1 0.75 • Temperature correction factor 1.04
105 Loadings
Because the model describes water quality of the system once the reservoir is created, it
was necessary to consider two conditions to define loading conditions:
o Loadings provided by land use. In this case, values recommended by Acevedo (2002)
were included, considering that there have been few or no changes in land use of the
watershed since then, and
o Loadings resulting by flooding of vegetation as a consequence of the elevation of
water level for creating the new reservoir.
For the first condition, values considered for BOD loadings of each watershed around the future reservoir are presented in Table 26. For those loads considered as a consequence of the flooding process, a time function was developed to adjust the simulation process to the flooding plan of the reservoir. According to this, the reservoir will be filled up in two main phases: from elevation 94 to elevation 98 m.a.s.l., in a period of time of time of 16 days; and from this level to elevation 127 m.a.s.l., in 45 days since the 35th day of the filling period. A representation of the
time function developed for two of these segments is presented in Figure 50.
Table 26. BOD loading conditions for adjacent watershed of the Tocoma Reservoir Watershed WASP Segment ID Time period of loading Value receiving the load(1) application (Kg/day)(2) El Merey Creek 5 700 Maria Luisa Creek 26 200 Cunaguaro River 17 03/01/2012 to 07/01/2012 1,800 Nameless 30 140 Guarapo Creek 9 150 (1) = Numbers correspond to water quality network under reservoir condition
(2) = Values according to Acevedo and Sherman (2002)
Loading estimations included during the filling up process for the new reservoir required
consideration of the area and types of vegetation that will be flooded during this process.
106 Specific degradation factors for each type of vegetation were also considered for establishing the amount of BOD that will be incorporated into the system as a consequence of the flooding process.
BOD loadings resulting from this estimation were spatially distributed in the water quality network considering the types of vegetation present in each segment, based on the vegetation map developed for that area (CVG EDELCA 2005b). Consequently, the time function for BOD loadings was practically developed for most of the segments of the water quality network while simulating the filling process of the new reservoir, as presented in Figure 51.
90.0
80.0
70.0
60.0
50.0
40.0
BOD Loading (Kg/day) Loading BOD 30.0
20.0
10.0
0.0
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /2 /1 /8 5 2 9 /5 2 9 6 /3 0 7 4 1 /7 4 1 8 3 3 /1 /2 /2 4 /1 /1 /2 5 /1 /1 /2 /3 6 /1 /2 /2 3 3 3 4 4 4 5 5 5 5 6 6 6
Seg. 5 - El Merey Creek Seg. 17 - Cunaguaro River
Figure 50. Temporal variation of BOD loadings for two segments of the water quality network, while simulating the filling up process for Tocoma Reservoir
107 Model Runs
Once calibrated, the EUTRO sub-model of WASP was used for simulating expected DO
and BOD concentrations of the new reservoir during the filling up and operation time.
Simulation period was of 120 days, corresponding to the period between March and July/2012, considering that the new reservoir will be created during those dates.
Table 27. Phytomass and degradation factors used for BOD loading estimates during the creation of Tocoma Reservoir Type of Vegetation Phytomass Phytomass Degradable Floodable Content Leaf Fraction 2 Area (m ) (Ton/Ha) Fraction (%) (%) Gallery Forest 3,700,000 80 River Forest 9,200,000 95 10 40 Medium Forest 2,260,000 80 Savanna with woody elements 20,070,000 22 Savanna w/woody elements and 630,000 38 bushes Savanna w/woody elements and 80,000 48.6 medium forest Low savanna 2,670,000 1.7 Low savanna with bushes 5,340,000 22 Low savanna with thicket 2,510,000 1.7 Low savanna with floodable 80,000 22 100 100 grassland Low savanna and savanna with 2,180,000 22 woody elements Floodable grassland 60,000 1.7 Bushes 1,270,000 48.6 Bushes with low forest 70,000 48.6 Thicket 430,000 38 Thicket with low forest 40,000 48.6 Thicket with low savanna 20,000 22
Water Quality Modeling Results for Tocoma Reservoir
The simulation of potential future dissolved oxygen and biochemical oxygen demand for
Tocoma Reservoir using WASP provided a better comprehension and understanding of the role
of hydrodynamic and aerodynamic components of the reaeration process for this water body. In
108 addition, the role of BOD loadings on water quality conditions for the new reservoir, resulting of vegetation decomposition, was also stated.
Dissolved Oxygen and Biochemical Oxygen Demand
Model runs included testing of the calibrated model for river condition of this portion of
the Caroni River, as proposed by Acevedo (2002), and two different scenarios for the reservoir
condition. The first simulation under reservoir condition considered absence of autochthonous
BOD loadings, with oxygen demand provided only by external loadings resulting from land use
around in the immediate watershed and oxygen demand from sediments; the second scenario for
simulating reservoir condition included the effect of vegetation decomposition assuming no
removal of existing vegetation in the area before filling up the reservoir.
WASP was previously calibrated for simulating river conditions of the Caroni River for
the portion that entails Guri Dam discharge channels and Tocoma Dam site. By that time, there
were no DO measurements for this entire portion of the river and boundary conditions for
calibration were based on estimations of expected DO concentrations at Guri Dam discharge
channels assuming DO concentrations similar to those measured in the surface layer of Guri
Reservoir (Acevedo 2002). Under this assumption, DO concentrations at the boundary of the
system were set at 6 mg/l and it did not account for the effect of turbinated water coming from deeper layers of Guri Reservoir, with low oxygen content, and its effect on the resulting DO concentration at Guri Dam discharge channels.
Consequently, for understanding the DO dynamic at this portion of the river and under future reservoir conditions, water quality measurements, including DO monitoring, were performed at Guri Dam discharge channels and at Tocoma Dam site between January and
August, 2005. These values were used for estimating model boundary conditions using the mass
109 balance approach described earlier in this chapter and for validation of this extended DO simulation for this portion of the Caroni River. DO data collected in a monthly interval during this period of time are presented in Table 28.
Islands Segments with BOD loading
Figure 51. Spatial distribution of water quality segments with autochthonous BOD loadings Table 28. DO values for Guri Dam discharge channels and Tocoma Dam site measured during January and August, 2005 Dissolved Oxygen (mg/l) Date Guri Dam discharge Guri Dam discharge Tocoma Dam site channel 1 channel 2 01//2005 1.4 4.0 7.2 02//2005 1.5 4.1 6.5 03//2005 1.8 4.4 6.8 04//2005 NM NM 7.7 05//2005 3.8 4.7 7.2 06//2005 3.1 4.9 7.8 07//2005 2.7 4.1 7.4 08//2005 (1) 7.4 4.6 NM NM = Not measured (1) = By the time of the sampling, the spillway was open.
110 Using DO data measured from January to July, 2005, DO concentration at the boundary
resulted in approximately 3.3 mg/l. The model accounted for spatial variations of average
concentration of DO in the river. Simulation of the system under river condition reflected a low increment of DO in the water while it flows through the Necuima Canyon and an important gain
of oxygen by the water once it reaches the widest portion of the river, as presented in Figure 52.
DO (mg/L) Min=3.3 Max=6.51
Tocoma Dam site
Guri Dam
Figure 52. Spatial distribution of DO concentrations for Caroni River between Guri Dam discharge channels and Tocoma Dam site
111 When analyzing the reaeration process that explains this pattern of DO spatial
distribution, it was found that DO gain at Necuima Canyon is low, even though reaeration coefficients at that portion of the river are high as a result of high water velocities (3 to 5 m/s approximately) and low depth (4 m average depth). The hydraulic travel time of the water at this
portion of the river explains this behavior.
As presented in Table 29, travel time at those segments at Necuima Cayon is short.
Consequently, a rapid flush of the water occurs and there is not enough time for diffusion of atmospheric oxygen through the deeper layers of the water column. Then, low oxygenated waters coming from Guri Reservoir, as indicated by the boundary condition, do not have enough time for reaeration and the resulting average of DO concentration for segments at this portion of the network remains low.
Table 29. Estimated values used for DO balance at water quality segments at Necuima Canyon under river condition WASP Length Flow Water Hydraulic Depth Reareation DO DO Segment (m) (m3/s) Velocity Travel (m) Coefficient Upstream Downstream ID (1) (m/s) Time at 27oC (mg/l) (mg/l) (day) (day-1) 1 1570 1333 2.96 0.0061 4 1.51 2.0 2.06 2 1570 2666 2.33 0.0078 3 1.79 4.0 4.06 3 1500 4000 5.23 0.0033 3 3.86 3.39 3.45 4 1200 4000 4.89 0.0028 4 2.80 3.45 3.48 5 1500 4000 5.17 0.0034 4 3.64 3.48 3.54 Total 0.0234 (1) as developed by Acevedo (2002) for Tocoma under river condition
Those segments located in the wider portion of the network also had high reaeration
coefficients, and because these coefficients are the most important hydrodynamic components
for determining the reaeration process simulated by the model (hk = 3 to 9 day-1, wk = 0.15 to
0.30 day-1), the DO concentration increased. However, their surface areas were high enough and
the hydraulic travel time got smaller in comparison with values obtained for segments at the
Necuima Canyon portion of the network. Therefore, oxygen diffusion to deeper layers of the
112 water column was more efficient and the resulting average DO concentration was higher, and
even close to saturation condition. The RMSE of errors for DO estimations for water quality
segment 28 was 11.7%, when comparing calibrated model results for river condition with DO
measured values at Tocoma Dam site in the Caroni River. WASP tended to under-predict DO
concentrations probably because those effects associated with turbulent flow were not included
in the simulation.
Simulations for Tocoma under reservoir condition showed not only spatial but also
temporal differences for expected DO concentrations in the reservoir. For both cases, with and
without BOD loadings, the model accounted for spatial variations of DO concentrations, as presented in Figure 53. Expected DO concentrations ranged between 3.3 and 6.5 mg/l when simulating the new reservoir without BOD loadings. The highest values corresponded to those segments closer to boundary segments with high DO initial concentrations. Simulated DO concentrations ranked between 0 and 4.5 mg/l when BOD loadings resulting from decomposition of flooded vegetation were included, with the exception for those segments closer to boundary segments with high initial DO concentrations.
The model reflected that the lowest DO concentrations for the new reservoir will be distributed along the central axis of this water body when simulating the system without loadings.
This pattern responds to hydrodynamic components of the simulation. The largest portion of practically deoxygenated waters coming from Guri Reservoir will flow through this area of the new reservoir. However, littoral zones will be shallower and oxygen diffusion will be easier at these portions of the reservoir.
This spatial pattern substantially changes when BOD loadings are incorporated for the simulation. An important DO depletion occurs along the entire future reservoir, with the
113 exception of those arms where boundary conditions allow the incorporation of oxygenated water, reducing the effect of the load in terms of dissolved oxygen concentrations for these areas.
Deeper areas closer to the dam receive not only deoxygenated water from Guri Reservoir, but also loading effects will be importantly affected considering that oxygen diffusion toward deeper layers of these water quality segments will be difficult.
a) DO (mg/l) with BOD loadings b) DO (mg/l) without BOD loadings
TocomaDam Tocoma Dam
Guri Dam Guri Dam
Figure 53. Expected spatial distribution of DO for Tocoma Reservoir when simulating the reservoir a) with, and b) without BOD loadings
114 CHAPTER 4
ROLE OF BOUNDARY CONDITIONS AND LOADINGS IN WATER QUALITY
SIMULATIONS
A sensitivity analysis was performed for testing the role of upstream boundary conditions and loadings for simulating some water quality variables in both of these systems, Lake Texoma and Tocoma Reservoir. A sensitivity analysis of a model allows evaluating the response of the model to variations in input values. If a change in an input value causes significant changes in output values (i.e. constituent concentrations), the model is said to be sensitive to that input. The sensitivity analysis was conducted by changing the magnitude of boundary conditions and loadings within reasonable limits while keeping all other parameters unchanged. All combinations of changes of boundary conditions and loadings used in Lake Texoma and Tocoma
Reservoir models were not evaluated because of the large number of simulations required. Only single changes of the boundary condition for chlorides, TSS, and DO; and simultaneous changes for chlorides and TSS loadings were tested for the Lake Texoma model. DO and BOD for
Tocoma Reservoir were considered for each of the tests described below. The model output for each of these simulations was compared to the original model output presented in the previous chapter, and the percent change was calculated. This analysis provided important insights of the role of these two elements of water quality simulations for determining the behavior of water quality conditions in both systems, Lake Texoma and Tocoma Reservoir.
Sensitivity Analysis to Upstream Boundary Conditions
Simulation of water quality for Lake Texoma and Tocoma Reservoir using WASP involved the definition of inflow boundaries. Chemical and physical boundary conditions
115 included in the Lake Texoma model were chlorides, total suspended solids, dissolved oxygen,
and carbonaceous biochemical oxygen demand. For the Tocoma Reservoir, dissolved oxygen
and carbonaceous biochemical oxygen demand were considered for defining chemical inflow
boundaries of the system.
Consequently, the sensitivity analyses of boundary conditions for Lake Texoma
comprised simultaneous changes of all boundary values for each simulated variable: chlorides,
TSS, and DO. This is to say, when evaluating sensitivity to changes of boundary conditions for
chlorides, daily functions defined for all the boundaries of the system were simultaneously
changed, and similarly was done for the rest of the variables.
For the Tocoma Reservoir, isolated changes of each boundary defined for the system
were tested for DO and BOD considering not only the differences found in the hydrodynamics of
the central body of the future reservoir with respect to the water movement at Cunaguaro River
arm, but also the role that DO and BOD exert on the water quality of this reservoir.
Sensitivity to Boundary Conditions for Lake Texoma
Three upstream boundaries were defined for Lake Texoma: Red River, Washita River,
and Big Mineral Creek boundaries. The sensitivity analysis for upstream boundaries of the Lake
Texoma model was performed by changing more or less certain percent of values used as boundary conditions for chlorides, DO, and TSS. The percent difference between the original predicted values and resulting values after changing the boundary condition was calculated, and the average of these values during the simulated period of time was used to establish the sensitivity of the model to the applied change.
Boundary conditions for chlorides in Lake Texoma model involved the development of a statistical regression to transform daily conductivity measurements made at Gainesville gaging
116 station in the Red River, into daily chloride concentrations, using the entire historical data available. For the Washita River boundary a numerical regression between chloride concentrations measured at Station 25 of Atkinson et al. study (1999) and conductivity values
measured at the gage station near Dickson, OK, was performed to develop a daily function for
chlorides. For Big Mineral Creek, an exponential function was found between chlorides
concentrations and flow, using historical data from USGS station 07316200, at a location near
Saddler, TX. Then, flow estimates using HEC-HMS simulations were converted into chlorides
concentrations using the function previously mentioned. Changes of +/- 10% were applied to
these three boundary conditions to test sensitivity of model results to boundary condition
changes.
Three different types of solids were considered for defining the boundary condition for
simulating TSS for Lake Texoma. After defining the function for developing daily
concentrations for TSS concentrations at the boundary, an additional function was used for
defining the fraction that each type of solids represented at the boundary. The sensitivity test was
performed by simultaneously changing by +/- 10% the concentration of each type of solids used
in the original simulation.
For testing sensitivity to changes in DO boundary conditions, those values estimated for
the original simulation were changed in more or less 30%, but taking care of not exceeding the
saturation concentration based on water temperature.
All values used for these sensitivity tests as boundaries are included in Appendix H, and
the results of the tests in terms of percent of changes are shown in Table 30 and Table 31.
117 The simulation results reflected that the system is very sensitive to changes of the
boundary condition for chlorides, but returned less intensity to changes in the boundary condition
of solids type 1 (clay), dissolved oxygen, and solids type 2 (silt) and 3 (sand), respectively.
However a spatial pattern of sensitivity was identified in this system. As reflected in
Figure 54, the Red River arm is the most sensitive to changes to the boundary condition,
especially for chlorides. This response becomes less significant as the water flows through the
system. Transition zones for the Red and Washita River arms reflected a lower sensitivity to this
change, and in the Main Lake zone, the effect was almost negligible. It can be considered that the
travel time could have affected these results, attenuating the effect of boundary condition
changes in the middle and terminal segments of the reservoir. Another factor to be considered is
the diffusion process, which was included in the simulation through a function defined for the
Washita River arm. The magnitude ascribed to this function could have overlapped the effect of
boundary condition changes at this portion of the reservoir. Consequently, the value of the
diffusion rate is another important factor that deserves to be tested for sensitive of the system to
changes in the inflow of chlorides.
When evaluating sensitivity for the Lake Texoma model to changes in the boundary
condition of DO, other factors, such as: reaeration processes and biochemical oxygen demand might have a definite effect on sensitivity to these changes. As shown in Table 30, changes in the boundary condition of dissolved oxygen did not represent a significant change in most of the areas of the reservoir. Only a relatively significant effect was detected at the proximities of the boundaries, and this effect was more remarkable for the Red River arm than for the Washita
River arm. At the Big Mineral Creek area the effect was negligible. The model reflected that vertical diffusion of oxygen is not affected at the shallower areas of the reservoir, and the
118 resulting average concentration for those water quality segments located at the arms of the
reservoir is relatively high, no matter the condition of the boundary. In future research, further
analysis of this situation can be performed by setting up the model in two-dimensions.
Figure 54. Spatial distribution of percent change of chloride concentrations in Lake Texoma when an increment of +10% of boundary condition for chlorides was applied In terms of TSS, the model showed an important response to changes in the boundary
condition for solids type 1, but with less intensity to changes in boundary condition for solids type 2 and 3, as can be understood from values presented in Table 31. Factors such as settling
velocities and resuspension velocities associated to each type of solids could be overwhelming
the observed results. However, when changes of the amount of each type of solid at the boundary
condition were applied, a change in the composition of the total amount of solids incorporating
into the system also occurred. Consequently, a closer look at the effect that these changes in the
119 proportions of each type of solid at the boundaries of the system exert on TSS results is
recommended.
Table 30. Results of the sensitivity test for chlorides and DO showing percent change in daily, volume-averaged concentrations for some WASP segments of the Lake Texoma model compare to the base case in March-May, 1997 Chlorides Dissolved Oxygen WASP Segment Input (Percent Output (Percent Output (Percent Output (Percent change) Change) Change) Change) RR4 + 10 10 + 30 8 - 10 -10 - 30 -8 RR8 + 10 9 + 30 5 - 10 -8 - 30 -5 RRT11 + 10 7 + 30 2 - 10 -6 - 30 -2 RRT14 + 10 6 + 30 1 - 10 -5 - 30 -1 WR20 + 10 5 + 30 9 - 10 -4 - 30 -9 WR22 + 10 3 + 30 3 - 10 -2 - 30 -3 WRT27 + 10 2 + 30 1 - 10 -1 - 30 -1 ML29 + 10 2 + 30 0.4 - 10 -1 - 30 -0.4 ML31 + 10 1 + 30 0.2 - 10 0 - 30 -0.2 BMC34 + 10 3 + 30 0.4 - 10 -3 - 30 -0.4
Sensitivity to Boundary Conditions for the Tocoma Reservoir
For the Tocoma Reservoir, three boundaries were set for DO-BOD simulations of the reservoir based on water quality measurements made at those sites. However, the sensitivity analysis of upstream boundary conditions was performed by changing ≤50% the average concentration of DO measured at Guri Dam discharge channels. The value originally simulated was 3.3 mg/l. The model output of each of these simulations was compared to the original model output and the relative difference was calculated. The spatial distribution of these differences
120 was graphically displayed, as shown in Figure 55, to establish whether there was a pattern for the
calculated differences.
Table 31. Results of the sensitivity test for solids 1, 2, and 3 showing percent change in daily, volume-averaged concentrations of TSS for some WASP segments of the Lake Texoma model compare to the base case in March- May, 1997 WASP Solids 1 Solids 2 Solids 3 Segment Input Output TSS Input Output TSS Input Output TSS (Percent concentration (Percent concentration (Percent concentration change) (Percent change) (Percent change) (Percent change) change) change) RR4 + 10 5.1 + 10 4.5 + 10 0.4 - 10 -4.9 - 10 -4.4 - 10 -0.3 RR8 + 10 5.4 + 10 3.6 + 10 0.4 - 10 -4.7 - 10 -3.0 - 10 0.2 RRT11 + 10 5.1 + 10 2.2 + 10 0.2 - 10 -4.7 - 10 -1.8 - 10 0.2 RRT14 + 10 4.6 + 10 1.0 + 10 4.6 - 10 -4.5 - 10 -0.9 - 10 0.1 WR20 + 10 2.0 + 10 4.9 + 10 1.0 - 10 -2.0 - 10 -4.9 - 10 -1.0 WRT22 + 10 2.6 + 10 3.1 + 10 0.2 - 10 -2.6 - 10 -3.1 - 10 -0.2 WRT27 + 10 2.6 + 10 1.4 + 10 0.01 - 10 -2.6 - 10 -1.4 - 10 -0.01 WRT29 + 10 2.0 + 10 1.1 + 10 0 - 10 -2.0 - 10 -1.1 - 10 0 ML31 + 10 1.5 + 10 1.1 + 10 0 - 10 -1.5 - 10 -1.1 - 10 0 BMC34 + 10 9.4 + 10 0.02 + 10 0 - 10 -9.4 - 10 -0.02 - 10 0
The effect of changing the boundary for DO was remarkable in the central portion of the
future reservoir. For those segments located in the path of the main flow of turbinated waters
coming from Guri Dam and flowing through the future Tocoma reservoir, DO concentrations
practically responded in the same magnitude of the change. The average DO concentration for
the entire reservoir increased in approximately 25% by increasing the boundary condition by
50%, and decreased in the same magnitude when the value at the boundary was reduced. Only
121 for those segments located at the confluence the tributaries to the main water body, the effect of
changes in the boundary condition was negligible. With this simulation, the role of vertical diffusion of oxygen to the deeper layers of the reservoir cannot be appreciated because the model estimated an average of the concentration for each water segment and as a consequence of the important effect that the boundary condition exerts over the entire system.
Figure 55. Spatial distribution of percent change obtained during the sensitivity test to changes in the boundary condition of DO for Tocoma Reservoir
Sensitivity Analysis for Loading Conditions
Loading conditions for TSS and chlorides were tested in the Lake Texoma model. For the
Tocoma Reservoir model, estimated allochthonous and autochthonous loadings of BOD were
tested. For both models, Lake Texoma and Tocoma Reservoir, loadings were varied by ≤50%
and the concentrations results were compared to original output of the models to quantify the
numerical difference after incorporating these changes. All values used for setting new loading
conditions are included in Appendix I.
122 Sensitivity to Loadings for Lake Texoma
All defined loadings for TSS were simultaneously changed by ≤50% for testing the sensitivity of the model to these changes. According to the simulation process, explained in the previous chapter of this dissertation, loadings for TSS were estimated based on simulation results of all sub-watersheds around Lake Texoma. TSS loadings were defined for 16 segments around the lake. Two runs were performed for testing the effect of increasing or decreasing in 50% inputs of solids type 2 (silt) from the surrounding watersheds on the resulting concentrations. For chlorides, the same procedure was applied to the unique load defined for the system, located at the Washita River Arm. An additional run was performed in this case, by eliminating the load, considering the magnitude of the results obtained during the sensitivity test. All results are presented in Table 32 and Table 33.
These results showed that the concentration of TSS estimated by Lake Texoma model is practically unaffected by important changes in loadings for solids type 2. Even though the magnitude of change for values of daily loading was remarkable and it represented significant amounts of solids entering into the system, the resulting TSS concentrations estimated for most of the water quality segments did not change. This might be the result of the application of high settling velocities in the model, which were related to low water velocities through the water body.
The evaluation of changes in the amount of solids accumulated in benthic layers defined in the calibration of WASP for Lake Texoma might provide important insights of the overall effect that these changes in TSS loadings produce in this reservoir.
Sensitivity to changes in loadings for chlorides was tested only at the Washita River Arm.
In this case, the results showed the effect to a lower extent. Percent changes of simulated
123 chloride concentrations were close to zero for almost all segments tested. However, an additional test was performed by totally removing the load at segment WRT22 and there was a noticeable effect in temporal trend of chlorides concentrations described by the model for the segment receiving the load (WRT22), as presented in Figure 56, and those downstream of the loading.
Table 32 . Output percent change for TSS concentrations at several segments of the water quality network of Lake Texoma during the sensitivity test to Solid 2 loading conditions WASP Input (Percent Output (Percent Segment change of boundary) change) +50 0.1 RR4 -50 - 0.1 +50 0.4 RR8 -50 - 0.4 +50 0.2 RRT11 -50 - 0.2 +50 0.1 RRT14 -50 - 0.1 +50 0.4 WR20 -50 - 0.4 +50 0.5 WRT22 -50 - 0.5 +50 0.4 WRT27 -50 - 0.4 +50 0.5 WRT29 -50 - 0.5 +50 1.7 ML31 -50 - 1.7 +50 0.03 BMC34 -50 - 0.03
These results showed that the permanent incorporation of a dissolved conservative substance, like chlorides, in this portion of the system, as described by the model for Lake
Texoma, is responsible of the temporal trend of the concentration of that substance at the
Washita River arm, and a lower influence is produced by the magnitude of the load. Diffusion processes defined in the model could be influencing these results and a deeper analysis of this process at this portion of the reservoir is recommended.
124 Table 33. Output percent change for chloride concentrations at several segments of the water quality network for Lake Texoma during the sensitivity test of the model to chloride loading conditions WASP Input (Percent Output (Percent Segment change of boundary) change) +50 5.6 x 10-4 WR20 -50 - 6.6 x 10-4 +50 6.8 x 10-4 WR21 -50 - 7.6 x 10-4 +50 9.8 x 10-4 WRT22 -50 - 9.3 x 10-4 +50 9.1 x 10-4 WRT23 -50 - 8.6 x 10-4 +50 7.4 x 10-4 WRT24 -50 - 6.3 x 10-4 +50 6.2 x 10-4 WRT25 -50 - 5.6 x 10-4 +50 4.8 x 10-4 WRT26 -50 - 4.6 x 10-4 +50 3.7 x 10-4 WRT27 -50 - 3.5 x 10-4 +50 2.3 x 10-4 WRT28 -50 - 2 x 10-4 +50 2.2 x 10-4 ML29 -50 - 1.4 x 10-4 +50 1.7 x 10-4 ML30 -50 - 1.2 x 10-4 +50 1.6 x 10-4 ML31 -50 - 1.2 x 10-4
Sensitivity to Loadings for Tocoma Reservoir
Testing sensitivity to loadings for Tocoma Reservoir included the analysis of
simultaneous changes to non-point BOD loadings defined for some water quality segments, as a
consequence of land use around the future reservoir, based on estimations found in the literature;
and changes of autochthonous BOD loadings as a consequence of decomposition of the
vegetation present in the reservoir area by flooding. Changes of ≤100% of loadings were applied
for testing allochthonous BOD loadings and changes of ≤50% of loadings were applied for testing changes of autochthonous BOD loadings. The metric used to measure sensitivity in this
125 case was the average of percent change of DO concentrations when comparing calibrated model results with those results obtained after applying the sensitivity test.
WASP Segment WRT22
180
160
140
120
100
80
Chlorides (mg/l) 60
40
20
0
7 7 7 7 99 /1997 /1 199 199 8/1997 /7 /2/1997 /9/1997 6 3/ 0/ 2 3 /28/1997 4/4/1997 /11/1997 /18/1997 /25/199 5 5 /1 /2 /3 2/ 3/14/1997 3/21/1997 3 4 4 4 5 5 5
Predicted w/o Load Predicted with Load
Figure 56. Chloride concentrations simulated for WASP segment WRT22 while testing sensitivity of Lake Texoma model to changes in chloride loadings The Tocoma Reservoir model did not reflect any change when increments or reductions
of non-point allochthonous BOD loadings were applied. The magnitude of the load and hydrodynamic conditions contributed to dissipate the effect of these loadings.
However, when changing autochthonous loading conditions, the results were different.
This difference was particularly remarkable when analyzing those values associated with the potentially most critical condition of DO concentration in this system. This is to say, when the progress of the reservoir flooding was simulated.
According to these simulation results, the potential reduction of DO expected for the
Tocoma Reservoir would occur in a lower overall intensity by removing 50% of the
126 autochthonous BOD provided by vegetation decomposition (See Figure 57). However, the final
condition expected for this reservoir would not significantly change by removing the load resulting of the mentioned decomposition (See Figure 58). This can be explained by the important effect that the DO boundary condition exerts over the entire system, as mentioned in a previous section of this chapter.
These results would be useful for managing the resulting environmental impact that the flooding plan could produce not only at the flooding area of the new reservoir, but also over the system of reservoirs located downstream of this new project.
127
Figure 57. Comparative spatial variation of potential DO change during the flooding phase of Tocoma Reservoir by reducing autochthonous BOD loading in 50% in comparison with the effect of the total BOD load
128 Figure 58. Comparative expected spatial distribution of DO for Tocoma reservoir after flooding with 50% and 100% of autochthonous BOD loadings
129 CHAPTER 5
CONCLUSIONS
DYNHYD, the hydrodynamic model coupled to the WASP system was calibrated for two
different systems: Lake Texoma and Tocoma Reservoir. For Lake Texoma, results were
validated with observed water elevation values of the reservoir and they showed a good correspondence between predicted and measured pool water elevation for the reservoir.
Validation of DYNHYD results for Tocoma Reservoir is still pending, considering that the
reservoir is still in project status. However, the performance of the model for describing river
conditions where Tocoma Reservoir is planned to be built was tested, showing a good agreement
between predicted and observed average water velocities when describing the river system.
This simulation also provided some insights of the different hydrodynamic conditions
developed in these reservoirs as a consequence of the hydrology of the main rivers that feed the
reservoirs and the morphology of the reservoirs. A calibrated model for each reservoir was used
to estimate the residence time of a simulated tracer in these reservoirs. A period of 78 days was
estimated as the longest residence time for Lake Texoma, considering the time that it took the
tracer to flow through the Red River Arm. On the other hand, 23 days was the longest travel time
estimated for Tocoma Reservoir, specifically at the Cunaguaro River arm. In the central portion
of this projected reservoir, the travel time was estimated to be 4.8 days.
Five different hydrodynamic areas were identified for Lake Texoma based on the
hydrodynamic results of DYNHYD. One of these areas is acting as an enlarged river (Red River
Zone), and the rest of the areas present clear lentic characteristics (Washita River Transition
130 Zone, Red River Transition Zone, Main Lake Zone, and Big Mineral Creek Zone). Back flow conditions and the largest residence time occur in the Big Mineral Creek Zone.
For Tocoma Reservoir, hydrodynamic results of DYNHYD showed that a clear central
flow is expected to be developed in the future reservoir, with several areas with practically
standing water, as a consequence of the shape of the reservoir and the bed channel.
These differences in the hydrodynamic of the evaluated water bodies were coupled to
water quality observations, and provided additional explanations of the functioning of physical
and chemical processes to be expected in these reservoirs. Also, this information can be used as
an input for the estimation of the assimilative capacity and as a tool for the explanation of the
biological productivity of these reservoirs in a future research.
The WASP model was capable of simulating longitudinal gradients in the physical
processes in these reservoirs. In Lake Texoma, the average of temporal and spatial variability of
chlorides in the reservoir was clearly depicted by the simulation, especially during those periods
of the year when there was not a marked stratification of chloride concentrations. DO
concentrations were used for understanding of temporal and spatial variability of a chemical compound in Tocoma Reservoir. Reaeration was the physical process driving DO concentrations in this reservoir
Transport and dilution processes determine the variations of chloride concentrations in
Lake Texoma. In the Red River arm, advection is the primary driver in the transport of chlorides, whereas in the Washita River arm advection and dispersion processes are both involved in the transport of chlorides. The main source of chlorides in the reservoir is the Red River. However, in the Washita River arm it was necessary to consider additional loads of chlorides for the area around Glasses Creek, in order to adjust predicted values. This was justified based on field
131 measurements of increased conductivity in this area of the reservoir that was not caused by the
Washita River (as measured at the gage station near Dickson) or by the Main Lake zone, which has lower conductivity. Additional studies are required to identify the sources of this additional load of chlorides in this portion of the reservoir, but oil and gas production appear to be likely sources.
The model predicts short-term variability of chloride concentrations in both the river and transition zones for Lake Texoma. However, monthly measurements do not account for this variability. On the other hand, chemical stratification reported for the reservoir is not accounted for in this depth-averaged arrangement of the model. In those areas where a clear difference between concentrations at the top and at the bottom of the water column was reported, predictions of the model underestimated the observed conditions because they are expected to represent an average of the concentration.
In terms of DO, physical processes considered in the Streeter-Phelps approach for DO simulations seems to be good enough for describing the average concentration of DO in Lake
Texoma, and for understanding the potential distribution of DO in the Tocoma Reservoir.
The model fits pretty well the temporal variation of the average concentration for this variable in the different zones of Lake Texoma, even though the patterns in these zones are different. However, some discrepancies were found. These differences between observed and predicted values of DO mainly occur in two ways: