A data input program (MFI) for the U.S. Geological Survey modular finite-difference ground-water flow model
Explore the source record for details and available documents.
SEARCH · USGS Science
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Model assumptions and parameters used in an earlier study of the vulnerability of ground water in Kent County, Michigan, to atrazine contamination were reviewed and compared with other studies. The review indicated that model assumptions are consistent with those used in other models and that the parameters assigned in the Kent County model are within the broad range commonly found in the literature. Model assumptions most likely to limit the accuracy of the previous study include those of uniform transport and steady-state flow. Simulation results are more sensitive to parameter estimates for atrazine half life, organic-carbon content, and organic-carbon partition coefficient than to other model parameters. Potential atrazine detection probabilities and concentrations of atrazine reported in the previous Kent County Study were compared with data from other studies. Detection probabilities of 14 measurements of atrazine concentrations in ground water from shallow wells in the southern Lower Peninsula of Michigan were compared with detection probabilities based on potential atrazine concentrations computed for Kent County. Results indicate that the distribution of detection probabilities based on measured concentrations is similar to that based on adjusted potential concentrations. Potential concentrations were adjusted for effects of differences between sampling and modeling depths and for differences between the uniform application rate used for potential concentrations and the percentage of Kent County that is likely to be treated with atrazine. Potential concentrations of atrazine in the Kent County Study were within the wide range of concentrations measured in other states. A stratified random sampling strategy was developed to verify expected atrazine concentrations in ground water within Kent County. The strategy helps identify strata, determine the optimum allocation of ground-water samples within defined strata, and project the sampling error. Implementation of the strategy was illustrated by use of potential atrazine concentrations computed in the previous Kent County Study. Once detailed information on historical application rates of atrazine is developed, expected atrazine concentrations can be computed by use of the vulnerability model and used to implement the sampling strategy. Sampling results may be used to verify the effectiveness of the vulnerability model and local estimates of historical atrazine application rates by use of analysis of variance.
The U.S. Geological Survey is developing nationally consistent water-use modeling approaches to replace previous methods relying on locally specific reported and estimated data. These national assessments require datasets that incorporate water withdrawal variability across the United States and over long periods. However, source data often have unclear definitions, missing or varied units, differing temporal resolutions, varied data quality, and inconsistent formats, which hinder automation and require individualized processing. The public-supply datasets described in this paper were used in machine learning models to estimate annual and monthly public-supply water use for 2000–2020 for the conterminous United States (CONUS) and in a model to estimate public-supply deliveries. Public-supply withdrawal data were acquired for the CONUS and the District of Columbia; however, 11 states had annual data for only 1 year, and 10 states had no monthly data. Annual withdrawal data were acquired for 81% of public-supply water service areas, and monthly withdrawal data were acquired for 47% for at least 1 year from 2000 to 2020. These datasets and methods provide the most comprehensive collection of reported public-supply withdrawals to date and can be used by water-use managers, the scientific community, and the broader public. The extensive data processing described herein can be applicable to datasets representing other categories of water use.
An interactive computer simulation program, based on physical models of system sensors, platform geometry, Earth environment, and spheroidal magnetically-permeable targets, was developed to generate synthetic magnetic field data from a conceptual tensor magnetic gradiometer system equipped with an active primary field generator. The system sensors emulate the prototype tensor magnetic gradiometer system (TMGS) developed under a separate contract for unexploded ordnance (UXO) detection and classification. Time-series data from different simulation scenarios were analyzed to recover physical dimensions of the target source. Helbig-Euler simulations were run with rectangular and rod-like source bodies to determine whether such a system could separate the induced component of the magnetization from the remanent component for each target. This report concludes with an engineering assessment of a practical system design.
The Great Basin Integrated Landscape Monitoring Pilot Project was developed in response to the need for a monitoring and predictive capability that addresses changes in broad landscapes and waterscapes. Human communities and needs are nested within landscapes formed by interactions among the hydrosphere, geosphere, and biosphere. Understanding the complex processes that shape landscapes and deriving ways to manage them sustainably while meeting human needs require sophisticated modeling and monitoring. This document summarizes current understanding of ecosystem structure and function for many of the ecosystems within the Great Basin using conceptual models. The conceptual ecosystem models identify key ecological components and processes, identify external drivers, develop a hierarchical set of models that address both site and landscape attributes, inform regional monitoring strategy, and identify critical gaps in our knowledge of ecosystem function. The report also illustrates an approach for temporal and spatial scaling from site-specific models to landscape models and for understanding cumulative effects. Eventually, conceptual models can provide a structure for designing monitoring programs, interpreting monitoring and other data, and assessing the accuracy of our understanding of ecosystem functions and processes.
Introduction: Magnetotelluric data were acquired during October 2001 by the U.S. Geological Survey (USGS) as part of a study to examine the structural nature of basins in the transition zone between the Sierra Nevada Mountains of California and the Basin and Range province of Nevada. Magnetotelluric (MT) geophysical studies assist the mapping of geologic structure and the inference of lithologic packages that are concealed beneath the Earth's surface. The Basin and Range province has a complicated geologic history, which includes extension and compression of the Earth's crust to form the basins and ranges that blanket much of Nevada. The basins and ranges in the vicinity of this study trend northeastward and are bounded by steeply dipping strike slip faults. Interestingly, deep east-west magnetic trends occur in the aeromagnetic data of this study area indicating that the northeast-trending basins and ranges represent only thin-skinned deformation at the surface with an underlying east-west structure. To investigate this issue, MT data were acquired at seven stations in eastern California, 20 km east of Mono Lake. The purpose of this report is to present a two-dimensional apparent resistivity model of the MT data acquired for this study.
The amount and quality of natural resources available for terrestrial and aquatic wildlife habitats are expected to decrease throughout the world in areas that are intensively managed for urban and agricultural uses. Changes in climate and management of increasingly limited water supplies may further impact water resources essential for sustaining habitats. In this report, we document adapting a Water Evaluation and Planning (WEAP) system model for the Central Valley of California. We demonstrate using this adapted model (WEAP-CV wh ) to evaluate impacts produced from plausible future scenarios on agricultural and wetland habitats used by waterbirds and other wildlife. Processed output from WEAP-CV wh indicated varying levels of impact caused by projected climate, urbanization, and water supply management in scenarios used to exemplify this approach. Among scenarios, the NCAR-CCSM3 A2 climate projection had a greater impact than the CNRM-CM3 B1 climate projection, whereas expansive urbanization had a greater impact than strategic urbanization, on annual availability of waterbird habitat. Scenarios including extensive rice-idling or substantial instream flow requirements on important water supply sources produced large impacts on annual availability of waterbird habitat. In the year corresponding with the greatest habitat reduction for each scenario, the scenario including instream flow requirements resulted in the greatest decrease in habitats throughout all months of the wintering period relative to other scenarios. This approach provides a new and useful tool for habitat conservation planning in the Central Valley and a model to guide similar research investigations aiming to inform conservation, management, and restoration of important wildlife habitats.
There are concerns that removal and trimming of vegetation during expansion of the Hayward Airport in Sawyer County, Wisconsin, could appreciably change the character of a nearby cold-water stream and its adjacent environs. In cooperation with the Wisconsin Department of Transportation, a two-dimensional, steady-state groundwater-flow model of the shallow groundwater-flow system near the Hayward Airport was refined from a regional model of the area. The parameter-estimation code PEST was used to obtain a best fit of the model to additional field data collected in February 2007 as part of this study. The additional data were collected during an extended period of low runoff and consisted of water levels and streamflows near the Hayward Airport. Refinements to the regional model included one additional hydraulic-conductivity zone for the airport area, and three additional parameters for streambed resistance in a northern tributary to the Namekagon River and in the main stem of the Namekagon River. In the refined Hayward Airport area model, the calibrated hydraulic conductivity was 11.2 feet per day, which is within the 58.2 to 7.9 feet per day range reported for the regional glacial and sandstone aquifer, and is consistent with a silty soil texture for the area. The calibrated refined model had a best fit of 8.6 days for the streambed resistance of the Namekagon River and between 0.6 and 1.6 days for the northern tributary stream. The previously reported regional groundwater-recharge rate of 10.1 inches per year was adjusted during calibration of the refined model in order to match streamflows measured during the period of extended low runoff; this resulted in an optimal groundwater-recharge rate of 7.1 inches per year during this period. The refined model was then used to simulate the capture zone of the northern tributary to the Namekagon River.
For over 28 years, the U.S. Geological Survey (USGS) has been acquiring seismic velocity and geologic data at a number of locations in California, many of which were chosen because strong ground motions from earthquakes were recorded at the sites. The method for all measurements involves picking first arrivals of P- and S-waves from a surface source recorded at various depths in a borehole (as opposed to noninvasive methods, such as the SASW method [e.g., Brown et al., 2002]). The results from most of the sites are contained in a series of U.S. Geological Survey Open-File Reports (see References). Until now, none of the results have been available as computer files, and before 1992 the interpretation of the arrival times was in terms of piecemeal interval velocities, with no attempt to derive a layered model that would fit the travel times in an overall sense (the one exception is Porcella, 1984). In this report I reanalyze all of the arrival times in terms of layered models for P- and for S-wave velocities at each site, and I provide the results as computer files. In addition to the measurements reported in the open-file reports, I also include some borehole results from other reports, as well as some results never before published. I include data for 277 boreholes (at the time of this writing; more will be added to the web site as they are obtained), all in California (I have data from boreholes in Washington and Utah, but these will be published separately). I am also in the process of interpreting travel time data obtained using a seismic cone penetrometer at hundreds of sites; these data can be interpreted in the same way of those obtained from surface-to-borehole logging. When available, the data will be added to the web site (see below for information on obtaining data from the World Wide Web (WWW)). In addition to the basic borehole data and results, I provide information concerning strong-motion stations that I judge to be close enough to the boreholes that the borehole velocity models can be used as the velocity models beneath the stations.
Thermoelectric power generation accounts for over 41% of total U.S. freshwater withdrawals, making understanding the determinants of power plants’ water withdrawals (WW) and consumption (WC) critical for reducing the sector’s reliance on increasingly scarce water resources. However, reported data inconsistencies and incomplete analysis of potential determinants of thermoelectric water use hinder such understanding. We address these challenges by introducing a novel data filtering method and a more complete assessment of water use determinants. First, we applied a power-cooling ratio as an operations-based data filter that removed operationally implausible records while retaining more original data, outperforming previous statistical filtering methods. Second, we found that different water use reporting methods (WURMs) provided statistically significantly different WW and WC values, revealing the importance of this previously unrecognized feature in reported water use records. Third, our data-driven approach showed that traditionally emphasized features─such as cooling technology and gross generation─are of primary importance but can be surpassed by other, often overlooked, features when modeling WW or WC individually. The plant configuration, cooling technology, and gross generation were the most important features of WW, whereas WURM, cooling technology, and reporting month were the most important for WC. These findings can improve thermoelectric power plant management, water use reporting accuracy, and water use modeling.
Survival is a key life history parameter that can inform management decisions and basic life history research. Because true survival is often confounded with emigration from the study area, many studies are forced to estimate apparent survival (i.e., probability of surviving and remaining inside the study area), which can be much lower than true survival for highly mobile species. One method for estimating true survival is the Barker joint live-recapture/live-resight (JLRLR) model, which combines capture data from a study area (hereafter the ‘capture site’) with resighting data from a broader geographic area. This model assumes that live resights occur throughout the entire area where animals can disperse to and this assumption is often not met in practice. Here we use simulation to evaluate survival bias from a JLRLR model under study design scenarios that differ in the site selection for resights: global, random, fixed including the capture site, and fixed excluding the capture site. Simulation results indicate that fixed designs that included the capture site showed negative survival bias, whereas fixed designs that excluded the capture site exhibited positive survival bias. The magnitude of the bias was dependent on movement and survival, where scenarios with high survival and frequent movement had minimal bias. In an effort to help minimize bias, we developed a multistate version of the JLRLR and demonstrated reductions in survival bias compared to the single-state version for most designs. Our results suggest minimizing bias can be accomplished by: 1) using a random resight design when feasible if global sampling is not possible, 2) using the multistate JLRLR model when appropriate, 3) including the capture site in the resight sampling frame when possible, and 4) reporting survival as apparent survival if fixed sites are used for resight with the single state JLRLR model.
U.S. Geological Survey (USGS) scientists have provided a model-based assessment of transport and deposition of residual Deepwater Horizon oil along the shoreline within the northern Gulf of Mexico in the form of mixtures of sand and weathered oil, known as surface residual balls (SRBs). The results of this USGS research, in combination with results from other components of the overall study, will inform operational decisionmaking. The results will provide guidance for response activities and data collection needs during future oil spills. In May 2012 the U.S. Coast Guard, acting as the Deepwater Horizon Federal on-scene coordinator, chartered an operational science advisory team to provide a science-based review of data collected and to conduct additional directed studies and sampling. The goal was to characterize typical shoreline profiles and morphology in the northern Gulf of Mexico to identify likely sources of residual oil and to evaluate mechanisms whereby reoiling phenomena may be occurring (for example, burial and exhumation and alongshore transport). A steering committee cochaired by British Petroleum Corporation (BP) and the National Oceanic and Atmospheric Administration (NOAA) is overseeing the project and includes State on-scene coordinators from four States (Alabama, Florida, Louisiana, and Mississippi), trustees of the U.S. Department of the Interior (DOI), and representatives from the U.S. Coast Guard. This report presents the results of hydrodynamic and sediment transport models and developed techniques for analyzing potential SRB movement and burial and exhumation along the coastline of Alabama and Florida. Results from these modeling efforts are being used to explain the complexity of reoiling in the nearshore environment and to broaden consideration of the different scenarios and difficulties that are being faced in identifying and removing residual oil. For instance, modeling results suggest that larger SRBs are not, under the most commonly observed low-energy wave conditions, likely to move very far alongshore. This finding suggests that SRBs from one source location may not (outside of storm conditions) be redistributed to other up or down coast locations. This information can guide operational response decisions. In addition, because SRBs are less mobile compared with sand, they are likely to become buried and unburied under normal sand transport processes thereby lengthening the time SRBs may take to move onshore. The rate of onshore movement was not specifically addressed by this study, yet the results resolve the cross-shore domain and cross-shore variations in alongshore transport that are relevant to achieving the primary objectives. Furthermore, during infrequent events (for example, winter storms and severe meteorological events such as Hurricane Isaac of August 2012), energy is shown to be sufficient to move a greater range of SRB sizes and potentially expose and break up submerged oil mats. When SRBs do move alongshore, the models indicate that there are regions that are more conducive to accumulation of SRB material than others. Accumulation can occur where there are reversals and decelerations in alongshore currents and where forces created by shear stress drops below critical thresholds to maintain or initiate SRB movement. In addition, flow and SRB mobility patterns around inlets indicate patterns in hydrodynamic forces that influence redistribution of SRBs and the surface oil that mixed with sediment to form oil mats in the first place.
Executive Summary The Hiram M. Chittenden (Ballard) Locks and Lake Washington Ship Canal connect freshwater Lake Washington and saline Shilshole Bay of Puget Sound in Seattle, Washington. The locks and canal allow for ships to traverse this reach. Anadromous salmonids also migrate through, transitioning between saline and freshwater environments, and making use of a fish ladder at the locks when traveling upstream. WEST Consultants, Inc., constructed a two-dimensional hydrodynamic and water-quality model (CE-QUAL-W2) simulating flow, water temperature, and salinity for the Ballard Locks and the Lake Washington Ship Canal. An initial model was built for calendar years 2014–15, and the model was updated using a more recent and modern dataset for calendar years 2016–20. The U.S. Army Corps of Engineers requested that the U.S. Geological Survey review this model and its documentation to evaluate the technical aspects of its development and calibration. Findings from this review include the following: Overall, the Lake Washington Ship Canal CE-QUAL-W2 model was well-documented and constructed largely following typical model-development methods. The Lake Washington Ship Canal model was built with CE-QUAL-W2 model version 4.5, compiled and released by Portland State University in April 2021. CE-QUAL-W2 updates and improvements are regularly released with bug fixes and new features, so any model updates would benefit from the use of the most-recent software release. The model grid that represents the Lake Washington Ship Canal bathymetry was 9.2 kilometers (5.7 miles) long, matching the expected length of the waterway. The deepest model segments were near sampling site LLLW (Large Locks site) near the locks. Lake Union is reported to constitute most of the volume of the Lake Washington Ship Canal and is depicted as such in the model grid. The model includes several water outflows at Ballard Locks, including the large and small locks, a saltwater drain, a spillway, smolt flumes, and a fish ladder. Flows from the spillway, smolt flumes, and fish ladder were combined into one structure outflow in the model and assigned one withdrawal elevation from the Lake Washington Ship Canal. The smolt flume and spillway withdraw from the same elevation, but the fish ladder flow withdraws from a higher elevation in Lake Washington Ship Canal, and that flow could be separated into its own withdrawal. The model input files were created using the Coordinated Universal Time standard instead of the more typical choice of using local standard time. This is not incorrect, but sub-daily results would need to be converted to local time for science-communication purposes. The meteorological dataset had some unexpected anomalies, such as a baseline shift in the wind-speed dataset. Other nearby meteorological datasets could be used instead or used to correct the current meteorological inputs. The upstream boundary was configured with water-temperature data from a continuous monitor buoy in Lake Washington. The boundary salinity was set at 0 parts per thousand for the duration of the model simulation. A more realistic estimate of salinity at the upstream boundary could be constructed using data from the same buoy. Saline inflow at the downstream boundary of the Lake Washington Ship Canal model through lock exchanges at the large lock was included as a tributary in the model. Salinity and temperature inputs in this tributary at the large locks were set as constant values for the entire simulation. Saline inflow through the small lock was not included in the model because few data were available, and the input was likely to be small because of the smaller surface area and volume of the small lock relative to the large lock. The model did not include any flow, water temperature, or salinity inputs to the Lake Washington Ship Canal other than at the locks and at the upstream boundary. Any point sources, small tributaries, or stormwater inputs were omitted from the model. It is unclear whether this is a substantial omission relative to model results. Most model parameters were set as defaults or to reasonable values. However, the value of the WINDH parameter, the height of the wind speed measurement, was different than the height of the meteorological site. Compared to measured data, the model simulated water-surface elevations and water temperatures with reasonable accuracy. Differences in the modeled and measured salinities revealed some opportunities to improve the simulation of salinity, both baseline salinity and the salinity maxima in summer and autumn.
This map portrays the geodynamics of Northeast Asia at a scale of 1:5,000,000 using the concepts of plate tectonics and analysis of terranes and overlap assemblages. The map is the result of a detailed compilation and synthesis at 5 million scale and is part of a major international collaborative study of the Mineral Resources, Metallogenesis, and Tectonics of Northeast Asia conducted from 1997 through 2002 by geologists from earth science agencies and universities in Russia, Mongolia, Northeastern China, South Korea, Japan, and the USA. This map is the result of extensive geologic mapping and associated tectonic studies in Northeast Asia in the last few decades and is the first collaborative compilation of the geology of the region at a scale of 1:5,000,000 by geologists from Russia, Mongolia, Northeastern China, South Korea, Japan, and the USA. The map was compiled by a large group of international geologists using the below concepts and definitions during collaborative workshops over a six-year period. The map is a major new compilation and re-interpretation of pre-existing geologic maps of the region. The map is designed to be used for several purposes, including regional tectonic analyses, mineral resource and metallogenic analysis, petroleum resource analysis, neotectonic analysis, and analysis of seismic hazards and volcanic hazards. The map consists of two sheets. Sheet 1 displays the map at a scale of 1:5,000,000, explanation. Sheet 2 displays the introduction, list of map units, and source references. Detailed descriptions of map units and stratigraphic columns are being published separately. This map is one of a series of publications on the mineral resources, metallogenesis, and geodynamics,of Northeast Asia. Companion studies and other articles and maps , and various detailed reports are: (1) a compilation of major mineral deposit models (Rodionov and Nokleberg, 2000; Rodionov and others, 2000; Obolenskiy and others, in press a); (2) a series of metallogenic belt maps (Obolenskiy and others, 2001; in press b); (3) a lode mineral deposits and placer districts location map for Northeast Asia (Ariunbileg and others, in press b); (4) descriptions of metallogenic belts (Rodionov and others, in press); and (5) a database on significant metalliferous and selected nonmetalliferous lode deposits, and selected placer districts (Ariunbileg and others, in press a).