Search USGSSearch

SEARCH · Search USGS

Results for “Modeling 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 145 records · Page 8Linked to original sources

Uniform California earthquake rupture forecast, version 3 (UCERF3): the time-independent model

In this report we present the time-independent component of the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3), which provides authoritative estimates of the magnitude, location, and time-averaged frequency of potentially damaging earthquakes in California. The primary achievements have been to relax fault segmentation assumptions and to include multifault ruptures, both limitations of the previous model (UCERF2). The rates of all earthquakes are solved for simultaneously, and from a broader range of data, using a system-level "grand inversion" that is both conceptually simple and extensible. The inverse problem is large and underdetermined, so a range of models is sampled using an efficient simulated annealing algorithm. The approach is more derivative than prescriptive (for example, magnitude-frequency distributions are no longer assumed), so new analysis tools were developed for exploring solutions. Epistemic uncertainties were also accounted for using 1,440 alternative logic tree branches, necessitating access to supercomputers. The most influential uncertainties include alternative deformation models (fault slip rates), a new smoothed seismicity algorithm, alternative values for the total rate of M≥5 events, and different scaling relationships, virtually all of which are new. As a notable first, three deformation models are based on kinematically consistent inversions of geodetic and geologic data, also providing slip-rate constraints on faults previously excluded because of lack of geologic data. The grand inversion constitutes a system-level framework for testing hypotheses and balancing the influence of different experts. For example, we demonstrate serious challenges with the Gutenberg-Richter hypothesis for individual faults. UCERF3 is still an approximation of the system, however, and the range of models is limited (for example, constrained to stay close to UCERF2). Nevertheless, UCERF3 removes the apparent UCERF2 overprediction of M6.5–7 earthquake rates and also includes types of multifault ruptures seen in nature. Although UCERF3 fits the data better than UCERF2 overall, there may be areas that warrant further site-specific investigation. Supporting products may be of general interest, and we list key assumptions and avenues for future model improvements.

California

Habitat Suitability Index Models: Gray partridge

This report is part of the Habitat Suitability Index model series which provides habitat information useful for impact assessment and habitat management. The Habitat Use Information section is largely constrained to those data that can be used to derive quantitative relationships between key environmental variables and habitat suitability. The habitat use information provides the foundation for the HSI model that follows. In addition, this same information may be useful in the development of other models more appropriate to specific assessment or evaluation needs for the gray partridge (Perdix perdix). The HSI model section documents a habitat model and information pertinent to its application. The model synthesizes the habitat use information into a framework appropriate for field application and is scaled to produce an index value between 0.0 (unsuitable habitat) and 1.0 (optimum habitat).

FWS/OBS

Estimated Perennial Streams of Idaho and Related Geospatial Datasets

The perennial or intermittent status of a stream has bearing on many regulatory requirements. Because of changing technologies over time, cartographic representation of perennial/intermittent status of streams on U.S. Geological Survey (USGS) topographic maps is not always accurate and (or) consistent from one map sheet to another. Idaho Administrative Code defines an intermittent stream as one having a 7-day, 2-year low flow (7Q2) less than 0.1 cubic feet per second. To establish consistency with the Idaho Administrative Code, the USGS developed regional regression equations for Idaho streams for several low-flow statistics, including 7Q2. Using these regression equations, the 7Q2 streamflow may be estimated for naturally flowing streams anywhere in Idaho to help determine perennial/intermittent status of streams. Using these equations in conjunction with a Geographic Information System (GIS) technique known as weighted flow accumulation allows for an automated and continuous estimation of 7Q2 streamflow at all points along a stream, which in turn can be used to determine if a stream is intermittent or perennial according to the Idaho Administrative Code operational definition. The selected regression equations were applied to create continuous grids of 7Q2 estimates for the eight low-flow regression regions of Idaho. By applying the 0.1 ft3/s criterion, the perennial streams have been estimated in each low-flow region. Uncertainty in the estimates is shown by identifying a 'transitional' zone, corresponding to flow estimates of 0.1 ft3/s plus and minus one standard error. Considerable additional uncertainty exists in the model of perennial streams presented in this report. The regression models provide overall estimates based on general trends within each regression region. These models do not include local factors such as a large spring or a losing reach that may greatly affect flows at any given point. Site-specific flow data, assuming a sufficient period of record, generally would be considered to represent flow conditions better at a given site than flow estimates based on regionalized regression models. The geospatial datasets of modeled perennial streams are considered a first-cut estimate, and should not be construed to override site-specific flow data.

Data Series

The Lake Tahoe Basin Land Use Simulation Model

This U.S. Geological Survey Open-File Report describes the final modeling product for the Tahoe Decision Support System project for the Lake Tahoe Basin funded by the Southern Nevada Public Land Management Act and the U.S. Geological Survey's Geographic Analysis and Monitoring Program. This research was conducted by the U.S. Geological Survey Western Geographic Science Center. The purpose of this report is to describe the basic elements of the novel Lake Tahoe Basin Land Use Simulation Model, publish samples of the data inputs, basic outputs of the model, and the details of the Python code. The results of this report include a basic description of the Land Use Simulation Model, descriptions and summary statistics of model inputs, two figures showing the graphical user interface from the web-based tool, samples of the two input files, seven tables of basic output results from the web-based tool and descriptions of their parameters, and the fully functional Python code.

Open-File Report

Simulation of contributing areas and surface-water leakage to potential replacement wells near the community of New Post, Sawyer County, Wisconsin, by means of a two-dimensional ground-water flow model

A two-dimensional, steady-state ground-water-flow model of the shallow ground-water-flow system near the community of New Post, Sawyer County, Wis., was refined from an existing model of the area. Hydraulic-conductivity and recharge values were not changed from the existing model for the scenario simulations described in this report. Rather, the model was refined by adding detail along the Chippewa Flowage and then was used to simulate contributing areas for three potential replacement wells pumping 30,000 gallons per day. The model also was used to simulate potential surface-water leakage out of the Chippewa Flowage captured by replacement-well pumping. A range in resistance to vertical ground-water flow was simulated along the Chippewa Flowage for each potential replacement-well location to bound the potential effects of representing three-dimensional flow with a two-dimensional model. Results indicate that pumping from a replacement well sited about 130 feet from the Chippewa Flowage could capture as much as 39 percent of the total pumping from the flowage. Pumping from either of two potential replacement wells sited at least 400 feet from the Chippewa Flowage did not induce surface-water leakage out of the flowage regardless of the resistance applied along the flowage for simulations described in this report.

Wisconsin

Additional descriptive models of industrial mineral deposits

This report includes previously written, but unpublished descriptive industrial mineral deposit models not released in earlier reports (Cox and Singer, 1986; Orris, 1992; Orris and Bliss, 1991, 1992). Although some of the models are preliminary, recent and repeated requests for several of the models has led to their release in this compilation. Initial drafts of the clay models were written for individual U.S. deposits or districts by J.W. Hosterman2 of the U.S. Geological Survey (USGS) in the late 1980/s. All of the models in this report have been modified from their original unpublished versions; for most models, the modifications have consisted of additional geologic detail, economic information, and (or) health and environmental considerations related to the deposit type. More extensive alterations have been made to the clay models written by Hosterman, which were put into standard format and extensively revised by G.J. Orris. As much as two-thirds of the information in the current versions of the models post-dates Hosterman's contributions and these models are shown as jointly authored.

Open-File Report

Upper Mississippi River system weighted wind fetch analysis (1989, 2000, 2010/2011)

Wind fetch is defined as the unobstructed distance that wind can travel over water in a constant direction. Fetches are limited by landforms surrounding the body of water. Fetch is an important characteristic of open water because longer fetches can result in larger wind-generated waves. The larger waves, in turn, can increase shoreline erosion and sediment resuspension (Rohweder and others 2012). Increases in sediment resuspension lead to increases in water turbidity, which in turn decreases light penetration and, therefore, create conditions less conducive to aquatic plant growth (Giblin and others 2010). A wind fetch model was developed by David Finlayson, U. S. Geological Survey, Pacific Science Center, while he was a Ph.D. student at the University of Washington (Finlayson 2005). This method calculates effective fetch using the recommended procedure of the Shore Protection Manual (USACE 1984). Scientists at the United States Geological Survey, Upper Midwest Environmental Sciences Center (UMESC) and the United States Army Corps of Engineers (USACE) further refined this model (Rohweder and others 2012) and structured it to operate using the most recent version of the ArcMap Geographic Information System platform (Esri, 2019). At the time the analysis was performed, the version of ArcMap used was 10.7.1. The model refined in 2012 was used for the analyses described in this report. Using this model, UMESC performed an analysis to model weighted wind fetch for the Upper Mississippi River System (UMRS) corresponding to three separate time periods of land cover spatial data acquisition (1989, 2000, and 2010/2011). The purpose of the analysis was to examine how fetch varies over time and space within the UMRS for potential management applications. For more detailed information on the wind fetch model, examine the USGS Open-File Report by Rohweder and others (2012).

Illinois, Iowa, Minnesota, Missouri, Wisconsin

Volcanic Aquifers of Hawai‘i—Construction and calibration of numerical models for assessing groundwater availability on Kaua‘i, O‘ahu, and Maui

Steady-state numerical groundwater-flow models were constructed for the islands of Kaua‘i, O‘ahu, and Maui to enable quantification of the hydrologic consequences of withdrawals and other stresses that can place limits on groundwater availability. The volcanic aquifers of Hawai‘i supply nearly all drinking water for the islands’ residents, freshwater for diverse industries, and natural discharge to springs, streams, and nearshore areas that support ecosystems, cultural practices, aesthetics, and recreation. Increases in groundwater withdrawal and changes in climate can cause water-table depression, saltwater rise, and reduction of natural groundwater discharge—all of which can limit fresh groundwater availability. The numerical models described in this report are designed to quantify these consequences. Separate models were created for each island using MODFLOW-2005 with the Seawater Intrusion package, which allows simulation of freshwater and saltwater in ocean-island aquifers. Calibration resulted in models that generally replicate observed water-level, stream base-flow, and spring-flow data, and simulate groundwater-flow directions and fresh groundwater thicknesses that are consistent with conceptual models. The calibrated models use hydraulic properties that are consistent with the ranges reported in previous studies. The models show that the relative distribution of fresh groundwater discharge to the ocean, streams, and springs and withdrawals for human use differ substantially among the three islands studied here. These differences indicate that consequences that limit the availability of fresh groundwater for human use are likely to differ among the three islands.

Hawaii

Ground-Water Flow Model for the Spokane Valley-Rathdrum Prairie Aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

This report presents a computer model of ground-water flow in the Spokane Valley-Rathdrum Prairie (SVRP) aquifer in Spokane County, Washington, and Bonner and Kootenai Counties, Idaho. The aquifer is the sole source of drinking water for more than 500,000 residents in the area. In response to the concerns about the impacts of increased ground-water withdrawals resulting from recent and projected urban growth, a comprehensive study was initiated by the Idaho Department of Water Resources, the Washington Department of Ecology, and the U.S. Geological Survey to improve the understanding of ground-water flow in the aquifer and of the interaction between ground water and surface water. The ground-water flow model presented in this report is one component of this comprehensive study. The primary purpose of the model is to serve as a tool for analyzing aquifer inflows and outflows, simulating the effects of future changes in ground-water withdrawals from the aquifer, and evaluating aquifer management strategies. The scale of the model and the level of detail are intended for analysis of aquifer-wide water-supply issues. The SVRP aquifer model was developed by the Modeling Team formed within the comprehensive study. The Modeling Team consisted of staff and personnel working under contract with the Idaho Department of Water Resources, personnel working under contract with the Washington Department of Ecology, and staff of the U.S. Geological Survey. To arrive at a final model that has the endorsement of all team members, decisions on modeling approach, methodology, assumptions, and interpretations were reached by consensus. The ground-water flow model MODFLOW-2000 was used to simulate ground-water flow in the SVPR aquifer. The finite-difference model grid consists of 172 rows, 256 columns, and 3 layers. Ground-water flow was simulated from September 1990 through September 2005 using 181 stress periods of 1 month each. The areal extent of the model encompasses an area of approximately 326 square miles. For the most part, the model extent coincides with the 2005 revised extent of the Spokane Valley-Rathdrum Prairie aquifer as defined in a previous report. However, the model excludes Spirit and Hoodoo Valleys because of uncertainties about the ground-water flow directions in those valleys and the degree of hydraulic connection between the valleys and northern Rathdrum Prairie. The SVRP aquifer is considered to be a single hydrogeologic unit except in Hillyard Trough and the Little Spokane River Arm. In those areas, a continuous clay layer divides the aquifer into an upper, unconfined unit and a lower, confined unit. The model includes all known components of inflows to and outflows from the aquifer. Inflows to the SVRP aquifer include (1) recharge from precipitation, (2) inflows from tributary basins and adjacent uplands, (3) subsurface seepage and surface overflows from lakes that border the aquifer, (4) flow from losing segments of the Spokane River to the aquifer, (5) return percolation from irrigation, and (6) effluent from septic systems. Outflows from the SVRP aquifer include (1) ground-water withdrawals from wells, (2) flow from the aquifer to gaining segments of the Spokane River, (3) aquifer discharge to the Little Spokane River, and (4) subsurface outflow from the lower unit at the western limit of the model area near Long Lake. These inflow and outflow components are represented in the model by using MODFLOW-2000 packages. The parameter-estimation program PEST was used to calibrate the SVRP aquifer model. PEST implements a nonlinear least-squares regression method to estimate model parameters so that the differences between measured and simulated quantities are minimized with respect to an optimal criterion. Calibration data include 1,573 measurements of water levels and 313 measurements of streamflow gains and losses along segments of the Spokane and Little Spokane Rivers. Model parameters estimated during calib

Scientific Investigations Report

Description of the National Hydrologic Model for use with the Precipitation-Runoff Modeling System (PRMS)

This report documents several components of the U.S. Geological Survey National Hydrologic Model of the conterminous United States for use with the Precipitation-Runoff Modeling System (PRMS). It provides descriptions of the (1) National Hydrologic Model, (2) Geospatial Fabric for National Hydrologic Modeling, (3) PRMS hydrologic simulation code, (4) parameters and estimation methods used to compute spatially and temporally distributed default values as required by PRMS, (5) National Hydrologic Model Parameter Database, and (6) model extraction tool named Bandit. The National Hydrologic Model Parameter Database contains values for all PRMS parameters used in the National Hydrologic Model. The methods and national datasets used to estimate all the PRMS parameters are described. Some parameter values are derived from characteristics of topography, land cover, soils, geology, and hydrography using traditional Geographic Information System methods. Other parameters are set to long-established default values and computation of initial values. Additionally, methods (statistical, sensitivity, calibration, and algebraic) were developed to compute parameter values on the basis of a variety of nationally-consistent datasets. Values in the National Hydrologic Model Parameter Database can periodically be updated on the basis of new parameter estimation methods and as additional national datasets become available. A companion ScienceBase resource provides a set of static parameter values as well as images of spatially-distributed parameters associated with PRMS states and fluxes for each Hydrologic Response Unit across the conterminuous United States.

Techniques and Methods

An Adaptive Management Approach for Summer Water Level Reductions on the Upper Mississippi River System

The primary purpose of this report is to provide an adaptive management approach for learning more about summer water level reductions (drawdowns) as a management tool, including where and how drawdowns can be applied most effectively within the Upper Mississippi River System. The report reviews previous drawdowns conducted within the system and provides specific recommendations for learning more about the lesser known effects of drawdowns and how the outcomes can be influenced by different implementation strategies and local conditions. The knowledge gained can be used by managers to determine how best to implement drawdowns in different parts of the UMRS to help achieve management goals. The information and recommendations contained in the report are derived from results of previous drawdown projects, insights from regional disciplinary experts, and the experience of the authors in experimental design, modeling, and monitoring. Modeling is a critical part of adaptive management and can involve conceptual models, simulation models, and empirical models. In this report we present conceptual models that express current understanding regarding functioning of the UMRS as related to drawdowns and highlight interactions among key ecological components of the system. The models were developed within the constraints of drawdown timing, magnitude (depth), and spatial differences in effects (longitudinal and lateral) with attention to ecological processes affected by drawdowns. With input from regional experts we focused on the responses of vegetation, fish, mussels, other invertebrates, and birds. The conceptual models reflect current understanding about relations and interactions among system components, the expected strength of those interactions, potential responses of system components to drawdowns, likelihood of the response occurring, and key uncertainties that limit our ability to make accurate predictions of effects (Table 1, Fig. 4-10). Based on this current understanding, the main questions still associated with drawdowns include (1) the effects of frequency of drawdowns (from once every few years to multiple years in succession); (2) timing of the beginning of drawdowns (follow the descending arm of the flood pulse versus always beginning in early summer); (3) long-term benefits (greater than 5-6 years), especially as compared to known short-term loses (e.g., mortality of mussels in exposed areas, loss of submersed vegetation in exposed areas, cost of advanced dredging); and (4) the effects in northern (above pool 14) versus southern pools (pool 14 and below, and the Illinois River). An adaptive management design should address these questions to reduce uncertainty in predictions of drawdown effects and help determine if different implementation strategies are needed in different parts of the system. Given that drawdowns will continue to be used as a management tool on the UMRS, we suggest that some drawdowns be conducted in an adaptive management context that helps meet management objectives, but also provides efficient learning about the questions listed above. We propose two different, but interrelated, experimental designs to address these questions. Both designs call for conducting multiple drawdowns in multiple pools (2-4 pools) to allow direct comparison of results and produce rapid learning. However, the report does not provide a detailed scope of work for carrying out the designs. If managers choose to implement one of the experimental designs, specifics of choosing appropriate pools and developing a monitoring plan will need to be determined through collaboration among managers, researchers, and statisticians. We suggest characteristics to consider in selecting treatment and reference pools (study sites) and also provide guidance for developing a monitoring plan. Some aspects of these two designs could be implemented individually, but by implementing individual elements, direct comparisons of some design features

Illinois;Iowa;Minnesota;Missouri;Wisconsin

Geospatial data for developing nutrient SPARROW models for the Midcontinental region of Canada and the United States

Through the International Watersheds Initiative of the International Joint Commission (IJC), the SpatiallyReferenced Regressions on Watershed attributes (SPARROW) model developed by the U.S. Geological Survey (USGS) is being applied to the Great Lakes, Rainy River – Lake of the Woods and Red-Assiniboine basins. The objective of this binational application of the SPARROW model is to better understand and quantify the sources of phosphorus (P) and nitrogen (N) that contribute to regional water-quality issues like algal blooms and eutrophication in Lake Erie and other parts of the Great Lakes, as well as Lake of the Woods. Led by the IJC, a team of researchers from the National Research Council of Canada – Ocean, Coastal and River Engineering Research Centre, USGS, and IJC are extending the SPARROW modelling work previously completed for the Red-Assiniboine basin and the U.S. portions of the Great Lakes, Ohio, Upper Mississippi, and Souris-Red-Rainy river basins to cover all of the Great Lakes, Rainy River – Lake of the Woods and RedAssiniboine basins. The current effort is termed the Midcontinent SPARROW modelling study. This report describes the data used to develop the Midcontinent SPARROW models, specifically the sources of original data, assembling the data, and the processing and harmonization required between the U.S. and Canada data needed to produce these models. Details provided include the: development of a digital stream network and related catchments – most significantly in the Canadian regions of the Great Lakes and Rainy River – Lake of the Woods basins where these data were not available to create a seamless binational network across the model domain; calculation of variables to aid in the determination of in-stream and in-reservoir decay of P and N; quantification of Canadian diversions within the Midcontinental region (i.e., Lake St. Joseph, Long Lake and the Ogoki Reservoir); development of binational input nutrient sources considered for model development (i.e., land cover, inorganic farm fertilizer, manure, atmospheric deposition, point-source pollution from wastewatertreatment plants and contribution from non-modelled watersheds); and development of delivery variables considered to be most predominant (i.e., temperature, precipitation and ensuing runoff, soil permeability and clay content, slope of the catchments, and tile drainage). The majority of the geospatial data collection and processing was required for Canadian datasets because many of the U.S. datasets were already assembled for previous SPARROW model applications in the U.S. The task of harmonizing data between the U.S. and Canada was important to ensure consistency of the datasets used in the models. The harmonized digital stream network, delineated catchments and input data for each catchment (i.e., source and delivery variables), created for the Midcontinent SPARROW models, are available for download at url: https://doi.org/10.4224/300.0001.

Technical Report

Documentation of the dynamic parameter, water-use, stream and lake flow routing, and two summary output modules and updates to surface-depression storage simulation and initial conditions specification options with the Precipitation-Runoff Modeling System (PRMS)

This report documents seven enhancements to the U.S. Geological Survey (USGS) Precipitation-Runoff Modeling System (PRMS) hydrologic simulation code: two time-series input options, two new output options, and three updates of existing capabilities. The enhancements are (1) new dynamic parameter module, (2) new water-use module, (3) new Hydrologic Response Unit (HRU) summary output module, (4) new basin variables summary output module, (5) new stream and lake flow routing module, (6) update to surface-depression storage and flow simulation, and (7) update to the initial-conditions specification. This report relies heavily upon U.S. Geological Survey Techniques and Methods, book 6, chapter B7, which documents PRMS version 4 (PRMS-IV). A brief description of PRMS is included in this report.

Techniques and Methods

One-Water Hydrologic Flow Model: A MODFLOW based conjunctive-use simulation software

The U.S. Geological Survey’s (USGS) Modular Ground-Water Flow Model (MODFLOW-2005) is a computer program that simulates groundwater flow by using finite differences. The MODFLOW-2005 framework uses a modular design that allows for the easy development and incorporation of new features called processes and packages that work with or modify inputs to the groundwater-flow equation. A process solves a flow equation or set of equations. For example, the central part of MODFLOW is the groundwater-flow process that solves the groundwater-flow equation; the surface-water routing process is an additional process that solves the surface-water flow equation. Packages are code related to the groundwater-flow process. For example, the subsidence package modifies the groundwater-flow process by including aquifer compaction effects on flow. With the development of new packages and processes, the MODFLOW-2005 base framework diverged into multiple independent versions designed for specific simulation needs. This divergence limited each independent MODFLOW release to its specific purpose, so that there was no longer a single, comprehensive, general-purpose hydraulic-simulation framework. The MODFLOW One-Water Hydrologic Flow Model (MF-OWHM, also informally known as OneWater) is an integrated hydrologic flow model that combines multiple MODFLOW-2005 variants in one cohesive simulation software; changes were made to enable multiple capabilities in one code. This fusion of the MODFLOW-2005 versions resulted in a simulation software that can be used to address and analyze a wide class of conjunctive-use, water-management, water-food-security, and climate-crop-water scenarios. As a second core version of MODFLOW-2005, MF-OWHM maintains backward compatibility with existing MODFLOW-2005 versions, with features that include the following: Process-based simulation. Saturated groundwater flow (three-dimensional). Surface-water flow (one- and two-dimensional). Stream and river flow. Lake and reservoir storage. Landscape simulation and irrigated agriculture. Land-use and crop simulation. Root uptake of groundwater. Actual evapotranspiration. Estimated irrigation demand. Reservoir operations. Aquifer compaction and subsidence by vertical model-grid deformation. Seawater intrusion by a sharp-interface assumption. Karst-aquifer and fractured-bedrock flow. Turbulent and laminar-pipe network flow. Unsaturated groundwater flow (one-dimensional). Internal linkages among the processes that couple hydraulic head, flow, and deformation. Redesigned code for faster simulation, increased user-input options, easier model updates, and more robust error reporting than in previous models. MF-OWHM is a MODFLOW-2005 based integrated hydrologic model that can simulate and analyze varying environmental conditions to allow for the evaluation of management options from many components of human and natural water movement through a physically based, supply and demand framework. The term “integrated,” in the context of this report, refers to the tight coupling of groundwater flow, surface-water flow, landscape processes, aquifer compaction and subsidence, reservoir operations, and conduit (karst) flow. Another benefit of this integrated hydrologic model is that models developed to run by MODFLOW-2005, MODFLOW-NWT, MODFLOW-CFP, or MODFLOW-FMP can also be simulated with MF-OWHM. At the time of this report’s publication, MF-OWHM version 2 (MF-OWHM2) does not include a direct internal simulation of snowmelt, advanced mountainous watershed rainfall-runoff simulation, detailed shallow soil-moisture accounting, or atmospheric moisture content. Atmospheric moisture may be accounted for indirectly by, optionally, specifying a pan-evaporation rate, reference evapotranspiration, and precipitation. These features are not included to ensure that simulation runtime remains short enough to enable the use of automated methods of calibrating model parameters to field observations, which typically require many simulation model runs. The MF-OWHM approach is to include as much detail as possible to simulate hydrological processes, providing the simulation runtimes remain reasonable enough to allow for robust parameter estimation and model calibration. To represent both natural and human-influenced flow, MF-OWHM integrates physically based flow processes derived from MODFLOW-2005 in a supply and demand framework. From this integration, the physically based movement of groundwater, surface water, imported water, and precipitation serve as supply to meet consumptive demands associated with irrigated and non-irrigated agriculture, natural vegetation, and urban water uses. Water consumption is determined by balancing the available water supply with water demand, leading to the concept of a demand-driven, supply-constrained simulation. The MF-OWHM Supply-and-Demand Framework is especially useful for the analysis of agricultural water use, where there are often few data available to describe changes in land-use through time, such as crop type and distribution, and the associated changes in groundwater pumpage. This framework attempts to satisfy each land-use water demand with available water supplies—that is, groundwater uptake, precipitation, and irrigation. An option provided in MF-OWHM2 is to automatically increase groundwater pumping for irrigation, which often is unknown, by the calculated residual between demand and the other available sources of supply. From large- to small-scale applications, the physically based supply and demand framework provides key capabilities for simulating and analyzing historical, current, and future conjunctive-use of surface water and groundwater. To achieve the physically based supply and demand framework, the MODFLOW-2005 standard of no inter-package and -process communication was relaxed for MF-OWHM2. Traditional MODFLOW simulation models required that all packages and processes interact through the groundwater-flow equation or by removing the water flow from the simulation domain. For example, the MODFLOW-2005 representation of a groundwater well extracts water from the groundwater-flow equation (by subtraction) and removes it from the simulation domain. This feature is available in the MF-OWHM framework, but options have been added to allow the specification of a use or destination of pumped groundwater within the model domain, for example, it can be used for irrigation, managed aquifer recharge, or return-flow to streams.

Techniques and Methods

Simulation of streamflow using a multidimensional flow model for white sturgeon habitat, Kootenai River near Bonners Ferry, Idaho – Supplement to Scientific Investigations Report 2005-5230

During 2005, the U.S. Geological Survey (USGS) developed, calibrated, and validated a multidimensional flow model for simulating streamflow in the white sturgeon spawning habitat of the Kootenai River in Idaho. The model was developed as a tool to aid understanding of the physical factors affecting quality and quantity of spawning and rearing habitat used by the endangered white sturgeon (Acipenser transmontanus) and for assessing the feasibility of various habitat-enhancement scenarios to re-establish recruitment of white sturgeon. At the request of the Kootenai Tribe of Idaho, the USGS extended the two-dimensional flow model developed in 2005 into a braided reach upstream of the current white sturgeon spawning reach. Many scientists consider the braided reach a suitable substrate with adequate streamflow velocities for re-establishing recruitment of white sturgeon. The 2005 model was extended upstream to help assess the feasibility of various strategies to encourage white sturgeon to spawn in the reach. At the request of the Idaho Department of Fish and Game, the USGS also extended the two-dimensional flow model several kilometers downstream of the white sturgeon spawning reach. This modified model can quantify the physical characteristics of a reach that white sturgeon pass through as they swim upstream from Kootenay Lake to the spawning reach. The USGS Multi-Dimensional Surface-Water Modeling System was used for the 2005 modeling effort and for this subsequent modeling effort. This report describes the model applications and limitations, presents the results of a few simple simulations, and demonstrates how the model can be used to link physical characteristics of streamflow to the location of white sturgeon spawning events during 1994-2001. Model simulations also were used to report on the length and percentage of longitudinal profiles that met the minimum criteria during May and June 2006 and 2007 as stipulated in the U.S. Fish and Wildlife Biological Opinion.

Idaho

Continuous real-time water-quality monitoring and regression analysis to compute constituent concentrations and loads in the North Fork Ninnescah River upstream from Cheney Reservoir, south-central Kansas, 1999–2012

Cheney Reservoir, located in south-central Kansas, is the primary water supply for the city of Wichita. The U.S. Geological Survey has operated a continuous real-time water-quality monitoring station since 1998 on the North Fork Ninnescah River, the main source of inflow to Cheney Reservoir. Continuously measured water-quality physical properties include streamflow, specific conductance, pH, water temperature, dissolved oxygen, and turbidity. Discrete water-quality samples were collected during 1999 through 2009 and analyzed for sediment, nutrients, bacteria, and other water-quality constituents. Regression models were developed to establish relations between discretely sampled constituent concentrations and continuously measured physical properties to compute concentrations of those constituents of interest that are not easily measured in real time because of limitations in sensor technology and fiscal constraints. Regression models were published in 2006 that were based on data collected during 1997 through 2003. This report updates those models using discrete and continuous data collected during January 1999 through December 2009. Models also were developed for four new constituents, including additional nutrient species and indicator bacteria. In addition, a conversion factor of 0.68 was established to convert the Yellow Springs Instruments (YSI) model 6026 turbidity sensor measurements to the newer YSI model 6136 sensor at the North Ninnescah River upstream from Cheney Reservoir site. Newly developed models and 14 years of hourly continuously measured data were used to calculate selected constituent concentrations and loads during January 1999 through December 2012. The water-quality information in this report is important to the city of Wichita because it allows the concentrations of many potential pollutants of interest to Cheney Reservoir, including nutrients and sediment, to be estimated in real time and characterized over conditions and time scales that would not be possible otherwise. In general, model forms and the amount of variance explained by the models was similar between the original and updated models. The amount of variance explained by the updated models changed by 10 percent or less relative to the original models. Total nitrogen, nitrate, organic nitrogen, E. coli bacteria, and total organic carbon models were newly developed for this report. Additional data collection over a wider range of hydrological conditions facilitated the development of these models. The nitrate model is particularly important because it allows for comparison to Cheney Reservoir Task Force goals. Mean hourly computed total suspended solids concentration during 1999 through 2012 was 54 milligrams per liter (mg/L). The total suspended solids load during 1999 through 2012 was 174,031 tons. On an average annual basis, the Cheney Reservoir Task Force runoff (550 mg/L) and long-term (100 mg/L) total suspended solids goals were never exceeded, but the base flow goal was exceeded every year during 1999 through 2012. Mean hourly computed nitrate concentration was 1.08 mg/L during 1999 through 2012. The total nitrate load during 1999 through 2012 was 1,361 tons. On an annual average basis, the Cheney Reservoir Task Force runoff (6.60 mg/L) nitrate goal was never exceeded, the long-term goal (1.20 mg/L) was exceeded only in 2012, and the base flow goal of 0.25 mg/L was exceeded every year. Mean nitrate concentrations that were higher during base flow, rather than during runoff conditions, suggest that groundwater sources are the main contributors of nitrate to the North Fork Ninnescah River above Cheney Reservoir. Mean hourly computed phosphorus concentration was 0.14 mg/L during 1999 through 2012. The total phosphorus load during 1999 through 2012 was 328 tons. On an average annual basis, the Cheney Reservoir Task Force runoff goal of 0.40 mg/L for total phosphorus was exceeded in 2002, the year with the largest yearly mean turbidity, and the long-term goal (0.10 mg/L) was exceeded in every year except 2011 and 2012, the years with the smallest mean streamflows. The total phosphorus base flow goal of 0.05 mg/L was exceeded every year. Given that base flow goals for total suspended solids, nitrate, and total phosphorus were exceeded every year despite hydrologic conditions, the established base flow goals are either unattainable or substantially more best management practices will need to be implemented to attain them. On an annual average basis, no discernible patterns were evident in total suspended sediment, nitrate, and total phosphorus concentrations or loads over time, in large part because of hydrologic variability. However, more rigorous statistical analyses are required to evaluate temporal trends. A more rigorous analysis of temporal trends will allow evaluation of watershed investments in best management practices.

Kansas

Regression model development and computational procedures to support estimation of real-time concentrations and loads of selected constituents in two tributaries to Lake Houston near Houston, Texas, 2005-9

In December 2005, the U.S. Geological Survey (USGS), in cooperation with the City of Houston, Texas, began collecting discrete water-quality samples for nutrients, total organic carbon, bacteria ( Escherichia coli and total coliform), atrazine, and suspended sediment at two USGS streamflow-gaging stations that represent watersheds contributing to Lake Houston (08068500 Spring Creek near Spring, Tex., and 08070200 East Fork San Jacinto River near New Caney, Tex.). Data from the discrete water-quality samples collected during 2005–9, in conjunction with continuously monitored real-time data that included streamflow and other physical water-quality properties (specific conductance, pH, water temperature, turbidity, and dissolved oxygen), were used to develop regression models for the estimation of concentrations of water-quality constituents of substantial source watersheds to Lake Houston. The potential explanatory variables included discharge (streamflow), specific conductance, pH, water temperature, turbidity, dissolved oxygen, and time (to account for seasonal variations inherent in some water-quality data). The response variables (the selected constituents) at each site were nitrite plus nitrate nitrogen, total phosphorus, total organic carbon, E. coli , atrazine, and suspended sediment. The explanatory variables provide easily measured quantities to serve as potential surrogate variables to estimate concentrations of the selected constituents through statistical regression. Statistical regression also facilitates accompanying estimates of uncertainty in the form of prediction intervals. Each regression model potentially can be used to estimate concentrations of a given constituent in real time. Among other regression diagnostics, the diagnostics used as indicators of general model reliability and reported herein include the adjusted R-squared, the residual standard error, residual plots, and p-values. Adjusted R-squared values for the Spring Creek models ranged from .582–.922 (dimensionless). The residual standard errors ranged from .073–.447 (base-10 logarithm). Adjusted R-squared values for the East Fork San Jacinto River models ranged from .253–.853 (dimensionless). The residual standard errors ranged from .076–.388 (base-10 logarithm). In conjunction with estimated concentrations, constituent loads can be estimated by multiplying the estimated concentration by the corresponding streamflow and by applying the appropriate conversion factor. The regression models presented in this report are site specific, that is, they are specific to the Spring Creek and East Fork San Jacinto River streamflow-gaging stations; however, the general methods that were developed and documented could be applied to most perennial streams for the purpose of estimating real-time water quality data.

Texas