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 919 records · Page 51Linked to original sources
USGSDate not supplied

Relations between geological structure and aeromagnetic anomalies in central Nevada

An aeromagnetic survey of a 6000-km 2 area in central Nevada has been complied by the U.S. Geological Survey. An interpretation of the aeromagnetic data is developed to establish relations between near-surface geological features and structure extending to intermediate depths in the crust. A three-dimensional distribution of magnetized material is described which would reproduce approximately the observed aeromagnetic anomalies. Anomaly sources extend from the surface to depths exceeding 15 km and require values of magnetic moment per unit volume averaging more than 2000 × 10 −6 (cgs). Magnetic anomaly lineations were analyzed by two methods. The first involves digital fan filtering and subsequent autocovariance computations. This method is sensitive to amplitude of linear anomalies, and indicated a dominant north-northwest trend with approximately 80 percent of the anomaly trend power in an azimuth range of northwest to northeast. The second method involves preparation of rosettes from measured lengths and azimuths of lineations indicated by noses and flexures on the contoured map. This method is not sensitive to anomaly amplitude. The dominant N. 10°E. trend determined in this way is indicative of the effect of near-surface geological structure on the magnetic field. The dominant north-northwest-trending zone of magnetic anomalies is located along an alinement of windows through the Roberts Mountains thrust fault. This suggests that elongate doming of the Paleozoic sedimentary rocks associated with the thrust may be related to emplacement of magnetic anomaly source material, and implies that this has been a zone of recurring tectonic activity since late Paleozoic time. Analysis of stresses possibly associated with forceful injection of magnetic anomaly sources and possible shear deformation along this north-northwest-trending zone suggests that these factors may have had a secondary influence on azimuths of high-angle faults and dikes.

Nevada↗
USGS2018

A surrogate regression approach for computing continuous loads for the tributary nutrient and sediment monitoring program on the Great Lakes

Water quality (WQ) in many Great Lake tributaries has been degraded (increased nutrient and sediment concentrations) due to changes in their watersheds, resulting in downstream eutrophication. As part of the Great Lakes Water Quality Agreement, specific goals were established for loading of specific constituents (e.g., phosphorus). In 2010, the Great Lakes Restoration Initiative was launched to identify problem areas, accelerate restoration efforts, and track their progress. In 2011, the U.S. Geological Survey established a monitoring program on 30 tributaries to the lakes, representing ~ 46% of the U.S. draining area and the spectrum of land uses. Discrete measurements of nutrients and suspended sediment, and continuous measurements of flow and WQ surrogates (turbidity, temperature, specific conductance, pH, and dissolved oxygen) are being collected in these tributaries to document their WQ and estimate continuous (5-min) loading. To estimate loadings, two regression models were developed for each constituent for each site: one using continuous flow and a seasonality factor; and one using flow, seasonality, and continuous surrogates. Variables included in the final models for each constituent were chosen from the explanatory variables that worked “best” for all sites. In computing loads, when continuous surrogate data were unavailable for short periods, loads were computed using the flow and seasonality models. Prediction intervals for all loads were calculated using results from both models. These results provide a better understanding of short-term variability and long-term changes in loading affecting the environmental health of the Great Lakes than traditional regression techniques that employ only flow and seasonality parameters.

Great Lake↗
USGS2009

Low-flow frequency and flow duration of selected South Carolina streams in the Pee Dee River basin through March 2007

Part of the mission of the South Carolina Department of Health and Environmental Control and the South Carolina Department of Natural Resources is to protect and preserve South Carolina's water resources. Doing so requires an ongoing understanding of streamflow characteristics of the rivers and streams in South Carolina. A particular need is information concerning the low-flow characteristics of streams; this information is especially important for effectively managing the State's water resources during critical flow periods such as the severe drought that occurred between 1998 and 2002 and the most recent drought that occurred between 2006 and 2009. In 2008, the U.S. Geological Survey, in cooperation with the South Carolina Department of Health and Environmental Control, initiated a study to update low-flow statistics at continuous-record streamgaging stations operated by the U.S. Geological Survey in South Carolina. Under this agreement, the low-flow characteristics at continuous-record streamgaging stations will be updated in a systematic manner during the monitoring and assessment of the eight major basins in South Carolina as defined and grouped according to the South Carolina Department of Health and Environmental Control's Watershed Water Quality Management Strategy. Depending on the length of record available at the continuous-record streamgaging stations, low-flow frequency characteristics are estimated for annual minimum 1-, 3-, 7-, 14-, 30-, 60-, and 90-day average flows with recurrence intervals of 2, 5, 10, 20, 30, and 50 years. Low-flow statistics are presented for 18 streamgaging stations in the Pee Dee River basin. In addition, daily flow durations for the 5-, 10-, 25-, 50-, 75-, 90-, and 95-percent probability of exceedance also are presented for the stations. The low-flow characteristics were computed from records available through March 31, 2007. The last systematic update of low-flow characteristics in South Carolina occurred more than 20 years ago and included data through March 1987. Of the 17 streamgaging stations included in this study, 15 had low-flow characteristics that were published in previous U.S. Geological Survey reports. A comparison of the low-flow characteristic for the minimum average flow for a 7-consecutive-day period with a 10-year recurrence interval from this study with the most recently published values indicated that 10 of the 15 streamgaging stations had values that were within ±25 percent of each other. Nine of the 15 streamgaging stations had negative percentage differences indicating the low-flow statistic had decreased since the previous study, 4 streamgaging stations had positive percent differences indicating that the low-flow statistic had increased since the previous study, and 2 streamgaging stations had a zero percent difference indicating no change since the previous study. The low-flow characteristics are influenced by length of record, hydrologic regime under which the record was collected, techniques used to do the analysis, and other changes that may have occurred in the watershed.

North Carolina, South Carolina↗
USGS2012

Archive of digital chirp subbottom profile data collected during USGS cruise 10BIM04 offshore Cat Island, Mississippi, September 2010

In September of 2010, the U.S. Geological Survey (USGS), in cooperation with the U.S. Army Corps of Engineers (USACE), conducted a geophysical survey to investigate the geologic controls on barrier island framework of Cat Island, Miss., as part of a broader USGS study on Barrier Island Mapping (BIM). These surveys were funded through the Mississippi Coastal Improvements Program (MsCIP) and the Northern Gulf of Mexico (NGOM) Ecosystem Change and Hazard Susceptibility Project as part of the Holocene Coastal Evolution of the Mississippi-Alabama Region Subtask. This report serves as an archive of unprocessed digital chirp subbottom data, trackline maps, navigation files, GIS files, Field Activity Collection System (FACS) logs, and formal FGDC metadata. Gained (showing a relative increase in signal amplitude) digital images of the seismic profiles are also provided. Refer to the Acronyms page for expansions of acronyms and abbreviations used in this report. The USGS Saint Petersburg Coastal and Marine Science Center (SPCMSC) assigns a unique identifier to each cruise or field activity. For example, 10BIM04 tells us the data were collected in 2010 during the fourth field activity for that project in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the field activity identification (ID). All chirp systems use a signal of continuously varying frequency; the EdgeTech SB-512i system used during this survey produces high-resolution, shallow-penetration (typically less than 50 milliseconds (ms)) profile images of sub-seafloor stratigraphy. The towfish contains a transducer that transmits and receives acoustic energy; it was housed within a float system (built at the SPCMSC), which allows the towfish to be towed at a constant depth of 1.07 meters (m) below the sea surface. As transmitted acoustic energy intersects density boundaries, such as the seafloor or sub-surface sediment layers, some energy is reflected back toward the transducer, received, and recorded by a Personal Computer (PC)-based seismic acquisition system. This process is repeated at regular time intervals (for example, 0.125 seconds (s)), and returned energy is recorded for a specific duration (for example, 50 ms). In this way, a two-dimensional (2-D) vertical image of the shallow geologic structure beneath the ship track is produced. Figure 1 displays the acquisition geometry. Refer to table 1 for a summary of acquisition parameters and table 2 for trackline statistics. The archived trace data are in standard Society of Exploration Geophysicists (SEG) SEG Y rev. 0 format (Barry and others, 1975); the first 3,200 bytes of the card image header are in American Standard Code for Information Interchange (ASCII) format instead of Extended Binary Coded Decimal Interchange Code (EBCDIC) format. The SEG Y files may be downloaded and processed with commercial or public domain software such as Seismic Unix (SU) (Cohen and Stockwell, 2010). See the How To Download SEG Y Data page for download instructions. The printable profiles provided here are GIF images that were processed and gained using SU software, and they can be viewed from the Profiles page or from links located on the trackline maps; refer to the Software page for links to example SU processing scripts. The SEG Y files are available on the DVD version of this report or on the Web, downloadable via the USGS Coastal and Marine Geoscience Data System (http://cmgds.marine.usgs.gov). The data are also available for viewing using GeoMapApp (http://www.geomapapp.org) and Virtual Ocean (http://www.virtualocean.org) multi-platform open source software.

Mississippi↗
USGS2021

Continuous turbidity data used to compute constituent concentrations in the South Loup River, Nebraska, 2017–18

The South Loup River in central Nebraska has been impaired by bacteria since at least 2004, which has resulted in the river not meeting its intended use as a recreational waterway. As part of a strategy for reducing the bacterial load in the river, the U.S. Geological Survey, in cooperation with the Lower Loup Natural Resources District, made continuous estimates of Escherichia coli ( E. coli ) and nutrient concentrations during seasonal monitoring at the South Loup River at Saint Michael, Nebraska, during 2017–18. Continuous turbidity data were collected from mid-April through October in 2017 and 2018 and were paired with 35 co-occurring discrete water samples that were analyzed for E. coli , nutrients, and suspended solids. Surrogate models relating the discrete concentrations to the continuous turbidity data were developed using ordinary-least-squares regression and were evaluated for model performance and uncertainty. Although the model assumptions were met for E. coli , the imprecision of the E. coli model was considerably higher than the other constituents, probably because of measurement imprecision and greater sensitivity to environmental factors. Once the models were developed, the turbidity data were used to predict continuous constituent concentrations and corresponding prediction intervals, which were made available online as part of the U.S. Geological Survey National Water Information System database. It is expected that results from these models will provide stakeholders with an understanding of constituent concentrations during the 2017–18 monitoring period and the results will also provide a good reference point for any future comparisons.

Nebraska↗
USGS2014

Surface wave site characterization at 27 locations near Boston, Massachusetts, including 2 strong-motion stations

The geotechnical properties of the soils in and around Boston, Massachusetts, have been extensively studied. This is partly due to the importance of the Boston Blue Clay and the extent of landfill in the Boston area. Although New England is not a region that is typically associated with seismic hazards, there have been several historical earthquakes that have caused significant ground shaking (for example, see Street and Lacroix, 1979; Ebel, 1996; Ebel, 2006). The possibility of strong ground shaking, along with heightened vulnerability from unreinforced masonry buildings, motivates further investigation of seismic hazards throughout New England. Important studies that are pertinent to seismic hazards in New England include source-parameter studies (Somerville and others, 1987; Boore and others, 2010), wave-propagation studies (Frankel, 1991; Viegas and others, 2010), empirical ground-motion prediction equations (GMPE) for computing ground-motion intensity (Tavakoli and Pezeshk, 2005; Atkinson and Boore, 2006), site-response studies (Hayles and others, 2001; Ebel and Kim, 2006), and liquefaction studies (Brankman and Baise, 2008). The shear-wave velocity (VS) profiles collected for this report are pertinent to the GMPE, site response, and liquefaction aspects of seismic hazards in the greater Boston area. Besides the application of these data for the Boston region, the data may be applicable throughout New England, through correlations with geologic units (similar to Ebel and Kim, 2006) or correlations with topographic slope (Wald and Allen, 2007), because few VS measurements are available in stable tectonic regions. Ebel and Hart (2001) used felt earthquake reports to infer amplification patterns throughout the greater Boston region and noted spatial correspondence with the dominant period and amplification factors obtained from ambient noise (horizontal-to-vertical ratios) by Kummer (1998). Britton (2003) compiled geotechnical borings in the area and produced a microzonation map based on generalized velocity profiles, where the amplifications were computed using Shake (Schnable and others, 1972), along with an assumed input ground motion. The velocities were constrained by only a few local measurements associated with the Central Artery/Tunnel project. The additional VS measurements presented in this report provide a number of benefits. First, these measurements provide improved spatial coverage. Second, the larger sample size provides better constraints on the mean and variance of the VS distribution for each layer, which may be paired with a three-dimensional (3D) model of the stratigraphy to generate one-dimensional (1D) profiles for use in a standard site-response analysis (for example, Britton, 2003). Third, the velocity profiles may also be used, along with a 3D model of the stratigraphy, as input into a 3D simulation of the ground motion to investigate the effects of basin-generated surface waves and the potential focusing of seismic waves. This report begins with a short review of the geology of the study area and the field methods that we used to estimate the velocity profiles. The raw data, processed data, and the interpreted VS profiles are given in appendix 1. Photographs and descriptions of the sites are provided in appendix 2.

Massachusetts↗
USGS2020

Conjoint use of hydraulic head and groundwater age data to detect hydrogeologic barriers

Hydraulic head and groundwater age data are effective in building understanding of groundwater systems. Yet their joint role in detecting and characterising low-permeability geological structures, i.e. hydrogeologic barriers such as faults and dykes, has not been widely studied. Here, numerical flow and transport models, using MODFLOW-NWT and MT3D-USGS, were developed with different hydrogeologic barrier configurations in a hypothetical aquifer. Computed hydraulic head and groundwater age distributions were compared to those without a barrier. The conjoint use of these datasets helps in detecting vertically-oriented barriers. Two forms of recharge were compared: (1) applied across the entire aquifer surface (uniform), and (2) applied to the upstream part of the aquifer (upgradient). The hydraulic head distribution is significantly impacted by a barrier that penetrates the aquifer’s full vertical thickness. This barrier also perturbs the groundwater age distribution when upgradient recharge prevails; however, with uniform recharge, groundwater age is not successful in detecting the barrier. When a barrier is buried, such as by younger sediment, hydraulic head data also do not clearly identify the barrier. Groundwater age data could, on the other hand, prove to be useful if sampled at depth-specific intervals. These results are important for the detection and characterisation of hydrogeologic barriers, which may play a significant role in the compartmentalisation of groundwater flow, spring dynamics, and drawdown and recovery associated with groundwater extraction.

Hydrogeology Journal↗
USGS2020

Automated location correction and spot height generation for named summits in the coterminous United States

Spot elevations published on historical U.S. Geological Survey topographic maps were established as needed to enhance information imparted by the quadrangle’s contours. In addition to other features, labels were routinely placed on mountain summits. While some elevations were established through field survey triangulation, many were computed during photogrammetric stereo-compilation. Today, Global Navigation Satellite System (GNSS) receivers have replaced expensive triangulation methods. However, since GNSS measurements require visiting the feature location, a national dataset containing high-accuracy spot elevations has not yet been created. Consequently, modern U.S. Topo maps are devoid of mountain peak or other spot elevations. Still, topographic map users continue to demand the display of spot heights. Therefore, a pilot study was conducted to evaluate the feasibility of automatically generating elevation values at named U.S. summits using available elevation data. The devised method uses an uphill stepping technique to find the most likely highest point in subsequently higher-resolution elevation models. Resulting elevation values are compared to other published sources. Results from 196 summits indicate that values derived from lidar are generally higher, whereas those populated from the one-third arc-second USGS Seamless 3DEP elevation dataset are generally lower. A thorough understanding of these relationships require the evaluation of more points.

International Journal of Digital Earth↗
USGS2001

Dynamic computer model for the metallogenesis and tectonics of the Circum-North Pacific

The digital files on this report consist of a dynamic computer model of the metallogenesis and tectonics of the Circum-North Pacific, and background articles, figures, and maps. The tectonic part of the dynamic computer model is derived from a major analysis of the tectonic evolution of the Circum-North Pacific which is also contained in directory tectevol. The dynamic computer model and associated materials on this CD-ROM are part of a project on the major mineral deposits, metallogenesis, and tectonics of the Russian Far East, Alaska, and the Canadian Cordillera. The project provides critical information on bedrock geology and geophysics, tectonics, major metalliferous mineral resources, metallogenic patterns, and crustal origin and evolution of mineralizing systems for this region. The major scientific goals and benefits of the project are to: (1) provide a comprehensive international data base on the mineral resources of the region that is the first, extensive knowledge available in English; (2) provide major new interpretations of the origin and crustal evolution of mineralizing systems and their host rocks, thereby enabling enhanced, broad-scale tectonic reconstructions and interpretations; and (3) promote trade and scientific and technical exchanges between North America and Eastern Asia.

Open-File Report↗
USGS2005

Variability of differences between two approaches for determining ground-water discharge and pumpage, including effects of time trends, Lower Arkansas River Basin, southeastern Colorado, 1998-2002

In the mid-1990s, the Colorado Division of Water Resources (CDWR) adopted rules governing measurement of tributary ground-water pumpage for the Arkansas River Basin. The rules allowed ground-water pumpage to be determined using one of two approaches?power conversion coefficient (PCC) or totalizing flowmeters (TFM). In addition, the rules allowed a PCC to be applied to the electrical power usage up to 4 years in the future to estimate ground-water pumpage. As a result of concerns about potential errors in applying the PCC approach forward in time, a study was done by the U.S. Geological Survey, in cooperation with CDWR and Colorado Water Conservation Board, to evaluate the variability in differences in pumpage between the two approaches, including the effects of time trends. This report compared measured ground-water pumpage using TFMs to computed ground-water pumpage using PCCs by developing statistical models of relations between explanatory variables, such as site, time, and pumping water level, and dependent variables, which are based on discharge, PCC, and pumpage. When differences in pumpage (diffP) were computed using PCC measurements and power consumption for the same year (1998-2002), the median diffP, depending on the year, ranged from +0.1 to -2.9 percent; the median diffP for the entire period was -1.5 percent. However, when diffP was computed using PCC measurements applied to the next year's power consumption, the median diffP was -0.3 percent; and when PCC measurements were applied 2, 3, or 4 years into the future, median diffPs were +1.8 percent for a 2-year forward lag and +5.3 percent for a 4-year forward lag, indicating that pumpage computed with the PCC approach, as generally applied under the ground-water pumpage measurement rules by CDWR, tended to overestimate pumpage as compared to pumpage using TFMs when PCC measurement was applied to future years of measured power consumption. Analyses were done to better understand the causes of the time trend; an estimate of the overall trend with time (uncorrected for pumping water-level changes) yielded a trend of about 2.2 percent per lag year for diffP. A separate analysis that incorporated a surface-water diversion term in the statistical model rendered the time-trend term insignificant, indicating that the time trend in the models served as a surrogate for other variables, some of which reflect underlying hydrologic conditions. A more precise explanation of the potential causes of the time trend was not obtained with the available data. However, the model results with the surface-water diversion term indicate that much of the trend of 2.2 percent per lag year in diffP resulted from applying a PCC to estimate pumpage under hydrologic conditions different from those under which the PCC was measured. Although there is no evidence to conclude that the upward time trend determined in the data for this 5-year period would hold in the future, historical static ground-water levels in the study area generally have exhibited small variations over multidecadal time scales. Therefore, the approximately 2 percent per lag year trend determined in these data is expected to be a reasonable guideline for estimating potential errors in the PCC approach resulting from temporally varying hydrologic conditions between time of PCC measurement and pumpage estimation. Comparisons also were made between total, or aggregated, pumpage for a network of wells as computed by the PCC approach and the TFM approach. For 100 wells and a lag of 4 years between PCC measurement and pumpage estimation, there was a 95-percent probability that the difference between total network pumpage measured by the PCC approach and that measured using a TFM would be between 5.2 and 14.4 percent. These estimates were based on a bias of 2.2 percent per lag year estimated for the period 1998-2002 during which hydrologic conditions were known to have changed. Using the same assumptions, the estimated d

Scientific Investigations Report↗
USGS2014

HydroClimATe: hydrologic and climatic analysis toolkit

The potential consequences of climate variability and climate change have been identified as major issues for the sustainability and availability of the worldwide water resources. Unlike global climate change, climate variability represents deviations from the long-term state of the climate over periods of a few years to several decades. Currently, rich hydrologic time-series data are available, but the combination of data preparation and statistical methods developed by the U.S. Geological Survey as part of the Groundwater Resources Program is relatively unavailable to hydrologists and engineers who could benefit from estimates of climate variability and its effects on periodic recharge and water-resource availability. This report documents HydroClimATe, a computer program for assessing the relations between variable climatic and hydrologic time-series data. HydroClimATe was developed for a Windows operating system. The software includes statistical tools for (1) time-series preprocessing, (2) spectral analysis, (3) spatial and temporal analysis, (4) correlation analysis, and (5) projections. The time-series preprocessing tools include spline fitting, standardization using a normal or gamma distribution, and transformation by a cumulative departure. The spectral analysis tools include discrete Fourier transform, maximum entropy method, and singular spectrum analysis. The spatial and temporal analysis tool is empirical orthogonal function analysis. The correlation analysis tools are linear regression and lag correlation. The projection tools include autoregressive time-series modeling and generation of many realizations. These tools are demonstrated in four examples that use stream-flow discharge data, groundwater-level records, gridded time series of precipitation data, and the Multivariate ENSO Index.

Techniques and Methods↗
USGS2015

Updating the USGS seismic hazard maps for Alaska

The U.S. Geological Survey makes probabilistic seismic hazard maps and engineering design maps for building codes, emergency planning, risk management, and many other applications. The methodology considers all known earthquake sources with their associated magnitude and rate distributions. Specific faults can be modeled if slip-rate or recurrence information is available. Otherwise, areal sources are developed from earthquake catalogs or GPS data. Sources are combined with ground-motion estimates to compute the hazard. The current maps for Alaska were developed in 2007, and included modeled sources for the Alaska-Aleutian megathrust, a few crustal faults, and areal seismicity sources. The megathrust was modeled as a segmented dipping plane with segmentation largely derived from the slip patches of past earthquakes. Some megathrust deformation is aseismic, so recurrence was estimated from seismic history rather than plate rates. Crustal faults included the Fairweather-Queen Charlotte system, the Denali–Totschunda system, the Castle Mountain fault, two faults on Kodiak Island, and the Transition fault, with recurrence estimated from geologic data. Areal seismicity sources were developed for Benioff-zone earthquakes and for crustal earthquakes not associated with modeled faults. We review the current state of knowledge in Alaska from a seismic-hazard perspective, in anticipation of future updates of the maps. Updated source models will consider revised seismicity catalogs, new information on crustal faults, new GPS data, and new thinking on megathrust recurrence, segmentation, and geometry. Revised ground-motion models will provide up-to-date shaking estimates for crustal earthquakes and subduction earthquakes in Alaska.

Quaternary Science Reviews↗
USGS2020

A two-stage step-wise framework for fast optimization of well placement in coalbed methane reservoirs

Coalbed methane (CBM) has emerged as a clean energy resource in the global energy mix, especially in countries such as Australia, China, India and the USA. The economical and successful development of CBM requires a thorough evaluation and optimization of well placement prior to field-scale exploitation. This paper presents a two-stage, step-wise optimization framework to obtain the optimal placement of wells for large-scale development of CBM reservoirs. In the first stage, an optimal uniform well pattern is obtained by optimizing well pattern description parameters with the particle swarm optimization (PSO) algorithm. Subsequently, the location and status (active/inactive) of each well are perturbed and optimized within the patterns through the integration of the generalized pattern search (GPS) algorithm and a quality map (QM) representing the production potential. This framework was tested in a synthetic anthracite CBM reservoir in the Qinshui basin (with high gas content and low permeability) and a real field high volatile bituminous reservoir in the Illinois basin (with low gas content and high permeability). The results show that: (i) significant variations in the net present value (NPV) exist with respect to different uniform well patterns (even for cases where the total number of wells are identical), the optima of which can be efficiently determined by the PSO within 100 numerical simulation runs; (ii) the optimization of well perturbations by the GPS results in a more noticeable improvement in NPVs for the synthetic (12.3%) than for the real field model (4.6%); (iii) for the low permeable synthetic model with narrow optimal well spacings (320 m × 200 m), the contribution of the optimization of well perturbation to the NPV increment is heavily dependent on the uniform well placement solution; (iv) for the high permeable real field model with large optimal well spacings (1300 m × 1300 m), the initial uniform well placement has a very minor effect on the subsequent well perturbation solutions in terms of NPV; (v) the proposed framework significantly outperforms the conventional well-by-well concatenation procedure in terms of computational efficiency, robustness and optimal criteria set for production potential.

International Journal of Coal Geology↗
USGS2017

Methods for estimating selected low-flow frequency statistics and harmonic mean flows for streams in Iowa

A statewide study was conducted to develop regression equations for estimating six selected low-flow frequency statistics and harmonic mean flows for ungaged stream sites in Iowa. The estimation equations developed for the six low-flow frequency statistics include: the annual 1-, 7-, and 30-day mean low flows for a recurrence interval of 10 years, the annual 30-day mean low flow for a recurrence interval of 5 years, and the seasonal (October 1 through December 31) 1- and 7-day mean low flows for a recurrence interval of 10 years. Estimation equations also were developed for the harmonic-mean-flow statistic. Estimates of these seven selected statistics are provided for 208 U.S. Geological Survey continuous-record streamgages using data through September 30, 2006. The study area comprises streamgages located within Iowa and 50 miles beyond the State's borders. Because trend analyses indicated statistically significant positive trends when considering the entire period of record for the majority of the streamgages, the longest, most recent period of record without a significant trend was determined for each streamgage for use in the study. The median number of years of record used to compute each of these seven selected statistics was 35. Geographic information system software was used to measure 54 selected basin characteristics for each streamgage. Following the removal of two streamgages from the initial data set, data collected for 206 streamgages were compiled to investigate three approaches for regionalization of the seven selected statistics. Regionalization, a process using statistical regression analysis, provides a relation for efficiently transferring information from a group of streamgages in a region to ungaged sites in the region. The three regionalization approaches tested included statewide, regional, and region-of-influence regressions. For the regional regression, the study area was divided into three low-flow regions on the basis of hydrologic characteristics, landform regions, and soil regions. A comparison of root mean square errors and average standard errors of prediction for the statewide, regional, and region-of-influence regressions determined that the regional regression provided the best estimates of the seven selected statistics at ungaged sites in Iowa. Because a significant number of streams in Iowa reach zero flow as their minimum flow during low-flow years, four different types of regression analyses were used: left-censored, logistic, generalized-least-squares, and weighted-least-squares regression. A total of 192 streamgages were included in the development of 27 regression equations for the three low-flow regions. For the northeast and northwest regions, a censoring threshold was used to develop 12 left-censored regression equations to estimate the 6 low-flow frequency statistics for each region. For the southern region a total of 12 regression equations were developed; 6 logistic regression equations were developed to estimate the probability of zero flow for the 6 low-flow frequency statistics and 6 generalized least-squares regression equations were developed to estimate the 6 low-flow frequency statistics, if nonzero flow is estimated first by use of the logistic equations. A weighted-least-squares regression equation was developed for each region to estimate the harmonic-mean-flow statistic. Average standard errors of estimate for the left-censored equations for the northeast region range from 64.7 to 88.1 percent and for the northwest region range from 85.8 to 111.8 percent. Misclassification percentages for the logistic equations for the southern region range from 5.6 to 14.0 percent. Average standard errors of prediction for generalized least-squares equations for the southern region range from 71.7 to 98.9 percent and pseudo coefficients of determination for the generalized-least-squares equations range from 87.7 to 91.8 percent. Average standard errors of prediction for weighted-least-squares equations developed for estimating the harmonic-mean-flow statistic for each of the three regions range from 66.4 to 80.4 percent. The regression equations are applicable only to stream sites in Iowa with low flows not significantly affected by regulation, diversion, or urbanization and with basin characteristics within the range of those used to develop the equations. If the equations are used at ungaged sites on regulated streams, or on streams affected by water-supply and agricultural withdrawals, then the estimates will need to be adjusted by the amount of regulation or withdrawal to estimate the actual flow conditions if that is of interest. Caution is advised when applying the equations for basins with characteristics near the applicable limits of the equations and for basins located in karst topography. A test of two drainage-area ratio methods using 31 pairs of streamgages, for the annual 7-day mean low-flow statistic for a recurrence interval of 10 years, indicates a weighted drainage-area ratio method provides better estimates than regional regression equations for an ungaged site on a gaged stream in Iowa when the drainage-area ratio is between 0.5 and 1.4. These regression equations will be implemented within the U.S. Geological Survey StreamStats web-based geographic-information-system tool. StreamStats allows users to click on any ungaged site on a river and compute estimates of the seven selected statistics; in addition, 90-percent prediction intervals and the measured basin characteristics for the ungaged sites also are provided. StreamStats also allows users to click on any streamgage in Iowa and estimates computed for these seven selected statistics are provided for the streamgage.

Iowa↗
USGS1987

Producing Alaska interim land cover maps from Landsat digital and ancillary data

In 1985, the U.S. Geological Survey initiated a research program to produce 1:250,000-scale land cover maps of Alaska using digital Landsat multispectral scanner data and ancillary data and to evaluate the potential of establishing a statewide land cover mapping program using this approach. The geometrically corrected and resampled Landsat pixel data are registered to a Universal Transverse Mercator (UTM) projection, along with arc-second digital elevation model data used as an aid in the final computer classification. Areas summaries of the land cover classes are extracted by merging the Landsat digital classification files with the U.S. Bureau of Land Management's Public Land Survey digital file. Registration of the digital land cover data is verified and control points are identified so that a laser plotter can products screened film separate for printing the classification data at map scale directly from the digital file. The final land cover classification is retained both as a color map at 1:250,000 scale registered to the U.S. Geological Survey base map, with area summaries by township and range on the reverse, and as a digital file where it may be used as a category in a geographic information system.

Alaska↗
USGS2004

Geohydrology of the French Creek Basin and simulated effects of drought and ground-water withdrawals, Chester County, Pennsylvania

This report describes the results of a study by the U.S. Geological Survey, in cooperation with the Delaware River Basin Commission, to develop a regional ground-water-flow model of the French Creek Basin in Chester County, Pa. The model was used to assist water-resource managers by illustrating the interconnection between ground-water and surface-water systems. The 70.7-square mile French Creek Basin is in the Piedmont Physiographic Province and is underlain by crystalline and sedimentary fractured-rock aquifers. Annual water budgets were calculated for 1969-2001 for the French Creek Basin upstream of streamflow-measurement station French Creek near Phoenixville (01472157). Average annual precipitation was 46.28 in. (inches), average annual streamflow was 20.29 in., average annual base flow determined by hydrograph separation was 12.42 in., and estimated average annual ET (evapotranspiration) was 26.10 in. Estimated average annual recharge was 14.32 in. and is equal to 31 percent of the average annual precipitation. Base flow made up an average of 61 percent of streamflow. Ground-water flow in the French Creek Basin was simulated using the finite-difference MODFLOW-96 computer program. The model structure is based on a simplified two-dimensional conceptualization of the ground-water-flow system. The modeled area was extended outside the French Creek Basin to natural hydrologic boundaries; the modeled area includes 40 square miles of adjacent areas outside the basin. The hydraulic conductivity for each geologic unit was calculated from reported specific-capacity data determined from aquifer tests and was adjusted during model calibration. The model was calibrated for above-average conditions by simulating base-flow and water-level measurements made on May 1, 2001, using a recharge rate of 20 in/yr (inches per year). The model was calibrated for below-average conditions by simulating base-flow and water-level measurements made on September 11 and 17, 2001, using a recharge rate of 6.2 in/yr. Average conditions were simulated by adjusting the recharge rate until simulated streamflow at streamflow-measurement station 01472157 matched the long-term (1968-2001) average base flow of 54.1 cubic feet per second. The recharge rate used for average conditions was 15.7 in/yr. The effect of drought in the French Creek Basin was simulated using a drought year recharge rate of 8 in/yr for 3 months. After 3 months of drought, the simulated streamflow of French Creek at streamflow-measurement station 01472157 decreased 34 percent. The simulations show that after 6 months of average recharge (15.7 in/yr) following drought, streamflow and water levels recovered almost to pre-drought conditions. The effect of increased ground-water withdrawals on stream base flow in the South Branch French Creek Subbasin was simulated under average and drought conditions with pumping rates equal to 50, 75, and 100 percent of the Delaware River Basin Commission Ground Water Protected Area (GWPA) withdrawal limit (1,393 million gallons per year) with all pumped water removed from the basin. For average recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 18, 28, and 37 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. After 3 months of drought recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 27, 40, and 52 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. The effect of well location on base flow, water levels, and the sources of water to the well was simulated by locating a hypothetical well pumping 200 gallons per minute in different places in the Beaver Run Subbasin with all pumped water removed from the basin. The smallest reduction in the base flow of Beaver Run was from a well on the drainage divide between the French Creek Basin and the Marsh Creek Basin to the south; the simulated base flow of Beaver Run at the mouth was reduced 1 percent. The greatest reduction in the base flow of Beaver Creek was from a well close to Beaver Run; the simulated base flow of Beaver Run at the mouth was reduced 8 percent. The simulations showed that (1) if the contributing area of a well is in a basin, pumping will affect stream base flow and water levels in that basin whether the well is inside or outside that basin; (2) wells in different areas of a basin away from a divide produce a similar reduction in base flow; (3) a well within a basin will derive more water from diverted base flow and less water from storage than a well on or near a basin divide; and (4) the reduction in base flow at the mouth of the stream is the same for a well in the headwaters and a well downstream near the confluence. Model simulations illustrate some of the typical analyses and results that can be produced. The model was calibrated using annual values for recharge and ground-water ET and then was run using the annual values in a seasonally independent transient mode to show changes with time. The timing and relative magnitude of some of the changes simulated with the model when viewed in terms of a normal climatic year may be subject to considerable uncertainty because of the variability in seasonal recharge and ground-water ET rates. Transient model simulations for short-term periods are indicative of possible hydrologic system response and are considered an approximation.

Pennsyvlania↗
USGS2018

Minimum energy requirements for desalination of brackish groundwater in the United States with comparison to international datasets

This paper uses chemical and physical data from a large 2017 U.S. Geological Survey groundwater dataset with wells in the U.S. and three smaller international groundwater datasets with wells primarily in Australia and Spain to carry out a comprehensive investigation of brackish groundwater composition in relation to minimum desalination energy costs. First, we compute the site-specific least work required for groundwater desalination. Least work of separation represents a baseline for specific energy consumption of desalination systems. We develop simplified equations based on the U.S. data for least work as a function of water recovery ratio and a proxy variable for composition, either total dissolved solids, specific conductance, molality or ionic strength. We show that the U.S. correlations for total dissolved solids and molality may be applied to the international datasets. We find that total molality can be used to calculate the least work of dilute solutions with very high accuracy. Then, we examine the effects of groundwater solute composition on minimum energy requirements , showing that separation requirements increase from calcium to sodium for cations and from sulfate to bicarbonate to chloride for anions , for any given TDS concentration. We study the geographic distribution of least work, total dissolved solids, and major ions concentration across the U.S. We determine areas with both low least work and high water stress in order to highlight regions holding potential for desalination to decrease the disparity between high water demand and low water supply. Finally, we discuss the implications of the USGS results on water resource planning, by comparing least work to the specific energy consumption of brackish water reverse osmosis plants and showing the scaling propensity of major electrolytes and silica in the U.S. groundwater samples.

Water Research↗