Geology Reports⌕ Search

SEARCH · Geology Reports

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 433 records · Page 24Linked to original sources
USGS2023

Geology and undiscovered resource assessment of the potash-bearing, Middle Devonian (Givetian), Prairie Evaporite, Elk Point Basin, Canada and United States

The U.S. Geological Survey (USGS) assessed undiscovered potash resources in the Elk Point Basin in Canada and the United States as part of a global mineral resource assessment. The Elk Point Basin is a large, Middle Devonian (Givetian) intracratonic evaporite basin covering approximately 1,200,000 square kilometers (km 2 ) and filled mainly with marine evaporite and minor clastic sedimentary rocks that contain stratabound potash-bearing salt. The potash-bearing salt is concentrated in four stratigraphic members (Patience Lake, Belle Plaine, White Bear, and Esterhazy) in the upper 100 meters (m) of the Prairie Evaporite and are separated by beds of halite (NaCl) that contain lesser—presently non-economic—amounts of sylvite (KCl) and carnallite (KMgCl 3 ·6H 2 O). The principal ore-bearing salt contains mainly sylvite. Four permissive tracts were defined that permit the presence of undiscovered stratabound potash (both sylvite- and carnallite-bearing salt) using geological criteria. Permissive tracts are defined by the spatial extent of each stratigraphic member that is at least 1 m thick, are less than 3 kilometers (km) from the surface, contain at least 4 percent equivalent potassium oxide (K 2 O), and contain the currently known resources. The permissive tracts include known potash deposits and potash occurrences as wells or mines not in production and show where undiscovered potash resources may be present. Well data are used to define the extent, thickness, average K 2 O equivalent grades, and volumes of each member. Data were supplied by the Saskatchewan Geological Survey or were obtained from published National Instrument (NI) 43-101 technical reports and other published reports, such as annual 10-K reports or news releases. The Elk Point Basin is the world’s largest source of potash, producing 23.0 million metric tons (Mt) of potassium chloride (KCl) (the equivalent of about 14.4 Mt of K 2 O) in 2018. In terms of global importance, the Elk Point Basin may contain 40 to greater than 50 percent of the world’s potash resources. Since 1962, potash companies have mined more than 1.5 trillion metric tons of ore containing 605 Mt of KCl (the equivalent of about 380 Mt of K 2 O). The total value of the ore produced through 2018 is on the order of $70 trillion (CAD). Potash is currently produced from eight conventional and three underground solution mines at depths ranging from 900 m to nearly 1,800 m. Estimates of the amount of potash in the Elk Point Basin vary considerably and the data and methods used in those estimations are not well documented. Known potash resources are approximately 99 billion metric tons (Bt) of ore containing 22 Bt of K 2 O equivalent. As a result of new mine openings and increased production capacity at existing mines, the total production capacity of mines in the Elk Point Basin has increased significantly (to about 32.8 Mt of KCl or 22.8 Mt of K 2 O equivalent per year). Additional production capacity of about 31 Mt of KCl (or 17 Mt of K 2 O equivalent) per year could be realized over the next decade if several current (as of 2019) exploration and development projects reach production status. Stratabound potash-bearing salt of the Prairie Evaporite presently underlies a total area of about 188,000 km 2 and has a total volume of about 2,690 cubic kilometers (km 3 ). Post-depositional solution processes considerably modified the mineralogy and presence of the potash-bearing salt. These changes had a profound effect on the volume and grade of potash resources that remained in the Prairie Evaporite and are a major consideration of exploration and mining operations as well as in this assessment of undiscovered potash resources. This USGS assessment includes the locations and possible amounts of undiscovered potash resources in the Prairie Evaporite. Volumes for each stratigraphic member were computed using member thicknesses and areal extent modified by actual, estimated geologic loss owing to salt dissolution and extraction ratios, as well as estimated distribution of carnallite and sylvite. Both sylvite- and carnallite-bearing salts were assessed for potash in this study. The assessment uses modern published grade and tonnage data. The amount of undiscovered potash is estimated by using Monte Carlo simulations to combine volume estimates of the potash-bearing members with probability distributions for average grade and bulk density. Mean potash grades (expressed as percentage of K 2 O equivalent) calculated using drill core analyses are 17.76 for the Patience Lake Member, 15.98 for the Belle Plaine Member, 10.66 for the White Bear Member, and 15.30 for the Esterhazy Member. Geologic losses reported as extraction ratios during mining may range from 27.5 to 41.6 percent and are dependent on mining method and local geologic conditions. The assessment determined that mean estimated undiscovered K 2 O equivalent resources for the Patience Lake, Belle Plaine, White Bear, and Esterhazy Members are 340, 220, 34, and 190 Bt, respectively, and estimated a total mean of 790 Bt for the entire Prairie Evaporite above a depth of 3 km. The total mineralized rock tonnage is estimated to be about 5,000 Bt. Most of the assessed potash is located within Saskatchewan with lesser amounts in Alberta and Manitoba as well as Montana and North Dakota within the United States. Although carnallite is mined for potash in Europe, it has historically been avoided in mining plans for potash-producing companies in Saskatchewan because of mining, processing, and grade considerations. Carnallite-rich salt is locally present in concentrations and volumes that could be a significant resource of magnesium chloride (MgCl 2 ) obtained as a byproduct of processing the carnallite for potash. Previously estimated reserves (not NI 43-101 compliant) of mineralized material from 1955 to 2019 are 695 Mt at 22.1 percent MgCl 2 . The total amount of K 2 O equivalent as carnallite was estimated during this USGS assessment to be about 120 Bt (or 180 Bt KCl). With uncertainties in defining the areal extent of carnallite in each of the potash-bearing members, the amount of MgCl 2 as carnallite in the Elk Point Basin could be approximately 180 Bt.

Elk Point basin↗
USGS1980

Computer-composite mapping for geologists

A computer program for overlaying maps has been tested and evaluated as a means for producing geologic derivative maps. Four maps of the Sugar House Quadrangle, Utah, were combined, using the Multi-Scale Data Analysis and Mapping Program, in a single composite map that shows the relative stability of the land surface during earthquakes. Computer-composite mapping can provide geologists with a powerful analytical tool and a flexible graphic display technique. Digitized map units can be shown singly, grouped with different units from the same map, or combined with units from other source maps to produce composite maps. The mapping program permits the user to assign various values to the map units and to specify symbology for the final map. Because of its flexible storage, easy manipulation, and capabilities of graphic output, the composite-mapping technique can readily be applied to mapping projects in sedimentary and crystalline terranes, as well as to maps showing mineral resource potential. ?? 1980 Springer-Verlag New York Inc.

Environmental Geology↗
USGS2015

csa2sac—A program for computing discharge from continuous slope-area stage data

Introduction Continuous Slope-Area (CSA) gages were developed by the Arizona Water Science Center to enable the estimation of hydrographs when direct measurements of discharge cannot be made (Smith and others, 2010). CSA gages extend standard U.S. Geological Survey (USGS) methods for determining peak discharges to mid and high flows over a hydrograph computed at regular intervals with indirect measurement methods (Benson and Dalrymple, 1967; Dalrymple and Benson, 1967). CSA gages combine continuous stage records at two or more (typically three or four) cross sections with crosssection surveys and estimates of channel roughness to compute discharge over a range of flows. With standard indirect methods of determining peak discharge, water-surface elevation in the study reach at the peak flow is estimated from surveys of debris associated with the peak-flow water line. With CSA gages, stages are continuously measured at the cross sections, at regular and synchronized intervals (typically 5 minutes) over a flow event, and discharge can be calculated at each interval. Calculation of discharge using indirect methods has been automated with the slope-area computation (SAC) program (Fulford, 1994). SAC is a widely used program within the USGS; it is easily run and displays output in a clear and convenient format, which includes flags that alert the user to shortcomings in the calculation. Use of SAC has been facilitated by SACGUI (Bradley, 2012; SACGUI uses a version of SAC called SAC7), a user interface that directly reads and displays survey data, allows for specification of water-surface slope and channel roughness, writes the input file for SAC7, runs SAC7, and displays SAC7 output. csa2sac is a program (appendix 1) that repeatedly runs SAC7 using stage data and a SAC7 input template file to compute the discharge at CSA gages. It is written in the C programming language, and is compatible with 64-bit Windows operating systems. The program reads a SAC7 input file and a file containing stage-data time series. It writes a new version of the SAC7 input file with the stage data for one time step, runs SAC7, then extracts computed discharges from the SAC7 output file and collates the discharges and stages to a separate file. It repeats these steps for each time interval in the stage file to produce a discharge time series from the stage data. csa2sac has been tested with two, three, four, and six cross sections and found to operate successfully. By running SAC7, csa2sac maintains consistency and comparability of both discharges calculated from CSA gages and of standard USGS methods for computing discharges indirectly. Brown and Metcalfe (2014) have made available alternative software for producing CSA discharges. In addition to csa2sac, the SAC7 program is required. It is the same as the original SAC program, except that it is compiled for 64-bit Windows operating systems and has a slightly different command line input. It is available online (http://water.usgs.gov/software/SAC/) as part of the SACGUI installation program. The program name, “SAC7.exe,” is coded into csa2sac, and must not be changed.

Open-File Report↗
USGS2002

Streamflow and nutrient data for the Yazoo River below Steele Bayou near Long Lake, Mississippi, 1996-2000

Increased nutrient loading to the Gulf of Mexico from off-continent flux has been identified as contributing to the increase in the areal extent of the low dissolved-oxygen zone that develops annually off the Louisiana and Texas coast. The proximity of the Yazoo River Basin in northwestern Mississippi to the Gulf of Mexico, and the intensive agricultural activities in the basin have led to speculation that the Yazoo River Basin contributes a disproportionate amount of nitrogen and phosphorus to the Mississippi River and ultimately to the Gulf of Mexico. An empirical measurement of the flux of nitrogen and phosphorus from the Yazoo Basin has not been possible due to the hydrology of the lower Yazoo River Basin. Streamflow for the Yazoo River below Steele Bayou is affected by backwater from the Mississippi River. Flow at the gage is non-uniform and varying, with bi-directional and reverse flows possible. Streamflow was computed by using remote sensing and acoustic and conventional discharge and velocity measurement techniques. Streamflow from the Yazoo River for the 1996-2000 period accounted for 2.8 percent of the flow of the Mississippi River for the same period. Water samples from the Yazoo River were collected from February 1996 through December 2000 and were analyzed for total nitrogen, nitrate, total phosphorus, and orthophosphorus as part of the U.S. Geological Survey National Water-Quality Assessment Program. These data were used to compute annual loads of nitrogen and phosphorus discharged from the Yazoo River for the period 1996-2000. Annual loads of nitrogen and phosphorus were calculated by two methods. The first method used multivariate regression and the second method multiplied the mean annual concentration by the total annual flow. Load estimates based on the product of the mean annual concentration and the total annual flow were within the 95 percent confidence interval for the load calculated by multivariate regression in 10 of 20 cases. The Yazoo River loads, compared to average annual loads in the Mississippi River, indicated that the Yazoo River was contributing 1.4 percent of the total nitrogen load, 0.7 percent of the nitrate load, 3.4 percent of the total phosphorus load, and 1.6 percent of the orthophosphorus load during 1996 - 2000. The total nitrogen, nitrate, and orthophosphorus loads in the Yazoo River Basin were less than expected, whereas the total phosphorus load was slightly higher than expected based on discharge.

Water-Resources Investigations Report↗
USGS2003

MODFLOW-2000, the U.S. Geological Survey modular ground-water model -- Documentation of MOD-PREDICT for predictions, prediction sensitivity analysis, and evaluation of uncertainty

This document describes the MOD-PREDICT program, which helps evaluate userdefined sets of observations, prior information, and predictions, using the ground-water model MODFLOW-2000. MOD-PREDICT takes advantage of the existing Observation and Sensitivity Processes (Hill and others, 2000) by initiating runs of MODFLOW-2000 and using the output files produced. The names and formats of the MODFLOW-2000 input files are unchanged, such that full backward compatibility is maintained. A new name file and input files are required for MOD-PREDICT. The performance of MOD-PREDICT has been tested in a variety of applications. Future applications, however, might reveal errors that were not detected in the test simulations. Users are requested to notify the U.S. Geological Survey of any errors found in this document or the computer program using the email address available at the web address below. Updates might occasionally be made to this document, to the MOD-PREDICT program, and to MODFLOW- 2000. Users can check for updates on the Internet at URL http://water.usgs.gov/software/ground water.html/.

Open-File Report↗
USGS2014

Response of Global Navigation Satellite System receivers to known shaking between 0.2 and 20 Hertz

Over the past decade, several technological advances have allowed Global Navigation Satellite Systems (GNSS) receivers to have the capability to record displacements at high frequencies, with sampling rates approaching 100 samples per second (sps). In addition, communication and computer hardware and software have allowed various institutions, including the U.S. Geological Survey (USGS), to retrieve, process, and display position changes recorded by a network of GNSS sites with small, less than 1-s delays between the time that the GNSS receiver records signals from a constellation of satellites and the time that the position is estimated (a method known as “real-time”). These improvements in hardware and software have allowed the USGS to process GNSS (or a subset of the GNSS, the Global Positioning System, GPS) data in real-time at 1 sps with the goal of determining displacements from earthquakes and volcanoes in real-time. However, the current set of GNSS equipment can record at rates of 100 sps, which allows the possibility of using this equipment to record earthquake displacements over the full range of frequencies that typically are recorded by acceleration and velocity transducers. The advantage of using GNSS to record earthquakes is that the displacement, rather than acceleration or velocity, is recorded, and for large earthquakes, the GNSS sensor stays on scale and will not distort the observations due to clipping of the signal at its highest amplitude. The direct observation of displacement is advantageous in estimating the size and spatial extent of the earthquake rupture. Otherwise, when using velocity or acceleration sensors, the displacements are determined by numerical integration of the observations, which can introduce significant uncertainty in the estimated displacements. However, GNSS technology can, at best, resolve displacements of a few millimeters, and for most earthquakes, their displacements are less than 1 mm. Consequently, to be useful, GNSS data are only relevant for the large earthquakes with magnitudes (M) exceeding M5.5 at best. With the capability to record GNSS data at high-rate, at sampling rates typical for seismological applications, experiments are needed to quantify the response of GNSS to shaking from earthquakes. There have been a few studies that examine the response of GNSS to strong shaking. One of the first was Elosegui and others (2006), where they simulated surface waves from a distant earthquake and mechanically applied the shaking to a GPS antenna. They processed the 1 sps observations and compared the estimated displacements with the simulated displacements. They determined that the GPS could accurately track the simulated surface wave whose primary frequency spans from 0.01 to 0.1 Hertz (Hz), which spanned the frequency band of the simulation. To test GNSS equipment due to shaking from a large earthquake in the near-field, Wang and others (2012) used a mechanical simulator or shake table with 6 degrees of freedom and studied two different inputs to the simulator—(1) the accelerometer record from one station that was located near the 2010 M8.8 Maule, Chile earthquake, and (2) a 2-Hz sinusoid. Wang and others (2012) analyzed the 2-Hz data with spectral analysis and determined that the displacements observed by the GPS included higher harmonics along with the 2-Hz signal. In addition, the background spectral amplitude was greater during periods of 2-Hz shaking than when at rest. With the simulated M 8.8 earthquake, Wang and others (2012) observed decreased signal to noise for L1 and L2 carrier frequencies of the GPS signal, at times corresponding to high acceleration and jerk (first derivative of acceleration). One of the principal limitations of these experiments was that the displacements of the shake table itself could not be measured independently. Although with the 2-Hz sinusoidal measurements, the input displacements were purely translational, Wang and others (2012) analysis of the data showed that the shake table also included rotational motions which affect horizontal inertial sensors like accelerometers and seismometers at first order. More recently, Ebinuma and Kato (2012) used a GPS simulator to electronically test several GNSS receivers and obtain the receiver characteristics at three frequencies: 1, 2, and 5 Hz. The results showed that the amplitude of 5-Hz displacements recorded by the GPS was, depending on the receiver model, between 30 and 125 percent more than the displacement input to the simulator. At low frequencies, the GPS displacement was nearly equal to the input displacement. In addition, Ebinuma and Kato (2012) examined how each receiver model amplified an earthquake displacement record in the 2–8 Hz band. The simulated earthquake was the 2008 moment magnitude (Mw) 6.8 Iwate-Miyagi earthquake where, for the simulated record, acceleration peaked at 1 G. The study discussed here builds on the tests by Ebinuma and Kato (2012), but rather than using electronic simulation, the tests are setup outdoors and closer to actual field installations of GNSS equipment. We used a one-dimensional shake table capable of 400 mm of displacement and high acceleration; the shake table also is constrained by a precision linear slider to have very low tilt that would affect inertial sensors. In addition, the stage position can be accurately monitored independent of the GNSS hardware and, importantly, provides a reference to compare with the estimated displacements from the GNSS data. Our tests spanned a greater frequency range from 0.2 to 20 Hz and we used equipment from three different manufacturers covering five different combinations of receivers and antennas. In addition, we have been able to simulate the frequency response of the GNSS equipment using a simple, causal filter. The quality of the filter was tested using additional test data where a step function in displacement was applied to the shake table. The observed displacements from the GNSS data show an overshoot in displacement at the time of the step or transition of the stage. That overshoot was accurately predicted using the filter design derived from our sinusoidal displacement tests. Similar to Wang and others (2012), we also examined the GPS displacement records using standard spectral techniques. However, we extended their work by evaluating several models of GNSS receivers using a variety of input frequencies. Because our shake table was limited on acceleration and displacement, we did not attempt to duplicate the high shaking associated with high magnitude earthquakes. However, because our shake table could measure the table displacement, we could directly compare the measured GPS displacements with the true displacements.

Open-File Report↗
USGS1999

HYPOELLIPSE; a computer program for determining local earthquake hypocentral parameters, magnitude, and first-motion pattern

This report provides Fortran source code and program manuals for HYPOELLIPSE, a computer program for determining hypocenters and magnitudes of near regional earthquakes and the ellipsoids that enclose the 68-percent confidence volumes of the computed hypocenters. HYPOELLIPSE was developed to meet the needs of U.S. Geological Survey (USGS) scientists studying crustal and sub-crustal earthquakes recorded by a sparse regional seismograph network. The program was extended to locate hypocenters of volcanic earthquakes recorded by seismographs distributed on and around the volcanic edifice, at elevations above and below the hypocenter. HYPOELLIPSE was used to locate events recorded by the USGS southern Alaska seismograph network from October 1971 to the early 1990s. Both UNIX and PC/DOS versions of the source code of the program are provided along with sample runs.

Open-File Report↗
USGS2007

Evaluation of abutment scour prediction equations with field data

The U.S. Geological Survey, in cooperation with FHWA, compared predicted abutment scour depths, computed with selected predictive equations, with field observations collected at 144 bridges in South Carolina and at eight bridges from the National Bridge Scour Database. Predictive equations published in the 4th edition of Evaluating Scour at Bridges (Hydraulic Engineering Circular 18) were used in this comparison, including the original Froehlich, the modified Froehlich, the Sturm, the Maryland, and the HIRE equations. The comparisons showed that most equations tended to provide conservative estimates of scour that at times were excessive (as large as 158 ft). Equations also produced underpredictions of scour, but with less frequency. Although the equations provide an important resource for evaluating abutment scour at bridges, the results of this investigation show the importance of using engineering judgment in conjunction with these equations.

Transportation Research Record↗
USGS2002

Strategic plan for the Coordinated Intermountain Restoration Project

In 1982, the Bureau of Land Management’s (BLM) Idaho State Office began the Intermountain Greenstripping and Rehabilitation Research Project (IGRRP), or the “Greenstripping Program,” to investigate plant materials and technologies that can reduce wildfire incidence and improve rehabilitation practices. Rehabilitation is normally applied as a reactive process to wildfires, yet land managers in the Great Basin wish to become proactive by replacing fire-prone invasive annual grasses with native plants. The Coordinated Intermountain Restoration Project (CIRP) evolved from the Greenstripping Program to conduct research studies and provide technical assistance on restoration of native ecosystems on rangelands that are infested with invasive annual grasses or other invasive or noxious weeds. To accomplish this objective, the CIRP will promote the understanding of ecosystem disturbance dynamics as well as evaluate plant materials, site preparation techniques, weed control methods, seeding equipment, management methods, and monitoring techniques for restoration projects. The CIRP will not address the restoration of forested or woodland (juniper [Juniperus]) ecosystems. It will include a component on fuel management to reduce the impacts of wildfires on semiarid rangeland ecosystems where exotic annual grasses provide the fuel. The people who will benefit directly from this research include land managers and users of public and private lands in the northern Great Basin, the Columbia Plateau, and the Snake River Plain. The CIRP will provide an integration framework for a multidisciplinary approach to research with numerous opportunities for input and collaboration. The U.S. Geological Survey will initially dedicate approximately \$1 million over 5 years (about \$200,000 per year) to jump-start this effort. U.S. Geological Survey funds will establish a science advisory board to oversee the project. This board will contain members of Federal research and management agencies within the region. U.S. Geological Survey funds will support (1) continued development of VegSpec, a computer program that is a restoration expert system, (2) research to examine changes in ecosystem processes when native plant-dominated communities shift to communities dominated by exotic annual grasses, and (3) research to address mechanisms for establishing native plants in locations dominated by exotic annual grasses. Through these initial funds, USGS hopes to leverage additional research with other agencies (e.g., BLM’s Great Basin Restoration Initiative or the Native Plant Materials Development Project, which is an interagency program to supply and manage native plant materials for restoration and rehabilitation on Federal lands) or funding organizations (e.g., the U.S. Department of Agriculture’s [USDA] National Research Initiative Competitive Grants Program, or the USDA’s and U.S. Department of the Interior’s [USDOI] Joint Fire Science Program), and to obtain additional research partners (e.g., university or Federal scientists) willing to expand this effort to address all aspects of this strategic plan.

Information and Technology Report↗
USGS2023

Evaluating oil and gas industry two-dimensional multichannel seismic data for use in near-surface assessment of geologic framework and potential marine minerals resources

Marine seismic reflection data acquired across the Gulf of Mexico during oil and gas exploration are available to the public through an online database archive. The data are archived as two-dimensional multichannel seismic data in two digital formats. The formats include image files in portable document format (PDF), and binary files in industry standard Society for Exploration Geophysicists revision Y (SEG-Y) format. Also included in the database are navigation files and acquisition information associated with the collection of the data. This study examines the data acquired within two geographic areas in the northern Gulf of Mexico. Although the seismic reflection data are acquired for oil and gas exploration many kilometers below the seafloor, this study focuses on the feasibility of using the data for near-surface geologic and seafloor morphologic studies (<100 meters below the seafloor). The report outlines the methodologies used to recover and process the data, including computer processing steps to convert the PDF imagery into SEG-Y format. The report includes two-dimensional profiles of the data to demonstrate the efficacy of the data in near-surface geologic studies. The study found that, for the two areas of interest, the seafloor reflectors in most of the available data are not resolvable. Although the data are readily available and computer processing can adequately image the uppermost reflectors of the seismic profiles, the resolution of the data in most cases are not suitable for near-surface geologic evaluations.

Alabama, Florida, Louisiana, Mississippi, Texas↗
USGS1997

An interactive program for computer-aided map design, display, and query: EMAPKGS2

EMAPKGS2 is a user-friendly, PC-based electronic mapping tool for use in hydrogeologic exploration and appraisal. EMAPKGS2 allows the analyst to construct maps interactively from data stored in a relational database, perform point-oriented spatial queries such as locating all wells within a specified radius, perform geographic overlays, and export the data to other programs for further analysis. EMAPKGS2 runs under Microsoft?? Windows??? 3.1 and compatible operating systems. EMAPKGS2 is a public domain program available from the Kansas Geological Survey. EMAPKGS2 is the centerpiece of WHEAT, the Windows-based Hydrogeologic Exploration and Appraisal Toolkit, a suite of user-friendly Microsoft?? Windows??? programs for natural resource exploration and management. The principal goals in development of WHEAT have been ease of use, hardware independence, low cost, and end-user extensibility. WHEAT'S native data format is a Microsoft?? Access?? database. WHEAT stores a feature's geographic coordinates as attributes so they can be accessed easily by the user. The WHEAT programs are designed to be used in conjunction with other Microsoft?? Windows??? software to allow the natural resource scientist to perform work easily and effectively. WHEAT and EMAPKGS have been used at several of Kansas' Groundwater Management Districts and the Kansas Geological Survey on groundwater management operations, groundwater modeling projects, and geologic exploration projects. ?? 1997 Elsevier Science Ltd.

Computers & Geosciences↗
USGS2007

The National Streamflow Statistics Program: A Computer Program for Estimating Streamflow Statistics for Ungaged Sites

The National Streamflow Statistics (NSS) Program is a computer program that should be useful to engineers, hydrologists, and others for planning, management, and design applications. NSS compiles all current U.S. Geological Survey (USGS) regional regression equations for estimating streamflow statistics at ungaged sites in an easy-to-use interface that operates on computers with Microsoft Windows operating systems. NSS expands on the functionality of the USGS National Flood Frequency Program, and replaces it. The regression equations included in NSS are used to transfer streamflow statistics from gaged to ungaged sites through the use of watershed and climatic characteristics as explanatory or predictor variables. Generally, the equations were developed on a statewide or metropolitan-area basis as part of cooperative study programs. Equations are available for estimating rural and urban flood-frequency statistics, such as the 1 00-year flood, for every state, for Puerto Rico, and for the island of Tutuila, American Samoa. Equations are available for estimating other statistics, such as the mean annual flow, monthly mean flows, flow-duration percentiles, and low-flow frequencies (such as the 7-day, 0-year low flow) for less than half of the states. All equations available for estimating streamflow statistics other than flood-frequency statistics assume rural (non-regulated, non-urbanized) conditions. The NSS output provides indicators of the accuracy of the estimated streamflow statistics. The indicators may include any combination of the standard error of estimate, the standard error of prediction, the equivalent years of record, or 90 percent prediction intervals, depending on what was provided by the authors of the equations. The program includes several other features that can be used only for flood-frequency estimation. These include the ability to generate flood-frequency plots, and plots of typical flood hydrographs for selected recurrence intervals, estimates of the probable maximum flood, extrapolation of the 500-year flood when an equation for estimating it is not available, and weighting techniques to improve flood-frequency estimates for gaging stations and ungaged sites on gaged streams. This report describes the regionalization techniques used to develop the equations in NSS and provides guidance on the applicability and limitations of the techniques. The report also includes a users manual and a summary of equations available for estimating basin lagtime, which is needed by the program to generate flood hydrographs. The NSS software and accompanying database, and the documentation for the regression equations included in NSS, are available on the Web at http://water.usgs.gov/software/.

Techniques and Methods↗
USGS1973

Stress changes ahead of an advancing tunnel

Instrumentation placed ahead of three model tunnels in the laboratory and ahead of a crosscut driven in a metamorphic rock mass detected stress changes several tunnel diameters ahead of the tunnel face. Stress changes were detected 4 diameters ahead of a model tunnel drilled into nearly elastic acrylic, 2·50 diameters ahead of a model tunnel drilled into concrete, and 2 diameters ahead of a model tunnel drilled into Silver Plume Granite. Stress changes were detected 7·50 diameters ahead of a crosscut driven in jointed, closely foliated gneisses and gneissic granites in an experimental mine at Idaho Springs, Colorado. These results contrast markedly with a theoretical elastic estimate of the onset of detectable stress changes at 1 tunnel diameter ahead of the tunnel face. A small compressive stress concentration was detected 2 diameters ahead of the model tunnel in acrylic, 1.25 diameters ahead of the model tunnel in concrete, and 1 diameter ahead of the model tunnel in granite. A similar stress peak was detected about 6 diameters ahead of the crosscut. No such stress peak is predicted from elastic theory. The 3-dimensional in situ stress determined in the field demonstrate that geologic structure controls stress orientations in the metamorphic rock mass. Two of the computed principal stresses are parallel to the foliation and the other principal stress is normal to it. The principal stress orientations vary approximately as the foliation attitude varies. The average horizontal stress components and the average vertical stress component are three times and twice as large, respectively, as those predicted from the overburden load. An understanding of the measured stress field appears to require the application of either tectonic or residual stress components, or both. Laboratory studies indicate the presence of proportionately large residual stresses. Mining may have triggered the release of strain energy, which is controlled by geologic structure.

Colorado↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS2025

Computing discharge using the entropy-based probability concept

This report describes the techniques and methods for computing the mean-channel velocity and discharge using the entropy-based probability concept (probability concept). The method is an alternative to or augments standard streamgaging methods adopted by the U.S. Geological Survey (USGS). Although sensor technology for measuring the mean velocity and discharge has advanced, standard streamgaging and computational methods have remained relatively unchanged since the USGS established its first streamgage at the Rio Grande at Embudo, New Mexico in 1889. Standard streamgaging methods rely on integrating velocities and depths measured at multiple verticals at a channel cross section (standard cross section) to compute a discharge. The probability concept computes discharge at a single vertical (y-axis) using the ratio of the mean-channel velocity (mean velocity) and maximum velocity, the measured maximum velocity, and the area as a function of stage at the standard cross section. Proper siting and operation and maintenance are required. If siting is conducted appropriately, the probability concept parameters and the y-axis stationing will be similar for different streamflow conditions. The timing of operation and maintenance visits should be based on hydrologic and meteorologic occurrences and seasonality and should capture low, medium, high, and opportunistic streamflow conditions. Advantages of the probability concept are the capacity to (1) compute discharge time series immediately after streamgage siting, (2) compute discharge for complex streamflow conditions that cannot be quantified by stage-discharge methods, (3) augment time-series data where gaps exist, and (4) integrate with surface velocity sensors such as Doppler velocity radars and cameras, which are not subject to damage caused by ice, debris, and flood flows. Potential sources of bias in discharge derived from the probability concept include (1) rain, (2) wind, and (3) geomorphologic and hydraulic instabilities. Recommendations to address these biases are provided. This report guides users through the steps to parameterize the probability concept, process field data, and compute the mean velocity and discharge using the probability concept.

Techniques and Methods↗
USGS2006

Computer program for the Kendall family of trend tests

The Seasonal Kendall (SK) test for trend was developed by the U.S. Geological Survey and has become the most frequently used test for trend in the environmental sciences. Recently the test was modified to form the Regional Kendall (RK) test for trend. In this form, trends at numerous locations within a region are tested to determine whether the direction of trend is consistent across the entire region. Computer code developed at the USGS in the 1980s to perform the SK test is no longer widely available. Other versions written by other scientists may or may not be easily available, and may require commercial software in order to be run. These other versions do not explicitly compute the RK test. Therefore, the original code for computing the SK test has been repackaged into a program that runs under the Windows operating system. This program may be used to verify that other implementations of the test give the same results as the original. The program also provides a means for computing the RK test and the simpler Mann-Kendall test for trend.

Scientific Investigations Report↗
USGS2004

NADM Conceptual Model 1.0 -- A Conceptual Model for Geologic Map Information

Executive Summary -- The NADM Data Model Design Team was established in 1999 by the North American Geologic Map Data Model Steering Committee (NADMSC) with the purpose of drafting a geologic map data model for consideration as a standard for developing interoperable geologic map-centered databases by state, provincial, and federal geological surveys. The model is designed to be a technology-neutral conceptual model that can form the basis for a web-based interchange format using evolving information technology (e.g., XML, RDF, OWL), and guide implementation of geoscience databases in a common conceptual framework. The intended purpose is to allow geologic information sharing between geologic map data providers and users, independent of local information system implementation. The model emphasizes geoscience concepts and relationships related to information presented on geologic maps. Design has been guided by an informal requirements analysis, documentation of existing databases, technology developments, and other standardization efforts in the geoscience and computer-science communities. A key aspect of the model is the notion that representation of the conceptual framework (ontology) that underlies geologic map data must be part of the model, because this framework changes with time and understanding, and varies between information providers. The top level of the model distinguishes geologic concepts, geologic representation concepts, and metadata. The geologic representation part of the model provides a framework for representing the ontology that underlies geologic map data through a controlled vocabulary, and for establishing the relationships between this vocabulary and a geologic map visualization or portrayal. Top-level geologic classes in the model are Earth material (substance), geologic unit (parts of the Earth), geologic age, geologic structure, fossil, geologic process, geologic relation, and geologic event.

Open-File Report↗