SEARCH · USGS Science

Results for “Rivers”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59Linked to original sources

USGS2009

Guide to Surficial Geology and River-Bluff Exposures, Noatak National Preserve, Northwestern Alaska

From its origin in rugged granitic highlands of the central Brooks Range, the Noatak River flows westward between the De Long Mountains and the Baird Mountains before turning south to enter Kotzebue Sound. Glaciers of middle and late Pleistocene age entered the Noatak River valley from the east, north, and south. Glaciers flowed down the upper Noatak River valley from the rugged peaks at its head, merging with tributary glaciers that issued from cirque-headed valleys along its south flank. Farther downvalley, small glaciers flowed northward from the Baird Mountains and much larger glaciers issued from the De Long Mountains. The De Long Mountains glaciers expanded southward to cover parts of the Noatak valley floor; they dammed the Noatak River during successive advances, creating a series of glacial lakes. The more extensive glacial advances dammed huge lakes that filled the Aniuk Lowland to overflowing. At various times, overflow waters spilled northward through Howard Pass, southward via Hunt River into the Kobuk River system, and westward down a series of channelways that skirted south of the glacier margins. Prominent bluffs along the Noatak River and its principal tributaries reveal glacial, glaciolacustrine, fluvial, and eolian sediments. More than 120 measured bluff exposures are described and illustrated in this report. These are dated by 92 radiocarbon age determinations and by the presence of the old Crow tephra, which was deposited about 130,000-140,000 years ago. Six geologic base maps, which cover sections of the Noatak River valley from east to west, show the locations of the river bluffs in relation to the glacial, glaciolacustrine, and fluvial deposits that cover the valley floor. The upper Noatak River valley is dominated by a bulky end moraine near Douglas Creek that was deposited during the last glacial maximum about 25,000-15,000 14C yr BP (termed the Itkillik II phase in the central Brooks Range glacial succession). Bluffs along this section of the Noatak River reveal thick till that underlies the moraine and interfingers downvalley with outwash and upvalley with moraine-dammed lake deposits. Remnants of older river gravels that underlie the set of glacial deposits contain wood fragments that are dated at about 35,000-30,000 14C yr BP. The Aniuk River area, which includes much of the eastern Aniuk Lowland, contains older moraines derived from headwaters of the Noatak valley that lie downvalley from the Douglas Creek moraine. These older moraines are assigned to the Itkillik IA and IB advances of the central Brooks Range glacial succession. Their deposits are seldom visible in river bluffs, but associated outwash and glaciolacustrine sediments are commonly exposed. More ancient end moraines farther downvalley are buried beneath lake deposits of the Aniuk Lowland, but are traceable as subdued arcuate drainage divides and as boulder concentrations in river bluffs or along their bases. The Cutler River area was occupied by glacial lakes assignable to three separate glacial phases. The oldest of these was probably dammed by the Cutler moraine, which crosses the Noatak valley floor near the mouth of Cutler River. The younger two are correlated with Itkillik-age deposits in the Aniuk Lowland. Glaciers of Itkillik age also flowed northwestward down the Cutler and Imelyak valley systems from cirques along the Noatak-Kobuk divide, but they did not reach the Noatak valley floor. The western Aniuk Lowland, which extends westward from the Cutler River mouth to the lower course of Nimiuktuk River, is dominated by a series of large end moraines deposited by glaciers from the De Long Mountains and that flowed southeastward down the Nimiuktuk valley system and then up the Noatak River valley. The Cutler moraine is the most extensive of these deposits. Following the Cutler glaciation, less extensive glacial advances built end moraines near the present-day mouths of Makpik Creek and Anisak River during in

Scientific Investigations Report↗

USGS2013

Groundwater and surface-water interaction within the upper Smith River Watershed, Montana 2006-2010

The 125-mile long Smith River, a tributary of the Missouri River, is highly valued as an agricultural resource and for its many recreational uses. During a drought starting in about 1999, streamflow was insufficient to meet all of the irrigation demands, much less maintain streamflow needed for boating and viable fish habitat. In 2006, the U.S. Geological Survey, in cooperation with the Meagher County Conservation District, initiated a multi-year hydrologic investigation of the Smith River watershed. This investigation was designed to increase understanding of the water resources of the upper Smith River watershed and develop a detailed description of groundwater and surface-water interactions. A combination of methods, including miscellaneous and continuous groundwater-level, stream-stage, water-temperature, and streamflow monitoring was used to assess the hydrologic system and the spatial and temporal variability of groundwater and surface-water interactions. Collectively, data are in agreement and show: (1) the hydraulic connectedness of groundwater and surface water, (2) the presence of both losing and gaining stream reaches, (3) dynamic changes in direction and magnitude of water flow between the stream and groundwater with time, (4) the effects of local flood irrigation on groundwater levels and gradients in the watershed, and (5) evidence and timing of irrigation return flows to area streams. Groundwater flow within the alluvium and older (Tertiary) basin-fill sediments generally followed land-surface topography from the uplands to the axis of alluvial valleys of the Smith River and its tributaries. Groundwater levels were typically highest in the monitoring wells located within and adjacent to streams in late spring or early summer, likely affected by recharge from snowmelt and local precipitation, leakage from losing streams and canals, and recharge from local flood irrigation. The effects of flood irrigation resulted in increased hydraulic gradients (increased groundwater levels relative to stream stage) or even reversed gradient direction at several monitoring sites coincident with the onset of nearby flood irrigation. Groundwater-level declines in mid-summer were due to groundwater withdrawals and reduced recharge from decreased precipitation, increased evapotranspiration, and reduced leakage in some area streams during periods of low flow. Groundwater levels typically rebounded in late summer, a result of decreased evapotranspiration, decreased groundwater use for irrigation, increased flow in losing streams, and the onset of late-season flood irrigation at some sites. The effect of groundwater and surface-water interactions is most apparent along the North and South Forks of the Smith River where the magnitude of streamflow losses and gains can be greater than the magnitude of flow within the stream. Net gains consistently occurred over the lower 15 miles of the South Fork Smith River. A monitoring site near the mouth of the South Fork Smith River gained (flow from the groundwater to the stream) during all seasons, with head gradients towards the stream. Two upstream sites on the South Fork Smith River exhibited variable conditions that ranged from gaining during the spring, losing (flowing from the stream to the groundwater) during most of the summer as groundwater levels declined, and then approached or returned to gaining conditions in late summer. Parts of the South Fork Smith River became dry during periods of losing conditions, thus classifying this tributary as intermittent. The North Fork Smith River is highly managed at times through reservoir releases. The North Fork Smith River was perennial throughout the study period although irrigation diversions removed a large percentage of streamflow at times and losing conditions persisted along a lower reach. The lowermost reach of the North Fork Smith River near its mouth transitioned from a losing reach to a gaining reach throughout the study period. Groundwater and surface-water interactions occur downstream from the confluence of the North and South Fork Smith Rivers, but are less discernible compared to the overall magnitude of the main-stem streamflow. The Smith River was perennial throughout the study. Monitoring sites along the Smith River generally displayed small head gradients between the stream and the groundwater, while one site consistently showed strongly gaining conditions. Synoptic streamflow measurements during periods of limited irrigation diversion in 2007 and 2008 consistently showed gains over the upper 41.4 river miles of the main stem Smith River where net gains ranged from 13.0 to 28.9 cubic feet per second. Continuous streamflow data indicated net groundwater discharge and small-scale tributary inflow contributions of around 25 cubic feet per second along the upper 10-mile reach of the Smith River for most of the 2010 record. A period of intense irrigation withdrawal during the last two weeks in May was followed by a period (early June 2010 to mid-July 2010) with the largest net increase (an average of 71.1 cubic feet per second) in streamflow along this reach of the Smith River. This observation is likely due to increased groundwater discharge to the Smith River resulting from irrigation return flow. By late July, the apparent effects of return flows receded, and the net increase in streamflow returned to about 25 cubic feet per second. Two-dimensional heat and solute transport VS2DH models representing selected stream cross sections were used to constrain the hydraulic properties of the Quaternary alluvium and estimate temporal water-flux values through model boundaries. Hydraulic conductivity of the Quaternary alluvium of the modeled sections ranged from 3x10-6 to 4x10-5 feet per second. The models showed reasonable approximations of the streambed and shallow aquifer environment, and the dynamic changes in water flux between the stream and the groundwater through different model boundaries.

Montana↗

USGS2014

Continuous water-quality monitoring and regression analysis to estimate constituent concentrations and loads in the Red River of the North at Fargo and Grand Forks, North Dakota, 2003-12

The Red River of the North (hereafter referred to as “Red River”) Basin is an important hydrologic region where water is a valuable resource for the region’s economy. Continuous water-quality monitors have been operated by the U.S. Geological Survey, in cooperation with the North Dakota Department of Health, Minnesota Pollution Control Agency, City of Fargo, City of Moorhead, City of Grand Forks, and City of East Grand Forks at the Red River at Fargo, North Dakota, from 2003 through 2012 and at Grand Forks, N.Dak., from 2007 through 2012. The purpose of the monitoring was to provide a better understanding of the water-quality dynamics of the Red River and provide a way to track changes in water quality. Regression equations were developed that can be used to estimate concentrations and loads for dissolved solids, sulfate, chloride, nitrate plus nitrite, total phosphorus, and suspended sediment using explanatory variables such as streamflow, specific conductance, and turbidity. Specific conductance was determined to be a significant explanatory variable for estimating dissolved solids concentrations at the Red River at Fargo and Grand Forks. The regression equations provided good relations between dissolved solid concentrations and specific conductance for the Red River at Fargo and at Grand Forks, with adjusted coefficients of determination of 0.99 and 0.98, respectively. Specific conductance, log-transformed streamflow, and a seasonal component were statistically significant explanatory variables for estimating sulfate in the Red River at Fargo and Grand Forks. Regression equations provided good relations between sulfate concentrations and the explanatory variables, with adjusted coefficients of determination of 0.94 and 0.89, respectively. For the Red River at Fargo and Grand Forks, specific conductance, streamflow, and a seasonal component were statistically significant explanatory variables for estimating chloride. For the Red River at Grand Forks, a time component also was a statistically significant explanatory variable for estimating chloride. The regression equations for chloride at the Red River at Fargo provided a fair relation between chloride concentrations and the explanatory variables, with an adjusted coefficient of determination of 0.66 and the equation for the Red River at Grand Forks provided a relatively good relation between chloride concentrations and the explanatory variables, with an adjusted coefficient of determination of 0.77. Turbidity and streamflow were statistically significant explanatory variables for estimating nitrate plus nitrite concentrations at the Red River at Fargo and turbidity was the only statistically significant explanatory variable for estimating nitrate plus nitrite concentrations at Grand Forks. The regression equation for the Red River at Fargo provided a relatively poor relation between nitrate plus nitrite concentrations, turbidity, and streamflow, with an adjusted coefficient of determination of 0.46. The regression equation for the Red River at Grand Forks provided a fair relation between nitrate plus nitrite concentrations and turbidity, with an adjusted coefficient of determination of 0.73. Some of the variability that was not explained by the equations might be attributed to different sources contributing nitrates to the stream at different times. Turbidity, streamflow, and a seasonal component were statistically significant explanatory variables for estimating total phosphorus at the Red River at Fargo and Grand Forks. The regression equation for the Red River at Fargo provided a relatively fair relation between total phosphorus concentrations, turbidity, streamflow, and season, with an adjusted coefficient of determination of 0.74. The regression equation for the Red River at Grand Forks provided a good relation between total phosphorus concentrations, turbidity, streamflow, and season, with an adjusted coefficient of determination of 0.87. For the Red River at Fargo, turbidity and streamflow were statistically significant explanatory variables for estimating suspended-sediment concentrations. For the Red River at Grand Forks, turbidity was the only statistically significant explanatory variable for estimating suspended-sediment concentration. The regression equation at the Red River at Fargo provided a good relation between suspended-sediment concentration, turbidity, and streamflow, with an adjusted coefficient of determination of 0.95. The regression equation for the Red River at Grand Forks provided a good relation between suspended-sediment concentration and turbidity, with an adjusted coefficient of determination of 0.96.

North Dakota↗

USGS2016

Flow characteristics and salinity patterns of tidal rivers within the northern Ten Thousand Islands, southwest Florida, water years 2007–14

Freshwater flow to the Ten Thousand Islands estuary has been altered by the construction of the Tamiami Trail and the Southern Golden Gate Estates. The Picayune Strand Restoration Project, which is associated with the Comprehensive Everglades Restoration Plan, has been implemented to improve freshwater delivery to the Ten Thousand Islands estuary by removing hundreds of miles of roads, emplacing hundreds of canal plugs, removing exotic vegetation, and constructing three pump stations. Quantifying the tributary flows and salinity patterns prior to, during, and after the restoration is essential to assessing the effectiveness of upstream restoration efforts. Tributary flow and salinity patterns during preliminary restoration efforts and prior to the installation of pump stations were analyzed to provide baseline data and preliminary analysis of changes due to restoration efforts. The study assessed streamflow and salinity data for water years 1 2007–2014 for the Faka Union River (canal flow included), East River, Little Wood River, Pumpkin River, and Blackwater River. Salinity data from the Palm River and Faka Union Boundary water-quality stations were also assessed. Faka Union River was the dominant contributor of freshwater during water years 2007–14 to the Ten Thousand Islands estuary, followed by Little Wood and East Rivers. Pumpkin River and Blackwater River were the least substantial contributors of freshwater flow. The lowest annual flow volumes, the highest annual mean salinities, and the highest percentage of salinity values greater than 35 parts per thousand (ppt) occurred in water year 2011 at all sites with available data, corresponding with the lowest annual rainfall during the study. The highest annual flow volumes and the lowest percentage of salinities greater than 35 ppt occurred in water year 2013 for all sites with available data, corresponding with the highest rainfall during the study. In water year 2014, the percentage of monitored annual flow contributed by East River increased and the percentage of flow contributed by Faka Union River decreased, compared to the earlier years. No changes in annual flow occurred at any sites west of Faka Union River. No changes in the relative flow contributions were observed during the wet season; however, the relative amounts of streamflow increased during the dry season at East River in 2014. East River had only 1 month of negative flow in 2014 compared to 6 months in 2011 and 7 months in 2008. Higher dry season flows in East River may be in response to restoration efforts. The sites to the west of Faka Union River had higher salinities on average than Faka Union River and East River. Faka Union River had the highest range in salinities, and Faka Union Boundary had the lowest range in salinities. Pumpkin River was the tributary with the lowest range in salinities. 1 Water year is defined as the 12-month period from October 1, for any given year, through September 30 of the following year.

Florida↗

USGS2018

Preliminary synthesis and assessment of environmental flows in the middle Verde River watershed, Arizona

A 3-year study was undertaken to evaluate the suitability of the available modeling tools for characterizing environmental flows in the middle Verde River watershed of central Arizona, describe riparian vegetation throughout the watershed, and estimate sediment mobilization in the river. Existing data on fish and macroinvertebrates were analyzed in relation to basin characteristics, flow regimes, and microhabitat, and a pilot study was conducted that sampled fish and macroinvertebrates and the microhabitats in which they were found. The sampling for the pilot study took place at five different locations in the middle Verde River watershed. This report presents the results of this 3-year study. The Northern Arizona Groundwater Flow Model (NARGFM) was found to be capable of predicting long-term changes caused by alteration of regional recharge (such as may result from climate variability) and groundwater pumping in gaining, losing, and dry reaches of the major streams in the middle Verde River watershed. Over the period 1910 to 2006, the model simulated an increase in dry reaches, a small increase in reaches losing discharge to the groundwater aquifer, and a concurrent decrease in reaches gaining discharge from groundwater. Although evaluations of the suitability of using the NARGFM and Basin Characteristic Model to characterize various streamflow intervals showed that smallerscale basin monthly runoff could be estimated adequately at locations of interest, monthly stream-flow estimates were found unsatisfactory for determining environmental flows. Orthoimagery and Moderate Resolution Imaging Spectroradiometer data were used to quantify stream and riparian vegetation properties related to biotic habitat. The relative abundance of riparian vegetation varied along the main channel of the Verde River. As would be expected, more upland plant species and fewer lowland species were found in the upper-middle section compared to the lower-middle section, and vice-versa. Vegetation changes within the upper-middle and lower-middle reaches are related to differences in climate and hydrology. In general, the riparian vegetation of the middle Verde River watershed is that of a healthy ecosystem’s mixed age, mixed patch structure, likely a result of the mostly unaltered disturbance regime. The frequency of in-river hydrogeomorphic features (pool, riffle, run) varied along the middle Verde River channel. There was a greater abundance of riffle habitat in the upper-middle reach; the lower-middle reach included more pool habitat. The Oak Creek tributary was more homogenous in geomorphic stream habitat composition than West Clear Creek, where runs dominated the upper reaches and pools dominated many of the lower reaches. On the basis of the period of record and discharges recorded at 15-minute intervals, five flows were found to reach the gravel-transport threshold. Sediment mobilization computed with flows averaged over daily time steps yielded just three flows that reached the gravel-transport threshold, and monthly averaged flows yielded none. In the middle Verde River watershed, 15-minute data should be used when possible to evaluate sediment transport in the river system. Data from more than 300 fish surveys conducted from 1992 to 2011 were analyzed using two schemes, one that divided the river into five reaches based on basin characteristics, and a second that divided the river into five reaches based on degree of flow alteration (specifically, diversions). Fish community metrics and assemblage data were used to analyze patterns of species composition and abundance in the two approaches. Overall, native and non-native species were regularly interacting and probably competing for similar resources. Fish abundances were also analyzed in response to floods and other flow metrics. Although the data are limited, native fish abundances increased more rapidly than non-native fish abundances in response to large floods. The basin-characteristic reach analysis showed native fish in greater abundance in the upper-middle reaches of the Verde River watershed and generally decreasing with downstream distance. The median relative abundance of native fish decreased by 50 percent from reach 1 to reach 5. Using the reach scheme based on degree of flow alteration, nondiverted reaches were found to have a greater abundance of native fish than diverted reaches. In heavily diverted reaches, non-native species outnumbered native species. Fish metrics and stream-flow metrics for the 30, 90, and 365-day periods before collection were computed and the results analyzed statistically. Only abundance of all fish species was associated with the 30-day flow metrics. The 90-day flow metrics were generally positively associated with fish metrics, whereas the 365-day flow metrics had more negative correlations. In particular, significant relations were found between fish metrics and the magnitude and frequency of high flows, including maximum monthly flow, median annual number of high-flow events, and median annual maximum streamflow. Native sucker (Catostomidae) populations tended to decrease in periods of extended base flow, and fish in the non-native sunfish family (Centrarchidae) decreased in periods of flashy, high magnitude flows. A pilot study surveyed fish at five locations in the upper part of the middle Verde River watershed as a means to measure microhabitat availability and quantify native and non-native fish use of that available microhabitat. Results indicated that native and non-native species exhibit some clear differences in microhabitat use. Although at least some native and non-native fish were found in each velocity, depth, and substrate category, preferential microhabitat use was common. On a percentage basis, non-native species had a strong preference for slow-moving and deeper water with silt and sand substrate, with a secondary preference for faster moving and very shallow water and a coarse gravel substrate. Native species showed a general preference for somewhat faster, moderate depth water over coarse gravel and had no clear secondary preference. Macroinvertebrate-variables index period, high-flow year, and collection location (upper-middle Verde River, lowermiddle Verde River, or Verde River tributaries) were found to be important explanatory variables in differentiating among community metrics. Overall richness (number of unique taxa), Shannon’s diversity index, and the percent of the most dominant taxa were all highly correlated, but their response to each macroinvertebrate variable was different. The percentage of mayfly (order Ephemeroptera) taxa was significantly higher in Oak Creek and the upper-middle and lower-middle Verde River reaches, locations which have higher flows and more urbanization than other reaches. When community metrics were related to hydrologic metrics, caddisfly (order Trichoptera) populations appeared to increase and mayfly populations to decrease in response to less flashy and more stable streamflows. Conversely, caddisfly populations appeared to decrease and mayfly populations to increase in response to greater flow variability. Six locations along the Verde River were sampled for macroinvertebrates as part of a pilot study associated with this report—(1) below Granite Creek, (2) near Campbell Ranch, (3) at the U.S. Geological Survey Paulden gage, (4) at the Perkinsville Bridge, (5) at the USGS Clarkdale gage, and (6) near the Reitz Ranch property. A nonmetric multidimensional scaling ordination of macroinvertebrate assemblages showed that the Verde River below Granite Creek site was different from the five other sites and that the Perkinsville Bridge and near Reitz Ranch samples had similar community structure. The near Campbell Ranch and Paulden gage locations had similar microhabitat characteristics, with the exception of riparian cover, yet the assemblage structure was very different. The different community composition at Verde River below Granite Creek was likely due to it having the smallest substrate sizes, lowest velocities, shallowest depths, and most riparian cover of the six sites.

Arizona↗

USGS2001

Simulation of ground-water flow in the Mojave River basin, California

The proximity of the Mojave River ground-water basin to the highly urbanized Los Angeles region has led to rapid growth in population and, consequently, to an increase in the demand for water. The Mojave River, the primary source of surface water for the region, normally is dry-except for a small stretch of perennial flow and periods of flow after intense storms. Thus, the region relies almost entirely on ground water to meet its agricultural and municipal needs. Ground-water withdrawal since the late 1800's has resulted in discharge, primarily from pumping wells, that exceeds natural recharge. To better understand the relation between the regional and the floodplain aquifer systems and to develop a management tool that could be used to estimate the effects that future stresses may have on the ground-water system, a numerical ground-water flow model of the Mojave River ground-water basin was developed, in part, on the basis of a previously developed analog model. The ground-water flow model has two horizontal layers; the top layer (layer 1) corresponds to the floodplain aquifer and the bottom layer (layer 2) corresponds to the regional aquifer. There are 161 rows and 200 columns with a horizontal grid spacing of 2,000 by 2,000 feet. Two stress periods (wet and dry) per year are used where the duration of each stress period is a function of the occurrence, quantity of discharge, and length of stormflow from the headwaters each year. A steady-state model provided initial conditions for the transient-state simulation. The model was calibrated to transient-state conditions (1931-94) using a trial-and-error approach. The transient-state simulation results are in good agreement with measured data. Under transient-state conditions, the simulated floodplain aquifer and regional aquifer hydrographs matched the general trends observed for the measured water levels. The simulated streamflow hydrographs matched wet stress period average flow rates and times of no flow at the Barstow and Afton Canyon gages. Steady-state particle-tracking was used to estimate travel times for mountain-front and streamflow recharge. The simulated travel times for mountain-front recharge to reach the area west of Victorville were about 5,000 to 6,000 years; this result is in reasonable agreement with published results. Steady-state particle-tracking results for streamflow recharge indicate that in most subareas along the river, the particles quickly leave and reenter the river. The complaint that resulted in the adjudication of the Mojave River ground-water basin alleged that the cumulative water production upstream of the city of Barstow had overdrafted the ground-water basin. In order to ascertain the effect of pumping on ground-water and surface-water relations along the Mojave River, two pumping simulations were compared with the 1931-90 transient-state simulation (base case). The first simulation assumed 1931-90 pumping in the upper region (Este, Oeste, Alto, and Transition zone model subareas) but with no pumping in the remainder of the basin, and the second assumed 1931-90 pumping in the lower region (Centro, Harper Lake, Baja, Coyote Lake, and Afton Canyon model subareas) but with no pumping in remainder of the basin. In the upper region, assuming pumping only in the upper region, there was no change in storage, recharge from the Mojave River, ground-water discharge to the Mojave River, or evapotranspiration when compared with the base case. In the lower region, assuming pumping only in the upper region, there was storage accretion, decreased recharge from the Mojave River, increased ground-water discharge to the Mojave River, and increased evapotranspiration when compared with the base case. In the upper region, assuming pumping only in the lower region, there was storage accretion, decreased recharge from the Mojave River, increased ground-water discharge to the Mojave River, and increased evapotranspiration when compared with the base case. In the lower region, assuming pumping only in the lower region, there was less storage depletion, increased recharge from the Mojave River, increased ground-water discharge to the Mojave River, and increased evapotranspiration when compared with the base case. Overall, pumping in the lower region does not negatively affect the upper region; however, pumping in the upper region negatively affects the lower region by decreasing recharge from the Mojave River. Streamflow, pumpage, and water-level data from calendar years 1995-99 were used to validate the calibrated ground-water flow model, that is, to test that the ground-water flow model will duplicate measured data for a noncalibration period without modification of the model parameters. In general, the simulated results are in good agreement with the measured data, and the simulated hydrographs for wells in the floodplain and regional aquifers follow the measured water-level trends. Simulated streamflow data for the 1995-99 wet and dry stress periods at the Lower Narrows, Barstow, and Afton Canyon were compared with the measured data for average streamflow for the same periods; in general, the model reflects 1995-99 streamflow conditions. The simulation results also indicate that the streambed conductance values calibrated to the 1931-94 conditions reasonably simulate the 1995-99 conditions and therefore can be used for predictive purposes. To visualize the magnitude, spatial distribution, and timing of water-level changes in the basin through time, simulated hydraulic heads for 1932-99 were compared with simulated hydraulic heads for 1931. Greater than average annual inflows to the Mojave River from the headwaters during the late 1930's and throughout much of the 1940's resulted in simulated hydraulic heads that were higher than the 1931 hydraulic heads along the Mojave River in most model subareas. Parts of the Baja and Harper Lake model subareas had declines in the simulated hydraulic head because of the increase in agricultural pumpage. By 1960, the simulated hydraulic heads were lower than the simulated hydraulic heads for 1931 in all model subareas of the floodplain and the regional aquifers because of pumpage. After 1960, the size and the magnitude of the areas of the regional aquifer for which simulated hydraulic heads were lower than those for 1931 continued to increase until the end of the simulation (1999). Along the Mojave River, hydraulic heads fluctuated in the floodplain aquifer in response to recharge during years with large inflows with little apparent effect on the simulated hydraulic heads in the regional aquifer. Three water-management alternatives were evaluated to determine their effect on ground-water resources using the calibrated ground-water flow model. The water-management alternatives consider the artificial recharge of imported water allocated to the Mojave Water Agency (MWA): the first assumes that zero percent of the MWA allocation is available (alternative 1), the second assumes that 50 percent of the MWA allocation is available (alternative 2), and the third assumes that 100 percent of the MWA allocation is available (alternative 3). Each of the three water-management alternatives were evaluated for a 20-year drought. Streamflow conditions were simulated using the 20-year drought of 1945-64 with associated calibrated stream parameters. Management alternative 1 results in a reduction in ground-water recharge from the Mojave River compared with average recharge for 1995-99; this reduction is reflected in simulated hydraulic-head declines between 1999 and 2019 of as much as 45 feet. Management alternatives 2 and 3 result in no change in recharge from the Mojave River for management alternative 2 and a small increase for management alternative 3 when compared with recharge for management alternative 1. The artificial recharge of imported water causes increases in simulated hydraulic head for both management alternatives at each of the artificial-recharge sites. Some of the increases are related to water that recharges into areas of low transmissivity which implies that the recharge operations may benefit from being distributed over a larger area.

California↗

USGS2000

Flow and salinity characteristics of the upper Suwannee River Estuary, Florida

Continuous stage and salinity data were recorded from August 1995 to December 1997 at four gages located in the upper Suwannee River Estuary. Continuous velocity data were recorded at two of the four gages and continuous discharge data were computed for these two gages. Additional salinity data were collected at 15 monitoring sites from November 1992 to October 1997. Wind-speed data collected at Cedar Key, Florida, during the study period were utilized in the regression analysis. Correlations were developed to describe the longitudinal extent of the saltwater/freshwater interface (defined as 0.5 parts per thousand (ppt) salinity) and salinity distribution in the upper Suwannee River Estuary. On East Pass, the median of difference between daily maximum and daily minimum stage ranged from 2.92 feet for a gage at river mile 3.8 to 3.33 feet for a gage at river mile 1.2. Velocities tended to be unidirectional with some instances of bilateral flow. Reversal in flow direction was common and coincided with rising tides. Monthly mean discharges for the Suwannee River near Wilcox, Florida, during the study period typically were lower than the average for the period of record (1931.97). Discharge near Wilcox averaged 4,000 cubic feet per second (ft3/s) lower than the long-term average from June to September 1996. An El Niï¿½o event induced precipitation that was responsible for higher than average monthly mean discharge measured near Wilcox during November and December 1997. The maximum observed salinity concentrations for the study period ranged from 28.20 ppt at river mile 3.8 to 31.00 ppt at river mile 1.9. Median daily fluctuations of salinity at river miles 3.8 and 1.2 were 0.12 and 11.31 ppt, respectively. The maximum daily upstream extent of the saltwater/freshwater interface was at or upstream from river mile 4.0 for about 50 percent of the study period. The interface was at or upstream from river mile 3.8 and river mile 2.8 40 and 57 percent of the time. The interface was downstream from river mile 1.2 and river mile 1.9 11 and 21 percent of the time, respectively. The median daily maximum salinity for the four gages ranged from 0.22 ppt at river mile 3.8 to 11.50 ppt at river mile 1.2. Multiple linear-regression models were developed to determine the isohaline location for 0.5, 2, 5, 10, 15, and 20 ppt salinity, and to predict the maximum daily salinity concentrations at gages as a function of stage, river discharge, and wind. The salinity at a location was inversely proportional to the daily mean discharge at the Suwannee River near Wilcox. Under extreme low-flow conditions (3,500 ft3/s), the regression models predicted that the interface would occur at river mile 7.2, upstream from the Gopher River confluence with the Suwannee River. Wind speed did not have a substantial influence on model predictions. The period of record for the Suwannee River at Wilcox was applied to appropriate regression models to produce a synthetic record of historical salinity distributions. Two withdrawal scenarios, a 10-percent diversion and a 1,000 ft3/s diversion, were evaluated relative to high-, medium-, and low-flow conditions and compared to actual salinity distributions. The 10-percent and 1,000 ft3/s withdrawals scenarios resulted in the isohaline of 0.5 ppt migrating 0.6 and 1.58 miles upstream from the actual isohaline location for a low-flow condition of 4,500 ft3/s, and migrating 0.14 and 0.65 miles upstream from the actual isohaline location for a high-flow conditions of 20,300 ft3/s for Wadley Pass.

Florida↗

USGS2004

Pesticide concentrations in water and in suspended and bottom sediments in the New and Alamo rivers, Salton Sea Watershed, California, April 2003

This report contains pesticide concentration data for water, and suspended and bed sediment samples collected in April 2003 from twelve sites along the New and Alamo Rivers in the Salton Sea watershed, in southeastern California. The study was done in collaboration with the California State Regional Water Quality Control Board, Colorado River Region, to assess inputs of current-use pesticides associated with water and sediment into the New and Alamo Rivers. Five sites along the New River and seven sites along the Alamo River, downstream of major agricultural drains, were selected and covered the lengths of the rivers from the international boundary to approximately 1.5 km from the river mouths. Sampling from bridges occurred at seven of the twelve sites. At these sites, streamflow measurements were taken. These same sites were also characterized for cross-stream homogeneity by measuring dissolved oxygen, pH, specific conductance, temperature, and suspended solids concentration at several vertical (depths) and horizontal (cross-stream) points across the river. Large volume water samples (200?300 L) were collected for isolation of suspended sediments by flow-through centrifugation. Water from the outflow of the flow-through centrifuge was sampled for the determination of aqueous pesticide concentrations. In addition, bottom sediments were sampled at each site. Current-use pesticides and legacy organochlorine compounds (p,p'-DDT, p,p'-DDE and p,p'-DDD) were extracted from sediments and measured via gas chromatography/mass spectrometry (GC/MS). Organic carbon and percentage of fines were also determined for suspended and bottom sediments. Cross-stream transects of dissolved constituents and suspended sediments showed that the rivers were fairly homogeneous at the sites sampled. Streamflow was higher at the outlet sites, with the Alamo River having higher flow (1,240 cfs) than the New River (798 cfs). Twelve current-use pesticides, one legacy organochlorine compound (p,p'-DDE), and the additive piperonyl butoxide were detected in water samples. Trifluralin was found in the highest concentration of all detected compounds (68.5?599 ng/L) at all sites in both rivers, except for the international boundary sites. Atrazine was also detected in high concentration (51.0?285 ng/L) at several sites. The outlet sites had among the highest numbers of pesticides detected and the international boundary sites had the lowest numbers of pesticides detected for both rivers. The numbers of pesticides detected were greater for the Alamo River than for the New River. Six current-use pesticides and two legacy organochlorines (p,p'-DDE and p,p'-DDD) were found associated with suspended and bed sediments. The DDT metabolite p,p'-DDE was detected in all suspended and bed sediments from the Alamo River, but only at two sites in the New River. Dacthal, chlorpyrifos, pendimethalin, and trifluralin were the most commonly detected current-use pesticides. Trifluralin was the compound found in the highest concentrations in suspended (14.5?120 ng/g) and bed (1.9?9.0 ng/g) sediments. The sites along the Alamo River had more frequent detections of pesticides in suspended and bed sediments when compared with the New River sites. The greatest number of pesticides that were detected in suspended sediments (seven) were in the samples from the Sinclair Road and Harris Road sites. For bottom sediments, the Alamo River outlet site had the greatest number of pesticide detections (eight).

Data Series↗

USGS2013

Behavior and movement of adult summer steelhead following collection and release, lower Cowlitz River, Washington, 2012--2013

Executive Summary Historically, adult summer steelhead Oncorhynchus mykiss returning to hatcheries on the lower Cowlitz River were sometimes transported and released in the river (recycled) to provide additional angling opportunity for the popular sport fishery in the basin. However, this practice has not been used in recent years because of concerns associated with interactions between hatchery fish and wild fish. Fishery managers were interested in resuming recycling but lacked information regarding effects of this practice on wild steelhead so we conducted a study during 2012–2013 to: (1) enumerate recycled steelhead that returned to the hatchery or were removed by anglers; and (2) determine if steelhead that were not removed from the river remained in the system where they could interact with wild fish. During June–August 2012, a total of 549 summer steelhead were captured at the Cowlitz Salmon Hatchery, tagged, and released downstream near the Interstate 5 Bridge. All recycled steelhead were tagged with a white Floy ® tag and opercle-punched; 109 (20 percent) of these fish also were radio-tagged. All adult steelhead that return to the hatchery were handled by hatchery staff so recycled steelhead that returned to the hatchery were enumerated daily. A creel survey and voluntary angler reports were used to determine the number of recycled steelhead that were caught by anglers. We established three fixed telemetry monitoring sites on the mainstem Cowlitz River and eight additional sites were deployed on tributaries to the lower Cowlitz River where wild winter steelhead are known to spawn. We also conducted mobile tracking from a boat during October 2012, November 2012, and January 2013 to locate radio-tagged fish. A total of 10,722 summer steelhead were captured at the Cowlitz Salmon Hatchery in 2012, which was the largest return since 2008. River flows during much of the study period were similar to 2008–2011 average flows, however, high-flow periods in July and November 2012 were nearly twice as high as the 2008–2011 average flows. We determined that 50 percent (273 fish) of the recycled steelhead returned to the hatchery and 18 percent (102 fish) of the recycled steelhead were caught by anglers. Most (243 fish; 89 percent) of the recycled steelhead that returned to the hatchery were recollected during July–August. The average elapsed time from release to recapture at the hatchery was 9 days (d) and 72 percent (182 fish) of the fish returned to the hatchery within 14 d of release. These trends were similar for recycled steelhead that were caught by anglers. Most fish were caught during July–August and the median time from release to capture was 10 d. We determined that 65 percent (70 fish) of the angler-caught fish returned to the hatchery within 14 d of release. River flows appeared to affect both hatchery returns and angler catch. The daily number of recycled steelhead that were recollected at the hatchery were low during periods when river flows were decreasing and high during periods when river flows were increasing. Conversely, daily angler catch of recycled steelhead generally was low when flows were increasing and high when flows were decreasing. We determined that 32 percent of the recycled steelhead (174 fish) were not removed from the lower Cowlitz River, based on observations from hatchery returns and angler reports, but results from the radio-tagged fish were insightful for understanding what may have happened to these fish. By comparison, we determined that 24 percent of the radio-tagged fish were not known to have been removed from the river. We determined that 12 percent of these fish were actively moving in the lower Cowlitz River during October 2012–January 2013. None of the radio-tagged fish were detected in tributaries during the study period except for a single fish that spent approximately 7 d in the Toutle River during early September. During October 2012–January 2013, 10 percent of the radio-tags from recycled steelhead were detected near popular fishing areas, and 2 percent of the radio-tagged steelhead were never detected during the study period. We suspect that a large proportion of these fish may have been harvested and not reported, or died. Detection patterns of radio-tagged steelhead showed that most fish (82 percent) moved upstream from the release site and were detected at the Trout Hatchery and the Barrier Dam sites. The median time from release to detection at these sites was 3.7 d and many of these fish made multiple trips between the two sites. Nearly one-third (29 percent) of the recycled steelhead that were detected at the Trout Hatchery and the Barrier Dam made at least two trips between the sites and some fish made as many as six trips. Radio-tagged fish that remained in the lower Cowlitz River during the spawning period (December 2012–January 2013) were observed in the river reach between the mouth of Ostrander Creek (river mile 10) and the Trout Hatchery (river mile 44). During this study, we collected data on opercle punch regrowth rates to understand the temporal effectiveness of this marking technique. We took opercle measurements for a total of 190 fish during the study. Fresh opercle punches were measured for 63 fish at the time of marking, and the remaining 127 fish were measured when fish returned to the hatchery. We determined that opercle punches remained open for about 30 d. The holes appeared to regrow slowly in the first 20 d after marking, but regrowth accelerated during the 20–30 d post-marking period. After 30 d, all opercle punches that we observed had completely closed due to tissue regrowth. Our study showed that a large proportion (68 percent) of the recycled steelhead were removed from the lower Cowlitz River. These fish primarily entered the hatchery or were caught by anglers within 14 d of release, which suggests that they present minimal risk to wild fish in the system. However, the remaining fish (32 percent) could not be accounted for, which may complicate fisheries management decisions associated with recycling summer steelhead. Findings from the radiotelemetry study suggest that unreported harvest or mortality could explain a large proportion of those fish that were not reported as having been removed from the river. Furthermore, intensive monitoring of the key spawning tributaries failed to detect a single fish during the spawning period. These findings were supported by observations from weir traps operated by the Washington Department of Fish and Wildlife. Our findings indicate that additional research may be warranted to further examine the effects of recycling hatchery summer steelhead in the lower Cowlitz River.

Washington↗

USGS2016

Behavior patterns and fates of adult steelhead, Chinook salmon, and coho salmon released into the upper Cowlitz River Basin, 2005–09 and 2012, Washington

A multiyear radiotelemetry evaluation was conducted to monitor adult steelhead ( Oncorhynchus mykiss ), Chinook salmon ( O. tshawytscha ), and coho salmon ( O. kisutch ) behavior and movement patterns in the upper Cowlitz River Basin. Volitional passage to this area was eliminated by dam construction in the mid-1960s, and a reintroduction program began in the mid-1990s. Fish are transported around the dams using a trap-and-haul program, and adult release sites are located in Lake Scanewa, the uppermost reservoir in the system, and in the Cowlitz and Cispus Rivers. Our goal was to estimate the proportion of tagged fish that fell back downstream of Cowlitz Falls Dam before the spawning period and to determine the proportion that were present in the Cowlitz and Cispus Rivers during the spawning period. Fallback is important because Cowlitz Falls Dam does not have upstream fish passage, so fish that pass the dam are unable to move back upstream and spawn. A total of 2,051 steelhead and salmon were tagged for the study, which was conducted during 2005–09 and 2012, and 173 (8.4 percent) of these regurgitated their transmitter prior to, or shortly after release. Once these fish were removed from the dataset, the final number of fish that was monitored totaled 1,878 fish, including 647 steelhead, 770 Chinook salmon, and 461 coho salmon. Hatchery-origin (HOR) and natural-origin (NOR) steelhead, Chinook salmon, and coho salmon behaved differently following release into Lake Scanewa. Detection records showed that the percentage of HOR fish that moved upstream and entered the Cowlitz River or Cispus River after release was relatively low (steelhead = 38 percent; Chinook salmon = 67 percent; coho salmon = 41 percent) compared to NOR fish (steelhead = 84 percent; Chinook salmon = 82 percent; coho salmon = 76 percent). The elapsed time from release to river entry was significantly lower for NOR fish than for HOR fish for all three species. Tagged fish entered the Cowlitz River in greater proportions than the Cispus River, regardless of origin. We found that 23–47 percent of the HOR fish entered the Cowlitz River and 12–38 percent entered the Cispus River. Similarly, 67–70 percent of the NOR fish entered the Cowlitz River and 38–66 percent entered the Cispus River. These behavioral differences translated into similar differences in fates during the spawning periods as higher percentages of tagged fish were assigned Cowlitz River fates than Cispus River fates. Fallback rates were affected by fish origin and release site. Overall, 12 percent of steelhead, 19 percent of Chinook salmon, and 8 percent of coho salmon fell back downstream of Cowlitz Falls Dam prior to spawning. Fallback rates were lower for fish that were released in the Cowlitz River or the Cispus River than for reservoir-released fish, but statistical comparisons were not robust because of small sample sizes at the river release sites. Fallback rates for fish released at the river release sites were 10 percent lower for steelhead, 4 percent lower for Chinook salmon, and 9 percent lower for coho salmon than for reservoir-released fish. However, fallback rates also were different between HOR and NOR fish. Fallback rates were significantly higher for HOR reservoir-released fish than for NOR reservoir-released fish. This study provided data that were insightful for understanding behavior and movement patterns in the upper Cowlitz River Basin and yielded estimates of fallback rates and fish fates that may be useful for fishery managers in the years to come. Studies from other systems have shown that factors such as prespawn mortality and fallback have resulted in substantial losses to spawning populations where trap-and-haul programs are being used as a restoration tool. Future research in the upper Cowlitz River Basin may use additional telemetry studies, genetic analyses, and spawning ground surveys to provide answers for new questions and to continue to monitor the progress of the reintroduction effort.

Washington↗

USGS2004

Hydrologic and water-quality characteristics for Bear Creek near Silver Hill, Arkansas, and selected Buffalo River sites, 1999-2004

The purpose of this report is to describe and compare the hydrologic and water-quality characteristics of Bear Creek near Silver Hill, Arkansas, to two sites on the Buffalo River upstream from the confluence of Bear Creek, to a site on Calf Creek, a smaller tributary to the Buffalo River, to selected undeveloped sites across the Nation, and to a developed site in Arkansas. A better understanding of the hydrology and water quality of Bear Creek is of interest to many, including the National Park Service, which administers the Buffalo National River, to evaluate its effects on the hydrology and water quality of the Buffalo River. The streamflow at Bear Creek near Silver Hill varied seasonally and annually from January 1999 to March 2004. The mean annual streamflow at Bear Creek for calendar years 1999 to 2003 was 86.0 cubic feet per second. The highest annual mean streamflow occurred in 2002 (158 cubic feet per second) and the lowest annual mean streamflow occurred in 1999 (56.4 cubic feet per second). The mean annual streamflow for calendar years 1999 to 2003 at the Buffalo River near Boxley and Buffalo River near St. Joe was 102 and 881 cubic feet per second, respectively. Concentrations of nitrogen measured for Bear Creek generally were greater than concentrations measured at the two Buffalo River sites and were similar to concentrations measured at Calf Creek. Concentrations of phosphorus measured at Bear Creek generally were greater than concentrations measured at the two Buffalo River sites and were similar to concentrations measured at Calf Creek. Fecal indicator bacteria concentrations generally were greater at Bear Creek than concentrations measured at the Buffalo River and similar to concentrations at Calf Creek. Bear Creek had significantly greater suspended-sediment concentrations than the Buffalo River near Boxley and the Buffalo River near St. Joe and similar concentrations to Calf Creek. Nutrients, dissolved organic carbon, and suspended-sediment loads at Bear Creek and two Buffalo River sites varied because of differences in land use and contributing drainage area for each site. In general, the Buffalo River near St. Joe had the greatest annual loads of nutrients, dissolved organic carbon, and suspended sediment. The Buffalo River near Boxley had the least annual nutrient and suspended-sediment loads amongthe three sites. Buffalo River near Boxley had lesser annual loads than the other two sites probably because of the higher percentage of forested land in the basin and smaller contributing drainage area. Mean annual nutrient, dissolved organic carbon, and suspended- sediment yields computed for Bear Creek were greater than yields computed for both of the Buffalo River sites. Bear Creek had greater median annual nutrient yields than selected undeveloped basins across the Nation and less median annual nutrient yields than the Illinois River south of Siloam Springs, Arkansas, which is representative of a developed basin. Bear Creek had greater median annual flow-weighted nutrient concentrations than the Buffalo River near St. Joe, the Buffalo River near Boxley, and selected undeveloped sites across the Nation. Bear Creek had less median flow-weighted nutrient concentrations than the Illinois River.

Scientific Investigations Report↗

USGS2006

Biomonitoring of Environmental Status and Trends (BEST) Program: Environmental contaminants, health indicators, and reproductive biomarkers in fish from the Colorado River basin

Seven fish species were collected from 14 sites on rivers in the Colorado River Basin (CDRB) from August to October 2003. Spatial trends in the concentrations of accumulative contaminants were documented and contaminant effects on the fish were assessed. Sites were located on the mainstem of the Colorado River and on the Yampa, Green, Gunnison, San Juan, and Gila Rivers. Common carp (Cyprinus carpio), black bass (Micropterus sp.), and channel catfish (Ictalurus punctatus) were the targeted species. Fish were field-examined for external and internal anomalies, selected organs were weighed to compute somatic indices, and tissue and fluid samples were preserved for fish health and reproductive biomarker analyses. Composite samples of whole fish, grouped by species and gender, from each site were analyzed for organochlorine and elemental contaminants using performance-based and instrumental methods. 2,3,7,8-tetrachlorodibenzo-p-dioxin-like activity (TCDD-EQ) was measured using the H4IIE rat hepatoma cell bioassay. Selenium (Se) and mercury (Hg) concentrations were elevated throughout the CDRB, and pesticides concentrations were greatest in fish from agricultural areas in the Lower Colorado River and Gila River. Selenium concentrations exceeded toxicity thresholds for fish (>1.0 ?g/g ww) at all sites except from the Gila River at Hayden, Arizona. Mercury concentrations were elevated (>0.1 ?g/g ww) in fish from the Yampa River at Lay, Colorado; the Green River at Ouray National Wildlife Refuge (NWR), Utah and San Rafael, Utah; the San Juan River at Hogback Diversion, New Mexico; and the Colorado River at Gold Bar Canyon, Utah, Needles, California, and Imperial Dam, Arizona. Concentrations of p,p'-DDE were relatively high in fish from Arlington, Arizona (>1.0 ?g/g ww) and Phoenix, Arizona (>0.5 ?g/g ww). Concentrations of other banned pesticides including toxaphene, total chlordanes, and dieldrin were also greatest at these two sites but did not exceed toxicity thresholds. Current-use or unlisted pesticides such as dacthal, endosulfan, '-HCH, and methoxychlor were also greatest in fish from Gila River. Total polychlorinated biphenyls (PCBs; >0.11 ?g/g ww) and TCDD-EQs (>5 pg/g ww) exceeded wildlife guidelines in fish from the Gila River at Phoenix, Arizona. Hepatic ethoxyresorufin O-deethylase (EROD) activity was also relatively high in carp from the Gila River at Phoenix, Arizona and in bass from the Green River at Ouray NWR, Utah. Altered biomarkers were noted in fish throughout the CDRB. Fish from some stations responded to chronic contaminant exposure as indicated by fish health and reproductive biomarker results. Multiple fish health indicators including altered body and organ weights and high health assessment index scores may be associated with elevated Se concentrations in fish from the Colorado River at Loma, Colorado and Needles, California. Although grossly visible external or internal lesions were found on most fish from some sites, histopathological analysis determined many of these to be inflammatory responses associated with parasites. Edema, exophthalmos, and cataracts were noted in fish from sites with elevated Se concentrations. Reproductive biomarkers including gonad development and maturation, vitellogenin concentrations, and steroid hormone concentrations were anomalous in fish from the Gila River at Hayden and Phoenix, Arizona. In addition, intersex fish were found at seven of 14 sites. The intersex condition was identified in smallmouth bass (M. dolomieu), largemouth bass (M. salmoides), channel catfish, and carp and may indicate exposure to endocrine disrupting compounds. Seven of ten male smallmouth bass from the Yampa River at Lay, Colorado were intersex. Male carp, bass, and channel catfish with low concentrations of vitellogenin were common in the CDRB. Comparatively high vitellogenin concentrations (>0.2 mg/mL) were measured in male fish from the Green River at Ouray NWR, Utah and the Colorado River at Im

Scientific Investigations Report↗

USGS2012

Surface-water salinity in the Gunnison River Basin, Colorado, water years 1989 through 2007

Elevated levels of dissolved solids in water (salinity) can result in numerous and costly issues for agricultural, industrial, and municipal water users. The Colorado River Basin Salinity Control Act of 1974 (Public Law 93-320) authorized planning and construction of salinity-control projects in the Colorado River Basin. One of the first projects was the Lower Gunnison Unit, a project to mitigate salinity in the Lower Gunnison and Uncompahgre River Basins. In cooperation with the Bureau of Reclamation (USBR), the U.S. Geological Survey conducted a study to quantify changes in salinity in the Gunnison River Basin. Trends in salinity concentration and load during the period water years (WY) 1989 through 2004 (1989-2004) were determined for 15 selected streamflow-gaging stations in the Gunnison River Basin. Additionally, trends in salinity concentration and load during the period WY1989 through 2007 (1989-2007) were determined for 5 of the 15 sites for which sufficient data were available. Trend results also were used to identify regions in the Lower Gunnison River Basin (downstream from the Gunnison Tunnel) where the largest changes in salinity loads occur. Additional sources of salinity, including residential development (urbanization), changes in land cover, and natural sources, were estimated within the context of the trend results. The trend results and salinity loads estimated from trends testing also were compared to USBR and Natural Resources Conservation Service (NRCS) estimates of off-farm and on-farm salinity reduction from salinity-control projects in the basin. Finally, salinity from six additional sites in basins that are not affected by irrigated agriculture or urbanization was monitored from WY 2008 to 2010 to quantify what portion of salinity may be from nonagricultural or natural sources. In the Upper Gunnison area, which refers to Gunnison River Basin above the site located on the Gunnison River below the Gunnison Tunnel, estimated mean annual salinity load was 110,000 tons during WY 1989-2004. Analysis of both study periods (WY 1989-2004 and WY 1989-2007) showed an initial decrease in salinity load with a minimum in 1997. The net change over either study period was only significant during WY 1989-2007. Salinity load significantly decreased at the Gunnison River near Delta by 179,000 tons during WY 1989-2004. Just downstream, the Uncompahgre River enters the Gunnison River where there also was a highly significant decrease in salinity load of 55,500 tons. The site that is located at the mouth of the study area is the Gunnison River near Grand Junction where the decrease was the largest. Salinity loads decreased by 247,000 tons during WY 1989-2004 at this site though the decrease attenuated by 2007 and the net change was a decrease of 207,000 tons. The trend results presented in this study indicate that the effect of urbanization on salinity loads is difficult to discern from the effects of irrigated agriculture and that natural sources contribute a fraction of the total salinity load for the entire basin. Based on the calculated yields and geology, 23-63 percent of the estimated annual salinity load was from natural sources at the Gunnison River near Grand Junction during WY 1989-2007. The largest changes in salinity load occurred at the Gunnison River near Grand Junction as well as the two sites located in Delta: the Gunnison River at Delta and the Uncompahgre River at Delta. Those three sites, especially the two sites at Delta, were the most affected by irrigated agriculture, which was observed in the estimated mean annual loads. Irrigated acreage, especially acreage underlain by Mancos Shale, is the target of salinity-control projects intended to decrease salinity loads. The NRCS and the USBR have done the majority of salinity control work in the Lower Gunnison area of the Gunnison River Basin, and the focus has been in the Uncompahgre River Basin and in portions of the Lower Gunnison River Basin (downstream from the Gunnison Tunnel). According to the estimates from the USBR and NRCS, salinity-control projects may be responsible for a reduction of 117,300 tons of salinity as of 2004 and 142,000 tons as of 2007 at the Gunnison River near Grand Junction, Colo. (streamflow-gaging station 09152500). USBR and NRCS estimates account for all but 130,000 tons in 2004 and 65,000 tons in 2007 of salinity load reduction. The additional reduction could be a reduction in natural salt loading to the streams because of land-cover changes during the study period. It is possible also that the USBR and NRCS have underestimated changes in salinity loads as a result of the implementation of salinity-control projects.

Colorado↗

USGS2013

Sources of suspended-sediment loads in the lower Nueces River watershed, downstream from Lake Corpus Christi to the Nueces Estuary, south Texas, 1958–2010

The U.S. Geological Survey (USGS), in cooperation with the U.S. Army Corps of Engineers, Fort Worth District; City of Corpus Christi; Guadalupe-Blanco River Authority; San Antonio River Authority; and San Antonio Water System, developed, calibrated, and tested a Hydrological Simulation Program-FORTRAN (HSPF) watershed model to simulate streamflow and suspended-sediment concentrations and loads during 1958-2010 in the lower Nueces River watershed, downstream from Lake Corpus Christi to the Nueces Estuary in south Texas. Data available to simulate suspended-sediment concentrations and loads consisted of historical sediment data collected during 1942-82 in the study area and suspended-sediment concentration data collected periodically by the USGS during 2006-7 and 2010 at three USGS streamflow-gaging stations (08211000 Nueces River near Mathis, Tex. [the Mathis gage], 08211200 Nueces River at Bluntzer, Tex. [the Bluntzer gage], and 08211500 Nueces River at Calallen, Tex. [the Calallen gage]), and at one ungaged location on a Nueces River tributary (USGS station 08211050 Bayou Creek at Farm Road 666 near Mathis, Tex.). The Mathis gage is downstream from Wesley E. Seale Dam, which was completed in 1958 to impound Lake Corpus Christi. Suspended-sediment data collected before and after completion of Wesley E. Seale Dam provide insights to the effects of the dam and reservoir on suspended-sediment loads transported by the lower Nueces River downstream from the dam to the Nueces Estuary. Annual suspended-sediment loads at the Nueces River near the Mathis, Tex., gage were considerably lower for a given annual mean discharge after the dam was completed than before the dam was completed. Most of the suspended sediment transported by the Nueces River downstream from Wesley E. Seale Dam occurred during high-flow releases from the dam or during floods. During October 1964-September 1971, about 536,000 tons of suspended sediment were transported by the Nueces River past the Mathis gage. Of this amount, about 473,000 tons, or about 88 percent, were transported by large runoff events (mean streamflow exceeding 1,000 cubic feet per second). To develop the watershed model to simulate suspended-sediment concentrations and loads in the lower Nueces River watershed during 1958-2010, streamflow simulations were calibrated and tested with available data for 2001-10 from the Bluntzer and Calallen gages. Streamflow data for the Nueces River obtained from the Mathis gage were used as input to the model at the upstream boundary of the model. Simulated streamflow volumes for the Bluntzer and Calallen gages showed good agreement with measured streamflow volumes. For 2001-10, simulated streamflow at the Calallen gage was within 3 percent of measured streamflow. The HSPF model was calibrated to simulate suspended sediment using suspended-sediment data collected at the Mathis, Bluntzer, and Calallen gages during 2006-7. Model simulated suspended-sediment loads at the Calallen gage were within 5 percent of loads that were estimated, by regression, from suspended-sediment sample analysis and measured streamflow. The calibrated watershed model was used to estimate streamflow and suspended-sediment loads for 1958-2010, including loads transported to the Nueces Estuary. During 1958-2010, on average, an estimated 288 tons per day (tons/d) of suspended sediment were delivered to the lower Nueces River; an estimated 278 tons/d were delivered to the estuary. The annual suspended-sediment load was highly variable, depending on the occurrence of runoff events and high streamflows. During 1958-2010, the annual total sediment loads to the estuary varied from an estimated 3.8 to 2,490 tons/d. On average, 113 tons/d, or about 39 percent of the estimated annual suspended-sediment contribution, originated from cropland in the study watershed. Releases from Lake Corpus Christi delivered an estimated 94 tons/d of suspended sediment or about 33 percent of the 288 tons/d estimated to have been delivered to the lower Nueces River. Erosion of stream-channel bed and banks accounted for 44 tons/d or about 15 percent of the estimated total suspended-sediment load. All other land categories, except cropland, accounted for an estimated 36 tons/d, or about 12 percent of the total. An estimated 10 tons/d of suspended sediment or about 3 percent of the suspended-sediment load delivered to the lower Nueces River were removed by water withdrawals before reaching the Nueces Estuary. During 2010, additional suspended-sediment data were collected during selected runoff events to provide new data for model testing and to help better understand the sources of suspended-sediment loads. The model was updated and used to estimate and compare sediment yields from each of 64 subwatersheds comprising the lower Nueces River watershed study area for three selected runoff events: November 20-21, 2009, September 7-8, 2010, and September 20-21, 2010. These three runoff events were characterized by heavy rainfall centered near the study area and during which minimal streamflow and suspended-sediment load entered the lower Nueces River upstream from Wesley E. Seale Dam. During all three runoff events, model simulations showed that the greatest sediment yields originated from the subwatersheds, which were largely cropland. In particular, the Bayou Creek subwatersheds were major contributors of suspended-sediment load to the lower Nueces River during the selected runoff events. During the November 2009 runoff event, high suspended-sediment concentrations in the Nueces River water withdrawn for the City of Corpus Christi public-water supply caused problems during the water-treatment process, resulting in failure to meet State water-treatment standards for turbidity in drinking water. Model simulations of the November 2009 runoff event showed that the Bayou Creek subwatersheds were the primary source of suspended-sediment loads during that runoff event.

Texas↗

USGS2013

A preliminary assessment of streamflow gains and losses for selected stream reaches in the lower Guadalupe River Basin, Texas, 2010-12

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers–Fort Worth District, the Texas Water Development Board, the Guadalupe-Blanco River Authority, and the Edwards Aquifer Authority, investigated streamflow gains and losses in the lower Guadalupe River Basin during four selected base-flow periods in March 2010, April 2011, August 2011, and, for a stream reach between Seguin, Tex., and Gonzales, Tex., in September 2012. Major sources of streamflow in this basin include releases from Canyon Lake, inflow from major springs (Comal Springs, San Marcos Springs, and Hueco Springs), and base flow (groundwater seeping to streams). Streamflow and spring-flow data were collected at 35 streamflow-gaging stations (including 6 deployed for this study) during the base-flow periods. This report describes streamflow in the lower Guadalupe River Basin, which consists of the Guadalupe River drainage basin downstream from Canyon Lake to the Guadalupe River near Tivoli, Tex. Streamflow conditions in the lower Guadalupe River Basin were analyzed by computing surface-water budgets for reaches of the lower Guadalupe River and tributary streams. Streamflow gains and losses were mapped for reaches where the computed gain or loss was greater than the uncertainty in the computed streamflow at the upstream and downstream ends of the reach. During the March 15–21, 2010, base-flow period, five reaches had gains greater than the uncertainty in the computed streamflow, including reach 1 on the Guadalupe River, which gained 130 cubic feet per second (ft 3 /s), and reach 3 on the Comal River, which gained 359 ft 3 /s. Streamflow gains during March 2010 primarily were derived from (1) inflow from the Edwards aquifer outcrop, including Hueco Springs and Comal Springs; (2) flow conveyed through the alluvium of the streambed; (3) inflows from the Carrizo-Wilcox aquifer and the Yegua Jackson aquifer; and (4) groundwater inflows from the Gulf Coast aquifer, which are enhanced by seepage losses from Coleto Creek Reservoir. During this base-flow period, none of the reaches had a loss greater in magnitude than the uncertainty in the computed streamflow. During the April 10–16, 2011, base-flow period, three reaches had gains greater than the uncertainty in the computed streamflow. Among these three reaches were reach 1 on the Guadalupe River, which gained 40.7 ft 3 /s, and reach 3 on the Comal River, which gained 271 ft 3 /s—reaches where streamflow gains were also measured in March 2010. Streamflow gains during April 2011 primarily were derived from (1) inflow from the Edwards aquifer outcrop, including Hueco Springs and Comal Springs; and (2) inflows from the Carrizo-Wilcox aquifer. During this base-flow period, three reaches had losses greater in magnitude than the uncertainty in the computed streamflow. A reach of the Blanco River near Kyle, Tex. (reach 10), lost 18.7 cubic feet per second (ft3/s). Much of this loss likely entered the groundwater system through the numerous faults that intersect the stream channel northwest of Kyle. The reach that included the confluence of the Guadalupe and San Marcos Rivers (reach 17) lost 155 ft 3 /s, likely as recharge to the Sparta and Queen City aquifers. During the August 19–25, 2011, base-flow period, three reaches had gains greater than the uncertainty in the computed streamflow, including reach 3 on the Comal River (168 ft 3 /s gain), which was one of the reaches where gains in streamflow also were measured in March 2010 and April 2011. Streamflow gains in August 2011 were primarily from (1) inflows from Comal Springs, (2) inflows from the Yegua Jackson aquifer, and (3) groundwater inflows from the Gulf Coast aquifer, which are enhanced by seepage losses from Coleto Creek Reservoir. During this base-flow period, five reaches had losses greater in magnitude than the uncertainty in the computed streamflow. The reach including the confluence of the Guadalupe and Comal Rivers lost 82.8 ft 3 /s. Much of that loss likely seeped into the local groundwater system. The reach of the Guadalupe River south of New Braunfels, Tex., to Seguin, Tex., lost 53.5 ft 3 /s. Part of that loss may have been from seepage through streambed alluvium. Reaches 9 and 10 of the Blanco River near Kyle lost 2.20 and 6.60 ft 3 /s, respectively, likely as infiltration through numerous faults intersecting the stream channel northwest of Kyle. Plum Creek between Lockhart, Tex., and Luling, Tex., lost 2.11 ft 3 /s, likely as recharge to the Carrizo-Wilcox aquifer. A base-flow period during September 22–28, 2012, was studied for the reach of the Guadalupe River between Seguin and Gonzalez, including flows from San Marcos River and Plum Creek. During this period, for the Guadalupe River reach between Seguin and Oak Forest, no computed gains or losses were greater in magnitude than the uncertainty in the computed streamflow.

Texas↗

USGS2017

Hydrogeology and simulation of groundwater flow and analysis of projected water use for the Canadian River alluvial aquifer, western and central Oklahoma

This report describes a study of the hydrogeology and simulation of groundwater flow for the Canadian River alluvial aquifer in western and central Oklahoma conducted by the U.S. Geological Survey in cooperation with the Oklahoma Water Resources Board. The report (1) quantifies the groundwater resources of the Canadian River alluvial aquifer by developing a conceptual model, (2) summarizes the general water quality of the Canadian River alluvial aquifer groundwater by using data collected during August and September 2013, (3) evaluates the effects of estimated equal proportionate share (EPS) on aquifer storage and streamflow for time periods of 20, 40, and 50 years into the future by using numerical groundwater-flow models, and (4) evaluates the effects of present-day groundwater pumping over a 50-year period and sustained hypothetical drought conditions over a 10-year period on stream base flow and groundwater in storage by using numerical flow models. The Canadian River alluvial aquifer is a Quaternary-age alluvial and terrace unit consisting of beds of clay, silt, sand, and fine gravel sediments unconformably overlying Tertiary-, Permian-, and Pennsylvanian-age sedimentary rocks. For groundwater-flow modeling purposes, the Canadian River was divided into Reach I, extending from the Texas border to the Canadian River at the Bridgeport, Okla., streamgage (07228500), and Reach II, extending downstream from the Canadian River at the Bridgeport, Okla., streamgage (07228500), to the confluence of the river with Eufaula Lake. The Canadian River alluvial aquifer spans multiple climate divisions, ranging from semiarid in the west to humid subtropical in the east. The average annual precipitation in the study area from 1896 to 2014 was 34.4 inches per year (in/yr). A hydrogeologic framework of the Canadian River alluvial aquifer was developed that includes the areal and vertical extent of the aquifer and the distribution, texture variability, and hydraulic properties of aquifer materials. The aquifer areal extent ranged from less than 0.2 to 8.5 miles wide. The maximum aquifer thickness was 120 feet (ft), and the average aquifer thickness was 50 ft. Average horizontal hydraulic conductivity for the Canadian River alluvial aquifer was calculated to be 39 feet per day, and the maximum horizontal hydraulic conductivity was calculated to be 100 feet per day. Recharge rates to the Canadian River alluvial aquifer were estimated by using a soil-water-balance code to estimate the spatial distribution of groundwater recharge and a water-table fluctuation method to estimate localized recharge rates. By using daily precipitation and temperature data from 39 climate stations, recharge was estimated to average 3.4 in/yr, which corresponds to 8.7 percent of precipitation as recharge for the Canadian River alluvial aquifer from 1981 to 2013. The water-table fluctuation method was used at one site where continuous water-level observation data were available to estimate the percentage of precipitation that becomes groundwater recharge. Estimated annual recharge at that site was 9.7 in/yr during 2014. Groundwater flow in the Canadian River alluvial aquifer was identified and quantified by a conceptual flow model for the period 1981–2013. Inflows to the Canadian River alluvial aquifer include recharge to the water table from precipitation, lateral flow from the surrounding bedrock, and flow from the Canadian River, whereas outflows include flow to the Canadian River (base-flow gain), evapotranspiration, and groundwater use. Total annual recharge inflows estimated by the soil-water-balance code were multiplied by the area of each reach and then averaged over the simulated period to produce an annual average of 28,919 acre-feet per year (acre-ft/yr) for Reach I and 82,006 acre-ft/yr for Reach II. Stream base flow to the Canadian River was estimated to be the largest outflow of groundwater from the aquifer, measured at four streamgages, along with evapotranspiration and groundwater use, which were relatively minor discharge components. Objectives for the numerical groundwater-flow models included simulating groundwater flow in the Canadian River alluvial aquifer from 1981 to 2013 to address groundwater use and drought scenarios, including calculation of the EPS pumping rates. The EPS for the alluvial and terrace aquifers is defined by the Oklahoma Water Resources Board as the amount of fresh water that each landowner is allowed per year per acre of owned land to maintain a saturated thickness of at least 5 ft in at least 50 percent of the overlying land of the groundwater basin for a minimum of 20 years. The groundwater-flow models were calibrated to water-table altitude observations, streamgage base flows, and base-flow gain to the Canadian River. The Reach I water-table altitude observation root-mean-square error was 6.1 ft, and 75 percent of residuals were within ±6.7 ft of observed measurements. The average simulated stream base-flow residual at the Bridgeport streamgage (07228500) was 8.8 cubic feet per second (ft 3 /s), and 75 percent of residuals were within ±30 ft 3 /s of observed measurements. Simulated base-flow gain in Reach I was 8.8 ft 3 /s lower than estimated base-flow gain. The Reach II water-table altitude observation root-mean-square error was 4 ft, and 75 percent of residuals were within ±4.3 ft of the observations. The average simulated stream base-flow residual in Reach II was between 35 and 132 ft 3 /s. The average simulated base-flow gain residual in Reach II was between 11.3 and 61.1 ft 3 /s. Several future predictive scenarios were run, including estimating the EPS pumping rate for 20-, 40-, and 50-year life of basin scenarios, determining the effects of current groundwater use over a 50-year period into the future, and evaluating the effects of a sustained drought on water availability for both reaches. The EPS pumping rate was determined to be 1.35 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 3.08 (acre-ft/acre)/yr in Reach II for a 20-year period. For the 40- and 50-year periods, the EPS rate was determined to be 1.34 (acre-ft/acre)/yr in Reach I and 3.08 (acre-ft/acre)/yr in Reach II. Storage changes decreased in tandem with simulated groundwater pumping and were minimal after the first 15 simulated years for Reach I and the first 8 simulated years for Reach II. Groundwater pumping at year 2013 rates for a period of 50 years resulted in a 0.2-percent decrease in groundwater-storage volumes in Reach I and a 0.6-percent decrease in the groundwater-storage volumes in Reach II. The small changes in storage are due to groundwater use by pumping, which composes a small percentage of the total groundwater-flow model budgets for Reaches I and II. A sustained drought scenario was used to evaluate the effects of a hypothetical 10-year drought on water availability. A 10-year period was chosen where the effects of drought conditions would be simulated by decreasing recharge by 75 percent. In Reach I, average simulated stream base flow at the Bridgeport streamgage (07228500) decreased by 58 percent during the hypothetical 10-year drought compared to average simulated stream base flow during the nondrought period. In Reach II, average simulated stream base flows at the Purcell streamgage (07229200) and Calvin streamgage (07231500) decreased by 64 percent and 54 percent, respectively. In Reach I, the groundwater-storage drought scenario resulted in a storage decline of 30 thousand acre-feet, or an average decline in the water table of 1.2 ft. In Reach II, the groundwater-storage drought scenario resulted in a storage decline of 71 thousand acre-feet, or an average decline in the water table of 2.0 ft.

Oklahoma↗

USGS2017

Hydrogeologic framework and selected components of the groundwater budget for the upper Umatilla River Basin, Oregon

Executive Summary This report presents a summary of the hydrogeology of the upper Umatilla River Basin, Oregon, based on characterization of the hydrogeologic framework, horizontal and vertical directions of groundwater flow, trends in groundwater levels, and components of the groundwater budget. The conceptual model of the groundwater flow system integrates available data and information on the groundwater resources of the upper Umatilla River Basin and provides insights regarding key hydrologic processes, such as the interaction between the groundwater and surface water systems and the hydrologic budget. The conceptual groundwater model developed for the study area divides the groundwater flow system into five hydrogeologic units: a sedimentary unit, three Columbia River basalt units, and a basement rock unit. The sedimentary unit, which is not widely used as a source of groundwater in the upper basin, is present primarily in the lowlands and consists of conglomerate, loess, silt and sand deposits, and recent alluvium. The Columbia River Basalt Group is a series of Miocene flood basalts that are present throughout the study area. The basalt is uplifted in the southeastern half of the study area, and either underlies the sedimentary unit, or is exposed at the surface. The interflow zones of the flood basalts are the primary aquifers in the study area. Beneath the flood basalts are basement rocks composed of Paleogene to Pre-Tertiary sedimentary, volcanic, igneous, and metamorphic rocks that are not used as a source of groundwater in the upper Umatilla River Basin. The major components of the groundwater budget in the upper Umatilla River Basin are (1) groundwater recharge, (2) groundwater discharge to surface water and wells, (3) subsurface flow into and out of the basin, and (4) changes in groundwater storage. Recharge from precipitation occurs primarily in the upland areas of the Blue Mountains. Mean annual recharge from infiltration of precipitation for the upper Umatilla River Basin during 1951–2010 is about 9.6 inches per year (in/yr). Annual recharge from precipitation for water year 2010 ranged from 3 in. in the lowland area to about 30 in. in the Blue Mountains. Using Kahle and others (2011) data and methods from the Columbia Plateau regional model, average annual recharge from irrigation is estimated to be about 2.2 in/yr for the 13 square miles of irrigated land in the upper Umatilla River Basin. Groundwater discharges to streams throughout the year and is a large component of annual streamflow in the upper Umatilla River Basin. Upward vertical hydraulic gradients near the Umatilla River indicate the potential for groundwater discharge. Groundwater discharge to the Umatilla River generally occurs in the upper part of the basin, upstream from the main stem. Groundwater development in the upper Umatilla River Basin began sometime after 1950 (Davies-Smith and others, 1988; Gonthier and Bolke, 1991). By water year 2010, groundwater use in the upper Umatilla River Basin was approximately 11,214 acre-feet (acre-ft). Total groundwater withdrawals for the study area were estimated at 7,575 acre-ft for irrigation, 3,173 acre-ft for municipal use, and 466 acre-ft for domestic use. Total groundwater flow into or from the study area depends locally on geology and hydraulic head distribution. Estimates of subsurface flow were calculated using the U.S. Geological Survey Columbia Plateau regional groundwater flow model. Net flux values range from 25,000 to 27,700 acre-ft per year and indicate that groundwater is moving out of the upper Umatilla River Basin into the lower Umatilla River Basin. Water level changes depend on storage changes within an aquifer, and storage changes depend on the storage properties of the aquifer, as well as recharge to or discharge from the aquifer. Groundwater level data in the upper Umatilla River Basin are mostly available from wells in Columbia River basalt units, which indicate areas of long-term water level declines in the Grande Ronde basalt unit near Pendleton and Athena, Oregon. Groundwater levels in the Wanapum basalt unit do not show long-term declines in the upper Umatilla River Basin. Because of pumping, some areas in the upper Umatilla River Basin have shown a decrease, or reversal, in the upward vertical head gradient. Key data needs are improvement of the spatial and temporal distribution of water-level data collection and continued monitoring of streamflow gaging sites. Additionally, refinement of recharge estimates would enhance understanding of the processes that provide the groundwater resources in the upper Umatilla River Basin.

Oregon↗

USGS2002

Measured and simulated runoff to the lower Charles River, Massachusetts, October 1999–September 2000

The lower Charles River, the water body between the Watertown Dam and the New Charles River Dam, is an important recreational resource for the Boston, Massachusetts, metropolitan area, but impaired water quality has affected its use. The goal of making this resource fishable and swimmable requires a better understanding of combined-sewer-overflow discharges, non-combined-sewer-overflow stormwater runoff, and constituent loads. This report documents the modeling effort used to calculate non-combined-sewer-overflow runoff to the lower Charles River. During the 2000 water year, October 1, 1999–September 30, 2000, the U.S. Geological Survey collected precipitation data at Watertown Dam and compiled data from five other precipitation gages in or near the watershed. In addition, surface-water discharge data were collected at eight sites—three relatively homogenous land-use sites, four major tributary sites, and the Charles River at Watertown Dam, which is the divide between the upper and lower watersheds. The precipitation and discharge data were used to run and calibrate Stormwater Management Models developed for the three land-use subbasins (single-family, multi-family, and commercial), and the two tributary subbasins (Laundry and Faneuil Brooks). These calibrated models were used to develop a sixth model to simulate 54 ungaged outfalls to the lower Charles River. Models developed by the U.S. Geological Survey at gaged sites were calibrated with up to 24 storms. Each model was evaluated by comparing simulated discharge against measured discharge for all storms with appreciable precipitation and reliable discharge data. The model-fit statistics indicated that the models generally were well calibrated to peak discharge and runoff volumes. The model fit of the commercial land-use subbasin was not as well calibrated compared to the other models because the measured flows appear to be affected by variable conditions not represented in the model. A separate Stormwater Management Model of the Stony Brook Subbasin previously developed by others was evaluated with the newly collected data from this study; this model had a model fit comparable to the models developed by the U.S. Geological Survey. The total annual runoff to the lower Charles River during the 2000 water year, not including contributions from combined-sewer-overflows except from the Stony Brook Subbasin, was 16,500 million cubic feet; 92 percent of the inflow was from the Charles River above Watertown Dam, 3 percent was from the Stony Brook Subbasin, 2 percent was from the Muddy River Subbasin, and less than 1 percent was from the combined inflows of Laundry and Faneuil Brooks. The remaining ungaged drainage area contributed about 2 percent of the total annual inflow to the lower Charles River. Excluding discharge from the Charles River above Watertown Dam, total annual runoff to the lower Charles River was 1,240 million cubic feet; 39 percent was from the Stony Brook Subbasin, 27 percent was from the Muddy River, which includes runoff that drains to the Muddy River conduit, 7 percent was from the Laundry Brook Subbasin, and 4 percent was from the Faneuil Brook Subbasin. Flow from the ungaged areas composed about 23 percent of the total annual inflow to the lower Charles River, excluding discharge from the Charles River above Watertown Dam. Runoff to the lower Charles River was calculated for two design storms representing a 3-month and a 1-year event, 1.84 and 2.79 inches of total rainfall, respectively. These simulated discharges were provided to the Massachusetts Water Resources Authority for use in a receiving-water model of the lower Charles River. Total storm runoff to the lower Charles River was 111 and 257 million cubic feet for the 3-month and 1-year storms, respectively. Excluding discharge from the Charles River above Watertown Dam, total runoff to the lower Charles River was 30 and 53 million cubic feet for the 3-month and 1-year storms, respectively. Runoff from the various tributary areas for the design storms was about in the same proportion as that for the annual runoff.

Massachusetts↗