USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water-Resources Investigations Report”

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

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

At least 937 records · Page 52Linked to original sources

Measurement of unsaturated hydraulic properties and evaluation of property-transfer models for deep sedimentary interbeds, Idaho National Laboratory, Idaho

Operations at the Idaho National Laboratory (INL) have the potential to contaminate the underlying Eastern Snake River Plain (ESRP) aquifer. Methods to quantitatively characterize unsaturated flow and recharge to the ESRP aquifer are needed to inform water-resources management decisions at INL. In particular, hydraulic properties are needed to parameterize distributed hydrologic models of unsaturated flow and transport at INL, but these properties are often difficult and costly to obtain for large areas. The unsaturated zone overlying the ESRP aquifer consists of alternating sequences of thick fractured volcanic rocks that can rapidly transmit water flow and thinner sedimentary interbeds that transmit water much more slowly. Consequently, the sedimentary interbeds are of considerable interest because they primarily restrict the vertical movement of water through the unsaturated zone. Previous efforts by the U.S. Geological Survey (USGS) have included extensive laboratory characterization of the sedimentary interbeds and regression analyses to develop property-transfer models, which relate readily available physical properties of the sedimentary interbeds (bulk density, median particle diameter, and uniformity coefficient) to water retention and unsaturated hydraulic conductivity curves. During 2013–14, the USGS, in cooperation with the U.S. Department of Energy, focused on further characterization of the sedimentary interbeds below the future site of the proposed Remote Handled Low-Level Waste (RHLLW) facility, which is intended for the long-term storage of low-level radioactive waste. Twelve core samples from the sedimentary interbeds from a borehole near the proposed facility were collected for laboratory analysis of hydraulic properties, which also allowed further testing of the property-transfer modeling approach. For each core sample, the steady-state centrifuge method was used to measure relations between matric potential, saturation, and conductivity. These laboratory measurements were compared to water-retention and unsaturated hydraulic conductivity parameters estimated using the established property-transfer models. For each core sample obtained, the agreement between measured and estimated hydraulic parameters was evaluated quantitatively using the Pearson correlation coefficient (r). The highest correlation is for saturated hydraulic conductivity ( K sat ) with an r value of 0.922. The saturated water content (q sat ) also exhibits a strong linear correlation with an r value of 0.892. The curve shape parameter (λ) has a value of 0.731, whereas the curve scaling parameter (y o ) has the lowest r value of 0.528. The r values demonstrate that model predictions correspond well to the laboratory measured properties for most parameters, which supports the value of extending this approach for quantifying unsaturated hydraulic properties at various sites throughout INL.

Idaho↗

Hydrogeologic framework and groundwater/surface-water interactions of the South Fork Nooksack River Basin, northwestern Washington

A hydrogeologic framework of the South Fork (SF) Nooksack River Basin in northwestern Washington was developed and hydrologic data were collected to characterize the groundwater-flow system and its interaction with surface‑water features. In addition to domestic, agricultural, and commercial uses of groundwater within the SF Nooksack River Basin, groundwater has the potential to provide ecological benefits by maintaining late-summer streamflows and buffering stream temperatures. Cold-water refugia, created and maintained in part by groundwater, have been identified by water-resource managers as key elements to restore the health and viability of threatened salmonids in the SF Nooksack River. The SF Nooksack River drains a 183-square mile area of the North Cascades and the Puget Lowland underlain by unconsolidated glacial and alluvial sediments deposited over older sedimentary, metamorphic, and igneous bedrock. The primary aquifer that interacts with the SF Nooksack River was mapped within unconsolidated glacial outwash and alluvial sediment. The lower extent of this unit is bounded by bedrock and fine-grained, poorly sorted unconsolidated glaciomarine and glaciolacustrine sediments. In places, these deposits overlie and confine an aquifer within older glacial sediments. The extent and thickness of the hydrogeologic units were assembled from mapped geologic units and lithostratigraphic logs of field-inventoried wells. Generalized groundwater-flow directions within the surficial aquifer were interpreted from groundwater levels measured in August 2012; and groundwater seepage gains and losses to the SF Nooksack River were calculated from synoptic streamflow measurements made in the SF Nooksack River and its tributaries in September 2012. A subset of the field-inventoried wells was measured at a monthly interval to determine seasonal fluctuations in groundwater levels during water year 2013. Taken together, these data provide the foundation for a future groundwater-flow model of the SF Nooksack River Basin that may be used to investigate the potential effects of future climate change, land use, and groundwater pumping on water resources in the study area. Site-specific hydrologic data, including time series of longitudinal temperature profiles measured with a fiber-optic distributed temperature sensor and continuous monitoring of stream stage and water levels measured in wells in adjacent wetlands and aquifers, also were measured to characterize the interaction among the SF Nooksack River, surficial aquifers, and riparian wetlands.

Washington↗

Montana StreamStats—A method for retrieving basin and streamflow characteristics in Montana: Chapter A in Montana StreamStats

The U.S. Geological Survey (USGS) provides streamflow characteristics and other related information needed by water-resource managers to protect people and property from floods, plan and manage water-resource activities, and protect water quality. Streamflow characteristics provided by the USGS, such as peak-flow and low-flow frequencies for streamflow-gaging stations, are frequently used by engineers, flood forecasters, land managers, biologists, and others to guide their everyday decisions. In addition to providing streamflow characteristics at streamflow-gaging stations, the USGS also develops regional regression equations and drainage area-adjustment methods for estimating streamflow characteristics at locations on ungaged streams. Regional regression equations can be complex and often require users to determine several basin characteristics, which are physical and climatic characteristics of the stream and its drainage basin. Obtaining these basin characteristics for streamflow-gaging stations and ungaged sites traditionally has been time consuming and subjective, and led to inconsistent results. StreamStats is a Web-based geographic information system application that was created by the USGS to provide users with access to an assortment of analytical tools that are useful for water-resource planning and management. StreamStats allows users to easily obtain streamflow and basin characteristics for USGS streamflow-gaging stations and user-selected locations on ungaged streams. The USGS, in cooperation with Montana Department of Transportation, Montana Department of Environmental Quality, and Montana Department of Natural Resources and Conservation, completed a study to develop a StreamStats application for Montana, compute streamflow characteristics at streamflow-gaging stations, and develop regional regression equations to estimate streamflow characteristics at ungaged sites. Chapter A of this Scientific Investigations Report describes the Montana StreamStats application and the datasets, streamflow-gaging stations, streamflow characteristics, and regression equations (as described fully in Chapters B through G of this report) that are used for development of the StreamStats application for Montana.

Montana↗

Hydrogeology of Valley-Fill Aquifers and Adjacent Areas in Eastern Chemung County, New York

The extent, hydrogeologic framework, and potential well yields of valley-fill aquifers within a 151-square-mile area of eastern Chemung County, New York, were investigated, and the upland distribution of till thickness over bedrock was characterized. The hydrogeologic framework of these valleyfill aquifers was interpreted from multiple sources of surficial and subsurface data and an interpretation of the origin of the glacial deposits, particularly during retreat of glacial ice from the region. Potential yields of screened wells are based on the hydrogeologic framework interpretation and existing well-yield data, most of which are from wells finished with open-ended well casing. Water-resource potential is greatest within saturated sand and gravel in the Chemung River valley (nearly 1 mile wide), especially where induced infiltration of additional water from the Chemung River is possible. The second most favorable area is the Newtown Creek valley at the confluence of Newtown Creek with North Branch Newtown Creek east of Horseheads, N.Y. Extensive sand and gravel deposits within the Breesport, N.Y., area are largely unsaturated but may have greater saturation along the east side of Jackson Creek immediately north of Breesport. Till deposits confine sand and gravel along Newtown Creek at Erin, N.Y., and along much of the upper reach of North Branch Newtown Creek; this confining unit may limit recharge and potential well yield. The north-south oriented valleys of Baldwin and Wynkoop Creeks end at notched divides that imply input of glacial meltwater and limited sediment from outside of the present watersheds. These two valleys are relatively narrow but contain variably sorted sand and gravel, which, in places, may be capable of supplying modest-size community water systems.

New York↗

Simulation of daily streamflow for nine river basins in eastern Iowa using the Precipitation-Runoff Modeling System

The U.S. Geological Survey, in cooperation with the Iowa Department of Natural Resources, constructed Precipitation-Runoff Modeling System models to estimate daily streamflow for nine river basins in eastern Iowa that drain into the Mississippi River. The models are part of a suite of methods for estimating daily streamflow at ungaged sites. The Precipitation-Runoff Modeling System is a deterministic, distributed- parameter, physical-process-based modeling system developed to evaluate the response of streamflow and general drainage basin hydrology to various combinations of climate and land use. Calibration and validation periods used in each basin mostly were October 1, 2002, through September 30, 2012, but differed depending on the period of record available for daily mean streamflow measurements at U.S. Geological Survey streamflow-gaging stations. A geographic information system tool was used to delineate each basin and estimate values for model parameters based on basin physical and geographical features. A U.S. Geological Survey auto-calibration tool that uses a shuffled complex evolution algorithm was used for initial calibration, and then manual modifications were made to parameter values to complete the calibration of each basin model. The main objective of the calibration was to match daily discharge values of simulated streamflow to measured daily discharge values. The accuracy of Precipitation-Runoff Modeling System model streamflow estimates of nine river basins in eastern Iowa as compared to measured values at U.S. Geological Survey streamflow-gaging stations varied. The Precipitation-Runoff Modeling System models of nine river basins in eastern Iowa were satisfactory at estimating daily streamflow at 57 of the 79 calibration sites and 13 of the 14 validation sites based on statistical results. Unsatisfactory performance can be contributed to several factors: (1) low flow, no flow, and flashy flow conditions in headwater subbasins having a small drainage area; (2) poor representation of the groundwater and storage components of flow within a basin; (3) lack of accounting for basin withdrawals and water use; and (4) the availability and accuracy of meteorological input data. The Precipitation- Runoff Modeling System models of nine river basins in eastern Iowa will provide water-resource managers with a consistent and documented method for estimating streamflow at ungaged sites and aid in environmental studies, hydraulic design, water management, and water-quality projects.

Iowa↗

Hydrologic data and groundwater-flow simulations in the Brown Ditch Watershed, Indiana Dunes National Lakeshore, near Beverly Shores and Town of Pines, Indiana

The U.S. Geological Survey (USGS) collected data and simulated groundwater flow to increase understanding of the hydrology and the effects of drainage alterations on the water table in the vicinity of Great Marsh, near Beverly Shores and Town of Pines, Indiana. Prior land-management practices have modified drainage and caused changes in the distribution of open water, streams and ditches, and groundwater abundance and flow paths. Collected hydrologic data indicate that the majority of water entering Great Marsh flows from the southern dune ridge beneath Town of Pines, Indiana. Groundwater flow is intercepted by Brown Ditch in the eastern portion of the study area and Derby Ditch in the western portion of the study area. A smaller amount of groundwater from the northern dune ridge beneath Beverly Shores also contributed water to Great Marsh. Continuous groundwater-level data collected indicate that the predominant north-south groundwater-flow gradients vary during the course of the year due to increased levels of precipitation or during periods of drainage obstructions. Continuous surface-water discharge and surface-water elevation were measured at three USGS streamgages, one each on Brown, Kintzele and Derby Ditches. The monthly mean discharge statistics indicate that during the period of record— June 2012 to September 2013—streamflow in Kintzele Ditch was lowest during July 2012 and highest during April 2013. In Derby Ditch, streamflow also was lowest during July 2012 and highest during April 2013. Periods of relatively high and low groundwater levels during August 1982, March 2013, and April 2014 were examined and simulated by using MODFLOW and companion software. Results from the simulation of conditions during March 2013 include that nearly 100 percent of all water entering the area simulating Town of Pines is from recharge. Of all the water simulated to enter the eastern and western portions of Great Marsh, nearly 20 and 18 percent, respectively, flows from Town of Pines to the western and eastern portions of Great Marsh. The dune ridges beneath Town of Pines and to a lesser extent beneath Beverly Shores are a major source of recharge to the surficial aquifer and Great Marsh. Results from the simulation of the conditions of April 2014 include that, despite increases in the amount of water entering Great Marsh due to a beaver-dam-modified hydrologic condition, there is still virtually zero simulated groundwater flow from Great Marsh to Town of Pines. The volume of water simulated to be entering the zone representing Beverly Shores decreased by 0.43 cubic foot per second from the results of the March 2013 simulation. This simulated difference in water budgets can be attributed to increased simulated recharge in Great Marsh and Town of Pines. Effects of the inclusion of the beaver dam included the increase of the simulated water table and simulated inundated area upstream of the beaver dam due to the effects of ponding surface water. Results from the simulation scenario that includes six proposed pool-riffle control structures in Brown Ditch under the hydrologic conditions of March 2013 indicate areas inundated by water are larger, including areas just to the north of the entrance of Brown Ditch into Great Marsh, and areas north of the confluence of Brown and Kintzele Ditches. Results from the scenario simulating the increase of the Lake Michigan water level to the historical high of May 31, 1998, showed inundated areas of Great Marsh south of Beverly Shores enlarged on both sides of Lakeshore County Road with the greatest enlargement simulated to be southeast of the intersection of Lakeshore County Road and Beverly Drive. For the scenario simulating the decrease of the Lake Michigan water level to the historical low of December 23, 2007, results show little change from the original March 2013 inundated area. The results of this study can be used by water-resource managers to understand how surrounding ditches affect water levels in Great Marsh and other inland wetlands and residential areas. The groundwater model developed can be applied to answer questions about how alterations to the drainage system in the area affects water levels in the public and residential areas surrounding Great Marsh. The modeling methods developed in this study provide a template for other studies of groundwater flow and groundwater/surface-water interactions within the shallow surficial aquifer in northern Indiana, and in similar hydrologic settings that include surficial sand aquifers in coastal areas.

Indiana↗

Water quality and hydrology of Silver Lake, Oceana County, Michigan, with emphasis on lake response to nutrient loading

Executive Summary Silver Lake is a 672-acre inland lake located in Oceana County, Michigan, and is a major tourist destination due to its proximity to Lake Michigan and the surrounding outdoor recreational opportunities. In recent years, Silver Lake exhibited patterns of high phosphorus concentrations, elevated chlorophyll a concentrations, and nuisance algal blooms. The U.S. Geological Survey (USGS), in cooperation with the Silver Lake Improvement Board and in collaboration with the Annis Water Resources Institute (AWRI) of Grand Valley State University, designed a study to assess the hydrologic and nutrient inputs to Silver Lake in order to identify the events and conditions that affect the nutrient chemistry and production of algal blooms in the lake. This information can inform water-resource managers in developing various management strategies to prevent or reduce the occurrence of future algal blooms. USGS and AWRI scientists collected data from November 2012 to December 2014 to provide information for future management decisions for Silver Lake. Silver Lake can be classified as a polymictic lake and has a residence time of approximately 223 days. Based on the mean lake Secchi depth, total phosphorus, and total nitrogen concentrations, Silver Lake is classified as a eutrophic lake. In-situ bioassay results indicate that algal growth in Silver Lake is colimited by both nitrogen and phosphorus. The nutrient budget for Silver Lake was calculated using the BATHTUB model based on 2 years of water-quality data collection. The BATHTUB model, developed by the U.S. Army Corps of Engineers, treats the lake as a well-mixed system with multiple inputs and outlets for both water and dissolved constituents, such as nutrients. Based on results of the BATHTUB model, which were conditioned on observed concentrations and flows, the mean annual input of phosphorus to Silver Lake was approximately 1,342 pounds (lb); the mean annual input of nitrogen to Silver Lake was approximately 51,998 lb. The major measured sources of phosphorus loading to Silver Lake were groundwater and Hunter Creek, whereas the major measured sources of nitrogen to Silver Lake were Hunter Creek, groundwater, and atmospheric deposition. The largest loading of phosphorus and nitrogen to Silver Lake occurred during the spring. Minimal phosphorus deposition (if any) occurred in the lakebed sediment; however, of the nitrogen that entered Silver Lake, approximately 42.2 percent was deposited in the lakebed sediment as simulated by the BATHTUB model. In addition to measured sources, a septic load model was used to estimate the likely range of septic contribution to groundwater and adjacent surface waters. The likely septic loading scenario estimates that septic systems contribute 47.8 percent of the phosphorus to groundwater and 22.3 percent of phosphorus to Hunter Creek. These results indicate that septic systems are a major source of phosphorus loading to Silver Lake. The likely septic loading scenario indicated that septic systems account for 0.95 percent of the nitrogen load to Hunter Creek and 1.1 percent of the contribution of nitrogen to groundwater. The BATHTUB model was used to estimate future nutrient loading and eutrophication scenarios based on water-quality data collected from Silver Lake, groundwater, major tributaries, and atmospheric deposition. A separate septic load model was used to estimate the septic contribution to groundwater or directly to surface water, and the nutrient load estimates were modeled using the BATHTUB model to determine subsequent water-quality changes to Silver Lake. BATHTUB model scenarios based on measured data: The first BATHTUB scenario evaluated the condition of Silver Lake and the change to lake water quality (trophic status) as a result of changes in nutrient loading from different sources. Based on BATHTUB model simulations, if groundwater loading of phosphorus and nitrogen only were decreased by 75 percent, and all of the other nutrient inputs stayed the same, the future condition of Silver Lake would most likely remain highly mesotrophic to eutrophic (the current [2014] condition of Silver Lake). If nutrient loading continued to increase in groundwater, the lake would continue to remain eutrophic with more frequent algal blooms. If nutrient loading from Hunter Creek only decreased by 50–75 percent, and all of the other nutrient inputs stayed the same as the baseline dataset, Silver Lake would remain eutrophic to highly mesotrophic. By reducing the input of manageable nutrient sources (Hunter Creek, groundwater, and lawn runoff) by 75 percent, the BATHTUB model simulation indicates that Silver Lake would be classified as mesotrophic, which is indicative of improved water quality, water clarity, and reduced algal bloom frequency. Simulations also were run using the BATHTUB model to evaluate the number of days Silver Lake could experience algal blooms (algal blooms are defined as modeled chlorophyll a in excess of 10 micrograms per liter [µg/L]) as a result of an increase/decrease in phosphorus and nitrogen loading from groundwater, Hunter Creek, and (or) a combination of sources. If the phosphorus and nitrogen loading from Hunter Creek is decreased (and all other sources are not altered), Silver Lake will continue to experience algal blooms, but less frequently than what is currently experienced. The same scenario holds true if the nutrient loading from groundwater is decreased. Another scenario was simulated using a combination of sources, which includes increases and decreases in phosphorus and nitrogen loading from sources that are the most likely to be managed, and includes groundwater (as a result of conversion of household septic to sewers), Hunter Creek (conversion of household septic to sewers), and lawn runoff. Results of the BATHTUB model indicated that a 50-percent reduction of phosphorus and nitrogen from these sources would result in a considerable decrease in algal bloom frequency (from 231 to 132 days) and severity, and a 75-percent reduction would greatly reduce algal bloom occurrence on Silver Lake (from 231 to 57 days). BATHTUB model scenarios based on septic load model: A scenario also was conducted using the BATHTUB model to simulate the conversion of septic to sewer and included a low, high, and medium (likely) scenario of nutrient loading to Silver Lake. Simulations of the BATHTUB model indicated that, under the likely scenario, the conversion of all onsite septic treatment to sewers would result in an overall change in lake trophic status from eutrophic to mesotrophic, thereby reducing the frequency of algal blooms and algal bloom intensity on Silver Lake (chlorophyll a >10 µg/L, from 231 to 184 days per year, or chlorophyll a >20 µg/L, from 80 to 49 days per year).

Michigan↗

Water use in the Apalachicola-Chattahoochee-Flint River Basin, Alabama, Florida, and Georgia, 2010, and water-use trends, 1985-2010

The Apalachicola-Chattahoochee-Flint (ACF) River Basin encompasses about 20,230 square miles in parts of Alabama, Florida, and Georgia. Increasing population growth and agricultural production from the 1970s to 2010 has prompted increases in water-resources development and substantially increased water demand in the basin. Since the 1980s, Alabama, Florida, Georgia, and the U.S. Army Corps of Engineers are parties to litigation concerning water management in the ACF River Basin. Estimating the 2010 water use in the ACF River Basin is one aspect of a multipart water resources study on the ACF River Basin that began in 2011. This ACF River Basin study is one focus area of the U.S. Geological Survey’s National Water Census program. The 2010 water-use estimates for the ACF River Basin are presented in this report. These estimates include an inventory of the quantity and sources of water withdrawn by category of use and location (State and river basin), and the surface-water returns in the ACF River Basin during 2010. Water-use trends from 1985 to 2010 in the basin also are presented. Offstream water-withdrawal data in the ACF River Basin are presented for each of the following categories: public supply, self-supplied domestic, self-supplied commercial, industrial, mining, agricultural (including crop irrigation, livestock, and aquaculture uses), and thermoelectric-power generation. Water-use data are compiled for the 14 subbasins in the ACF River Basin. For the counties in Alabama, Florida, and Georgia that are partially within the ACF River Basin, data are presented for only that part of the county that lies within the basin. A variety of Federal, State, local, private, and online sources in Alabama, Florida, and Georgia were used to gather surface-water and groundwater withdrawal, surface-water discharges (return flows), and water-use data for the ACF River Basin in 2010. The population in the ACF River Basin was 3.835 million in 2010, a 45-percent increase from the 1990 population of nearly 2.636 million. About 92 percent of the 2010 ACF population resided in Georgia with nearly 75 percent living in the Atlanta metropolitan area. In 2010, 1,645 million gallons per day (Mgal/d) of water were withdrawn from groundwater (576 Mgal/d) and surface-water (1,069 Mgal/d) sources in the ACF River Basin. About 89 percent of the groundwater and 83 percent of the surface-water withdrawals were from Georgia. About 5.6 percent of the total groundwater and nearly 4 percent of the total surface-water withdrawals in the ACF River Basin were from Florida, whereas about 5.3 percent of groundwater and nearly 16 percent of surface water were withdrawn in Alabama. Total water use (withdrawals plus public-supplied deliveries) in the ACF River Basin was 1,593 Mgal/d in 2010. About 56 Mgal/d of water withdrawn in the ACF River Basin was delivered (interbasin transfer) to basins beyond the ACF River Basin. About 564 Mgal/d of water was returned to surface-water bodies in the ACF River Basin. Most of that amount, 63 percent, was treated wastewater discharged by public wastewater-treatment facilities. Water used for once-through cooling by thermoelectric-power facilities accounted for nearly 24 percent of the surface-water returns in the basin. About 70 percent of all water withdrawals in the ACF River Basin were by self-supplied agricultural water users and public water suppliers. Agricultural withdrawals were greatest in the Flint River Basin (501 Mgal/d) with ground-water representing 84 percent of the withdrawals from that basin. Within the Flint River Basin, agricultural withdrawals were greatest in the Lower Flint River and Spring Creek subbasins. About 3.52 million people were served by public water suppliers in the ACF River Basin during 2010, and 88 percent of that population used surface water. Georgia had the largest public-supplied population, representing nearly 93 percent (3.17 million) of the public-supplied population in the ACF River Basin. Public water suppliers served 193,700 people (5.7 percent) in Alabama and 31,880 people in Florida (1.3 percent). Public-supply losses were estimated at 101 Mgal/d. Withdrawals for public supply (483 Mgal/d) and self-supplied industry (141 Mgal/d) were greatest in the Chattahoochee River Basin. Surface water accounted for 96 percent of all withdrawals in the Chattahoochee River Basin. Withdrawals for public supply were greatest in the Upper Chattahoochee River subbasin (366 Mgal/d), whereas self-supplied industrial withdrawals were greatest in the Lower Chattahoochee River subbasin (110 Mgal/d). Water-use trends in the ACF River Basin have varied during the 25 years between 1985 and 2010. Surface-water withdrawals declined between 1985 and 2000, sharply increased in 2000, and declined again between 2000 and 2010. In contrast, groundwater withdrawals increased between 1985 and 2000, declined in 2005, and increased between 2005 and 2010.

Alabama, Florida, Georgia↗

Potential effects of alterations to the hydrologic system on the distribution of salinity in the Biscayne aquifer in Broward County, Florida

To address concerns about the effects of water-resource management practices and rising sea level on saltwater intrusion, the U.S. Geological Survey in cooperation with the Broward County Environmental Planning and Community Resilience Division, initiated a study to examine causes of saltwater intrusion and predict the effects of future alterations to the hydrologic system on salinity distribution in eastern Broward County, Florida. A three-dimensional, variable-density solute-transport model was calibrated to conditions from 1970 to 2012, the period for which data are most complete and reliable, and was used to simulate historical conditions from 1950 to 2012. These types of models are typically difficult to calibrate by matching to observed groundwater salinities because of spatial variability in aquifer properties that are unknown, and natural and anthropogenic processes that are complex and unknown; therefore, the primary goal was to reproduce major trends and locally generalized distributions of salinity in the Biscayne aquifer. The methods used in this study are relatively new, and results will provide transferable techniques for protecting groundwater resources and maximizing groundwater availability in coastal areas. The model was used to (1) evaluate the sensitivity of the salinity distribution in groundwater to sea-level rise and groundwater pumping, and (2) simulate the potential effects of increases in pumping, variable rates of sea-level rise, movement of a salinity control structure, and use of drainage recharge wells on the future distribution of salinity in the aquifer. Results from the simulation of historical conditions indicate that the model generally represents the observed greater westward extent of elevated salinity in the central part of the intruded area relative to the northern and southernmost parts of the intruded area. Results of sensitivity testing indicate that the extent of elevated salinity is most sensitive to pumping in areas where the source of saltwater is largely offshore, from the Atlantic Ocean, and is most sensitive to sea-level rise in areas where the source of salinity is downward leakage of brackish water from canals. Simulations of future scenarios indicate that increases in pumping near the existing interface may cause the interface to advance and decreases in pumping may cause it to retreat. Climatic effects, such as periods of prolonged drought or high precipitation, may augment or counteract long-term effects of changes in pumping on aquifer salinity at well fields. With increasing rates of sea-level rise, the freshwater-saltwater interface advances progressively inland, and flow-averaged salinities at well fields near the existing interface increase commensurately. Hypothetical southeastward (downstream) re-positioning of the existing G–54 salinity-control structure may prevent the interface from moving northwestward along and near the North New River canal, but beneficial effects are localized. Implementation of freshwater recharge wells in the city of Hallandale Beach may also have only a localized freshening effect in the aquifer and little appreciable effect on the freshwater-saltwater interface or on concentrations of salinity at well fields. Model accuracy and use are limited by uncertainty in the physical properties and boundary conditions of the system, uncertainty in historical and future conditions, and generalizations made in the mathematical relationships used to describe the physical processes of groundwater flow and transport. Because of these limitations, model results should be considered in relative rather than absolute terms. Nonetheless, model results do provide useful information on the relative scale of response of the system to changes in pumping distribution, sea-level rise, and mitigation activities.

Florida↗

Estimating flood magnitude and frequency at gaged and ungaged sites on streams in Alaska and conterminous basins in Canada, based on data through water year 2012

Estimates of the magnitude and frequency of floods are needed across Alaska for engineering design of transportation and water-conveyance structures, flood-insurance studies, flood-plain management, and other water-resource purposes. This report updates methods for estimating flood magnitude and frequency in Alaska and conterminous basins in Canada. Annual peak-flow data through water year 2012 were compiled from 387 streamgages on unregulated streams with at least 10 years of record. Flood-frequency estimates were computed for each streamgage using the Expected Moments Algorithm to fit a Pearson Type III distribution to the logarithms of annual peak flows. A multiple Grubbs-Beck test was used to identify potentially influential low floods in the time series of peak flows for censoring in the flood frequency analysis. For two new regional skew areas, flood-frequency estimates using station skew were computed for stations with at least 25 years of record for use in a Bayesian least-squares regression analysis to determine a regional skew value. The consideration of basin characteristics as explanatory variables for regional skew resulted in improvements in precision too small to warrant the additional model complexity, and a constant model was adopted. Regional Skew Area 1 in eastern-central Alaska had a regional skew of 0.54 and an average variance of prediction of 0.45, corresponding to an effective record length of 22 years. Regional Skew Area 2, encompassing coastal areas bordering the Gulf of Alaska, had a regional skew of 0.18 and an average variance of prediction of 0.12, corresponding to an effective record length of 59 years. Station flood-frequency estimates for study sites in regional skew areas were then recomputed using a weighted skew incorporating the station skew and regional skew. In a new regional skew exclusion area outside the regional skew areas, the density of long-record streamgages was too sparse for regional analysis and station skew was used for all estimates. Final station flood frequency estimates for all study streamgages are presented for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities. Regional multiple-regression analysis was used to produce equations for estimating flood frequency statistics from explanatory basin characteristics. Basin characteristics, including physical and climatic variables, were updated for all study streamgages using a geographical information system and geospatial source data. Screening for similar-sized nested basins eliminated hydrologically redundant sites, and screening for eligibility for analysis of explanatory variables eliminated regulated peaks, outburst peaks, and sites with indeterminate basin characteristics. An ordinary least‑squares regression used flood-frequency statistics and basin characteristics for 341 streamgages (284 in Alaska and 57 in Canada) to determine the most suitable combination of basin characteristics for a flood-frequency regression model and to explore regional grouping of streamgages for explaining variability in flood-frequency statistics across the study area. The most suitable model for explaining flood frequency used drainage area and mean annual precipitation as explanatory variables for the entire study area as a region. Final regression equations for estimating the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probability discharge in Alaska and conterminous basins in Canada were developed using a generalized least-squares regression. The average standard error of prediction for the regression equations for the various annual exceedance probabilities ranged from 69 to 82 percent, and the pseudo-coefficient of determination (pseudo-R 2 ) ranged from 85 to 91 percent. The regional regression equations from this study were incorporated into the U.S. Geological Survey StreamStats program for a limited area of the State—the Cook Inlet Basin. StreamStats is a national web-based geographic information system application that facilitates retrieval of streamflow statistics and associated information. StreamStats retrieves published data for gaged sites and, for user-selected ungaged sites, delineates drainage areas from topographic and hydrographic data, computes basin characteristics, and computes flood frequency estimates using the regional regression equations.

Alaska↗

Groundwater-flow model for the Wood River Valley aquifer system, south-central Idaho

A three-dimensional numerical model of groundwater flow was developed for the Wood River Valley (WRV) aquifer system, Idaho, to evaluate groundwater and surface-water availability at the regional scale. This mountain valley is located in Blaine County and has a drainage area of about 2,300 square kilometers (888 square miles). The model described in this report can serve as a tool for water-rights administration and water-resource management and planning. The model was completed with support from the Idaho Department of Water Resources, and is part of an ongoing U.S. Geological Survey effort to characterize the groundwater resources of the WRV. A highly reproducible approach was taken for constructing the WRV groundwater-flow model. The collection of datasets, source code, and processing instructions used to construct and analyze the model was distributed as an R statistical-computing and graphics package. Flow in the WRV aquifer was simulated using the MODFLOW-USG groundwater flow model. The transient flow model simulates groundwater flow between 1995 and 2010. The model uses a 100-meter (328-feet) uniform grid spacing with 54,922 active model cells distributed over three model layers. A confining unit in the south-central part of the Bellevue fan necessitated the use of a multi-layer model. Specified-flow boundaries were used to simulate the groundwater inflows from each of the major tributary basins (also known as tributary basin underflow) and the areal recharge of precipitation and applied irrigation. Head‑dependent flow boundaries were used to simulate the stream-aquifer flow exchange in river reaches and the groundwater discharge at the outlet boundaries of Stanton Crossing and Silver Creek. The model was calibrated by adjusting aquifer hydraulic properties to match simulated and measured water levels and stream-aquifer flow exchange, using the parameter-estimation program PEST. The model reasonably simulated the measured water-table elevation, orientation, and gradients. Stream-aquifer flow exchange along river reaches also was reasonably simulated by the model. Inflow into the WRV aquifer system originates from three sources (from largest to smallest): Streamflow loss from the Big Wood River and Silver Creek, Areal recharge of precipitation and applied irrigation, and Tributary basin underflow. Outflow from the WRV aquifer system originates from five sources (from largest to smallest): Aquifer discharge into the Big Wood River and Silver Creek, Production-well pumping, Discharge across the water table into the vadose zone, Subsurface outflow beneath Silver Creek near Picabo, and Subsurface outflow beneath the Big Wood River near Stanton Crossing. Temporal changes in aquifer storage are most affected by areal recharge and groundwater pumping, and also contribute to changes in streamflow gains.

Idaho↗

Effects of water-supply reservoirs on streamflow in Massachusetts

State and local water-resource managers need modeling tools to help them manage and protect water-supply resources for both human consumption and ecological needs. The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, has developed a decision-support tool to estimate the effects of reservoirs on natural streamflow. The Massachusetts Reservoir Simulation Tool is a model that simulates the daily water balance of a reservoir. The reservoir simulation tool provides estimates of daily outflows from reservoirs and compares the frequency, duration, and magnitude of the volume of outflows from reservoirs with estimates of the unaltered streamflow that would occur if no dam were present. This tool will help environmental managers understand the complex interactions and tradeoffs between water withdrawals, reservoir operational practices, and reservoir outflows needed for aquatic habitats. A sensitivity analysis of the daily water balance equation was performed to identify physical and operational features of reservoirs that could have the greatest effect on reservoir outflows. For the purpose of this report, uncontrolled releases of water (spills or spillage) over the reservoir spillway were considered to be a proxy for reservoir outflows directly below the dam. The ratio of average withdrawals to the average inflows had the largest effect on spillage patterns, with the highest withdrawals leading to the lowest spillage. The size of the surface area relative to the drainage area of the reservoir also had an effect on spillage; reservoirs with large surface areas have high evaporation rates during the summer, which can contribute to frequent and long periods without spillage, even in the absence of water withdrawals. Other reservoir characteristics, such as variability of inflows, groundwater interactions, and seasonal demand patterns, had low to moderate effects on the frequency, duration, and magnitude of spillage. The reservoir simulation tool was used to simulate 35 single- and multiple-reservoir systems in Massachusetts over a 44-year period (water years 1961 to 2004) under two water-use scenarios. The no-pumping scenario assumes no water withdrawal pumping, and the pumping scenario incorporates average annual pumping rates from 2000 to 2004. By comparing the results of the two scenarios, the total streamflow alteration can be parsed into the portion of streamflow alteration caused by the presence of a reservoir and the additional streamflow alteration caused by the level of water use of the system. For each reservoir system, the following metrics were computed to characterize the frequency, duration, and magnitude of reservoir outflow volumes compared with unaltered streamflow conditions: (1) the median number of days per year in which the reservoir did not spill, (2) the median duration of the longest consecutive period of no-spill days per year, and (3) the lowest annual flow duration exceedance probability at which the outflows are significantly different from estimated unaltered streamflow at the 95-percent confidence level. Most reservoirs in the study do not spill during the summer months even under no-pumping conditions. The median number of days during which there was no spillage was less than 365 for all reservoirs in the study, indicating that, even under reported pumping conditions, the reservoirs refill to full volume and spill at least once during nondrought years, typically in the spring. Thirteen multiple-reservoir systems consisting of two or three hydrologically connected reservoirs were included in the study. Because operating rules used to manage multiple-reservoir systems are not available, these systems were simulated under two pumping scenarios, one in which water transfers between reservoirs are minimal and one in which reservoirs continually transferred water to intermediate or terminal reservoirs. These two scenarios provided upper and lower estimates of spillage under average pumping conditions from 2000 to 2004. For sites with insufficient data to simulate daily water balances, a proxy method to estimate the three spillage metrics was developed. A series of 4,000 Monte Carlo simulations of the reservoir water balance were run. In each simulation, streamflow, physical reservoir characteristics, and daily climate inputs were randomly varied. Tobit regression equations that quantify the relation between streamflow alteration and physical and operational characteristics of reservoirs were developed from the results of the Monte Carlo simulations and can be used to estimate each of the three spillage metrics using only the withdrawal ratio and the ratio of the surface area to the drainage area, which are available statewide for all reservoirs. A graphical user-interface for the Massachusetts Reservoir Simulation Tool was developed in a Microsoft Access environment. The simulation tool contains information for 70 reservoirs in Massachusetts and allows for simulation of additional scenarios than the ones considered in this report, including controlled releases, dam seepage and leakage, demand management plans, and alternative water withdrawal and transfer rules.

Massachusetts↗

Simulation and assessment of groundwater flow and groundwater and surface-water exchanges in lakes of the northeast Twin Cities Metropolitan Area, Minnesota, 2003 through 2013: Chapter B of Water levels and groundwater and surface-water exchanges in lakes of the northeast Twin Cities Metropolitan Area, Minnesota, 2002 through 2015

Water levels during 2003 through 2013 were less than mean water levels for the period 1925–2013 for several lakes in the northeast Twin Cities Metropolitan Area in Minnesota. Previous periods of low lake-water levels generally were correlated with periods with less than mean precipitation. Increases in groundwater withdrawals and land-use changes have brought into question whether or not recent (2003–13) lake-water-level declines are solely caused by decreases in precipitation. A thorough understanding of groundwater and surface-water exchanges was needed to assess the effect of water-management decisions on lake-water levels. To address this need, the U.S. Geological Survey, in cooperation with the Metropolitan Council and the Minnesota Department of Health, developed and calibrated a three-dimensional, steady-state groundwater-flow model representing 2003–13 mean hydrologic conditions to assess groundwater and lake-water exchanges, and the effects of groundwater withdrawals and precipitation on water levels of 96 lakes in the northeast Twin Cities Metropolitan Area. Lake-water budgets for the calibrated groundwater-flow model indicated that groundwater is flowing into lakes in the northeast Twin Cities Metropolitan Area and lakes are providing water to underlying aquifers. Lake-water outflow to the simulated groundwater system was a major outflow component for Big Marine Lake, Lake Elmo, Snail Lake, and White Bear Lake, accounting for 45 to 64 percent of the total outflows from the lakes. Evaporation and transpiration from the lake surface ranged from 19 to 52 percent of the total outflow from the four lakes. Groundwater withdrawals and precipitation were varied from the 2003‒13 mean values used in the calibrated model (30-percent changes in groundwater withdrawals and 5-percent changes in precipitation) for hypothetical scenarios to assess the effects of groundwater withdrawals and precipitation on water budgets and levels in Big Marine Lake, Snail Lake, and White Bear Lake. Simulated lake-water levels and budgets for Snail Lake and White Bear Lake were affected by 30-percent changes in groundwater withdrawals and 5-percent changes in precipitation in the area, whereas the water level in Big Marine Lake was mainly affected by 5-percent precipitation changes. The effects of groundwater withdrawals on the lake-water levels depend on the number of wells and amount of withdrawals from wells near the lakes. Although lake-water levels are sensitive to precipitation changes, increases in groundwater withdrawals during dry periods exacerbate lake-water level declines. The calibrated, groundwater-flow model is a tool that water-resources managers can use to address future water management issues in the northeast Twin Cities Metropolitan Area.

Minnesota↗

Groundwater-flow model of the northern High Plains aquifer in Colorado, Kansas, Nebraska, South Dakota, and Wyoming

The High Plains aquifer is a nationally important water resource underlying about 175,000 square miles in parts of eight states: Colorado, Kansas, Oklahoma, Nebraska, New Mexico, South Dakota, Texas, and Wyoming. Droughts across much of the Northern High Plains from 2001 to 2007 have combined with recent (2004) legislative mandates to elevate concerns regarding future availability of groundwater and the need for additional information to support science-based water-resource management. To address these needs, the U.S. Geological Survey began the High Plains Groundwater Availability Study to provide a tool for water-resource managers and other stakeholders to assess the status and availability of groundwater resources. A transient groundwater-flow model was constructed using the U.S. Geological Survey modular three-dimensional finite-difference groundwater-flow model with Newton-Rhapson solver (MODFLOW–NWT). The model uses an orthogonal grid of 565 rows and 795 columns, and each grid cell measures 3,281 feet per side, with one variably thick vertical layer, simulated as unconfined. Groundwater flow was simulated for two distinct periods: (1) the period before substantial groundwater withdrawals, or before about 1940, and (2) the period of increasing groundwater withdrawals from May 1940 through April 2009. A soil-water-balance model was used to estimate recharge from precipitation and groundwater withdrawals for irrigation. The soil-water-balance model uses spatially distributed soil and landscape properties with daily weather data and estimated historical land-cover maps to calculate spatial and temporal variations in potential recharge. Mean annual recharge estimated for 1940–49, early in the history of groundwater development, and 2000–2009, late in the history of groundwater development, was 3.3 and 3.5 inches per year, respectively. Primary model calibration was completed using statistical techniques through parameter estimation using the parameter estimation suite of software with Tikhonov regularization. Calibration targets for the groundwater model included 343,067 groundwater levels measured in wells and 10,820 estimated monthly stream base flows at streamgages. A total of 1,312 parameters were adjusted during calibration to improve the match between calibration targets and simulated equivalents. Comparison of calibration targets to simulated equivalents indicated that, at the regional scale, the model correctly reproduced groundwater levels and stream base flows for 1940–2009. This comparison indicates that the model can be used to examine the likely response of the aquifer system to potential future stresses. Mean calibrated recharge for 1940–49 and 2000–2009 was smaller than that estimated with the soil-water-balance model. This indicated that although the general spatial patterns of recharge estimated with the soil-water-balance model were approximately correct at the regional scale of the Northern High Plains aquifer, the soil-water-balance model had overestimated recharge, and adjustments were needed to decrease recharge to improve the match of the groundwater model to calibration targets. The largest components of the simulated groundwater budgets were recharge from precipitation, recharge from canal seepage, outflows to evapotranspiration, and outflows to stream base flow. Simulated outflows to irrigation wells increased from 7 percent of total outflows in 1940–49 to 38 percent of 1970–79 total outflows and 49 percent of 2000–2009 total outflows.

Colorado, Kansas, Nebraska, South Dakota, Wyoming↗

Development and evaluation of a reservoir model for the Chain of Lakes in Illinois

Forecasts of flows entering and leaving the Chain of Lakes reservoir on the Fox River in northeastern Illinois are critical information to water-resource managers who determine the optimal operation of the dam at McHenry, Illinois, to help minimize damages to property and loss of life because of flooding on the Fox River. In 2014, the U.S. Geological Survey; the Illinois Department of Natural Resources, Office of Water Resources; and National Weather Service, North Central River Forecast Center began a cooperative study to develop a system to enable engineers and planners to simulate and communicate flows and to prepare proactively for precipitation events in near real time in the upper Fox River watershed. The purpose of this report is to document the development and evaluation of the Chain of Lakes reservoir model developed in this study. The reservoir model for the Chain of Lakes was developed using the Hydrologic Engineering Center–Reservoir System Simulation program. Because of the complex relation between the dam headwater and reservoir pool elevations, the reservoir model uses a linear regression model that relates dam headwater elevation to reservoir pool elevation. The linear regression model was developed using 17 U.S. Geological Survey streamflow measurements, along with the gage height in the reservoir pool and the gage height at the dam headwater. The Nash-Sutcliffe model efficiency coefficients for all three linear regression model variables ranged from 0.90 to 0.98. The reservoir model performance was evaluated by graphically comparing simulated and observed reservoir pool elevation time series during nine periods of high pool elevation. In addition, the peak elevations during these time periods were graphically compared to the closest-in-time observed pool elevation peak. The mean difference in the simulated and observed peak elevations was -0.03 feet, with a standard deviation of 0.19 feet. The Nash-Sutcliffe coefficient for peak prediction was calculated as 0.94. Evaluation of the model based on accuracy of peak prediction and the ability to simulate an elevation time series showed the performance of the model was satisfactory.

Illinois↗

Hydrogeologic framework and hydrologic conditions of the Piney Point aquifer in Virginia

The Piney Point aquifer in Virginia is newly described and delineated as being composed of six geologic units, in a study conducted by the U.S. Geological Survey in cooperation with the Virginia Department of Environmental Quality (VA DEQ). The eastward-dipping geologic units include, in stratigraphically ascending order, the sand of the Nanjemoy Formation Woodstock Member, interbedded limestone and sand of the Piney Point Formation, silty and clayey sand of the Gosport Formation equivalent sediments, silty sand of the Oligocene-age sediments, silty fine-grained sand of the Old Church Formation, and silty sand of the Calvert Formation, Newport News unit and basal Plum Point Member. Identification of geologic units is based on typical sediment lithologies of geologic formations. Fine-grained sediments that compose confining units positioned immediately above and below the Piney Point aquifer are also described. The Piney Point aquifer is one of several confined aquifers within the Virginia Coastal Plain and includes a highly porous and solution-channeled indurated limestone within the Piney Point Formation from which withdrawals are made. The limestone is relatively continuous laterally across central parts of the Northern Neck, Middle Peninsula, and York-James Peninsula. Other geologic units are of variable extent. The configurations of most of the geologic units are further affected by newly identified faults that are aligned radially from the Chesapeake Bay impact crater and create constrictions or barriers to groundwater flow. Some geologic units are also truncated beneath the lower Rappahannock River by a resurge channel associated with the impact crater. Groundwater withdrawals from the Piney Point aquifer increased from approximately 1 million gallons per day (Mgal/d) during 1900 to 7.35 Mgal/d during 2004. As a result, a water-level cone of depression in James City and northern York Counties was estimated to be as low as 70 feet (ft) below the National Geodetic Vertical Datum of 1929 (NGVD 29) by 2005. Withdrawals decreased to 5.01 Mgal/d by 2009 as withdrawals were shifted toward other sources, and by 2015 water levels had recovered to approximately 50 ft below NGVD 29. The mean estimated transmissivity of the Piney Point aquifer in York and James City Counties is 16,300 feet squared per day (ft 2 /d), but farther north it is only 925 ft 2 /d. The mean well specific capacity in York and James City Counties is 11.4 gallons per minute per foot (gal/min/ft). Farther north in Virginia, mean specific capacity is only 2.26 gal/min/ft, and in Maryland it is 0.99 gal/min/ft. The northward decrease in specific capacity probably reflects the northward decrease in transmissivity, which results from poor development of the solution-channeled limestone. An aquifer test in northern York County induced vertical leakage to the solution-channeled limestone from overlying silty sand and a change in response of the aquifer to pumping from a single layer to two layers. Transmissivity of the limestone of approximately 19,800 ft 2 /d was distinguished from the silty sand of approximately 2,500 ft 2 /d. Most of the water in the Piney Point aquifer is slightly alkaline with moderate concentrations primarily of sodium and bicarbonate that are slightly undersaturated with respect to calcite. Iron concentrations are generally less than 0.3 milligrams per liter (mg/L). Mixing of freshwater with seawater has elevated chloride concentrations to the southeast to as much as 7,120 mg/L. Information on the Piney Point aquifer can benefit water-resource management in siting production wells, predicting likely well yield, and anticipating water-level response to withdrawals. Models that vertically discretize individual geologic units can potentially be used to evaluate groundwater flow in greater detail by representing lateral flow and vertical leakage among the geologic units. Because groundwater withdrawals are made primarily from the limestone and sand of the Piney Point Formation, the VA DEQ has considered regarding the limestone and sand singly as a regulated aquifer apart from the other geologic units. Under current policy in Virginia, if only the limestone and sand were regarded as a regulated aquifer, a greater amount of drawdown would be allowed than is allowed for the Piney Point aquifer consisting of six geologic units. Some production wells intercept multiple geologic units, and the units can undergo water-level decline and vertical leakage induced by pumping from the limestone and sand. Whether the other geologic units are to be regarded as regulated aquifers is an additional consideration for the VA DEQ.

Virginia↗

Low-flow frequency and flow-duration characteristics of selected streams in Alabama through March 2014

Low-flow statistics are needed by water-resource engineers, planners, and managers to protect and manage the water resources of Alabama. The accuracy of these statistics is influenced by such factors as length of record and specific hydrologic conditions measured in those records. As such, it is generally recommended that flow statistics be updated about every 10 years to provide improved and representative low-flow characteristics. The previous investigation of low-flow characteristics for Alabama included data through September 1990. Since that time, Alabama has experienced several historic droughts highlighting the need to update the low-flow characteristics at U.S. Geological Survey streamgaging stations. Consequently, this investigation was undertaken in cooperation with a number of State and local agencies to update low-flow frequency and flow-duration statistics at 210 continuous-record streamgaging stations in Alabama and 67 stations from basins that are shared with surrounding States. The flow characteristics were computed on the basis of available data through March 2014.

Alabama↗

Simulation of daily streamflow for 12 river basins in western Iowa using the Precipitation-Runoff Modeling System

The U.S. Geological Survey, in cooperation with the Iowa Department of Natural Resources, constructed Precipitation-Runoff Modeling System models to estimate daily streamflow for 12 river basins in western Iowa that drain into the Missouri River. The Precipitation-Runoff Modeling System is a deterministic, distributed-parameter, physical-process-based modeling system developed to evaluate the response of streamflow and general drainage basin hydrology to various combinations of climate and land use. Calibration periods for each basin varied depending on the period of record available for daily mean streamflow measurements at U.S. Geological Survey streamflow-gaging stations. A geographic information system tool was used to delineate each basin and estimate initial values for model parameters based on basin physical and geographical features. A U.S. Geological Survey automatic calibration tool that uses a shuffled complex evolution algorithm was used for initial calibration, and then manual modifications were made to parameter values to complete the calibration of each basin model. The main objective of the calibration was to match daily discharge values of simulated streamflow to measured daily discharge values. The Precipitation-Runoff Modeling System model was calibrated at 42 sites located in the 12 river basins in western Iowa. The accuracy of the simulated daily streamflow values at the 42 calibration sites varied by river and by site. The models were satisfactory at 36 of the sites based on statistical results. Unsatisfactory performance at the six other sites can be attributed to several factors: (1) low flow, no flow, and flashy flow conditions in headwater subbasins having a small drainage area; (2) poor representation of the groundwater and storage components of flow within a basin; (3) lack of accounting for basin withdrawals and water use; and (4) limited availability and accuracy of meteorological input data. The Precipitation-Runoff Modeling System models of 12 river basins in western Iowa will provide water-resource managers with a consistent and documented method for estimating streamflow at ungaged sites and aid in environmental studies, hydraulic design, water management, and water-quality projects.

Iowa↗