SEARCH · USGS Science

Results for “Scientific Monthly”

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.

778 records · Page 44Linked to original sources

USGS2018

Santa Barbara and Foothill groundwater basins Geohydrology and optimal water resources management—Developed using density dependent solute transport and optimization models

Groundwater has been a part of the city of Santa Barbara’s water-supply portfolio since the 1800s; however, since the 1960s, the majority of the city’s water has come from local surface water, and the remainder has come from groundwater, State Water Project, recycled water, increased water conservation, and as needed, seawater desalination. Although groundwater from the Santa Barbara and Foothill groundwater basins only accounts for a small percentage of the long-term supply, it is an important source of supplemental water during times of surface-water shortages. During the late 1980s and early 1990s, production wells extracted additional groundwater to compensate for drought related water-delivery shortfalls from other sources; in response, water levels declined substantially in the Santa Barbara and Foothill groundwater basins (below sea level in the Santa Barbara groundwater basin). In coastal basins that have groundwater extraction near shore, seawater intrusion is often a problem. Seawater intrusion in the Santa Barbara groundwater basin is thought to be more limited than in other coastal basins because of an offshore fault that acts as a partial barrier to groundwater flow. During the late 1980s and early 1990s, seawater intrusion was observed in the Santa Barbara groundwater basin, as indicated by increased chloride concentrations at several monitoring wells that ranged from 200 ft to 1,300 ft from the ocean and as close as 2,900 ft to the nearest pumping well. This demonstrated that seawater can intrude into the Santa Barbara groundwater basin when groundwater levels fall below sea level near the coast. The city of Santa Barbara is interested in developing a better understanding of the sustainability of its groundwater supplies. In 2014, California adopted historic legislation to manage its groundwater: the Sustainable Groundwater Management Act (SGMA). The SGMA requires the development and implementation of “Groundwater Sustainability Plans” in 127 priority groundwater basins; although Santa Barbara was not a designated priority basin, the city is taking steps to achieve sustainability. Sustainability was defined in the SGMA in terms of avoiding undesirable results: significant and unreasonable groundwater-level declines, reduction in groundwater storage, seawater intrusion, water-quality degradation, land subsidence, and surface-water depletion. In this project, a cooperative study between the U.S. Geological Survey (USGS) and the city of Santa Barbara, sustainable yield is defined as the volume of groundwater that can be pumped from storage without causing water-level drawdowns and the associated increases in seawater intrusion (as indicated by increases in measured chloride concentrations) at selected wells. In order to estimate the sustainability of Santa Barbara’s groundwater basins, a three-dimensional density-dependent groundwater-flow and solute-transport model (the Santa Barbara Flow and Transport Model, or SBFTM) was developed on the basis of an existing groundwater-flow model. To simulate seawater intrusion to the Santa Barbara Basin under various management strategies, the SBFTM uses the USGS code SEAWAT to simulate salinity transport and variable-density flow. The completed SBFTM was coupled with a management optimization tool, in this case a multi-objective evolutionary algorithm, to determine optimal pumping strategies that maximize the sustainable yield and at the same time satisfy user-defined drawdown and chloride-concentration constraints. As part of this study, a three-dimensional hydrogeologic framework model was developed to quantify the extent and hydrogeologic characteristics of the Santa Barbara and Foothill groundwater basins and to help define the discretization and hydraulic properties used in the SBFTM. The development of the hydrogeologic framework model required the collection and reconciliation of geologic and geophysical data from existing maps, reports, and databases, along with geologic and hydrologic data from recently drilled wells. These data were integrated into a three-dimensional hydrogeologic framework model that defines the stratigraphy and geometry of the aquifer zones and the major geologic structures in the basin. The hydrogeologic framework model also quantifies the variation in sediment grain size within each aquifer zone as the percentage of coarse-grained sediment. Previous studies indicated that there are two principal water-producing zones in the Santa Barbara groundwater basin, the upper and lower producing zones; an additional thin, productive zone was identified as part of this study. This “middle producing zone” is not as areally extensive as the upper and lower producing zones and only exists in the coastal part of Storage Unit I. These producing zones are bounded at depth by less productive shallow, middle, and deep zones. Two versions of the SBFTM were constructed: an initial-condition model and a modern transient model. The initial-condition model is a long-term transient model that simulates flow and solute-transport conditions during a period with limited anthropogenic influences preceeding the modern transient model. The simulation-transient model simulates flow and transport conditions from 1929 through 2013; however, because of data availability, the focus of the model calibration was 1972–2013. The SBFTM was calibrated to measured groundwater levels and drawdown, as well as measured chloride concentrations and change in concentrations, using a combination of automated and trial-and-error parameter-estimation techniques. A sensitivity analysis indicated that, in general, the SBFTM was most sensitive to recharge- and pumping-distribution parameters, specifically those controlling the amount of small-catchment recharge and the distribution of water extraction by hydrogeologic layer for production wells. The model was also sensitive to parameters controlling stream-recharge rates, horizontal and vertical hydraulic conductivity, and porosity. From 1929 to 1971, most of the water entering the area represented by the SBFTM was from creek and small-catchment recharge, and the majority of water leaving the SBFTM area was from pumping, discharge to creeks, and drains. In addition, about 37 percent of the total pumpage came from a net reduction in groundwater storage. From 1972 to 2013, the amount of water entering and leaving the SBFTM was fairly similar as that from 1929 to 1971, except the reduction in pumpage added about 17,000 acre-ft of water to storage. During this later period, there were also times of storage loss. For example, during July 1990, a month when approximately 705 acre-ft of groundwater was pumped in the study area, the pumpage was much greater than all sources of recharge combined, and about 382 acre-ft of water was removed from groundwater storage. Simulated hydraulic heads replicated the observed data to an acceptable matching of the measured water-level, flow direction, and vertical gradients. Simulated hydrographs for selected wells were in good agreement with the measured data, with an average residual of -2.7 ft and a standard deviation of 14.5 ft, indicating that the simulated heads, on average, underestimated the observed water levels. An examination of the model fit indicated that most of the discrepancies were lower simulated heads at wells proximal to production well sites. The simulated chloride concentrations reasonably matched the rising limbs of the measured breakthrough curves in terms of timing and magnitude; however, the simulation overestimated the chloride concentrations on the falling limbs. The overestimation of low chloride concentrations was attributed to the model overestimating the advance of the chloride front during periods of heavy pumping and underestimating the retreat of the chloride front during periods of low pumping. These simulation errors would result in a conservative response by local water managers to seawater intrusion. The SBFTM was used to develop a collection of predictive simulations optimized to produce pumping schedules that maximize yield, subject to a set of constraints and competing objectives. The simulations were grouped as scenarios that differed in their time horizon, initial conditions for groundwater levels and chloride concentrations, as well as precipitation, which was incorporated into the model through simulated recharge. Overall, five scenarios were developed in a multi-objective framework to obtain optimal pumping rates for all of the wells managed by the city, while minimizing excessive drawdown and seawater intrusion. For the current study, complexities in the simulation model and the optimization formulation required additional considerations. Incorporating the solute-transport equations to simulate chloride transport added a highly nonlinear process that is solved iteratively in each time step of the groundwater-flow model. These nonlinearities, coupled with the highly refined grid in the current model, creates challenges for many traditional optimization methods. Therefore, an optimization method was needed that could address nonlinear relationships as well as a very large problem size. Lastly, the optimization problem was reformulated to include multiple objectives without requiring convergence to a single solution. This approach, guided by the city’s objectives, allowed the maximum extraction of information from the complex simulation. Borg, a multi-objective evolutionary algorithm, was chosen as the optimization algorithm for this study for several reasons: (1) it is very computationally efficient; (2) it can run in parallel; (3) it requires little user input; and (4) it can solve for multiple competing objectives. The first three points allow the algorithm to proceed toward the optimal solutions at the fastest possible rate. The fourth point is advantageous for large, complex optimization problems because it is difficult to formulate the optimization problem in a way that produces only one optimal solution. The problem formulation consisted of four competing objectives and a constraint set in accordance with the main concerns of the city. The objectives were maximizing total pumpage, minimizing seawater intrusion, minimizing total drawdown in production wells, and minimizing the maximum drawdown. The constraints were pump capacity, meeting drinking-water standards for chloride, maintaining a specified minimum flowrate to a groundwater treatment plant, and maintaining minimum water levels in pumping wells. The decision variables either were quarterly pumpage by well or total pumpage by basin. Five optimization scenarios were developed that allow the decision makers to evaluate a range of optimal solutions for a variety of water levels and chloride concentrations as well as potential future climatic conditions. Three scenarios (1, 2, and 5) were multi-objective optimization formulations that allowed for variations in management preferences and climatic conditions. The other two scenarios (3 and 4) were designed to examine the optimization results to answer specific questions. Scenario 1 described the best-case sustainable yield assuming a “full” basin (that is, high initial water levels) and typical climate conditions for 10 years. Scenario 2 also started with a “full” basin; however, this was followed by a 10-year drought. Scenario 3 determined if an “empty” basin (that is, low initial water levels) would recover to full conditions (1998 conditions) given climate assumptions and optimal pumping schedules from scenarios 1 and 2. Scenario 4 was designed to produce decision rules that can be used by water managers to help choose an optimal pumping schedule based on measured water-level or chloride data. Scenario 5 identified future pumping schedules based on short-term climate variations during a 2-year management horizon. The results from scenarios 1 and 2 described the differences in maximum pumpage in the basin under typical and dry long-term climate projections, respectively. The scenario 1 results indicated the maximum 10-year pumpage of the basin was about 31,300 acre-ft under typical conditions and controlling simulated seawater intrusion and drawdowns. For scenario 2, less recharge over the 10-year dry climate produced a maximum pumpage estimate of 30,000 acre-ft to control seawater intrusion and drawdowns. The larger pumpage for scenario 1 resulted in more seawater intrusion, but less total drawdown, compared to that of scenario 2. Results for scenarios 3 and 4 showed the basin’s response to management actions combined with climate projections. Both scenarios used the optimal pumping schedules and the 10-year climates from scenarios 1 and 2. The scenario 3 results showed that under minimal pumping, the basin did not fully recover to 1998 water levels within 10 years under either climate scenario. The relatively larger recharge from the typical climate resulted in less drawdown at coastal monitoring wells after the 10-year recovery period than that from the dry climate. The location of the seawater intrusion front was not appreciably different between the scenarios, however. Scenario 4 used the optimal results from scenarios 1 and 2 to produce decision-rule curves that illustrated the pumpage for each basin, given measured levels of chloride concentration or drawdown. This allowed the use of additional measurements at monitoring wells to assess future management decisions on the basis of the sensitivity of observations of drawdown and seawater intrusion to various pumping rates. Scenario 5 allowed managers to investigate the effects of short-term climate variations on optimal pumping schedules. Three specific 2-year simulations were optimized: typical-to-dry (scenario 5A), dry-to-typical (scenario 5B), and dry-to-dry (scenario 5C). The most noteable result from scenario 5 was the overall reduction in optimal pumpage for most schedules in scenario 5C, when the climate is simulated as dry-to-dry. There are also many optimal pumping schedules that produced an overall increase in waterlevels over the two-year simulation period, regardless of climatic condition. Similar to scenario 2, the scenario 5C results represents conservative yield estimates under a minimal-precipitation climatic condition.

California↗

USGS2013

Sources and characteristics of organic matter in the Clackamas River, Oregon, related to the formation of disinfection by-products in treated drinking water

This study characterized the amount and quality of organic matter in the Clackamas River, Oregon, to gain an understanding of sources that contribute to the formation of chlorinated and brominated disinfection by-products (DBPs), focusing on regulated DBPs in treated drinking water from two direct-filtration treatment plants that together serve approximately 100,000 customers. The central hypothesis guiding this study was that natural organic matter leaching out of the forested watershed, in-stream growth of benthic algae, and phytoplankton blooms in the reservoirs contribute different and varying proportions of organic carbon to the river. Differences in the amount and composition of carbon derived from each source affects the types and concentrations of DBP precursors entering the treatment plants and, as a result, yield varying DBP concentrations and species in finished water. The two classes of DBPs analyzed in this study-trihalomethanes (THMs) and haloacetic acids (HAAs)-form from precursors within the dissolved and particulate pools of organic matter present in source water. The five principal objectives of the study were to (1) describe the seasonal quantity and character of organic matter in the Clackamas River; (2) relate the amount and composition of organic matter to the formation of DBPs; (3) evaluate sources of DBP precursors in the watershed; (4) assess the use of optical measurements, including in-situ fluorescence, for estimating dissolved organic carbon (DOC) concentrations and DBP formation; and (5) assess the removal of DBP precursors during treatment by conducting treatability "jar-test" experiments at one of the treatment plants. Data collection consisted of (1) monthly sampling of source and finished water at two drinking-water treatment plants; (2) event-based sampling in the mainstem, tributaries, and North Fork Reservoir; and (3) in-situ continuous monitoring of fluorescent dissolved organic matter (FDOM), turbidity, chlorophyll- a , and other constituents to continuously track source-water conditions in near real-time. Treatability tests were conducted during the four event-based surveys to determine the effectiveness of coagulant and powdered activated carbon (PAC) on the removal of DBP precursors. Sample analyses included DOC, total particulate carbon (TPC), total and dissolved nutrients, absorbance and fluorescence spectroscopy, and, for regulated DBPs, concentrations of THMs and HAAs in finished water and laboratory-based THM and HAA formation potentials (THMFP and HAAFP, respectively) for source water and selected locations throughout the watershed. The results of this study may not be typical given the record and near record amounts of precipitation that occurred during spring that produced streamflow much higher than average in 2010-11. Although there were algal blooms, lower concentrations of chlorophyll- a were observed in the water column during the study period compared to historical data. Concentrations of DBPs in finished (treated) water averaged 0.024 milligrams per liter (mg/L) for THMs and 0.022 mg/L for HAAs; maximum values were about 0.040 mg/L for both classes of DBPs. Although DBP concentrations were somewhat higher within the distribution system, none of the samples collected for this study or for the quarterly compliance monitoring by the water utilities exceeded levels permissible under existing U.S. Environmental Protection Agency (USEPA) regulations: 0.080 mg/L for THMs and 0.060 mg/L for HAAs. DOC concentrations were generally low in the Clackamas River, typically about 1.0-1.5 mg/L. Concentrations in the mainstem occasionally increased to nearly 2.5 mg/L during storms; DOC concentrations in tributaries were sometimes much higher (up to 7.8 mg/L). The continuous in-situ FDOM measurements indicated sharp rises in DOC concentrations in the mainstem following rainfall events; concentrations were relatively stable during summer base flow. Even though the first autumn storm mobilized appreciable quantities of carbon, higher concentrations of DBPs in finished water were observed 3-weeks later, after the ground was saturated from additional rainfall. The majority of the DOC in the lower Clackamas River appears to originate from the upper basin, suggesting terrestrial carbon was commonly the dominant source. Lower-basin tributaries typically contained the highest concentrations of DOC and DBP precursors and contributed substantially to the overall loads in the mainstem during storms. During low-flow periods, tributaries were not major sources of DOC or DBP precursors to the Clackamas River. Although the dissolved fraction of organic carbon contributed the majority of DBP precursors, at times the particulate fraction (inorganic sediment and organic particles including detritus and algal material) contributed a substantial fraction of DBP precursors. Considering just the main-stem sites, on average, 10 percent of THMFP and 32 percent of HAAFP were attributed to particulate carbon. This finding suggests water-treatment methods that remove particles prior to chlorination would reduce finished-water DBP concentrations to some degree. Overall, concentrations of THM and HAA precursors were closely linked to DOC concentrations; laboratory DBP formation potentials (DBPFPs) clearly showed that THMFP and HAAFP were greatest in the downstream tributaries that contained elevated carbon concentrations. However, carbon-normalized "specific" formation potentials for THMs and HAAs (STHMFP and SHAAFP, respectively) revealed changes in carbon character over time that affected the two types of DBP classes differently. HAA precursors were elevated in waters containing aromatic-rich soil-derived material arising from forested areas. In contrast, THM precursors were associated with carbon having a lower aromatic content; highest STHMFP occurred in autumn 2011 in the mainstem from North Fork Reservoir downstream to LO DWTP. This pattern suggests the potential for a link between THM precursors and algal-derived carbon. The highest STHMFP value was measured within North Fork Reservoir, indicating reservoir derived carbon may be important for this class of DBPs. Weak correlations between STHMFP and SHAAFP emphasize that precursor sources for these types of DBPs may be different. This highlights not only that different locations within the watershed produce carbon with different reactivity (specific DBPFP), but also that different management approaches for each class of DBP precursors could be required for control. Treatability tests conducted on source water during four basin-wide surveys demonstrated that an average of about 40 percent of DOC can be removed by coagulation. While the decrease in THMFP following coagulation was similar to DOC, the decrease in HAAFP was much greater (approximately 70 percent), indicating coagulation is particularly effective at removing HAA precursors'likely because of the aromatic nature of the carbon associated with HAA precursors. Several findings from this study have direct implications for managing drinking-water resources and for providing useful information that may help improve treatment-plant operations. For example, the use of in-situ fluorometers that measure FDOM provided an excellent proxy for DOC concentration in this system and revealed short-term, rapid changes in DOC concentration during storm events. In addition, the strong correlation between FDOM values measured in-situ and HAA5 concentrations in finished water may permit estimation of continuous HAA concentrations, as was done here. As part of this study, multiple in-situ FDOM sensors were deployed continuously and in real-time to characterize the composition of dissolved organic matter. Although the initial results were promising, additional research and engineering developments will be needed to demonstrate the full utility of these sensors for this purpose. In conclusion, although DBPFPs were strongly correlated to DOC concentration, some DBPs formed from particulate carbon, including terrestrial leaf material and algal material such as planktonic species of blue-green algae and sloughed filaments, stalks, and cells of benthic algae. Different precursor sources in the watershed were evident from the data, suggesting specific actions may be available to address some of these sources. In-situ measurements of FDOM proved to be an excellent proxy for DOC concentration as well as HAA formation during treatment, which suggests further development and refinement of these sensors have the potential to provide real-time information about complex watershed processes to operators at the drinking-water treatment plants. Follow-up studies could examine the relative roles that terrestrial and algal sources have on the DBP precursor pool to better understand how watershed-management activities may be affecting the transport of these compounds to Clackamas River drinking-water intakes. Given the low concentrations of algae in the water column during this study, additional surveys during more typical river conditions could provide a more complete understanding of how algae contribute DBP precursors. Further development of FDOM-sensor technology can improve our understanding of carbon dynamics in the river and how concentrations may be trending over time. This study was conducted in collaboration with Clackamas River Water and the City of Lake Oswego water utilities. Other research partners included Oregon Health and Science University in Hillsboro, Oregon, Alexin Laboratory in Tigard, Oregon, U.S. Geological Survey National Research Program Laboratory in Denver, Colorado, and the U.S. Geological Survey Water Science Centers in Portland, Oregon, and Sacramento, California. This project was supported with funding from Clackamas River Water, City of Lake Oswego, the U.S. Geological Survey, and the Water Research Foundation.

Oregon↗

USGS2010

Effects of selected low-impact-development (LID) techniques on water quality and quantity in the Ipswich River Basin, Massachusetts: Field and modeling studies

During the months of August and September, flows in the Ipswich River, Massachusetts, dramatically decrease largely due to groundwater withdrawals needed to meet increased residential and commercial water demands. In the summer, rates of groundwater recharge are lower than during the rest of the year, and water demands are higher. From 2005 to 2008, the U.S. Geological Survey, in a cooperative funding agreement with the Massachusetts Department of Conservation and Recreation, monitored small-scale installations of low-impact-development (LID) enhancements designed to diminish the effects of storm runoff on the quantity and quality of surface water and groundwater. Funding for the studies also was contributed by the U.S. Environmental Protection Agency’s Targeted Watersheds Grant Program through a financial assistance agreement with Massachusetts Department of Conservation and Recreation. The monitoring studies examined the effects of (1) replacing an impervious parking lot surface with a porous surface on groundwater quality, (2) installing rain gardens and porous pavement in a neighborhood of 3 acres on the quantity and quality of stormwater runoff, and (3) installing a 3,000-square foot (ft2) green roof on the quantity and quality of stormwater runoff. In addition, the effects of broad-scale implementation of LID techniques, reduced water withdrawals, and water-conservation measures on streamflow in large areas of the basin were simulated using the U.S. Geological Survey’s Ipswich River Basin model. From June 2005 to 2007, groundwater quality was monitored at the Silver Lake town beach parking lot in Wilmington, MA, prior to and following the replacement of the conventional, impervious-asphalt surface with a porous surface consisting primarily of porous asphalt and porous pavers. Changes in the concentrations of the water-quality constituents, phosphorus, nitrogen, cadmium, chromium, copper, lead, nickel, zinc, and total petroleum hydrocarbons, were monitored. Increased infiltration of precipitation did not result in discernible increases in concentrations of these potential groundwater contaminants. Concentrations of dissolved oxygen increased slightly in groundwater profiles following the removal of the impervious asphalt parking lot surface. In Wilmington, MA, in a 3-acre neighborhood, stormwater runoff volume and quality were monitored to determine the ability of selected LID enhancements (rain gardens and porous paving stones) to reduce flows and loads of the above constituents to Silver Lake. Flow-proportional water-quality samples were analyzed for nutrients, metals, total petroleum hydrocarbons, and total-coliform and Escherichia coli bacteria. In general, when all storms were considered, no substantial decreases were observed in runoff volume as a result of installing LID enhancements. However, the relation between rainfall and runoff did provide some insight into how the LID enhancements affected the effective impervious area for the neighborhood. A decrease in runoff was observed for storms of 0.2 inches (in.) or less of precipitation, which indicated a reduction in effective impervious area from approximately 10 percent to about 4.5 percent for the 3-acre area. Water-quality-monitoring results were inconclusive; there were no statistically significant differences in concentrations or loads when the pre- and post-installation-period samples were compared. Three factors were probably most important in minimizing differences: (1) the small decrease in effective impervious area, (2) the differences in the size of storms sampled for water-quality constituents before and after installation of the infiltration enhancing measures, and (3) small sample sizes. In a third field study, the characteristics of runoff from a vegetated “green” roof and a conventional, rubber-membrane roof were compared. The amount of precipitation and the length of the antecedent dry period were the two primary factors affecting the green roof’s water-storage capacity. The green roof retained more than 50 percent of the precipitation from storms with 0.04 to 1.0 in. of rain. Approximately 95 percent of the precipitation from one storm of nearly 2 in. was retained by the green roof. On the rubber-membrane roof, only a small, shallow puddle of insubstantial volume ever remained after a storm. Bulk precipitation from 10 storms was monitored for the same constituents (nutrients, metals, and total petroleum hydrocarbons) as the roof runoff, and the results were compared with those for roof-runoff samples. The use of fertilizers to help establish the vegetation during the study probably distorted any effect the plants and growing medium may have had on the retention of target analytes. As a result of the fertilizer and growing medium chemistry, median concentrations of total nitrogen, total phosphorus, cadmium, copper, and nickel in runoff from the green roof were greater than in the runoff from the conventional roof or in bulk precipitation. Concentrations of lead and zinc were greater in runoff from the conventional roof, probably a result of passage through the old, metal drainpipes. Simulations of the effects of LID on streamflow in the Ipswich River Basin were conducted with a previously calibrated Hydrological Simulation Program-FORTRAN (HSPF) precipitation-runoff model. Simulations were conducted at multiple spatial scales to evaluate the effects of (1) updated water withdrawals for the towns of Reading and Wilmington; (2) potential land-use changes at buildout (potential future development); (3) effective impervious area reductions upstream from the South Middleton streamgage to represent the effects of widespread implementation of LID retrofit techniques; (4) basin-scale water withdrawal reductions scaled up (expanded to the town level) from water-conservation pilot programs conducted by the Massachusetts Department of Conservation and Recreation; and (5) land-use change and LID techniques at a local scale, which is smaller than the HSPF subbasin. Effects on streamflow generally were evaluated by comparing results of two or more related simulations for selected reaches in the basin; thus, relative rather than absolute changes in simulated flow were the focus of the assessment. Simulations indicated that reduced withdrawals for the towns of Reading and Wilmington led to substantially higher medium and low flows in most of the reaches upstream from the South Middleton streamgage. Simulations of water-conservation measures resulted in negligible effects on streamflow. Overall, simulations indicated that spatial scale is an important factor in determining the effects of land-use change and LID practices on streamflow. Potential land-use changes at buildout had modest (percent differences of less than 20 percent) effects on streamflow in most subbasins because relatively little land in the basin was available for development (about 17 percent); moreover, most of the available open land is zoned for low-density residential development, and this land-use category was simulated to contain relatively little effective impervious area and to be similar hydrologically to the forested land in place prior to development. Results of the simulations conducted to evaluate widespread effective impervious area reductions upstream from the South Middleton streamgage indicated that the percentage of urban land use and associated effective impervious area was too small for a 50-percent reduction of effective impervious area to appreciably affect streamflow (percent differences of less than 20 percent) in most subbasins. In contrast, the results of the hypothetical local-scale simulations indicated that for smaller streams, where the percentage of urban land use and associated effective impervious area in the drainage area may be substantially higher, land-use change, development patterns, and LID practices potentially have much greater effects on streamflow. Modeling results also indicated that LID was potentially most beneficial for minimizing streamflow alteration when applied to dense urban development, largely because larger tracts of effective impervious area were available for reduction than were available for other land-use categories. For example, commercial-industrial-transportation land use is composed of 37 percent pervious area and 63 percent effective impervious area in the HSPF model, whereas low-density residential area is composed of 97.5 percent pervious area and only 2.5 percent effective impervious area. Field and modeling studies concurred in the assessment that LID enhancements would likely have the greatest effect on decreasing stormwater runoff when broadly applied to highly impervious urban areas. A measurable effect for small rainfall events (less than 0.25 inch) was determined in the small, highly pervious area that was monitored in this study, but the volume difference was not great.

Massachusetts↗

USGS2016

Hydrogeological framework, numerical simulation of groundwater flow, and effects of projected water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma

This report describes a study of the hydrology, hydrogeological framework, numerical groundwater-flow models, and results of simulations of the effects of water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma. The purpose of the study was to provide analyses, including estimating equal-proportionate-share (EPS) groundwater-pumping rates and the effects of projected water use and droughts, pertinent to water management of the Beaver-North Canadian River alluvial aquifer for the Oklahoma Water Resources Board. The Beaver-North Canadian River alluvial aquifer consists of unconsolidated sand, gravel, silt, and clay in varying proportions that underlies the Beaver and North Canadian River Valleys for approximately 175 miles (mi) from the Oklahoma Panhandle to the western edge of Oklahoma City in central Oklahoma. The aquifer as delineated for this study varies from 4 to 12 mi wide and is as thick as 308 feet (ft) in the northwest where the aquifer includes the Ogallala Formation. There are two distinct but in most areas hydraulically connected alluvial units that compose the Beaver-North Canadian River alluvial aquifer: a Quaternary-age topographically higher terrace deposit and a topographically lower, younger alluvium along the active river channel that includes active and Quaternary-age alluvium. The Beaver River composes the headwaters of the North Canadian River, which begins at the confluence of the Beaver River and Wolf Creek. The aquifer is divided for water management into two geographic areas: Reach I upstream from Canton Dam and Reach II downstream from Canton Dam. Reach I covers an area of approximately 874 square miles (mi 2 ), and Reach II covers an area of approximately 371 mi 2 . The Beaver-North Canadian River alluvial aquifer crosses several climatic zones, from semiarid in the west to continental subhumid in the east. Mean annual precipitation varies from 23.5 inches (in.) in the western part of this aquifer to 35.7 in. in the east. Surface-water demands were met through numerous temporary and permanent surface-water diversions from the Beaver and North Canadian Rivers during the period of study. During the study period, seven diversions removed a mean annual 2,000 acre-feet (acre-ft) of water from Reach I. There were 14 diversions from Reach II with a mean annual permitted volume of approximately 81,000 acre-ft, including diversion into the Lake Hefner Canal for the Oklahoma City public water supply. During the period of this study, 17 temporary surface-water diversion permits were active in Reach I, with total permitted volumes of 2,000 acre-ft, and 41 diversions were active in Reach II, with total permitted volumes of 38,000 acre-ft. The total water use for each temporary permit was assumed to be taken over the 3-month period allotted to temporary withdrawal permits. The groundwater-use analysis full period of record, 1967–2011, was divided into two sub-intervals because of varying water use, 1970–80 and 1981–2011. Groundwater use in Reach I and Reach II was substantially greater from 1970 to 1980 compared to the rest of the period, and the sub-period 1981–2011 was used because this period includes recent population growth and modern irrigation methods. The total mean annual groundwater use in Reach I was 15,309 acre-feet per year (acre-ft/yr) during 1967–2011; 20,724 acre-ft/yr during 1970–80, and 13,739 acre-ft/yr during 1981–2011. Total mean annual groundwater use in Reach II was similar but slightly less than in Reach I, with 14,098 acre-ft/yr during 1967–2011; 19,963 acre-ft/yr during 1970–80; and 12,285 acre-ft/yr during 1981–2011. Irrigation composed 72 percent of groundwater use in Reach I and 48 percent of groundwater use in Reach II during the 1967–2011 period. Public water supply was a much smaller proportion of total groundwater use in Reach I (15 percent) than in Reach II (39 percent). The proportion of groundwater use for power was 10 percent in Reach I and 5.2 percent in Reach II. All other water-use categories in Reach I only composed 2.2 percent of groundwater use in Reach I. In Reach II, industrial, mining, and commercial categories combined accounted for 4.4 percent of groundwater use; recreation, fish, and wildlife groundwater use accounted for 2.3 percent; and nonirrigated agriculture accounted for 1.5 percent of groundwater use. Permian-age bedrock underlies the Beaver-North Canadian River alluvial aquifer. In the east, the Dog Creek Shale, the Duncan Sandstone, and the Blaine and Chickasha Formations, none of which are notable sources of groundwater in the study area, underlie the Beaver-North Canadian River alluvial aquifer. In the northwestern part of Reach I, bedrock is composed of the Rush Springs and Marlow Formations, which are productive aquifers in some areas. The Cloud Chief Formation is not a source of groundwater. One hydrogeological unit was delineated in the Beaver-North Canadian River alluvial aquifer, composed of the terrace deposits and alluvium, with limited flow between this unit and bedrock units. Groundwater in this aquifer generally flows from northwest to southeast and across the aquifer toward the Beaver and North Canadian Rivers. Groundwater recharge from precipitation was estimated for the entire Beaver-North Canadian River alluvial aquifer and then itemized for both reaches by using a soil-water-balance (SWB) model. At two locations in Reach I, a water-table fluctuation method was used to estimate local recharge. Total mean annual groundwater recharge from the soil-water-balance method was estimated to be approximately 136,400 acre-ft in Reach I and 82,400 acre-ft in Reach II; the mean annual recharge for both reaches combined was approximately 218,800 acre-ft. Two sites in Reach I located at observation wells with continuous water-level measurements and nearby streamflow-gaging stations with precipitation gages were used to estimate the percentage of precipitation that becomes groundwater recharge. The Woodward site was located at observation well OW-4 near the Woodward, Okla. (07237500), streamflow-gaging station. Total precipitation and recharge for the Woodward and Seiling sites were calculated for the water year 2013. The Woodward site had a total of 14.18 in. of precipitation and 6.3 in. of recharge was calculated, equaling 44 percent of precipitation. The mean percentage of precipitation that was estimated to become recharge in the SWB model for the period 1980–2011 at that location was 9.2 percent, although adjacent SWB-model cells were as high as 20 percent of precipitation. The Seiling site had a total of 26.84 in. of precipitation during the water year 2013, and a total of 6.9 in. of recharge was estimated, equaling 25.9 percent of precipitation. At the Seiling site, the mean percentage of precipitation that became recharge in the SWB model for the period 1980–2011 was 23.0 percent. The principal inflow to the Beaver-North Canadian River alluvial aquifer was estimated to be surface recharge from precipitation, and plant evapotranspiration was estimated to be the greatest discharge, followed by stream and lake base flow, groundwater pumping, and flow to seeps and springs along the eastern margin of the aquifer. Reach I also included inflow from the High Plains aquifer as lateral inflow of groundwater, though this flow was estimated to be a very minor component of the total water budget. Most of the Beaver and North Canadian Rivers were determined to be gaining streamflow from groundwater, but several reaches in Reach I upstream from Wolf Creek were determined to be losing streamflow through infiltration to the aquifer. Aquifer hydrogeologic characteristics were estimated from borehole lithologic logs, well-construction information, and published aquifer tests and during numerical model calibration. The maximum saturated aquifer thickness in Reach I was estimated to be 308 ft, and the mean thickness was estimated to be 36 ft. The maximum saturated thickness in Reach II was estimated to be 86 ft, and the mean thickness was estimated to be 29 ft. Mean hydraulic conductivity of Reach I was estimated to be 70 feet per day (ft/d) with a range of 7–279 ft/d. Mean hydraulic conductivity in Reach II was estimated to be 92 ft/d with a range of 4–279 ft/d. Both reach models were calibrated manually by using trial-and-error adjustment of recharge, hydraulic conductivity, specific yield, and conductance of boundary conditions. The Reach I model used 28 head observations during the steady-state period of 1980 and 487 head observations during the transient period of 1981–2011. The root-mean-square error of head residuals (observed minus simulated head) was 3.86 ft, and 83 percent of head residuals were between -5 and 5 ft. The Reach II model was calibrated to 75 steady-state head observations and 134 head observations during the transient period. The root-mean-square error of head residuals for that reach was 3.58 ft, and similar to Reach I, 85 percent of residuals were between -5 and 5 ft. Several analyses were performed by using the numeric groundwater-flow models as predictive tools, including estimating the EPS pumping rate for both reaches. The EPS is defined by the Oklahoma Water Resources Board as an annual per-acre groundwater-pumping rate that will reduce saturated thickness in half of the aquifer to 5 ft or less over a period of 20 years; additional estimates were made for periods of 40 and 50 years. Other analyses included using models to estimate the effects of groundwater pumping and a prolonged drought on groundwater in storage and streamflow and lake storage of water. The EPS pumping rate was found to be approximately 0.57 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 0.73 (acre-ft/acre)/yr in Reach II for a 20-year period. For a 40-year period, the annual EPS pumping rate was determined to be 0.54 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. For a 50-year period, the EPS pumping rate was determined to be 0.53 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. Groundwater pumping at the 2011 rate for 50 years resulted in a 3.6-percent decrease in the amount of water in groundwater storage in Reach I and a decrease of 2.5 percent in the amount of groundwater in storage in Reach II. A cumulative 32-percent increase in pumping greater than the 2011 rate over a period of 50 years caused a decrease in groundwater storage of 4.0 percent in Reach I and 3.3 percent in Reach II. A hypothetical severe drought was simulated by using aquifer recharge flow rates during the drought year of 2011 for a period of 10 years. All other flows including evapotranspiration and groundwater pumping were set at estimated 2011 rates. The hypothetical drought caused a decrease in water in aquifer storage by about 7 percent in Reach I and 7 percent in Reach II. Another analysis of the effects of hypothetical drought estimated the effects of drought on streamflow and lake storage. The hypothetical drought was simulated by decreasing recharge by 75 percent for a selected 10-year period (1994–2004) during the 1980–2011 simulation. In Reach I, the amounts of water stored in Canton Lake and streamflow at the Seiling, Okla., streamflow-gaging station were analyzed. Streamflow at the Seiling station decreased by a mean of 75 percent and was still diminished by 10 percent after 2011. In Reach II, the effect of drought on the streamflow at the Yukon, Okla., streamflow-gaging station was examined. The greatest mean streamflow decrease was approximately 60 percent during the simulated drought, and after 2011, the mean decrease in streamflow was still about 5 percent. Canton Lake storage decreased by as much as 83 percent during the simulated drought and did not recover by 2011.

Oklahoma↗