USGS ScienceSearch

SEARCH · USGS Science

Results for “Geochronology, Thermochronology and Time Scale Research”

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.

1,663 records · Page 25Linked to original sources

New developments at the Center for Engineering Strong-Motion Data (CESMD)

The Center for Engineering Strong-Motion Data (CESMD), an internationally utilized joint center of the U.S. Geological Survey (USGS) and the California Geological Survey (CGS), provides a single access point for earthquake strong-motion records and station metadata from the CGS California Strong-Motion Instrumentation Program (CSMIP), the USGS National Strong-Motion Project (NSMP), the USGS Advanced National Seismic System, and other affiliates. The CESMD has been continuously improving its webtools to facilitate the access of strong-motion data and metadata for use in post-earthquake response and for scientific and engineering research applications. The Center provides raw and processed strong-motion data via the Engineering Data Center (EDC) and the Virtual Data Center (VDC) web portals. This paper focuses on the strong-motion products provided by the EDC where more than 48,000 records with peak ground accelerations greater than 0.1% g from over 2400 earthquakes are currently hosted. and on the ongoing efforts to develop data access tools and applications. The new developments and ongoing efforts in the EDC include: 1) enhancements to the CESMD webservices to facilitate access to station metadata, earthquake information, and strong motion records 2) new features to the interactive map interface, improving the visualization and access to earthquake, station, and record information, 3) efforts to develop a new web application tool for data format conversion from a number of data formats, 4) efforts to unify varying waveform data formats into a consistent format, 5) ongoing efforts to compile seismic station site geology, measured or inferred Vs30 values, shear-wave profiles, NEHRP site class, and available structural instrument deployment schematics, and 6) a special studies pages for research topic-specific ground motion datasets that offer uniform processing of records from a variety of sources.

Conference Paper

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow flood frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota

Climate change and future water availability in the United States

The steady rise in global temperature as a result of human activity is causing changes in Earth’s water cycle. The balance of water stored within and moving between vapor, liquid, and frozen states in the water cycle is shifting, with consequences for water availability that include increases in drought, fire weather, flooding, and heavy precipitation, as well as cryosphere decline and sea-level rise. In this chapter of the U.S. Geological Survey Integrated Water Availability Assessment—2010–20, we provide an overview of climate-change observations and projections from Earth-system model simulations that relate to future water availability, from global and national climate assessments and from the published literature. Effects of climate change on primary water-cycle components are discussed in context of how global-scale hydroclimate drivers influence regional processes within the United States. Understanding the major climate drivers impacting the water cycle is crucial to predicting future changes in water availability and developing adaptation strategies to ensure human and ecosystem water supplies. First, we provide background information on the water cycle, the climate-model ensemble simulations developed to produce projections based on warming scenarios, and attribution and certainty levels. Tipping points, self-reinforcing feedbacks, cascading effects, and compound extremes are introduced. The framework of climatic impact drivers (CIDs) outlined in the Intergovernmental Panel on Climate Change Sixth Assessment Report (IPCC AR6) is used to show primary drivers of physical change to the water cycle and to understand and predict changes in future water availability. Specific climate-change related observations and projections are discussed for water cycle components of precipitation, evapotranspiration, soil moisture, streamflow, lakes and wetlands, ice and snow, and groundwater, as well as their implications for future water availability for humans and ecosystems. The chapter concludes with a synthesis discussion of three examples of complex regional-scale hydroclimate processes that influence water availability for populations in the United States, including (1) mountain and coastal precipitation, (2) aridification and drought, and (3) the influence of forest-cover change on terrestrial water-vapor recycling.

Professional Paper

Recent scientific contributions by the U.S. Geological Survey in the San Francisco Bay and Sacramento–San Joaquin Delta Estuary

Introduction The San Francisco Bay and Sacramento–San Joaquin Delta Estuary (hereafter, Bay-Delta) is the largest estuary on the West Coast of the United States. The Bay-Delta covers more than 1,600 square miles and drains a watershed of more than 75,000 square miles, which is greater than 40 percent of California. The region surrounding the Bay-Delta is home to about 10 million people, and its habitats (fig. 1) support more than 800 plant and animal species. The waterways of the Bay-Delta are the central hub of California’s extensive freshwater delivery system, supplying water to more than 27 million Californians and 4 million acres of farmland in the Central Valley. The U.S. Geological Survey (USGS) is the primary science agency of the U.S. Department of the Interior, providing a broad range of Earth, water, biological, and mapping data and expertise to inform natural resource management across the country. This fact sheet focuses on research conducted by the USGS in the Bay-Delta region, mostly within the past 5 years. The fact sheet is organized across five major themes: water flow and water quality, fish and wildlife, wetland restoration, invasive species, and hazards. In each of these areas, the USGS works closely with Federal, State, and local agencies; academia; and non-governmental organizations to provide objective, evidence-based science. The data and knowledge gained from USGS research and monitoring are publicly available and are routinely used by agencies, including the U.S. Bureau of Reclamation, U.S. Army Corps of Engineers, U.S. Fish and Wildlife Service, California Department of Water Resources, California Department of Fish and Wildlife, California State Water Resources Control Board, Delta Stewardship Council, and many additional organizations and stakeholders to inform management of the Bay-Delta region and its resources.

California

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra

Geologic map of the northwest flank of Mauna Loa volcano, Island of Hawai‘i, Hawaii

Mauna Loa, the largest active volcano on Earth, has erupted 34 times since written descriptions became available in A.D. 1832. The most recent eruption of Mauna Loa occurred on November 27, 2022, after a 38 year hiatus; it lasted for 12 days. Some eruptions began with only brief seismic unrest, whereas others followed several months to a year of increased seismicity. Once underway, Mauna Loa’s eruptions can produce lava flows that may reach the sea in less than 24 hours, severing roads and utilities. For example, lava flows that erupted from the Southwest Rift Zone in 1950 advanced at an average rate of 9.3 kilometers per hour (5.8 miles per hour); all three lobes reached the ocean within ~24 hours. Near the eruptive vents, the flows likely traveled even faster. In terms of eruption frequency, pre-eruption warning, and rapid flow emplacement, Mauna Loa has great volcanic-hazard potential for the Island of Hawai‘i. Volcanic hazards on Mauna Loa can be anticipated, and risk substantially mitigated, by documenting its past activity to refine our knowledge of the hazards, and by alerting the public and local government officials of our findings and their implications for hazards assessments and risk. The map of the north and west flanks of Mauna Loa shows the distribution and relation of volcanic and surficial sedimentary deposits. It incorporates previously reported work published as generalized small-scale maps and a more detailed map. Within the mapped area, lava has flowed from three different source regions: the Northeast Rift Zone (22 percent), the summit (64 percent), and radial vents (14 percent). All three have different points of origin which, in turn, affect the flow characteristics and periodicity of activity. The map area includes the uppermost part of the NERZ and extends from the highest elevation––13,040 feet at the south end of the Kokoolau quadrangle, just below the summit caldera––to the sea northwest and west of the summit. Lava that erupts from the north and west flanks typically flows to the west, northwest, or north, depending on the vent location. Both morphologic lava flow types—‘a‘ā and pāhoehoe—are present. Pāhoehoe units tend to spread out or widen in low-slope regions, such as in the saddle regions between Mauna Loa and Mauna Kea or between Mauna Loa and Hualālai. In comparison, ʻaʻā flows generally produce narrower flow lobes that have higher relief. This map is the fifth in a series of five maps that will cover Mauna Loa volcano. NOTE: Map sheet 1 contains lines and type with overprint. This feature may be turned on or off in the Adobe Acrobat page display preferences.

Hawaii

Impact of land subsidence on housing sale values: Evidence from the San Joaquin Valley, California

This study assesses the impact of land subsidence on housing sale values in the San Joaquin Valley, California. The study utilizes home sale transactions and vertical land-surface displacement data from Interferometric Synthetic Aperture Radar techniques. Using fine-scale fixed effects, matching, as well as a repeat-sales approach, our results indicate that land subsidence resulted in a 2.4% to 5.8% reduction in housing sale values, with the largest reductions occurring in areas where substantial subsidence occurred. Such findings may have implications for groundwater management and can potentially help inform policy design to help mitigate the causes and impacts of land subsidence.

California

Ground-motion characterization for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands

We develop the ground-motion characterization (GMC) for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands (NSHM-PRVI) for earthquakes in active crustal, subduction interface, and subduction intraslab regimes. Using ground-motion models (GMMs) from the Next-Generation Attenuation (NGA)-West2 and NGA-Subduction projects, the GMC is parameterized by scaled-backbone models for median ground motions and by independent logic trees of aleatory variability. We introduce several novel GMC features into the U.S. NSHM: (1) use of regional ground-motion data for modeling median ground motions; (2) development of scaled-backbone models for median ground motions; and (3) development of independent logic trees of aleatory variability from variance components of GMMs and computed from multiple ground-motion datasets, and incorporating regional ground-motion data effects on variability. We compute probabilistic seismic hazard curves and maps to evaluate the contributions from the GMC components and for comparison with the 2003 GMC. Contributions of the GMC to epistemic uncertainty in seismic hazard are evaluated through spatial variations in epistemic uncertainty in hazard maps, comparisons between mean hazard curves and fractiles, as well as investigations of the contributions of logic-tree branches to hazard maps and curves. Comparisons between seismic hazard from the 2025 and 2003 GMCs allow examination of the changes in hazard curves and mapped values with 2% and 10% probabilities of exceedance in 50 years. The 2025 GMC exhibits modest changes in median predictions, relative to the 2003 GMC; overall, values of aleatory variability are higher, except at long periods ( ), resulting in changes to probabilistic ground motions at low probabilities of exceedance (<10% probability of exceedance in 50 years). Changes in hazard at 2% and 10% probabilities of exceedance in 50 years are also relatively modest (within 20%) at most sites, with the impacts from the 2025 GMC exhibiting minor reductions and increases at 0.2 and 1.0 s periods, respectively, relative to the 2003 GMC.

Puerto Rico, U.S. Virgin Islands

Treatability study to evaluate bioremediation of trichloroethene at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, 2020–22

Executive Summary Chlorinated solvents, including trichloroethene (TCE) and other chlorinated volatile organic compounds (cVOCs), are widespread contaminants that can be treated by bioremediation approaches that enhance anaerobic reductive dechlorination. Reductive dechlorination can be enhanced either through the addition of an electron donor (biostimulation) or the addition of a known dechlorinating culture (bioaugmentation) along with an electron donor. Although bioremediation has been applied at many TCE-contaminated groundwater sites, application in source zones at sites where residual dense nonaqueous phase liquid (DNAPL) is present is more limited. In this study, laboratory and field treatability tests were completed to evaluate the potential application of anaerobic bioremediation for a shallow groundwater plume containing TCE in a perched alluvial aquifer at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, which was on the National Priorities List as the New Brighton/Arden Hills Superfund site until 2019. In addition to the presence of residual DNAPL at the site, temporal variability in groundwater flow directions and input of oxygenated recharge were possible complicating factors for the application of enhanced anaerobic biodegradation in the shallow plume. The Site K plume extends beneath the footprint of Building 103, which was demolished in 2006, and soil excavations to a maximum depth of 6 feet (ft) below ground surface in 2014 were known to leave some deeper contaminated soil in place in the TCE source area. Groundwater treatment at the site, formalized as part of the 1997 Record of Decision, has been in operation since 1986 and consists of an extraction trench at the downgradient edge of the plume to collect groundwater, which is then pumped to an on-site air stripper. Groundwater concentrations in the plume have been relatively stable since treatment began, indicating a continued source of TCE in the aquifer. The desire for a destructive remedy that would enhance the removal of cVOCs in the aquifer at Site K and shorten the remediation timeframe led the U.S. Army to request that the U.S. Geological Survey conduct a groundwater treatability study to assess bioremediation. This report describes the U.S. Geological Survey bioremediation treatability study conducted during 2020–22, including pre-design site characterization to assist in formulating the bioremediation approach, laboratory experiments to support the design of the field pilot test, and implementation and 1-year performance monitoring results for the pilot test. Pre-design site characterization included the collection of soil cores for cVOC analysis and lithologic descriptions and the re-installment of three wells to obtain hydrologic measurements and initial groundwater chemistry. Relatively flat head gradients were measured at the site, and substantial decreases in water-level elevations occurred from spring to summer (May–July 2021). Continuous water-level monitoring indicated a rapid response to precipitation. Groundwater flow velocities were consistently less than 0.5 foot per day, and the pilot bioremediation test was therefore designed with short lateral distances (about 5 ft) between injection and individual monitoring points. Soil analyses confirmed that high volatile organic compound contamination was left in place in the source area. The highest concentrations were near or in clay at the base of the perched aquifer. Concentrations of cVOCs measured in the replaced wells were consistent with historical data and had a maximum TCE concentration of 57,700 micrograms per liter (μg/L), indicative of nearby residual DNAPL based on the general rule of observed concentrations exceeding 1 percent of solubility. The primary TCE daughter product detected was 1,2-cis-dichloroethene (cisDCE), which indicated limited reductive dechlorination in the plume. Groundwater in both the source and downgradient areas was relatively reducing during the pre-design characterization, particularly in the source area where methane concentrations greater than 400 μg/L were measured. Initial laboratory tests conducted using native aquifer microorganisms from the three replacement wells showed that anaerobic TCE biodegradation rates were low when biostimulated with the addition of sodium lactate as an electron donor, also known as a carbon donor, and resulted in the production of only cisDCE. Addition of a known dechlorinating culture, WBC-2, however, resulted in rapid biodegradation and production of ethene, verifying complete reductive dechlorination of TCE. Microcosms constructed with aquifer soil collected from the site were used to evaluate other electron donors besides lactate to support reductive dechlorination by WBC-2, including corn syrup as an alternative fast-release compound and whey, soy-based vegetable oil, and 3-D Microemulsion (Regenesis, San Clemente, California) as slow-release compounds. First-order rate constants for total organic chlorine removal in these WBC-2 amended microcosms were greatest with either lactate or vegetable oil as the donor, ranging between 0.061 and 0.047 per day or corresponding half-lives of 11–15 days. Testing of commercial products in other WBC-2-bioaugmented microcosms led to selection for the field pilot test of an emulsified vegetable oil product that also contained some sodium lactate as a fast-release donor. Delaying the addition of WBC-2 relative to the donor in the microcosms resulted in the most rapid overall biodegradation rates. The selected design for the pilot test utilized three separate test plots, each about 30-ft wide and 60-ft long: plots GS1 and GS2 in the source area of the plume and plot GS3 in the downgradient area of the plume near the excavation trench. Each test plot had one injection well, one monitoring well upgradient from the injection point, and 12 surrounding monitoring wells in a grid to capture variable groundwater flow directions. Donor injections, which included a bromide tracer, were completed in October 2021, immediately following baseline sampling, and the WBC-2 culture was injected about 40 days later, between November 30 and December 2, 2021. Performance monitoring conducted until December 2022 included hydrologic measurements and analyses of cVOCs, redox-sensitive constituents, dissolved organic carbon, bromide, volatile fatty acids, compound-specific carbon isotopes, and microbial communities. The biogeochemical data collected during the pilot tests in the three treatment plots showed that enhanced, complete reductive dechlorination of cVOCs in the groundwater was achieved in the GS1 and GS3 plots. In contrast, evidence of distribution of the injected amendments and subsequent biodegradation was limited in GS2, which was in an area of more heterogeneous soil lithology and low water table elevations. The molar composition of volatile organic compounds in the GS1 and GS3 plots was dominated by ethene in wells that were reached by the injected amendments by the end of the monitoring period. In the GS1 and GS3 plots, similar patterns were observed of cVOC concentrations decreasing to near detection levels, or below, at some wells sampled in July and October 2022, whereas ethene became dominant and indicated sustained complete reductive dechlorination. Baseline cVOC concentrations were more than a factor of 10 higher in the groundwater in the GS1 plot than in GS3, but no apparent inhibition of complete dechlorination occurred. As expected from the initial pre-design site data and the laboratory experiments, enhanced dissolution of residual DNAPL coupled to biodegradation was evident in the GS1 plot, where a marked increase in dichloroethene (DCE) above the initial baseline and upgradient TCE and DCE concentrations occurred. DCE concentrations subsequently declined where DNAPL dissolution was evident, concurrent with production of vinyl chloride and then predominantly ethene. Thus, overall biodegradation rates outpaced the DNAPL dissolution and desorption and DCE production in the source area. This success in complete degradation to predominantly ethene was achieved even in areas where the DCE concentrations reached a maximum of about 30,000 μg/L. Compound specific isotope analysis of carbon in TCE, cisDCE, trans-1,2-dichloroethene, and vinyl chloride was conducted to provide another line of evidence of the occurrence and extent of anaerobic biodegradation. Along a flow path in each plot that was affected by the injected amendments, carbon isotopes in the TCE and daughter cVOCs in the groundwater became isotopically heavier, indicating biodegradation. Enhanced biodegradation rates calculated from the field tests in GS1 and GS3 showed half-lives of 36.9–75.3 days for DCE degradation and 9.48–38.5 days for ethene production. Notably, these ethene production rates calculated from the field tests are consistent with the results of WBC-2-bioaugmented microcosms amended with either lactate or vegetable oil, which had half-lives for total organic chlorine removal that ranged from 11 to 15 days. These rates indicated rapid enhanced biodegradation, which is promising for application of a full-scale bioremediation remedy. Ultimately, however, the mass of residual or sorbed TCE in the aquifer that remains accessible for dissolution and biodegradation would likely control the time required for a full-scale bioremediation effort to achieve performance goals for TCE and cisDCE specified in the Record of Decision for Site K. The field pilot tests showed that the relatively low hydraulic head gradients and temporal changes in groundwater flow directions in the shallow aquifer would add complexity to a full-scale bioremediation effort. The radius of influence (ROI) at GS1 and GS3 (16.3 ft and 12.7 ft, respectively) were close to the design ROI of 15 ft. The estimated ROI at GS2 was about four times the design ROI, but may be less reliable at this location owing to groundwater flow direction. In addition, the low temperatures following WBC-2 injection in late November to early December 2021, in combination with the low hydraulic head gradients, were probably major factors in the delay observed before the onset of enhanced biodegradation following injection of the culture. Additional test injections could be beneficial to optimize the timing of donor and culture injections with the variable temperatures and hydraulic head in the shallow aquifer.

Minnesota

A comparison of electronic particle counting and pipette techniques in routine mud analysis

Electronic particle counting and pipette size-analysis techniques are compared for their relative effectiveness in large-scale routine mud analysis. Optimum laboratory procedures for conducting electronic analyses with a Model TA Coulter Counter are presented. Comparative electronic and pipette analyses conducted on the same samples resulted in equivalent precision for both techniques. Results of electronic analyses consistently showed size distributions significantly coarser than those obtained by pipette analyses. The laboratory efficiency of the relatively rapid electronic technique is far superior to that of the pipette technique, although more extensive precautionary measures are required to maintain a high level of analytical accuracy and precision. The quality of electronic analyses largely depends upon the effectiveness of standardized operational procedures established within the laboratory.

Open-File Report

Large springs in the Valley and Ridge Province in Tennessee

Approximately 2,700 miscellaneous discharge measurements for 171 large springs in 28 counties of East Tennessee, predominantly within the Valley and Ridge physiographic province, were analyzed statistically and results tabulated to summarize data useful to the Appalachian Valleys-Piedmont Regional Aquifer System Analysis study. The number of measurements at each spring ranged from 1 to 65. Information from each spring includes station numbers, latitude, longitude, spring name and where data are sufficient, an estimate of the mean discharge and the discharge exceeded by 75%, 50%, and 25% of the miscellaneous measurements. Data are referenced to locations on a 1:1,000 ,000-scale map. The highest mean spring discharge was 32.2 cu ft/sec (14,5000 gallons/min).

Tennessee

Reconciling bias in moderate magnitude earthquake ground motions predicted by numerical simulations

Recent studies found a significant underprediction in ground motion intensity measures for finite-fault simulations of moderate magnitude events in southern California relative to established ground motion models. This study aims to understand the source(s) of this bias by evaluating ground motion residuals. For this, simulations have been performed for a total of 27 well-recorded earthquakes in southern California. Systematic efforts have been employed to identify the source(s) of bias by ruling out factors that are insignificant. Preliminary findings indicate that the magnitude-area scaling used in the simulations is the likely major cause of the observed bias. Adjustment in the source attributes on event-by-event basis is underway to study if the observed bias can be reconciled.

Conference Paper

A geospatial analysis of water-quality threats from orphan wells in principal and secondary aquifers of the United States

Throughout the history of oil and gas production in the United States, millions of wells have been drilled for exploration and energy production. Hundreds of thousands of unplugged wells are no longer actively producing and are currently under orphan status, with no responsible party obligated for plugging. Orphan wells can pose threats to water resources by providing pathways for contaminants such as hydrocarbons and brines to migrate into water-supply aquifers. In this study, we investigate the potential threats to groundwater resources posed by orphan wells at the national scale. Water-quality data is extremely sparse in relation to orphan wells nationally and may not be suitable for identifying contamination from oil and gas development. We used geospatial and statistical methods to evaluate which principal and secondary aquifer systems may be most susceptible to contamination from orphan wells. Analysis involved three sets of susceptibility factors including: 1) factors related to the number and density of orphan wells; 2) factors that can threaten well integrity and contribute to transport of contaminants; and 3) factors related to groundwater withdrawal rates and the affected populations/communities in the event of water quality disturbances. From a dataset of 117,672 documented orphan wells, 64,203 fall within a principal aquifer system, while the remainder fall within a secondary aquifer system. By assessing the combination of well integrity and hydrogeologic factors within these aquifer systems, five groupings of principal aquifers were identified, where groups ranged from aquifer systems with high numbers of orphan wells, multiple well integrity threats and high withdrawals, to aquifers with a relatively low number of orphan wells, limited well integrity threats and minimal water use. Three regions of the country emerge containing aquifers with higher susceptibility to contamination from orphan oil and gas wells. These regions include 1) The Appalachian Basin (including the Pennsylvanian Aquifer System), 2) The Gulf Coast Aquifers (including the Coastal Lowlands Aquifer system) and 3) The California Aquifers (including the California Coastal Basin Aquifer system). This work is the first multivariate geospatial investigation of orphan wells and groundwater resources on a national scale, and sheds light on which aquifers are most susceptible to groundwater contamination from orphan wells.

Science of the Total Environment

Niobium-rich minerals from the Sheep Creek carbonatite-related deposits, Montana, USA

The Sheep Creek area in southwest Montana contains a number of small, tabular, carbonatite-related deposits that are rich in Nb- and REE-minerals. Two of the “vein-dikes” were mined at a small scale for Nb in the late 1950s. The most abundant Nb-bearing minerals are columbite-(Fe), nioboaeschynite-(Ce), baotite, and two varieties of Ti-Nb-Fe-oxide. An early oxide with Ti/Nb (atomic) near unity later broke down to a mixture of columbite + TiO2 with distinctly lower Nb content. Zircon intergrown with the Nb-REE minerals was dated by U-Pb methods to 116±2 Ma. This age is in agreement with previous dating (Th-Pb on monazite, Re-Os on molybdenite), and indicates that the Sheep Creek deposits are early Cretaceous in age.

Montana

Bedrock geologic map of the Woodstock quadrangle, Grafton County, New Hampshire

The bedrock geology of the Woodstock 7.5-minute quadrangle consists of highly deformed metasedimentary rocks of the Central Maine trough, including the Silurian Rangeley and Perry Mountain Formations and the Devonian Littleton Formation. The central, northern, and eastern parts of the quadrangle are underlain by the oldest rocks in the area, the Rangeley Formation. In the southwest and south-central part of the quadrangle, metaturbidites of the Perry Mountain Formation and subsequent Littleton Formation overly the Rangeley Formation in a deformed F 1 synform, herein informally called the Bagley Brook basin. The metasedimentary rocks were intruded by widespread syn- to post-tectonic granitoids of the Devonian New Hampshire Plutonic Suite and minor post-metamorphic Jurassic-Cretaceous mafic dikes of the White Mountain Plutonic-Volcanic Suite. The metasedimentary rocks were affected by at least two episodes of deformation in the Devonian Acadian orogeny. The dominant regional foliation is second-generation (S 2 /D 2 ) and formed during the development of sillimanite-muscovite mineral assemblages. Large bodies of the Early Devonian Kinsman Granodiorite intruded the metasedimentary rocks semi-concordantly during D 2 deformation. Dikes of the Late Devonian Concord Granite cut the Kinsman Granodiorite and the metasedimentary rocks and were emplaced either syn- or post-D 2 . The map pattern in the Rangeley Formation is dominated by northeast to northwest trending, moderately to steeply north-dipping F 2 and F 3 folds. Map-scale F 1 folds are defined by the Bagley Brook basin. Previous division of Rangeley Formation stratigraphy in this region into “upper” and “lower” parts was not corroborated by 1:24,000-scale mapping of lithodemic units, and rocks previously mapped as part of the Smalls Falls and Madrid Formations are here reassigned to the Rangeley Formation. Some rocks previously mapped as the lower part of the Littleton Formation are now assigned to the Perry Mountain Formation. The Littleton Formation on this map is approximately equivalent to rocks previously mapped as the upper part of the same formation. Steeply dipping fractures in the quadrangle show a preferred northeast orientation, consistent with subsurface fracture orientations in the well fields near Mirror Lake. Jurassic-Cretaceous mafic dikes and normal faults show preferred northeast orientations, similar to the fractures, suggesting that the extensional stress field that controlled dike orientation during the Mesozoic also produced the dominant brittle fabrics in the area.

New Hampshire

Bathymetric and acoustic-backscatter mapping of Lake Sammamish, Washington, during USGS field activity 2021-656-FA

In 2021, scientists from the U.S. Geological Survey (USGS) acquired bathymetric and acoustic-backscatter data within Lake Sammamish, Washington. Mapping was completed to find evidence of past earthquakes such as underwater landslides and is part of a larger USGS project to understand the overall geologic hazards history of the Cascadia Margin region. The survey was conducted using the USGS research vessel Parke Snavely outfitted with an interferometric sidescan-sonar system for swath mapping. Three estimates of bathymetric uncertainty were calculated from the post-processed bathymetric data. This report provides a summary of the mapping mission.

Washington

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington