USGS ScienceSearch

Geology topics

Jerad D. Bales

Publications and source records attributed to Jerad D. Bales.

At least 19 recordsLinked to original sources

The water-energy nexus: An earth science perspective

Water availability and use are closely connected with energy development and use. Water cannot be delivered to homes, businesses, and industries without energy, and most forms of energy development require large amounts of water. The United States faces two significant and sometimes competing challenges: to provide sustainable supplies of freshwater for humans and ecosystems and to ensure adequate sources of energy for future generations. This report reviews the complex ways in which water and energy are interconnected and describes the earth science data collection and research that can help the Nation address these important challenges. The earth sciences have been a cornerstone in developing our current understanding of the water-energy nexus. A full understanding of the nexus, however, is limited by uncertainty in our knowledge of fundamental issues, such as the quantity of freshwater that is available, the amount of water that is used in energy development, the effects that emerging energy development technologies have on water quality and quantity, and the amount of energy required to treat and deliver freshwater. Enhanced data collection and research can improve our understanding of these important issues and thereby lay the groundwork for informed resource management. Relevant earth science issues analyzed and discussed herein include freshwater availability; water use; ecosystems health; assessment of saline water resources; assessment of fossil-fuel, uranium, and geothermal resources; subsurface injection of wastewater and carbon dioxide and related induced seismicity; climate change and its effect on water availability and energy production; byproducts and waste streams of energy development; emerging energy-development technologies; and energy for water treatment and delivery.

Circular

Progress in data collection and dissemination in water resources – 1974-2014

In the 50 years since the founding of the American Water Resources Association (AWRA), there has been tremendous and likely unforeseen progress in water-re- sources data collection and dissemination. Langford and Doyel (1974) (henceforth L&D) described progress during the decade following the founding of AWRA, and focused their description around seven topics. L&D described the changes as being “more philosophical than technical,” and noted the importance to the water-resources com- munity of the more than 30 Federal Acts or Amendments enacted in the decade. The purpose of this article is to provide an update to L&D by reviewing L&D’s predictions of anticipated changes in water resources data collection and dissemi-nation, providing an overview of some of the drivers of change in the water-resources community in the last 40 years, identifying some key advances in water-resources data collection and dissemination since 1974, and out-lining some important near-term challenges. The overview is necessarily incomplete, but represents one perspective based on years of collaboration throughout the water-resources community.

Water Resources Impact

Storm surge modeling and applications in coastal areas

This chapter introduces the reader to a wide spectrum of storm surge modeling systems used to assess the impact of tropical cyclones, covering a range of numerical methods, model domains, forcing and boundary conditions, and purposes. New technologies to obtain data such as deployment of temporary sensors and remote sensing practices to support modeling are also presented. Extensive storm surge modeling applications have been made with existing modeling systems and some of them are described in this chapter. The authors recognize the importance of evaluating river-ocean interactions in coastal environments during tropical cyclones. Therefore, the coupling of hydraulic (riverine) and storm surge models is discussed. In addition, results from studies performed in the coast of India are shown which generated maps to help emergency managers and reduce risk due to coastal inundation.

Book chapter

Use of a two-dimensional hydrodynamic model to evaluate extreme flooding and transport of dissolved solids through Devils Lake and Stump Lake, North Dakota, 2006

The U.S. Geological Survey in cooperation with the North Dakota Department of Transportation, North Dakota State Water Commission, and U.S. Army Corps of Engineers, developed a two-dimensional hydrodynamic model of Devils Lake and Stump Lake, North Dakota to be used as a hydrologic tool for evaluating the effects of different inflow scenarios on water levels, circulation, and the transport of dissolved solids through the lake. The numerical model, UnTRIM, and data primarily collected during 2006 were used to develop and calibrate the Devils Lake model. Performance of the Devils Lake model was tested using 2009 data. The Devils Lake model was applied to evaluate the effects of an extreme flooding event on water levels and hydrological modifications within the lake on the transport of dissolved solids through Devils Lake and Stump Lake. For the 2006 calibration, simulated water levels in Devils Lake compared well with measured water levels. The maximum simulated water level at site 1 was within 0.13 feet of the maximum measured water level in the calibration, which gives reasonable confidence that the Devils Lake model is able to accurately simulate the maximum water level at site 1 for the extreme flooding scenario. The timing and direction of winddriven fluctuations in water levels on a short time scale (a few hours to a day) were reproduced well by the Devils Lake model. For this application, the Devils Lake model was not optimized for simulation of the current speed through bridge openings. In future applications, simulation of current speed through bridge openings could be improved by more accurate definition of the bathymetry and geometry of select areas in the model grid. As a test of the performance of the Devils Lake model, a simulation of 2009 conditions from April 1 through September 30, 2009 was performed. Overall, errors in inflow estimates affected the results for the 2009 simulation; however, for the rising phase of the lakes, the Devils Lake model accurately simulated the faster rate of rise in Devils Lake than in Stump Lake, and timing and direction of wind-driven fluctuations in water levels on a short time scale were reproduced well. To help the U.S. Army Corps of Engineers determine the elevation to which the protective embankment for the city of Devils Lake should be raised, an extreme flooding scenario based on an inflow of one-half the probable maximum flood was simulated. Under the conditions and assumptions of the extreme flooding scenario, the water level for both lakes reached a maximum water level around 1,461.9 feet above the National Geodetic Vertical Datum of 1929. One factor limiting the extent of pumping from the Devils Lake State Outlet is sulfate concentrations in West Bay. If sulfate concentrations can be reduced in West Bay, pumping from the Devils Lake State Outlet potentially can increase. The Devils Lake model was used to simulate the transport of dissolved solids using specific conductance data as a surrogate for sulfate. Because the transport of dissolved solids was not calibrated, results from the simulations were not actual expected concentrations. However, the effects of hydrological modifications on the transport of dissolved solids could be evaluated by comparing the effects of hydrological modifications relative to a baseline scenario in which no hydrological modifications were made. Four scenarios were simulated: (1) baseline condition (no hydrological modification), (2) diversion of Channel A, (3) reduction of the area of water exchange between Main Bay and East Bay, and (4) combination of scenarios 2 and 3. Relative to scenario 1, mean concentrations in West Bay for scenarios 2 and 4 were reduced by approximately 9 percent. Given that there is no change in concentration for scenario 3, but about a 9-percent reduction in concentration for scenario 4, the diversion of Channel A was the only hydrologic modification that appeared to have the potential to reduce sulfate c

Scientific Investigations Report

Assessing the response of the Pamlico Sound, North Carolina, USA to human and climatic disturbances: Management implications

The Pamlico Sound (PS) with its sub-estuaries is the largest lagoonal ecosystem in the United States. It exhibits periodically strong salinity stratification and an average freshwater residence time of 1 year for the sound proper. This relatively long residence time promotes effective use and cycling of nutrients, allowing the system to support high rates of primary and secondary production, and serve as a vitally important fisheries nursery. This hydrologic characteristic also makes the system highly sensitive to nutrient over-enrichment and eutrophication. The PS is experiencing ecological change in response to increasing human activity and climatic perturbations. Human impacts include a rise in nutrient, sediment, and other pollutant loads that accompany urbanization and agricultural and industrial growth in its watersheds and airsheds. Since the mid-1990s, the PS has witnessed a sudden rise in tropical storm and hurricane impacts, with eight hurricanes and four tropical storms having made landfall in the PS watershed during the 1996 to 2007 period. Each of these storms had unique hydrologic, nutrient, and other pollutant loading effects. In addition, since the early 2000s, the region has experienced record droughts, which are continuing. Variable freshwater discharges from storms and droughts have caused large oscillations in nutrient enrichment, reflected ultimately in differential phytoplankton production, biomass, and community compositional responses. Floodwaters from the two wettest hurricanes, Fran (1996) and Floyd (1999), and from Tropical Storm Ernesto (2006) exerted long-term (months) effects on hydrology, nutrient loads, and algal production. Windy but relatively dry hurricanes, like Irene (1999) and Isabel (2003), caused strong vertical mixing, storm surges, but relatively minor changes in river flow, flushing, and nutrient loads. These contrasting effects are accompanied by biogeochemical (hypoxia, nutrient cycling) and habitat alterations, and associated food web disturbances. Each storm type influenced algal growth and compositional dynamics; however, their respective ecological impacts differed substantially. Changes in hydrologic and wind forcing resulting from changes in frequency and intensity of storms and droughts strongly influence water and habitat quality. These changes must be integrated with nutrient loading/dilution effects when assessing and predicting ecological responses to nutrient and hydrologic variability on this and other large lagoonal ecosystems.

North Carolina

Relation Between Flow and Dissolved Oxygen in the Roanoke River Between Roanoke Rapids and Jamesville, North Carolina, 1998-2005

Understanding the relation between dam release characteristics and downstream water quality in the lower Roanoke River, North Carolina, is important for natural-resource management and ecosystem protection. Data from four raingages, four water-quality monitoring sites, and one streamflow-measurement site were used to identify statistical relations and discernible quantitative or qualitative patterns linking Roanoke River instream dissolved-oxygen (DO) levels to releases at Roanoke Rapids Dam for the period 1998-2005. The time-series DO data, complicated by the occurrence of major hurricanes in the short period of hourly DO data collection at the dam, present a mixed picture of the effects of hydropower peaking (a technique used by hydropower dam operators to produce electricity when consumption is high by passing a large volume of water through the dam turbines, which dramatically increases the volume of flow below the dam) on downstream DO. Other than in 2003 when dissolved-oxygen concentrations in the Roanoke River were likely affected by runoff from Hurricane Isabel rains, there were not consistent, statistically significant differences detected in the annual medians of hourly and(or) daily DO values during peaking versus nonpeaking periods. Along the Roanoke River, downstream of Roanoke Rapids Dam at Oak City, North Carolina, using a 95-percent confidence interval, the median value of the May-November daily mean DO concentrations for each year was lower during peaking periods for 2 years, higher for 2 years, and not significantly different for 4 years. Downstream at Jamesville, North Carolina, also using a 95-percent confidence interval, the median value of the annual May-November daily mean DO concentrations during hydropower peaking was lower for 4 years, higher for 2 years, and not significantly different for 2 years. In summary, the effect of hydropower peaking on downstream DO was inconsistent. Conversely, large precipitation events downstream from the dam resulted in consistent, statistically significant decreases in DO in the mainstem of the Roanoke River at Oak City and Jamesville.

North Carolina

Simulation of Hydrodynamics and Water Quality in Pueblo Reservoir, Southeastern Colorado, for 1985 through 1987 and 1999 through 2002

Pueblo Reservoir is west of Pueblo, Colorado, and is an important water resource for southeastern Colorado. The reservoir provides irrigation, municipal, and industrial water to various entities throughout the region. In anticipation of increased population growth, the cities of Colorado Springs, Fountain, Security, and Pueblo West have proposed building a pipeline that would be capable of conveying 78 million gallons of raw water per day (240 acre-feet) from Pueblo Reservoir. The U.S. Geological Survey, in cooperation with Colorado Springs Utilities and the Bureau of Reclamation, developed, calibrated, and verified a hydrodynamic and water-quality model of Pueblo Reservoir to describe the hydrologic, chemical, and biological processes in Pueblo Reservoir that can be used to assess environmental effects in the reservoir. Hydrodynamics and water-quality characteristics in Pueblo Reservoir were simulated using a laterally averaged, two-dimensional model that was calibrated using data collected from October 1985 through September 1987. The Pueblo Reservoir model was calibrated based on vertical profiles of water temperature and dissolved-oxygen concentration, and water-quality constituent concentrations collected in the epilimnion and hypolimnion at four sites in the reservoir. The calibrated model was verified with data from October 1999 through September 2002, which included a relatively wet year (water year 2000), an average year (water year 2001), and a dry year (water year 2002). Simulated water temperatures compared well to measured water temperatures in Pueblo Reservoir from October 1985 through September 1987. Spatially, simulated water temperatures compared better to measured water temperatures in the downstream part of the reservoir than in the upstream part of the reservoir. Differences between simulated and measured water temperatures also varied through time. Simulated water temperatures were slightly less than measured water temperatures from March to May 1986 and 1987, and slightly greater than measured data in August and September 1987. Relative to the calibration period, simulated water temperatures during the verification period did not compare as well to measured water temperatures. In general, simulated dissolved-oxygen concentrations for the calibration period compared well to measured concentrations in Pueblo Reservoir. Spatially, simulated concentrations deviated more from the measured values at the downstream part of the reservoir than at other locations in the reservoir. Overall, the absolute mean error ranged from 1.05 (site 1B) to 1.42 milligrams per liter (site 7B), and the root mean square error ranged from 1.12 (site 1B) to 1.67 milligrams per liter (site 7B). Simulated dissolved oxygen in the verification period compared better to the measured concentrations than in the calibration period. The absolute mean error ranged from 0.91 (site 5C) to 1.28 milligrams per liter (site 7B), and the root mean square error ranged from 1.03 (site 5C) to 1.46 milligrams per liter (site 7B). Simulated total dissolved solids generally were less than measured total dissolved-solids concentrations in Pueblo Reservoir from October 1985 through September 1987. The largest differences between simulated and measured total dissolved solids were observed at the most downstream sites in Pueblo Reservoir during the second year of the calibration period. Total dissolved-solids data were not available from reservoir sites during the verification period, so in-reservoir specific-conductance data were compared to simulated total dissolved solids. Simulated total dissolved solids followed the same patterns through time as the measured specific conductance data during the verification period. Simulated total nitrogen concentrations compared relatively well to measured concentrations in the Pueblo Reservoir model. The absolute mean error ranged from 0.21 (site 1B) to 0.27 milligram per liter as nitrogen (sites 3B and 7

Scientific Investigations Report

Review of Selected Documents Related to Flooding at City of Salisbury Facilities on the Yadkin River Upstream from High Rock Dam, North Carolina, September 2007

This report documents a review of the hydraulic and sediment-transport models developed by the City of Salisbury, Alcoa Power Generating, Inc., and the Federal Energy Regulatory Commission to address issues of flooding and sedimentation in the vicinity of Salisbury's water-supply intake 19.4 miles upstream from High Rock Dam. The objective of the review was to determine if the modeling results submitted by Salisbury clearly demonstrate that the presence of High Rock Dam has led to an increase in water levels at Salisbury facilities or, conversely, if the documents of Alcoa Power Generating, Inc., demonstrate that High Rock Dam has not had an effect on water levels at Salisbury facilities. No new data were collected as a part of the review, and the models developed by involved parties were not tested during the review. Some historical discharge-measurement notes and previously published reports were checked as part of the review. The one-dimensional hydraulic modeling results submitted by Alcoa Power Generating, Inc., did not assess the effects of changes in bathymetry on changes in flood levels at Salisbury's facilities because pre-impoundment conditions were not simulated. Hydraulic modeling performed by consultants for the City of Salisbury seems to indicate that both the presence of the dam in the absence of any post-impoundment sedimentation and changes in bathymetry between pre-impoundment and 1997 conditions have resulted in increased water levels relative to pre-impoundment conditions at Salisbury facilities on the Yadkin River for a fairly wide range of flows. The degree to which the dam and the changes in bathymetry have affected flood levels at the Salisbury facilities relative to pre-impoundment conditions is open to discussion because of uncertainty in topographic/bathymetric data and the absence of calibration and sensitivity testing of the hydraulic models. None of the three hydraulic models appears to have been calibrated to or tested against measurements, and no sensitivity testing was reported. Sediment-transport modeling results submitted by the City of Salisbury were calibrated, well documented, and provide a good understanding of the expected growth of the sediment delta in the upper end of High Rock Lake. Simulations made using this model seem to have demonstrated that the presence of the dam and the growth of the delta have resulted in increases in water-surface elevations at Salisbury's facilities over a range of flows and that these increases are expected to increase through time if current conditions remain unchanged.

Open-File Report

LiDAR-Derived Flood-Inundation Maps for Real-Time Flood-Mapping Applications, Tar River Basin, North Carolina

Flood-inundation maps were created for selected streamgage sites in the North Carolina Tar River basin. Light detection and ranging (LiDAR) data with a vertical accuracy of about 20 centimeters, provided by the Floodplain Mapping Information System of the North Carolina Floodplain Mapping Program, were processed to produce topographic data for the inundation maps. Bare-earth mass point LiDAR data were reprocessed into a digital elevation model with regularly spaced 1.5-meter by 1.5-meter cells. A tool was developed as part of this project to connect flow paths, or streams, that were inappropriately disconnected in the digital elevation model by such features as a bridge or road crossing. The Hydraulic Engineering Center-River Analysis System (HEC-RAS) model, developed by the U.S. Army Corps of Engineers, was used for hydraulic modeling at each of the study sites. Eleven individual hydraulic models were developed for the Tar River basin sites. Seven models were developed for reaches with a single gage, and four models were developed for reaches of the Tar River main stem that receive flow from major gaged tributaries, or reaches in which multiple gages were near one another. Combined, the Tar River hydraulic models included 272 kilometers of streams in the basin, including about 162 kilometers on the Tar River main stem. The hydraulic models were calibrated to the most current stage-discharge relations at 11 long-term streamgages where rating curves were available. Medium- to high-flow discharge measurements were made at some of the sites without rating curves, and high-water marks from Hurricanes Fran and Floyd were available for high-stage calibration. Simulated rating curves matched measured curves over the full range of flows. Differences between measured and simulated water levels for a specified flow were no more than 0.44 meter and typically were less. The calibrated models were used to generate a set of water-surface profiles for each of the 11 modeled reaches at 0.305-meter increments for water levels ranging from bankfull to approximately the highest recorded water level at the downstream-most gage in each modeled reach. Inundated areas were identified by subtracting the water-surface elevation in each 1.5-meter by 1.5-meter grid cell from the land-surface elevation in the cell through an automated routine that was developed to identify all inundated cells hydraulically connected to the cell at the downstream-most gage in the model domain. Inundation maps showing transportation networks and orthoimagery were prepared for display on the Internet. These maps also are linked to the U.S. Geological Survey North Carolina Water Science Center real-time streamflow website. Hence, a user can determine the near real-time stage and water-surface elevation at a U.S. Geological Survey streamgage site in the Tar River basin and link directly to the flood-inundation maps for a depiction of the estimated inundated area at the current water level. Although the flood-inundation maps represent distinct boundaries of inundated areas, some uncertainties are associated with these maps. These are uncertainties in the topographic data for the hydraulic model computational grid and inundation maps, effective friction values (Manning's n), model-validation data, and forecast hydrographs, if used. The Tar River flood-inundation maps were developed by using a steady-flow hydraulic model. This assumption clearly has less of an effect on inundation maps produced for low flows than for high flows when it typically takes more time to inundate areas. A flood in which water levels peak and fall slowly most likely will result in more inundation than a similar flood in which water levels peak and fall quickly. Limitations associated with the steady-flow assumption for hydraulic modeling vary from site to site. The one-dimensional modeling approach used in this study resulted in good agreement between measurements and simulations. T

North Carolina

Automated Routines for Calculating Whole-Stream Metabolism: Theoretical Background and User's Guide

In order to standardize methods and facilitate rapid calculation and archival of stream-metabolism variables, the Stream Metabolism Program was developed to calculate gross primary production, net ecosystem production, respiration, and selected other variables from continuous measurements of dissolved-oxygen concentration, water temperature, and other user-supplied information. Methods for calculating metabolism from continuous measurements of dissolved-oxygen concentration and water temperature are fairly well known, but a standard set of procedures and computation software for all aspects of the calculations were not available previously. The Stream Metabolism Program addresses this deficiency with a stand-alone executable computer program written in Visual Basic.NET?, which runs in the Microsoft Windows? environment. All equations and assumptions used in the development of the software are documented in this report. Detailed guidance on application of the software is presented, along with a summary of the data required to use the software. Data from either a single station or paired (upstream, downstream) stations can be used with the software to calculate metabolism variables.

Techniques and Methods

Flow and salt transport in the Suwannee River estuary, Florida, 1999–2000: Analysis of data and three-dimensional simulations

A three-dimensional numerical model was developed to assist in the evaluation of the effects of changes in freshwater flow on the salinity regime of the lower Suwannee River, its estuary, and Suwannee Sound. Hydrodynamic and salt-transport modeling were supported by data from a comprehensive data-collection network operated in the lower Suwannee River during 1998-2000. The study area included all of the downstream-most 12 kilometers of the Suwannee River, Suwannee Sound, and part of the Gulf of Mexico. Development, calibration, and application of the hydrodynamic and salt-transport model were completed by using data primarily collected during October 1999-September 2000.

Florida

Simulation of constituent transport in the Red River of the North basin, North Dakota and Minnesota, during unsteady-flow conditions, 1977 and 2003-04

The Bureau of Reclamation identified eight water-supply alternatives for the Red River Valley Water Supply Project. Of those alternatives, six were considered for this study. Those six alternatives include a no-action alternative, two in-basin alternatives, and three interbasin alternatives. To address concerns of stakeholders and to provide information for an environmental impact statement, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, developed and applied a water-quality model to simulate the transport of total dissolved solids, sulfate, chloride, sodium, and total phosphorus during unsteady-flow conditions and to simulate the effects of the water-supply alternatives on water quality in the Red River and the Sheyenne River. The physical domain of the model, hereinafter referred to as the Red River model, includes the Red River from Wahpeton, North Dakota, to Emerson, Manitoba, and the Sheyenne River from below Baldhill Dam, North Dakota, to the confluence with the Red River. Boundary conditions were specified for May 15 through October 31, 2003, and January 15 through June 30, 2004. Measured streamflow data were available for August 1 through October 31, 2003, and April 1 through June 30, 2004, but water-quality data were available only for September 15 through 16, 2003, and May 10 through 13, 2004. The water-quality boundary conditions were assumed to be time invariant for the entire calibration period and to be equal to the measured value. The average difference between the measured and simulated streamflows was less than 4 percent for both calibration periods, and most differences were less than 2 percent. The average differences are considered to be acceptable because the differences are less than 5 percent, or the same as the error that would be expected in a typical streamflow measurement. Simulated total dissolved solids, sulfate, chloride, and sodium concentrations generally were less than measured concentrations for both calibration periods. The average absolute differences generally were less than 25 percent. Total phosphorus was simulated as a nonconservative constituent by assuming that concentrations change according to a first-order decay rate. The average difference between the measured and simulated total phosphorus concentrations was 6.2 percent for the 2003 calibration period and -24 percent for the 2004 calibration period. The Red River model demonstrates sensitivity to changes in boundary conditions so a reasonable assumption is that the model can be used to compare relative effects of the various water-supply alternatives. The calibrated Red River model was used to simulate the effects of the six water-supply alternatives by using measured streamflows for September 1, 1976, through August 31, 1977, when streamflows throughout the Red River Basin were relatively low. Streamflows for the Red River at Fargo, North Dakota, were less than 17.9 cubic feet per second on 159 days of that 12-month period, and monthly average streamflows for the Red River at Grand Forks, North Dakota, and the Red River at Emerson, Manitoba, were less than 30 percent of the respective long-term average monthly streamflows for 11 of the 12 months during September 1976 through August 1977. Water-quality boundary conditions were generated using a stochastic approach in which probability distributions derived from all available historical data on instream concentrations were used to produce daily concentrations at model boundaries. Return flow concentrations were estimated from source concentrations and current (2006) wastewater-treatment technology. Because no historical information on ungaged local inflow constituent concentrations is available to estimate those boundary conditions, time-invariant concentrations for the low-flow 2003 calibration period were used as the ungaged local inflow boundary conditions. The effects of the water-supply alternatives on water quality in the Red River and

Minnesota, North Dakota

Simulation of conservative-constituent transport in the Red River of the North Basin, North Dakota and Minnesota, 2003-04

Population growth along with possible future droughts in the Red River of the North (Red River) Basin will create an increasing need for reliable water supplies. Therefore, as a result of the Dakota Water Resources Act of 2000, the Bureau of Reclamation identified eight water-supply alternatives (including a no-action alternative) to meet future water needs in the basin. Because of concerns about the possible effects of the alternatives on water quality in the Red River and the Sheyenne River and in Lake Winnipeg, Manitoba, the Bureau of Reclamation needs to prepare an environmental impact statement that describes the specific environmental effects of each alternative. To provide information for the environmental impact statement, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, conducted a study to develop and apply a water-quality model, hereinafter referred to as the Red River water-quality model, to part of the Red River and the Sheyenne River to simulate conservative-constituent transport in the Red River Basin. The Red River water-quality model is a one-dimensional, steady-state flow and transport model for selected constituents in the Red River and the Sheyenne River. The model simulates the flow and transport of total dissolved solids, sulfate, and chloride during steady-state conditions. The physical model domain includes the Red River from the confluence of the Bois de Sioux and Otter Tail Rivers to the Red River at Emerson, Manitoba, and the Sheyenne River from above Harvey, N. Dak., to the confluence with the Red River. The Red River water-quality model was calibrated and tested using data collected at 34 sites from September 15 through 16, 2003, and from May 10 through 13, 2004. Water-quality samples were collected during low, steady-flow conditions from September 15 through 16, 2003, and during medium, unsteady-flow conditions from May 10 through 13, 2004. The simulated total dissolved-solids, sulfate, and chloride concentrations generally were within 5 percent of the measured concentrations. The Red River water-quality model was used to simulate conservative-constituent transport in the Red River and the Sheyenne River for the eight water-supply alternatives identified by the Bureau of Reclamation. For the first set of eight simulations, September 2003 streamflows were used with projected 2050 return flows and withdrawals. For the second set of eight simulations, the September 2003 streamflows were reduced by 25 percent. The simulated concentrations for three of the alternatives generally were lower than for the no-action alternative. Of those alternatives, one would result in a decrease in concentrations for two constituents, one would result in a decrease in concentrations for all three constituents, and one would result in a decrease in concentrations for one constituent and an increase in concentrations for another constituent. For four of the alternatives, the differences between the mean simulated concentrations were less than calibration errors, indicating the effects of those alternatives on water quality in the rivers is uncertain. The effects of reduced streamflow on simulated total dissolved-solids, sulfate, and chloride concentrations were greatest for alternative 2. Reduced streamflow probably has an effect on simulated total dissolved-solids concentrations for alternatives 2, 3, 5, and 7 and on simulated sulfate concentrations for alternatives 2 and 5. Except for alternative 2, reduced streamflow had little effect on simulated chloride concentrations.

Scientific Investigations Report

Two-dimensional hydrodynamic simulation of surface-water flow and transport to Florida Bay through the Southern Inland and Coastal Systems (SICS)

Successful restoration of the southern Florida ecosystem requires extensive knowledge of the physical characteristics and hydrologic processes controlling water flow and transport of constituents through extremely low-gradient freshwater marshes, shallow mangrove-fringed coastal creeks and tidal embayments, and near-shore marine waters. A sound, physically based numerical model can provide simulations of the differing hydrologic conditions that might result from various ecosystem restoration scenarios. Because hydrology and ecology are closely linked in southern Florida, hydrologic model results also can be used by ecologists to evaluate the degree of ecosystem restoration that could be achieved for various hydrologic conditions. A robust proven model, SWIFT2D, (Surface-Water Integrated Flow and Transport in Two Dimensions), was modified to simulate Southern Inland and Coastal Systems (SICS) hydrodynamics and transport conditions. Modifications include improvements to evapotranspiration and rainfall calculation and to the algorithms that describe flow through coastal creeks. Techniques used in this model should be applicable to other similar low-gradient marsh settings in southern Florida and elsewhere. Numerous investigations were conducted within the SICS area of southeastern Everglades National Park and northeastern Florida Bay to provide data and parameter values for model development and testing. The U.S. Geological Survey and the National Park Service supported investigations for quantification of evapotranspiration, vegetative resistance to flow, wind-induced flow, land elevations, vegetation classifications, salinity conditions, exchange of ground and surface waters, and flow and transport in coastal creeks and embayments. The good agreement that was achieved between measured and simulated water levels, flows, and salinities through minimal adjustment of empirical coefficients indicates that hydrologic processes within the SICS area are represented properly in the SWIFT2D model, and that the spatial and temporal resolution of these processes in the model is adequate. Sensitivity analyses were conducted to determine the effect of changes in boundary conditions and parameter values on simulation results, which aided in identifying areas of greatest uncertainty in the model. The parameter having the most uncertainty (most in need of further field study) was the flow coefficient for coastal creeks. Smaller uncertainties existed for wetlands frictional resistance and wind. Evapotranspiration and boundary inflows indicated the least uncertainty as determined by varying parameters used in their formulation and definition. Model results indicated that wind was important in reversing coastal creek flows. At Trout Creek (the major tributary connecting Taylor Slough wetlands with Florida Bay), flow in the landward direction was not simulated properly unless wind forcing was included in the simulation. Simulations also provided insight into the major influence that wind has on salinity mixing along the coast, the varying distribution of wetland flows at differing water levels, and the importance of topography in controlling flows to the coast. Slight topographic variations were shown to highly influence the routing of water. A multiple regression analysis was performed to relate inflows at the northern boundary of Taylor Slough bridge to a major pump station (S-332) north of the SICS model area. This analysis allows Taylor Slough bridge boundary conditions to be defined for the model from operating scenarios at S-332, which should facilitate use of the SICS model as an operational tool.

Florida

Relations among floodplain water levels, instream dissolved-oxygen conditions, and streamflow in the lower Roanoke River, North Carolina, 1997-2001

The lower Roanoke River corridor in North Carolina contains a floodplain of national significance. Data from a network of 1 streamflow-measurement site, 13 river-stage sites, 13 floodplain water-level sites located along 4 transects, and 5 in situ water-quality monitoring sites were used to characterize temporal and spatial variations of floodplain and river water levels during 1997-2000 and to describe dissolved-oxygen conditions in the lower Roanoke River for the period 1998-2001. Major differences in the relation of floodplain inundation to flow occurred both among sites at a given transect and among transects. Several floodplain sites were inundated for the full range of flow conditions measured during the study. These included one site on the Big Swash transect (at about river kilometer 119); one site on the Broadneck Swamp transect (river kilometer 97), which was inundated 91 percent of the time during the study; one site on the Devils Gut transect (river kilometer 44), which was inundated throughout the study; and three sites on the Cow Swamp transect (near river kilometer 10). The relation of floodplain inundation depth to Roanoke River flow was highly variable among sites. There was no relation between flow and inundation depth at one of the Big Swash sites or at any of the four Cow Swamp sites. At two of the Big Swash transect sites, there was some relation between inundation depth and 10-day mean flow for flows greater than 700 cubic meters per second. A relatively strong relation between inundation depth and 10-day mean flow occurred at two of the Broadneck Swamp sites and, to a lesser degree, at two of the Devils Gut transect sites. There was much greater interannual variability in floodplain water levels, as represented by the difference between the maximum and minimum daily water level for a given calendar date during January-May and September-October than during the summer and late fall months. If data from this study are representative of long-term conditions, then this means that there is less uncertainty about what future floodplain water levels will be during June-August and November-December than during other months. Rates of ground-water decline, primarily due to evapotranspiration, were fairly similar at all sites, ranging from about 3 to 4 centimeters per day. For a 10-day mean flow of 300 cubic meters per second, an evaporative loss of 2 centimeters per day is equal to about 56 cubic meters per second. Evapotranspiration rates are much lower during the fall and winter months, so losses of river flow to floodplain processes likely are much lower during those months. The ground-water gradient at most sites was from the floodplain to the river, indicating a potential for ground-water movement into the river from the floodplain. At two of the Devils Gut sites, however, the water level often was higher in the river than in the floodplain when floodplain sites were not inundated. This indicates that there is a potential for river water to move as ground water from the river into the floodplain. It seems likely that this feature observed at the Devils Gut transect occurs elsewhere in the lower Roanoke River corridor. Dissolved-oxygen concentrations typically decrease with increasing distance from Roanoke Rapids Dam. During the 1998-2001 study period, the median dissolved-oxygen concentration at Halifax (river kilometer 187), the upstream-most station, was 8.4 milligrams per liter, and the median concentration at the downstream-most station (NC-45, bottom sensor; river kilometer 2.6) was 6.6 milligrams per liter. Several synoptic measurements of dissolved-oxygen concentration down the river identified the presence of a dissolved-oxygen sag in the vicinity of Halifax, with some recovery of concentrations between Halifax and about Scotland Neck at river kilometer 156. Data from the synoptic measurements also indicated that the greatest rate of dissolved-oxygen change with distance along the riv

North Carolina

Simulation of Temperature, Nutrients, Biochemical Oxygen Demand, and Dissolved Oxygen in the Catawba River, South Carolina, 1996-97

Time-series plots of dissolved-oxygen concentrations were determined for various simulated hydrologic and point-source loading conditions along a free-flowing section of the Catawba River from Lake Wylie Dam to the headwaters of Fishing Creek Reservoir in South Carolina. The U.S. Geological Survey one-dimensional dynamic-flow model, BRANCH, was used to simulate hydrodynamic data for the Branched Lagrangian Transport Model. Waterquality data were used to calibrate the Branched Lagrangian Transport Model and included concentrations of nutrients, chlorophyll a, and biochemical oxygen demand in water samples collected during two synoptic sampling surveys at 10 sites along the main stem of the Catawba River and at 3 tributaries; and continuous water temperature and dissolved-oxygen concentrations measured at 5 locations along the main stem of the Catawba River. A sensitivity analysis of the simulated dissolved-oxygen concentrations to model coefficients and data inputs indicated that the simulated dissolved-oxygen concentrations were most sensitive to watertemperature boundary data due to the effect of temperature on reaction kinetics and the solubility of dissolved oxygen. Of the model coefficients, the simulated dissolved-oxygen concentration was most sensitive to the biological oxidation rate of nitrite to nitrate. To demonstrate the utility of the Branched Lagrangian Transport Model for the Catawba River, the model was used to simulate several water-quality scenarios to evaluate the effect on the 24-hour mean dissolved-oxygen concentrations at selected sites for August 24, 1996, as simulated during the model calibration period of August 23 27, 1996. The first scenario included three loading conditions of the major effluent discharges along the main stem of the Catawba River (1) current load (as sampled in August 1996); (2) no load (all point-source loads were removed from the main stem of the Catawba River; loads from the main tributaries were not removed); and (3) fully loaded (in accordance with South Carolina Department of Health and Environmental Control National Discharge Elimination System permits). Results indicate that the 24-hour mean and minimum dissolved-oxygen concentrations for August 24, 1996, changed from the no-load condition within a range of - 0.33 to 0.02 milligram per liter and - 0.48 to 0.00 milligram per liter, respectively. Fully permitted loading conditions changed the 24-hour mean and minimum dissolved-oxygen concentrations from - 0.88 to 0.04 milligram per liter and - 1.04 to 0.00 milligram per liter, respectively. A second scenario included the addition of a point-source discharge of 25 million gallons per day to the August 1996 calibration conditions. The discharge was added at S.C. Highway 5 or at a location near Culp Island (about 4 miles downstream from S.C. Highway 5) and had no significant effect on the daily mean and minimum dissolved-oxygen concentration. A third scenario evaluated the phosphorus loading into Fishing Creek Reservoir; four loading conditions of phosphorus into Catawba River were simulated. The four conditions included fully permitted and actual loading conditions, removal of all point sources from the Catawba River, and removal of all point and nonpoint sources from Sugar Creek. Removing the point-source inputs on the Catawba River and the point and nonpoint sources in Sugar Creek reduced the organic phosphorus and orthophosphate loadings to Fishing Creek Reservoir by 78 and 85 percent, respectively.

South Carolina

Mountain Island Lake, North Carolina: Analysis of ambient conditions and simulation of hydrodynamics, constituent transport, and water-quality characteristics, 1996–97

Mountain Island Lake is an impoundment of the Catawba River in North Carolina and supplies drinking water to more than 600,000 people in Charlotte, Gastonia, Mount Holly, and several other communities. The U.S. Geological Survey, in cooperation with the Charlotte-Mecklenburg Utilities, conducted an investigation of the reservoir to characterize hydrologic and water-quality conditions and to develop and apply a simulation model to predict the response of the reservoir to changes in constituent loadings or the flow regime. During 1996–97, flows into Mountain Island Lake were dominated by releases from Cowans Ford Dam on Lake Norman, with more than 85 percent of the total inflow to the reservoir coming from Lake Norman. Riverbend Steam Station discharges accounted for about 12 percent of the inflows to the reservoir, and inflows from tributary streams contributed less than 1.5 percent of the total inflows. Releases through Mountain Island Dam accounted for about 81 percent of outflows from the reservoir, while Riverbend Steam Station withdrawals, which were equal to discharge from the facility, constituted about 13 percent of the reservoir withdrawals. About 5.5 percent of the withdrawals from the reservoir were for water supply. Strong thermal stratification was seldom observed in Mountain Island Lake during April 1996-September 1997. As a result, dissolved-oxygen concentrations were only infrequently less than 4 milligrams per liter, and seldom less than 5 milligrams per liter throughout the entire reservoir, including the coves. The Riverbend Steam Station thermal discharge had a pronounced effect on surface-water temperatures near the outfall. McDowell Creek, which drains to McDowell Creek cove, receives treated wastewater from a large municipal facility and has exhibited signs of poor water-quality conditions in the past. During April 1996-September 1997, concentrations of nitrate, ammonia, total phosphorus, and chlorophyll a were higher in McDowell Creek cove than elsewhere throughout the reservoir. Nevertheless, the highest chlorophyll a concentration measured during the study was 13 micrograms per liter—well below the North Carolina ambient water-quality standard of 40 micrograms per liter. In the mainstem of the reservoir, near-bottom ammonia concentrations occasionally were greater than near-surface concentrations. However, the relatively large top-to-bottom differences in ammonia and phosphorus that have been observed in other Catawba River reservoirs were not present in Mountain Island Lake. External loadings of suspended solids, nitrogen, phosphorus, and biochemical oxygen demand were determined for May 1996-April 1997. Flows through Cowans Ford Dam contributed more than 80 percent of the biochemical oxygen demand and nitrogen load to the reservoir, with McDowell Creek contributing about 15 percent of the biochemical oxygen demand load. In contrast, McDowell Creek contributed about half of the phosphorus load to the reservoir, while inflows through Cowans Ford Dam contributed about one-fourth of the phosphorus load, and the McDowell Creek wastewater-treatment plant contributed about 15 percent of the total phosphorus load. The remainder of the phosphorus loadings came from Gar Creek and the discharge from the Riverbend ash settling pond. Mountain Island Lake is a relatively small (11.3-square-kilometer surface area) impoundment. An area of 181 square kilometers drains directly to the reservoir, but much of this area is undergoing development. In addition, the reservoir receives treated effluent from a municipal wastewater-treatment facility. The two-dimensional, laterally averaged model CE-QUAL-W2 was applied to Mountain Island Lake. The model was configured to simulate water level, water temperature, and 12 water-quality constituents. The model included the mainstem, four coves, three point-source discharges, and three withdrawals. Simulated water levels generally were within 10 centimeters of measured values, indicating a good calibration of the water balance for the reservoir. The root-mean-square difference between measured and simulated water temperatures was about 1 to 1.5 degrees Celsius, and vertical distributions of water temperature were accurately simulated in both the mainstem and coves. Seasonal and spatial patterns of nitrate, ammonia, orthophosphorus, and chlorophyll a were reasonably reproduced by the water-quality model. Because of the absence of the denitrification process in the model formulation, nitrate concentrations typically were overpredicted. Simulated and measured ammonia concentrations seldom differed by more than 0.01 milligram per liter, and simulations of seasonal fluctuations in chlorophyll a were representative of measured conditions. The root mean square of the difference between measured and simulated dissolved-oxygen concentrations was about 1 milligram per liter. The calibrated water-quality model was applied to evaluate (1) the movement of a conservative, neutrally buoyant material, or tracer, through the reservoir for several sets of conditions; (2) the effects of the Riverbend thermal discharge on water temperature in the reservoir; (3) the effects of changes in water-supply withdrawal rates on water-quality conditions; and (4) changes in reservoir water quality in response to changes in point- and nonpoint-source loadings. In general, dissolved material entering Mountain Island Lake from both Cowans Ford Dam and McDowell Creek during the summer moves along the bottom of the lake toward Mountain Island Dam, with little mixing of dissolved material into the surface layers. Simulations suggest that dissolved material can move upstream in the reservoir when flows from Cowans Ford Dam are near zero. Dissolved material can remain in Mountain Island Lake for a period far in excess of the theoretical retention time of 12 days. Simulations indicated that the Riverbend thermal discharge increases water temperature in the surface layers of the downstream part of the reservoir by as much as 5 degrees Celsius. However, the discharge has little effect on near-bottom water temperature. Based on model simulations, a proposed doubling of the water-supply withdrawals from Mountain Island Lake has no readily apparent effect on water quality in the reservoir. The increased withdrawal rate may have some localized effects on circulation in the reservoir, but a more detailed model of the intake zone would be required to identify those effects. The effects of a 20-percent increase in water-chemistry loadings through Cowans Ford Dam and from McDowell Creek were simulated separately. Increased loadings from Cowans Ford Dam had about the same effect on water-quality conditions near Mountain Island Dam as did increased loadings from McDowell Creek. Maintaining good water quality in Mountain Island Lake depends on maintaining good water quality in Lake Norman as well as in the inflows from the McDowell Creek watershed.

North Carolina