USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computers and Geology”

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 1,333 records · Page 74Linked to original sources
USGS2019

Isotopic ratios of Saturn's rings and satellites: Implications for the origin of water and Phoebe

Isotopic ratios have long been used to learn about physical processes acting over a wide range of geological environments, and in constraining the origin and/or evolution of planetary bodies. We report the spectroscopic detection of deuterium in Saturn's rings and satellites, and use these measurements to determine the (D/H) ratios in their near-surface regions. Saturn's moons, Phoebe and Iapetus, show a strong signature of CO2 and the 13C component of this molecule is detected and quantified. Large averages of spectra obtained by the Cassini Visual and Infrared Mapping Spectrometer, VIMS, were computed for the rings and icy satellites. The observed intensities of the infrared absorptions in H2O and CO2 and their isotopes were calibrated using laboratory data and radiative transfer models to derive the D/H and 13C/12C ratios. We find that the D/H in Saturn's rings and satellites is close to the Vienna Standard Mean Ocean Water (VSMOW) and bulk Earth (4% lower than VSMOW) value except for Phoebe, which is 8.3 times the VSMOW value. This is the highest value for any Solar-System surface yet measured, and suggests that Phoebe formed from material with a different D/H ratio than the other satellites in the Saturn system. Phoebe’s 13C/12C ratio is also unusual: 4.7 times greater than terrestrial, and greater than values measured for the interstellar medium and the galactic center. The high 13C abundance in the CO2 suggests that Phoebe was never warm enough for the large D/H ratio in its surface to have originated by evaporative fractionation of its waterice (e.g., from heating in the inner Solar System before its eventual capture by Saturn). We also report the detection of a probable O-D stretch absorption due to OD in minerals on Phoebe at 3.62 μm. This absorption is not detected on other Saturnian satellites. Stronger signatures of bound water absorptions are found in the dark material of Iapetus and we report a new detection of bound water at 1.9 μm. The position of this absorption matches that seen in spectra of hydrated iron oxides but does not match absorptions seen in spectra of tholins. Despite the strong bound water signature in the Iapetus dark material, no 3.62-μm OD absorption is seen in the spectra, further indicating the high deuterium level on Phoebe is unusual. As such, it is likely that Phoebe originated in a colder part of the outer Solar System, relative to the prevailing temperatures at Saturn’s distance from the Sun.

Icarus↗
USGS2005

An efficient algorithm for double-difference tomography and location in heterogeneous media, with an application to the Kilauea volcano

Improving our understanding of crustal processes requires a better knowledge of the geometry and the position of geological bodies. In this study we have designed a method based upon double-difference relocation and tomography to image, as accurately as possible, a heterogeneous medium containing seismogenic objects. Our approach consisted not only of incorporating double difference in tomography but also partly in revisiting tomographic schemes for choosing accurate and stable numerical strategies, adapted to the use of cross-spectral time delays. We used a finite difference solution to the eikonal equation for travel time computation and a Tarantola-Valette approach for both the classical and double-difference three-dimensional tomographic inversion to find accurate earthquake locations and seismic velocity estimates. We estimated efficiently the square root of the inverse model's covariance matrix in the case of a Gaussian correlation function. It allows the use of correlation length and a priori model variance criteria to determine the optimal solution. Double-difference relocation of similar earthquakes is performed in the optimal velocity model, making absolute and relative locations less biased by the velocity model. Double-difference tomography is achieved by using high-accuracy time delay measurements. These algorithms have been applied to earthquake data recorded in the vicinity of Kilauea and Mauna Loa volcanoes for imaging the volcanic structures. Stable and detailed velocity models are obtained: the regional tomography unambiguously highlights the structure of the island of Hawaii and the double-difference tomography shows a detailed image of the southern Kilauea caldera-upper east rift zone magmatic complex. Copyright 2005 by the American Geophysical Union.

Hawaii↗
USGS2011

Hawaii Volcano Observatory 75th anniversary

The 75th anniversary of the founding of the U.S. Geological Survey (USGS) Hawaiian Volcano Observatory (HVO) was celebrated during January 1987. The festivities began on January 9 with the opening in Hilo of a major exhibit at the Wailoa Center on the current work of HVO, its history, and its special relationship to Hawaii Volcanoes National Park. In addition to spectacular photographs of volcanic activity and HVO field work, the exhibit featured continuous showings of underwater footage taken during the recent eruptive activity (December 1986) and a computer monitor that displayed the activity of HVO seismic net in real time. The exhibit, which was coordinated by HVO librarian/photo archivist Taeko Jane Takahashi and which featured the professional work of HVO photographer J . D. Griggs, was well received and will probably form the core of a future traveling exhibit.

Hawaii↗
USGS2016

User guide for MODPATH Version 7—A particle-tracking model for MODFLOW

MODPATH is a particle-tracking post-processing program designed to work with MODFLOW, the U.S. Geological Survey (USGS) finite-difference groundwater flow model. MODPATH version 7 is the fourth major release since its original publication. Previous versions were documented in USGS Open-File Reports 89–381 and 94–464 and in USGS Techniques and Methods 6–A41. MODPATH version 7 works with MODFLOW-2005 and MODFLOW–USG. Support for unstructured grids in MODFLOW–USG is limited to smoothed, rectangular-based quadtree and quadpatch grids. A software distribution package containing the computer program and supporting documentation, such as input instructions, output file descriptions, and example problems, is available from the USGS over the Internet ( http://water.usgs.gov/ogw/modpath/) .

Open-File Report↗
USGS1975

Abundances of uranium, thorium, and potassium for some Australian crystalline rocks

This report contains a tabulation of the basic radioelement and radiogenic heat data obtained during an Australian National University (ANU) - United States Geological Survey (USGS) heat-flow project, directed jointly by J. C. Jaeger (ANU) and J. H. Sass (USGS). Most samples were collected during the periods June through September, 1971 and 1972. The measurements were made subsequently by two of us (C. M. Bunker and C. A. Bush) using the gamma-ray spec trometric techniques described by Bunker and Bush (1966, 1967). Interpreting the spectra for quantitative analyses of the radioelements was accomplished with an iterative leastsquares computer program modified from one by Schonfeld (1966). Uranium content determined by gamma-ray spectrometry is based on a measurement of the daughter products of 226 Ra. Equilibrium in the uranium-decay series was assumed for these analyses . Throughout the report, when U content is stated, radium-equivalent uranium is implied. The coefficient of variation for the accuracy of the radioelement data, when compared to ana lyses by isotope dilution and flame photometry is about 3 percent for radium-equivalent uranium and thorium and about 1 percent for potassium. These percentages are in addition to minimum standard deviations of about 0.05 ppm for U and Th, and about 0.03 percent for K.

Open-File Report↗
USGS2005

Implementation and use of direct-flow connections in a coupled ground-water and surface-water model

The U.S. Geological Survey's MODFLOW finite-difference ground-water flow model has been coupled with three surface-water packages - the MODBRANCH, River, and Stream packages - to simulate surface water and its interaction with ground water. Prior to the development of the coupling packages, the only interaction between these modeling packages was that leakage values could be passed between MODFLOW and the three surface-water packages. To facilitate wider and more flexible uses of the models, a computer program was developed and added to MODFLOW to allow direct flows or stages to be passed between any of the packages and MODFLOW. The flows or stages calculated in one package can be set as boundary discharges or stages to be used in another package. Several modeling packages can be used in the same simulation depending upon the level of sophistication needed in the various reaches being modeled. This computer program is especially useful when any of the River, Stream, or MODBRANCH packages are used to model a river flowing directly into or out of wetlands in direct connection with the aquifer and represented in the model as an aquifer block. A field case study is shown to illustrate an application.

Groundwater↗
USGS1979

Streamflow statistical summaries for Colorado streams through September 30, 1975; Volume 2: Colorado River basin above Gunnison River

This report contains statistical summaries of daily streamflow data as of September 30, 1975, for 189 stations west of the Continental Divide in Colorado for use by agencies and individuals engaged in water studies. To be included in this report, a station had to have a minimum of 3 complete water years (October 1 to September 30) of record and 3 complete climatic years (April 1 to March 31) of record having virtually constant effects of regulation and diversion, and to have had streamflow records previously reviewed and published by the U.S. Geological Survey. Duration tables show the distribution of daily discharges for each water year ending September 30. High-flow sequence tables show the highest mean discharge for 1, 3, 7, 15, 30, 60, 90, 120, and 183 consecutive days of each water year. Low-flow sequence tables show the lowest mean discharge for 1, 3, 7, 14, 30, 60, 90, 120, and 183 consecutive days of each climatic year. For monthly flows the following statistics have been computed for each station having 5 or more complete water years of record: (1) Mean, (2) variance, (3) standard deviation, (4) skewness, (5) coefficient of variation, and (6) percentage of average flow. The first five items were also computed for the water year annual flows along with the first order serial correlation coefficient. For stations which have had virtually instantaneous significant changes in the effects of regulation or diversion, a maximum of two sets of summaries are shown. Referenced reports describe procedures for processing and interpreting the streamflow statistics.

Colorado↗
USGS1982

Catalog of available clear mylar templates used to determine latitude and longitude of sample localities between the latitudes of 51 degrees 00' and 71 degrees 30' at the scale of 1:63,360 and between the latitude of 49 degrees 00'00" north or south at the scale of 1:250,000

The present use of electro-mechanical plotting of sample points can give rise to location errors, particularly when map scale changes are involved, unless careful control is maintained. The following catalog of computer-generated mylar templates, with plots of latitude and longitude (longitudes equidistant north to south), sufficient to cover the latitudes between 49?00'00'' and 71?30'00'' N. or S., will help facilitate the plotting of points when there is a map scale change for two different map scales (1:63,360 and 1:250,000). Using the template overlays (available as individual page numbers in Open-File Report 82-724--xerographic clear film overlays from U.S. Geological Survey Open-File Services Section, Branch of Distribution, USGS, Box 25425, Denver, Colorado 80225, telephone (303) 234-5888), the point plotted on a composite final copy may be checked against field maps and notes when only degrees latitude and longitude are given.

Open-File Report↗
USGS1979

Streamflow statistical summaries for Colorado streams through September 30, 1975; Volume 3, Colorado River basin from Gunnison River to San Juan River

This report contains statistical summaries of daily streamflow data as of September 30, 1975, for 234 stations west of the Continental Divide in Colorado and adjacent States for use by agencies and individuals engaged in water studies. To be included in this report, a station had to have a minimum of 3 complete water years (October 1 to September 30) of record and 3 complete climatic years (April I to March 31) of record having virtually constant effects of regulation and diversion, and to have had streamflow records previously reviewed and published by the U.S. Geological Survey. Duration tables show the distribution of daily discharges for each water year ending September 30. High-flow sequence tables show the highest mean discharge for 1, 3, 7, 15, 30, 60, 90, 120, and 183 consecutive days of each water year. Low-flow sequence tables show the lowest mean discharge for 1, 3, 7, 14, 30, 60, 90, 120, and 183 consecutive days of each climatic year. For monthly flows the following statistics have been computed for each station having 5 or more complete water years of record: (1) mean, (2) variance, (3) standard deviation, (4) skewness, (5) coefficient of variation, and (6) percentage of average flow. The first five items were also computed for the water year annual flows along with the first order serial correlation coefficient. For stations which have had virtually instantaneous significant changes in the effects of regulation and diversion, a maximum of two sets of summaries are shown. Referenced reports describe procedures for processing and interpreting the streamflow statistics.

Open-File Report↗
USGS2001

Concentrations of Escherichia coli in streams in the Kankakee and lower Wabash River watersheds in Indiana, June-September 1999

Water samples collected from 58 surface- water sites in the Kankakee and Lower Wabash River Watersheds from June through September 1999 were analyzed for concentrations of Escherichia coli bacteria. Each sitewas sampled five times in a 30-day period. Twentynine sites were sampled during June and July, and 29 different sites were sampled during August and September. A five-sample geometric mean of concentrations was computed for each site. Concentrations of Escherichia coli (E. coli) in 126 of the 289 samples exceeded the State of Indiana single-sample standard of 235 colonies per 100 milliliters for waters used for recreation. Concentrations in samples from 38 of the 58 sites exceeded the State of Indiana standard for a five-sample geometric mean of 125 colonies per 100 milliliters for waters used for recreation. Ten of the 58 sites were at or near U.S. Geological Survey streamflow-gaging stations. Based on records from the streamflowgaging stations, 18 percent of the samples collected at these sites were collected at streamflows above the median daily discharge for each station. E. coli concentrations and turbidity measurements collected during 1999 were analyzed in concert with similar concentration and turbidity data collected in 1998 at streams within the Upper Wabash River Watershed in Indiana to investigate the relation between concentrations of bacteria and turbidity. The analysis indicated a statistically significant correlation between concentrations of E. coli and turbidity. If the turbidity was greater than 83 nephelometric turbidity units, the E. coli concentration always exceeded the singlesample standard. If, however, the turbidity was less than 83 nephelometric turbidity units, concentrations of E. coli were not always below the single-sample standard.

Illinois, Indiana↗
USGS2011

Programming PHREEQC calculations with C++ and Python a comparative study

The new IPhreeqc module provides an application programming interface (API) to facilitate coupling of other codes with the U.S. Geological Survey geochemical model PHREEQC. Traditionally, loose coupling of PHREEQC with other applications required methods to create PHREEQC input files, start external PHREEQC processes, and process PHREEQC output files. IPhreeqc eliminates most of this effort by providing direct access to PHREEQC capabilities through a component object model (COM), a library, or a dynamically linked library (DLL). Input and calculations can be specified through internally programmed strings, and all data exchange between an application and the module can occur in computer memory. This study compares simulations programmed in C++ and Python that are tightly coupled with IPhreeqc modules to the traditional simulations that are loosely coupled to PHREEQC. The study compares performance, quantifies effort, and evaluates lines of code and the complexity of the design. The comparisons show that IPhreeqc offers a more powerful and simpler approach for incorporating PHREEQC calculations into transport models and other applications that need to perform PHREEQC calculations. The IPhreeqc module facilitates the design of coupled applications and significantly reduces run times. Even a moderate knowledge of one of the supported programming languages allows more efficient use of PHREEQC than the traditional loosely coupled approach.

Conference Paper↗
USGS1961

Drainage basins, channels, and flow characteristics of selected streams in central Pennsylvania

The hydraulic, basin, and geologic characteristics of 16 selected streams in central Pennsylvania were measured for the purpose of studying the relations among these general characteristics and their process of development. The basic parameters which were measured include bankfull width and depth, channel slope, bed material size and shape, length of stream from drainage divide, and size of drainage area. The kinds of bedrock over which the streams flow were noted. In these streams the bankfull channel is filled by flows approximating the 2.3-year flood. By measuring the breadth and mean depth of the channel, it was possible to compute the bankfull mean velocity for each of the 119 sampling stations. These data were then used to compute the downstream changes in hydraulic geometry of the streams studied. This method has been called an indirect computation of the hydraulic geometry. The results obtained by the indirect method are similar to those of the direct method of other workers. The basins were studied by examining the relations of drainage area, discharge, and length of stream from drainage divide. For the streams investigated, excellent correlations were found to exist between drainage area and the 2.3-year flood, as well as between length of stream from the basin divide and drainage area. From these correlations it is possible to predict the discharge for the 2.3-year flood at any arbitrary point along the length of the stream. The long, intermediate, and short axes of pebbles sampled from the bed of the stream were recorded to study both size and sphericity changes along individual streams and among the streams studied. No systematic downstream changes in sphericity were found. Particle size changes are erratic and show no consistent relation to channel slope. Particle size decreases downstream in many streams but remains constant or increases in others. Addition of material by tributaries is one factor affecting particle size and another is the parent material. Wear does not appear to account for some of the changes noted in particle size in a downstream direction. Comparison with laboratory studies indicates that at least in some streams the downstream decrease in size is much greater than would be expected from wear alone. The type of bedrock underlying the channels included in this study appears to affect both channel slope and particle size. For a given length of stream, a stream channel underlain by sandstone tends to have a steeper slope and larger bed material than channels underlain by shale or limestone. Hence, a stream which heads in sandstone and ends in limestone tends to have a more rapid decrease in slope and particle size than a stream heading in limestone and ending in sandstone. The association of steep slopes and small particles for limestone channels implies that slope and particle size may show a vague correlation between lithologic groups although no correlation may exist within a given lithologic type. In addition to the effect of bedrock on slope and particle size, there is some evidence that channels in limestone or dolomite have a slightly smaller cross section at bankfull stage than channels in shale or sandstone. Near the headwaters of many of these streams, a deposit of periglacial rubble affects the slope and bed material size. Some of the debris contains residual boulders which are too large to be moved by ordinary floods and, therefore, impose larger particle sizes in the bed of the stream. The addition of this very coarse debris to the bed material is another example of the influence of geologic factors on stream channels even though the channel consists of unconsolidated debris instead of bedrock. The influence of geologic factors noted in selected streams in central Pennsylvania may not be directly applicable to areas other than the Appalachian Mountains, but the general process is no doubt similar in most areas. In large alluvial valleys bedrock cannot be much of an influencing factor; yet large, thick alluvial deposits and terraces are in a sense "bedrock" materials upon which the stream works to form the landscape.

Pennsylvania↗
USGS2016

The global Landsat archive: Status, consolidation, and direction

New and previously unimaginable Landsat applications have been fostered by a policy change in 2008 that made analysis-ready Landsat data free and open access. Since 1972, Landsat has been collecting images of the Earth, with the early years of the program constrained by onboard satellite and ground systems, as well as limitations across the range of required computing, networking, and storage capabilities. Rather than robust on-satellite storage for transmission via high bandwidth downlink to a centralized storage and distribution facility as with Landsat-8, a network of receiving stations, one operated by the U.S. government, the other operated by a community of International Cooperators (ICs), were utilized. ICs paid a fee for the right to receive and distribute Landsat data and over time, more Landsat data was held outside the archive of the United State Geological Survey (USGS) than was held inside, much of it unique. Recognizing the critical value of these data, the USGS began a Landsat Global Archive Consolidation (LGAC) initiative in 2010 to bring these data into a single, universally accessible, centralized global archive, housed at the Earth Resources Observation and Science (EROS) Center in Sioux Falls, South Dakota. The primary LGAC goals are to inventory the data held by ICs, acquire the data, and ingest and apply standard ground station processing to generate an L1T analysis-ready product. As of January 1, 2015 there were 5,532,454 images in the USGS archive. LGAC has contributed approximately 3.2 million of those images, more than doubling the original USGS archive holdings. Moreover, an additional 2.3 million images have been identified to date through the LGAC initiative and are in the process of being added to the archive. The impact of LGAC is significant and, in terms of images in the collection, analogous to that of having had two additional Landsat-5 missions. As a result of LGAC, there are regions of the globe that now have markedly improved Landsat data coverage, resulting in an enhanced capacity for mapping, monitoring change, and capturing historic conditions. Although future missions can be planned and implemented, the past cannot be revisited, underscoring the value and enhanced significance of historical Landsat data and the LGAC initiative. The aim of this paper is to report the current status of the global USGS Landsat archive, document the existing and anticipated contributions of LGAC to the archive, and characterize the current acquisitions of Landsat-7 and Landsat-8. Landsat-8 is adding data to the archive at an unprecedented rate as nearly all terrestrial images are now collected. We also offer key lessons learned so far from the LGAC initiative, plus insights regarding other critical elements of the Landsat program looking forward, such as acquisition, continuity, temporal revisit, and the importance of continuing to operationalize the Landsat program.

Remote Sensing of Environment↗
USGS2020

Partly cloudy with a chance of lava flows: Forecasting volcanic eruptions in the 21st century

A primary goal of volcanology is forecasting hazardous eruptive activity. Despite much progress over the last century, however, volcanoes still erupt with no detected precursors, lives and livelihoods are lost to eruptive activity, and forecasting the onsets of eruptions remains fraught with uncertainty. Long‐term forecasts are generally derived from the geological and historical records, from which recurrence intervals and styles of activity can be inferred, while shorter‐term forecasts are derived from patterns in monitoring data. Information from geology and monitoring data can be evaluated and combined using statistical analysis, expert elicitation, and conceptual and or mathematical models. Integrative frameworks, such as event trees, combine this diversity of information to produce probabilistic forecasts that can inform the style and scale of the societal response to a potential future eruption. Several developments show promise to revolutionize the utility and accuracy of these forecasts. These include growth in the quantity and quality of multidisciplinary monitoring data, coupled with increases in computing power; machine learning algorithms, which will allow far better utilization of this growing volume of data; and new physiochemical volcano models and data assimilation algorithms, which take advantage of a wide range of monitoring data and realistic physics to better predict the evolution of a given physical state. Although eruption forecasts may never be as generally reliable as weather forecasts, and great caution must be exercised when attempting to predict highly complex volcanic behavior, these and other innovations—particularly when combined in integrative, fully probabilistic forecasting frameworks—should help volcanologists to better issue warnings of volcanic activity on societally relevant time frames.

Hawaii↗
USGS2021

Regression relations and long-term water-quality constituent concentrations, loads, yields, and trends in the North Fork Ninnescah River, south-central Kansas, 1999–2019

Cheney Reservoir, in south-central Kansas, is the primary water supply for the city of Wichita, Kansas. The North Fork Ninnescah River is the largest tributary to Cheney Reservoir and contributes about 70 percent of the inflow. The U.S. Geological Survey, in cooperation with the City of Wichita, has been continuously monitoring water quality (including water temperature, specific conductance, pH, dissolved oxygen, and turbidity) on the North Fork Ninnescah River upstream from Cheney Reservoir (U.S. Geological Survey site 07144780) since November 1998. Continued data collection would be beneficial to update and describe changing water-quality conditions in the drainage basin and in the reservoir over time. Regression models were developed to describe relations between discretely measured constituent concentrations and continuously measured physical properties. The models updated in this report include total suspended solids (TSS), suspended-sediment concentration (SSC), nitrate plus nitrite, nitrate, orthophosphate (OP), total phosphorus (TP), and total organic carbon (TOC). Daily computed concentrations for TSS, TP, and nitrate plus nitrite during 1999–2019 were compared with Cheney Reservoir Task Force (CRTF) goals for base-flow and runoff conditions. CRTF goals for base-flow concentrations were exceeded more frequently (70 to 99.9 percent of the time) than runoff goals (0 to 11 percent of the time). Except for 2012, annual mean TSS concentrations exceeded the base-flow goal every year. Nitrate plus nitrite and TP annual mean concentrations exceeded the base-flow goals every year. TSS and nitrate plus nitrite annual mean concentrations during runoff conditions never exceeded the CRTF runoff goal. TP annual mean concentrations during runoff conditions only exceeded the CRTF runoff goal during 2002. Sedimentation is progressively reducing the storage capacity of Cheney Reservoir. During 1999–2019, 55 percent of the computed suspended-sediment load was transported during the top 1 percent of loading days (76 days); 22 percent of the total load was transported in the top 10 loading days, indicating that substantial parts of suspended-sediment loads continue to be delivered during disproportionately small periods in Cheney Reservoir. Successful sediment management efforts necessitate reduction techniques that account for these large load events. Flow-normalized concentrations and fluxes were computed during 1999 through 2019 using Weighted Regressions on Time, Discharge, and Season (WRTDS) statistical models and WRTDS bootstrap tests. Flow-normalized concentrations of TSS, SSC, OP, TP, and TOC had upward trend probabilities; conversely, nitrate plus nitrite had a downward trend. Flow-normalized fluxes for OP, TP, and TOC had an upward trend. No discernible patterns were identified for flow-normalized flux of TSS or suspended sediment. Nitrate plus nitrite flow-normalized flux indicated a downward trend. Flow-normalized concentrations for TSS were less than the CRTF long-term goal of 100 milligrams per liter (mg/L), but the upward trend indicated the long-term goal may be exceeded if no changes are made. Flow-normalized TP concentrations exceeded the CRTF long-term goal (0.1 mg/L) and were assigned a very likely upward trend. Flow-normalized nitrate plus nitrite concentrations exceeded the CRTF long-term goal of 1.2 mg/L during the beginning of the study period, then were less than the CRTF goal for the remainder of the study; however, during 2010–19 flow-normalized concentrations increased by 6 percent. Linking water-quality changes to causal factors requires consistent monitoring before, during, and after changes; this presents challenges related to length and frequency of data collection and available concomitant land-use and conservation practice data. As such, attribution of water-quality trends to land-use changes or conservation practices was not possible for this study because of a lack of land-use and conservation practice data. Additionally, because precipitation frequency and intensity are projected to continue to increase in the Great Plains region, accounting for extreme episodic events may be an important consideration in future sediment and nutrient load reduction plans.

Kansas↗
USGS2023

Documentation of linear regression models for computing water-quality constituent concentrations using continuous real-time water-quality data for the North Fork Ninnescah River and Cheney Reservoir, Kansas, 2014–21

Cheney Reservoir, in south-central Kansas, was constructed to provide a reliable municipal water supply for the city of Wichita, Kansas, and to provide downstream flood control, wildlife habitat, and recreation. Cheney Reservoir will continue to be important for municipal water supply use as needs increase with ongoing population growth and urban development. Advanced notification of changing water-quality conditions near water-treatment facility intakes and in source waters allows water-treatment facilities and resource planning officials to proactively monitor changing conditions. The U.S. Geological Survey (USGS), in cooperation with the City of Wichita, collected water-quality data at the North Fork Ninnescah River above Cheney Reservoir (USGS station 07144780) and Cheney Reservoir near Cheney, Kans. (USGS station 07144790), monitoring sites to update and develop regression models relating continuous water-quality constituents, streamflow, reservoir storage, and seasonal components to discretely sampled water-quality constituent concentrations of interest. Linear regression analysis was used to update and develop models for alkalinity, major ions, nutrients (nitrogen and phosphorus species), total and dissolved organic carbon, total suspended solids, suspended sediment, fecal indicator bacteria, and atrazine at the North Fork Ninnescah River site and total dissolved solids, major ions, hardness as calcium carbonate, nutrients (nitrogen and phosphorus species), chlorophyll a , and suspended sediment at the Cheney Reservoir site. New and updated models for both sites are applicable to the period of YSI EXO water-quality monitor and sensor deployment (November 14, 2015, through September 30, 2021, at the North Fork Ninnescah River site; October 1, 2014, through September 30, 2021, at the Cheney Reservoir site). Models and resulting water-quality information included in this report can be used in real time, potentially as guidance for water-treatment processes, and can be used to characterize changes in water-quality conditions over time in Cheney Reservoir and its contributing drainage basin provided that the deployed equipment, sensors, and location do not change.

Kansas↗
USGS2023

Calibration of imperfect geophysical models by multiple satellite interferograms with measurement bias

Model calibration consists of using experimental or field data to estimate the unknown parameters of a mathematical model. The presence of model discrepancy and measurement bias in the data complicates this task. Satellite interferograms, for instance, are widely used for calibrating geophysical models in geological hazard quantification. In this work, we used satellite interferograms to relate ground deformation observations to the properties of the magma chamber at K i¯lauea Volcano in Hawai‘i. We derived closed-form marginal likelihoods and implemented posterior sampling procedures that simultaneously estimate the model discrepancy of physical models, and the measurement bias from the atmospheric error in satellite interferograms. We found that model calibration by aggregating multiple interferograms and downsampling the pixels in the interferograms can reduce the computation complexity compared to calibration approaches based on multiple data sets. The conditions that lead to no loss of information from data aggregation and downsampling are studied. Simulation illustrates that both discrepancy and measurement bias can be estimated, and real applications demonstrate that modeling both effects helps obtain a reliable estimation of a physical model’s unobserved parameters and enhance its predictive accuracy. We implement the computational tools in the RobustCalibration package available on CRAN.

Technometrics↗
USGS2004

SutraGUI, a graphical-user interface for SUTRA, a model for ground-water flow with solute or energy transport

This report describes SutraGUI, a flexible graphical user-interface (GUI) that supports two-dimensional (2D) and three-dimensional (3D) simulation with the U.S. Geological Survey (USGS) SUTRA ground-water-flow and transport model (Voss and Provost, 2002). SutraGUI allows the user to create SUTRA ground-water models graphically. SutraGUI provides all of the graphical functionality required for setting up and running SUTRA simulations that range from basic to sophisticated, but it is also possible for advanced users to apply programmable features within Argus ONE to meet the unique demands of particular ground-water modeling projects. SutraGUI is a public-domain computer program designed to run with the proprietary Argus ONE? package, which provides 2D Geographic Information System (GIS) and meshing support. For 3D simulation, GIS and meshing support is provided by programming contained within SutraGUI. When preparing a 3D SUTRA model, the model and all of its features are viewed within Argus 1 in 2D projection. For 2D models, SutraGUI is only slightly changed in functionality from the previous 2D-only version (Voss and others, 1997) and it provides visualization of simulation results. In 3D, only model preparation is supported by SutraGUI, and 3D simulation results may be viewed in SutraPlot (Souza, 1999) or Model Viewer (Hsieh and Winston, 2002). A comprehensive online Help system is included in SutraGUI. For 3D SUTRA models, the 3D model domain is conceptualized as bounded on the top and bottom by 2D surfaces. The 3D domain may also contain internal surfaces extending across the model that divide the domain into tabular units, which can represent hydrogeologic strata or other features intended by the user. These surfaces can be non-planar and non-horizontal. The 3D mesh is defined by one or more 2D meshes at different elevations that coincide with these surfaces. If the nodes in the 3D mesh are vertically aligned, only a single 2D mesh is needed. For nonaligned meshes, two or more 2D meshes of similar connectivity are used. Between each set of 2D meshes (and model surfaces), the vertical space in the 3D mesh is evenly divided into a user-specified number of layers of finite elements. Boundary conditions may be specified for 3D models in SutraGUI using a variety of geometric shapes that may be located freely within the 3D model domain. These shapes include points, lines, sheets, and solids. These are represented by 2D contours (within the vertically-projected Argus ONE view) with user-defined elevations. In addition, boundary conditions may be specified for 3D models as points, lines, and areas that are located exactly within the surfaces that define the model top and the bottoms of the tabular units. Aquifer properties may be specified separately for each tabular unit. If the aquifer properties vary vertically within a unit, SutraGUI provides the Sutra_Z function that can be used to specify such variation.

Open-File Report↗