USGS Science⌕ Search

SEARCH · USGS Science

Results for “New Zealand Journal of Geology and Geophysics”

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,689 records · Page 20Linked to original sources
USGS2024

Formation of vertical columnar seismic structures and seafloor depressions by groundwater discharge in the drowned Miami Terrace platform and overlying deep-water carbonates, southeastern Florida

The presence of vertical cross-formational fluid migration passageways within sedimentary basins can profoundly impact aquifer and reservoir fluid-flow and their identification is fundamental to informing management of subsurface fluid resources (groundwater, oil, gas). In an onshore and offshore southeastern part of Florida, 2D/3D seismic-reflection and bathymetry data document ∼153 vertical columnar structures composed of reflection disruptions up to 790 m in the height and averaging 360 m in diameter, and ∼219 subcircular to circular seafloor depressions up to 1334 m wide. Our study focuses on these features found within the offshore shallow-marine carbonate Miami Terrace platform, which drowned approximately at the end of the middle Miocene, and within overlying Plio-Quaternary deep-water carbonate slope and drift deposits. Most columnar structures are rooted in stratiform aquifers of the Miami Terrace platform and associated with faults or fault intersections produced by Eocene and circa late Miocene tectonics. The columns commonly terminate within the platform or as subcircular depressions along an amalgamated karstic and drowning unconformity at the platform top. The columns typically stretch upwards from a zone of deep karst cavity collapse through the Miami Terrace platform with upward decreasing sag on internal reflections. Following drowning and Plio-Quaternary partial burial of the Miami Terrace platform by deep-water deposits, the subcircular depressions and faults along the platform top were points of origin for a second phase of column growth upward into the deep-water deposits. The continuation of deep platform cavity collapse and column evolution produced pockmarks along paleo-seafloors within the deep-water deposits and at the present-day sea floor. The Plio-Quaternary pockmarks formed at water depths too deep to suggest an origin related to meteoric karst above or near sea level, but rather their formation is suggested to be related to cyclic sea level falls that drove increased groundwater head and density gradients, and seafloor discharge of offshore freshened groundwater sourced from the underlying platform. Plausibly, mixing of freshened groundwater and seawater at the seafloor discharge sites drove dissolution of the host deep-water deposits, which together with erosion by groundwater venting and current scouring formed the pockmarks. Seaward of the Plio-Quaternary seafloor pockmarks, at the late-middle Miocene upper slope of the Miami Terrace platform and along the regional karst/drowning unconformity is a slope-parallel band of ∼189 densely distributed subcircular seafloor depressions with diameters up to 1334 m at water depths up to ∼660 m. It is plausible that along the upper slope, faults and fractures produced by gravity-driven slope instability and possibly tectonics formed a dense network of fluid passageways that promoted upward artesian freshened groundwater flow to sites of discharge where mixing with seawater generated limestone dissolution and the depressions. But tectonic uplift may have forced emersion and initial meteoric sinkhole formation circa late Miocene with later enhancement by freshened groundwater discharge and bottom current erosion.

Florida↗
USGS2024

A journey to the center of the USGS National Strong-motion Project processing and beyond

The United States Geological Survey (USGS) National Strong Motion Project (NSMP) has the primary U.S. government responsibility to acquire, process, and disseminate significant strong-motion earthquake ground motion records measured at surficial free-field stations, structures (buildings, dams, and bridges, and geotechnical arrays to the earthquake engineering community. As a result of the deployment of modern seismic instrumentation and growth of tools such as web-services, earthquake data from U.S. and international seismic networks are more accessible than ever. Our mission is to provide raw and processed strong-motion waveforms with PGA values greater than 0.1%g for M3.0 earthquakes and larger in California and M4.0 and larger within the conterminous US, Hawaii, Puerto Rico, and Alaska. Datasets of interest to the engineering and geophysics communities, such as event sequences in areas of induced seismicity and significant global events, are also processed and posted at the Center for Engineering Strong Motion Data (CESMD) at strongmotioncenter.org when available through collaboration with the international strong-motion data community. Here we outline (1) the NSMP’s current workflow to acquire, process, and distribute data at CESMD; (2) our new endeavours and collaborations focusing on comparison and integration of waveform processing software, development of techniques for metadata quality checks before and after earthquakes, and construction of a dynamic site characterization repository; and (3) our topics for possible collaboration topics across the global strong-motion community.

Conference Paper↗
USGS2024

Hyperspectral (VNIR-SWIR) analysis of roll front uranium host rocks and industrial minerals from Karnes and Live Oak Counties, Texas Coastal Plain

VNIR-SWIR (400–2500 nm) reflectance measurements were made on the surfaces of various cores, cuttings and sample splits of sedimentary rocks from the Tertiary Jackson Group, and Catahoula, Oakville and Goliad Formations. These rocks vary in composition and texture from mudstone and claystone to sandstone and are known host rocks for roll front uranium occurrences in Karnes and Live Oak Counties, Texas. Spectral reflectance profiles, 569 in total, were reduced to 125 representative spectral signatures, which were analyzed using the U.S. Geological Survey's (USGS) Material Identification and Characterization Algorithm (MICA). MICA uses an automated continuum-removal procedure together with a least-squares linear regression to determine the fit of observed sample spectral absorption features to those of reference mineral standards in a spectral library. The reference minerals include various clay, mica, carbonate, ferric and ferrous iron minerals and their mixtures. In addition, absorption feature band-depth analysis was done to identify rock surfaces exhibiting absorption features related to uranium and zeolite minerals, which were not included in the command files used to execute MICA. Rocks from each of the four geologic units produced broadly similar spectral signatures as a result of comparable mineral compositions, but there were some notable differences. For example, Ca- and Na-montmorillonite was matched most frequently to the spectral absorption features in 2-μm (∼2000–2500 nm) wavelengths, while goethite occurred often at 1-μm (∼400–1000 nm) wavelengths. The latter is related to limonitic iron-staining in and around oxidized zones of the uranium roll front as described in previous papers. Rocks of the Jackson Group differed from those of the Catahoula, Oakville and Goliad units in that the former exhibited spectral features we interpret as being due to the presence of lignite-bearing mudstone layers. Goliad rocks exhibit spectral features related to dolomite, gypsum, anhydrite, and an unidentified green clay mineral that is possibly glauconite . Jackson Group rocks also exhibit weak but well-resolved absorption features at 964 and 1157 nm related to either or both zeolite minerals clinoptilolite and heulandite. These zeolite minerals and a few spectra exhibiting hydrous silica absorption features are indicative of alteration of volcanic glass in tuffaceous mudstone and claystone layers. A few sample spectra exhibited strong absorption features at around 1135 nm related to the uranium mineral coffinite. Both the 1135 nm coffinite and 1157 nm zeolite absorption features overlap somewhat, potentially making them difficult to distinguish without additional hyperspectral field, laboratory or remote sensing data. The results of this study were compared to mixtures of minerals described for ore, gangue and alteration minerals in deposit models for sandstone-hosted uranium, sedimentary bentonite and sedimentary zeolite. Use of these spectra can help facilitate mapping of both waste materials from the legacy mining of the above commodities, as well as future exploration and resource assessment activities.

Texas↗
USGS2025

Tracking diagenetic alteration of magnetic susceptibility in thrust ridge and slope basin sediments of the Cascadia margin (ODP Sites 1249 and 1252; IODP Site U1325)

We investigated sediment core records from the Cascadia Margin (Ocean Drilling Program Sites 1249 and 1252 at Hydrate Ridge; Integrated Ocean Drilling Program Site U1325 offshore Vancouver Island) using a Zr/Rb heavy mineral proxy from X-ray fluorescence (XRF) core scanning to identify intervals of primary detrital magnetic susceptibility (κ) and predict intervals where diagenesis caused magnetite dissolution by hydrogen sulfide. We also measured total sulfur (TS) content, grain size distributions, total organic carbon (TOC) content, and the magnetic mineral assemblage to further constrain the role of diagenesis on κ. Understanding how κ can be used to better characterize the varied effects of detrital and diagenetic signals in marine settings is important for understanding biogeochemical cycling and records of paleoenvironmental change. The upper 100 m of slope basin Site 1252 contains multiple intervals (> 90 m total) of decreased κ correlated with elevated TS content, consistent with dissolution of magnetite and precipitation of pyrite, iron monosulfides, and/or elemental sulfur. Similarly at the other slope basin site, Site U1325, κ is lower and TS is elevated in the interval between 24 and 51 mbsf, due to sulfide formation. At both slope basin sites, these low κ intervals correspond with high TOC, suggesting the possibility that organoclastic sulfate reduction (OSR) is likely a major driver of diagenetic alteration of κ at these sites. High TS:TOC ratios at Site U1325 suggests anaerobic oxidation of methane (AOM) during sulfate-methane transition zone (SMTZ) migration may have contributed to alteration of κ. In contrast, within the upper 90 m of Site 1249, a methane seep site at the summit of Hydrate Ridge, κ is almost entirely altered by diagenetic processes, with much of the low κ explained by a high degree of iron sulfide formation, while some intervals are affected by precipitation of magnetic iron sulfides that maintain or even increase κ. The presence of abundant methane seepage and gas hydrate as well as chemosynthetic seafloor fauna at this site, suggests that sulfide is released to the water column and AOM, rather than OSR, drives diagenetic alteration of κ at this site. Overall, the slope basin sites show episodic variation of κ that is influenced by TOC content, likely driven by changes in marine primary productivity and sedimentation rate, while the seep site shows consistently altered κ with lower TS content and no correlation with TOC. Methane seep environments likely experience loss of hydrogen sulfide to the water column and oxidation of hydrogen sulfide by seafloor seep fauna, which limits the amount of solid phase sulfur (pyrite, iron monosulfides, elemental sulfur) that can be precipitated within the sediments. In contrast, the migration of a buried SMTZ at slope basin sites results in enhanced sulfur precipitation within the sediments. This integrated magnetic and geochemical approach reveals the diagenetic production pathway and residence time of sulfide with the sediment column ultimately controls the style and degree of diagenetic loss of κ in marine sediments. This approach works best in environments with unaltered reference intervals/sites, consistent magnetic mineralogy, and clay-to-silt grain sizes.

Cascadia margin, Pacific Ocean↗
USGS2025

Rapid emplacement of the Keaiwa Lava Flow of 1823 from the Great Crack in the Southwest Rift Zone of Kilauea volcano

The Keaīwa Lava Flow of 1823 in the Southwest Rift Zone of Kīlauea volcano is unusual for its expansive pāhoehoe sheet flow morphology and lack of constructive vent topography, despite having a similar tholeiitic basalt composition to other lavas erupted from Kīlauea. This lava flow issued from a ∼10-km-long continuous fissure now known as the Great Crack, and has an unusually thin sheet flow morphology with margin thicknesses of ∼15–110 cm (average of 42 cm). Based on field observations of the lava flow at its fissure vent (e.g., drain-back features), we propose that the Great Crack formed, or at least significantly widened, just prior to and syn-eruptively with this 1823 eruption. The absence of pyroclastic cones or spatter ramparts indicates that the eruption consisted of a rapid outpouring of relatively degassed lava as the fissure unzipped. The rapidly moving lava flow overtopped pre-existing tumuli and scoria cones (e.g., Lava Plastered Cones) up to ∼10 m tall. Glass and whole-rock chemistry yield homogeneous compositions for the lavas erupted from the Great Crack, with glass compositions of 6.40 ± 0.10 wt% MgO and whole-rock compositions of 7.39 ± 0.07 wt% MgO. Lava pads erupted from a short western fissure system are richer in mafic minerals (e.g., olivine and clinopyroxene), and show slightly more MgO-rich whole-rock compositions (7.79 ± 0.05 wt%). MgO-in-glass thermometry on juvenile spatter yield eruption temperatures of 1153 ± 13°C that are typical of Kīlauea lavas. Thus, the extensive sheet-like lava flow morphology is not a direct consequence of unusual magmatic or rheological conditions (i.e., low viscosity). Instead, the flow morphology is associated with high effusion rates caused by sudden drainage of uprift magma as it erupted from the Great Crack. Lava flow modeling on a 2-m-resolution digital elevation model indicates that a minimum bulk effusion rate of ∼5800 m 3 /s (∼3500 m 3 /s dense rock equivalent) and a minimum flow velocity of ∼11 m/s are required for the lava flow to overcome the topography of the Lava Plastered Cones. This effusion rate is among the highest inferred for eruptions in Hawaiʻi and around the world. This study highlights a less frequent eruption style at Hawaiian volcanoes characterized by a sudden outpouring of lava from an unusual fissure system. Local eyewitness accounts indicate that the 1823 eruption was preceded by seismicity. Given the complex magmatic-volcanic-tectonic relations across Kīlauea, we speculate that the south flank could have slipped over one or more events that ultimately triggered unzipping of the Great Crack and passive release of briefly stored uprift magma. An eruption similar to 1823 at Kīlauea or Mauna Loa, with an eruptive timeframe that could be as short as an hour, with high effusion rates and rapid flow front velocities, would not easily allow for a timely response.

Hawaii↗
USGS2025

Utilizing downhole datasets for modelling the aeromagnetic signature of the Iron Creek Co-Cu deposit in the Idaho Cobalt Belt

The Idaho Cobalt Belt in east-central Idaho is host to some of the largest domestic Co resources, including the informal Iron Creek deposit. The two main ore zones of this deposit, the Iron Creek and the Ruby, are hosted in greenschist-grade interbedded argillite/siltstone and quartz-rich units of the Mesoproterozoic Apple Creek Formation of Lemhi Group. The primary ore mineral is cobalt-bearing pyrite, which occurs with pyrrhotite, chalcopyrite, and magnetite. This study integrates mineralogical and core-scale geophysical measurements with regional-scale aeromagnetic surveys. The high magnetite content within the Ruby zone produces elevated magnetic susceptibility, but the relatively limited spatial footprint of the ore zone would produce a small-scale anomaly that may be overlooked in regional surveys. The low magnetite content in the Iron Creek zone results in low magnetic susceptibility, creating a relatively low amplitude geophysical response. By characterizing the magnetic properties and mineralogy of these ore zones, this study enhances the interpretation of aeromagnetic data, enabling the identification of small or faint anomalies as potential Co targets. These findings improve can improve exploration strategies, both within the Idaho Cobalt Belt as well as for similar deposit types globally.

Idaho↗
USGS2024

On-orbit calibration and performance of the EMIT imaging spectrometer

The Earth surface Mineral dust source InvesTigation (EMIT) is a remote visible to shortwave infrared (VSWIR) imaging spectrometer that has been operating onboard the International Space Station since July 2022. This article describes EMIT's on-orbit spectroradiometric calibration and validation. Accurate spectroscopy is vital to achieve consistent mapping results with orbital imaging spectrometers. EMIT takes a unique approach to this challenge, with just six optical elements, no shutter, and no onboard calibration systems. Its simple design focuses on uniformity and stability to enable vicarious spectroradiometric calibration. Our experiments demonstrate that this approach is successful, approaching the fidelity of manual field spectroscopy in some cases, and enabling new and more accurate products across diverse Earth science disciplines. EMIT achieves several notable firsts for an instrument of its class. It demonstrates successful on-orbit adjustments of Focal Plane Array (FPA) alignment with sub-micron precision. It offers spectral uniformity better than 98%. Optical artifacts in the measurement channels are at least three orders of magnitude below the primary solar-reflected surface signals. Its noise performance enables percent-level discrimination in the depths of mineral absorption features. In these aspects, EMIT satisfies the stringent performance needs for the next generation of VSWIR imaging spectrometers to observe the Earth's ecosystems, geology , and water resources.

Remote Sensing of Environment↗
USGS2024

Limited evidence of late Quaternary tectonic surface deformation in the eastern Tennessee seismic zone, USA

The ~300-km-long eastern Tennessee seismic zone (ETSZ), USA, is the second-most seismically active region east of the Rocky Mountains. Seismicity generally occurs below the Paleozoic fold-and-thrust belt within the Mesoproterozoic basement, at depths of 5–26 km, and earthquake magnitudes during the instrumental record have been moment magnitude ( M w )≤4.8. Evidence of surface deformation may not exist or be difficult to detect because of the vegetated and soil-mantled landscape, landslides, locally steep topography, anthropogenic landscape modification, or long, irregular recurrence intervals between surface-rupturing earthquakes. Despite the deep seismicity, analog models indicate that accumulation of strike-slip or oblique-slip displacement at depth could be expected to propagate upward through the Paleozoic section, producing a detectable surficial signal of distributed faulting. To identify potential surface deformation, we interrogated the landscape at different spatial scales. We evaluated morphotectonic and channel metrics, such as channel sinuosity and catchment-scale hypsometry. Additionally, we mapped possible fault-related topographic features on 1-m lidar. Finally, we integrated our observations with available bedrock and Quaternary surficial mapping and subsurface geophysical data. At a regional scale, most morphotectonic and channel metrics have a strong lithologic control. Within smaller regions of similar lithology, we observe changes in landscape metrics like channel sinuosity and catchment-scale hypsometry that spatially correlate with new lineaments identified in this study and previously mapped east–west Cenozoic faults. These faults have apparent left-lateral offsets, are optimally oriented to slip in the current stress field, and match kinematics from recent focal mechanisms, but do not clearly preserve evidence of late Pleistocene or Holocene tectonic surface deformation. Most newly mapped lineaments might be explained by either tectonic or non-tectonic origins, such as fluvial or karst processes. We also re-evaluated a previously described paleoseismic site and interpret that the exposure does not record evidence of late Pleistocene faulting but instead is explained by fluvial stratigraphy.

Tennessee↗
USGS2025

Shallow geologic framework of the Mississippi Sound and the potential for sediment resources

The Mississippi Sound, an estuarine environment located between the mainland and barrier islands bordering the northern Gulf of America (formerly the Gulf of Mexico), serves as a vital ecosystem for the States of Mississippi and Alabama. Spanning approximately 100 kilometers from east to west and covering 1,400 square kilometers, the sound is home to marine industry and ports, and its shallow and brackish waters sustain a diverse array of marine life. Barrier islands along the southern edge of the sound separate the microtidal estuary from the Gulf of America. This protection from gulf wave action mediates current flow within the sound, resulting in predominantly fine-grained sediment deposition along the seafloor. This study, conducted by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers, provides insight on fluvial and tidal processes spanning the past 5,000 years. The report synthesizes existing research to provide a comprehensive overview of the sound geology, from Pleistocene origins to present-day morphology, and utilizes high-resolution single channel seismic profiles and sediment data to identify and map sedimentary deposits and morphologic features at and below the seafloor. Despite its ecological significance, the Mississippi Sound faces environmental challenges, including water-quality issues, habitat degradation, storm-induced erosion, and the ongoing threats of sea-level rise and environmental changes. This study uses the present-day understanding of the sound's geology to inform coastal management decisions, hazard assessment, and potential mineral resources.

Louisiana, Mississippi↗
USGS2025

Melt generation sources and conditions in the wake of a migrating slab window: Geochemistry and petrology of the million-year history of primitive volcanism at Clear Lake volcanic field, California

Clear Lake volcanic field (CLVF) is the northernmost and youngest (~2.2 Ma to 8 ka) of the volcanic centers distributed along the San Andreas transform fault in western California. The initial phase of CLVF volcanism (interval one) occurred between ~2.2 and 1.3 Ma and extends ~35 km southeast of Clear Lake, forming a semi-continuous upland plateau capped by lava flows, with isolated volcanic remnants on the periphery. This volcanism is broadly characterized by geochemically primitive compositions that reflect three source compositions and conditions of melt generation. (1) Partial melting of upwelling asthenospheric mantle lherzolite at moderate pressures (1.2–1.4 GPa) and temperatures (1297–1329 °C) produced high-CaO (9.8–11.3 wt %) basalts with high Al 2 O 3 (16.8–17.6 wt %), Mg#s (66–70), MgO (8–10 wt %), Ni (103–262 μg/g), and Cr (284–609 μg/g). These high-CaO basalts contain olivine (Fo 87–91 ) phenocrysts with Cr-spinel inclusions ± subordinate plagioclase and crop out only in the southern part of the CLVF. (2) Partial melting of depleted sub-continental lithospheric mantle harzburgite at variable pressures (0.7–1.5 GPa) and temperatures (1097–1299 °C) produced a compositional continuum of med-K 2 O, calc-alkaline, high-MgO basalts through high-MgO andesites with high Mg#s (67–77), MgO (8–14 wt %) and high Ni and Cr abundances (154–439 and 340–1124 μg/g, respectively). Mineral assemblages are olivine (Fo 88–93 ) with Cr-spinel inclusions ± subordinate clinopyroxene, orthopyroxene and plagioclase. Small (<2.5 cm) mantle harzburgite xenoliths and mantle olivine xenocrysts are also found in several of these samples. These high-MgO basalts through andesites represent the largest volume of primitive compositions and have erupted predominantly along the main, fault-controlled northwest-southeast trending axis of volcanism with peripheral outcrops to the north, west, and east. (3) Partial melting of the Gorda eclogite slab edge produced adakitic silicic slab melts with strong depletion in the heavy rare earth elements (Yb = 0.6 μg/g). Subsequent reaction of those melts with depleted ultramafic rocks during ascent imprinted the adakitic dacites with high Mg#s (65–78) and elevated Ni (117–210 μg/g) and Cr (191–283 μg/g). Phenocrysts of orthopyroxene (En 87–94 ) with spinel inclusions (Cr# = 80–88) and extremely Ni-rich (9483 μg/g) olivine cores (Fo 84–93 ) record those reactions. Small-volume outcrops of the adakites on the eastern periphery of the CLVF track the passing slab edge. The trio of melting sources recorded by early CLVF magmatism reflect the tectonically complex environment and the hot (1097–1329 °C), shallow (0.7–1.5 GPa) melting conditions for these primitive compositions and provide estimates of the heat delivered to the crust. Over time, this flux led to maturation of the CLVF magmatic system toward the more voluminous and silicic volcanism that characterizes the balance of its subsequent volcanic history and maintains the present-day anomalously high heat flow in the region. The current interval (interval four) of volcanic activity at CLVF is characterized by low-volume, fault-controlled eruptions of basaltic andesite and andesite suggestive of mantle magma and heat delivery to the crust, similar to interval one. This analogous activity provides motivation for the current study and begs the question of whether the system is undergoing thermal priming for renewed silicic volcanism.

California↗
USGS1997

Integrity of production wells and confining unit at the Naval Weapons Industrial Reserve Plant, Dallas, Texas, 1995

Ground water in the shallow alluvial aquifer is contaminated at the Naval Weapons Industrial Reserve Plant, Dallas, Texas. Five production wells at the site are cased through the alluvial aquifer and underlying units and are screened in either the Paluxy or Twin Mountains aquifers. Three abandoned wells, originally completed in the Twin Mountains aquifer but filled with drilling mud in 1958, also penetrate the alluvial aquifer. The Paluxy and Twin Mountains aquifers are used for drinking-water supplies in and around the Dallas-Fort Worth metroplex. Trichloroethylene and its degradation products, dichloroethylene and vinyl chloride, and the metal chromium previously have been detected in the shallow alluvial aquifer. Current (1995) analyses of water-quality samples taken from the static water column of the five production wells and one of the abandoned wells indicate no trichloroethylene, dichloroethylene, or vinyl chloride in the water column of these wells. Chromium was detected in all samples, but concentrations were less than the practical quantitation limit, which is the regulatory action level for this site. The results of borehole geophysical log analysis indicate that two of the production wells could have vertically connected intervals where cement bonding in the well annulus is poor. The other production wells have overall good bonding. Temperature logs do not indicate flow behind casing except in the screened interval of one well. Geophysical logs show the Eagle Ford Shale ranges from 147 to 185 feet thick at the site. The Eagle Ford Shale has low permeability and a high plasticity index. These physical characteristics make the Eagle Ford Shale an excellent confining unit.

Texas↗
USGS2024

Development of a hydrogeologic visualization model for western Sarpy County, Nebraska

Population in western Sarpy County, Nebraska, has steadily increased over the last several decades and has led to increased groundwater use for domestic purposes. To meet the increase in demand, the Papio-Missouri River Natural Resources District is seeking to use all available sources of groundwater in western Sarpy County. Additionally, elevated groundwater nitrate plus nitrite as nitrogen concentrations were detected, indicating the need to better understand the groundwater quality of the area. Although the general geology of the area is understood, the area does not have detailed information on the extent of the various aquifers, particularly the Dakota aquifer. To characterize these aquifers, the Papio-Missouri River Natural Resources District invested in airborne electromagnetic surveys of the area to better understand the subsurface geology. Although these surveys improved understanding of the groundwater systems in the area, the Papio-Missouri River Natural Resources District wanted to integrate the subsurface information with available water-quality and groundwater-level data. In response, the U.S. Geological Survey, in cooperation with the Papio-Missouri River Natural Resources District, the Nebraska Natural Resources Commission, and the Nebraska Department of Natural Resources, assembled geologic, hydrogeologic and nitrate plus nitrite as nitrogen information for the selected area into a three-dimensional visualization computer software package called GeoScene3D. The completed GeoScene3D project was assembled to provide a visualization of the groundwater systems and associated water-quality results in Sarpy County and to provide the Papio-Missouri River Natural Resources District managers with information that can be used to make more informed groundwater resource-planning decisions in the future. This report details the development of a three-dimensional model created within GeoScene3D to visualize the subsurface, particularly the Dakota Sandstone in western Sarpy County.

Nebraska↗
USGS2026

2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results

The U.S. Geological Survey recently updated the National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The first version of the PRVI NSHM was released in 2003, and therefore this 2025 update includes over 20 years of new geologic, geophysical, and engineering data, methods, and models. Updates follow similar efforts performed in the recent 2023 50-state NSHM. However, this is the first NSHM in which we: (1) apply an inversion methodology to subduction interface fault sources in the earthquake rupture forecast (ERF) model; (2) develop scaled backbone median ground-motion models and independent aleatory variability models that are applied in the ground-motion characterization (GMC) model; and (3) calculate epistemic uncertainty related to alternative scenarios in the ERF and GMC models for all grid points in the study region. Long-term time-independent mean hazard calculations were performed for peak ground acceleration and 5%-damped pseudospectral acceleration at 21 spectral periods from 0.01- to 10.0-s, for eight National Earthquake Hazards Reduction Program site conditions ranging from V S30 = 150 to 1500 m/s, and for 2%, 5%, and 10% in 50-year probabilities of exceedance (return periods of 2475, 975, and 475 years, respectively). Epistemic uncertainty, in the form of selected percentiles, is also provided for a suite of test sites and all grid points in the study region for limited periods, site conditions, and probabilities of exceedance. Selected results, including comparisons with the 2003 PRVI NSHM, are shown and discussed for selected periods, site conditions, and probabilities of exceedance. When comparing the 2025 PRVI NSHM with the 2003 PRVI NSHM, hazard is generally higher at shorter periods and lower at longer periods, as a result of updates in both ERF and GMC models. The 2025 PRVI NSHM is applicable for return periods greater than ∼475 or less than ∼10,000 years.

Puerto Rico, U.S. Virgin Islands↗
USGS2025

Update of the 2008 provisional Enhanced Geothermal Systems (EGS) assessment for the Great Basin, USA

In response to the Energy Act of 2020, the U.S. Geological Survey (USGS) is updating the Enhanced Geothermal Systems (EGS) resource assessment for the Great Basin, USA. The previous 2008 provisional assessment estimated how much electricity could be generated from EGS resources of the western United States using models of electric-grade heat, models of heat extraction over time, and estimates of how much rock might be stimulated to produce viable amounts of heat. Herein, a similar conceptual strategy is applied, using updated models of heat extraction as a function of fracture spacing and well distance. Previously used reservoir heat delivery models are updated to have a dependence on fracture and well spacing, potentially improving future estimates of EGS resources as ongoing research provides a better understanding about the success of reservoir stimulation as a function of geology and location. For a range of well distances (250-1000 m) and fracture spacings (1-50 m), heat extraction efficiency ranges from 25-62%, demonstrating the importance of accounting for the most likely results of proven viable fracturing technologies. Although fracturing is important, the biggest uncertainty by far in estimating the EGS resource for the Great Basin is estimating which geologic units at what depths can be stimulated sufficiently to produce geothermal energy economically and efficiently. Uncertainties in these factors yield estimates that range over two orders of magnitude with an upper limit of ~174 terawatts-thermal (TWth) produced for 30 years from the upper 7 km of the crust. This upper limit would require significant technological advances to access most of the electric-grade resource across the Great Basin. Assuming that 1% of this estimate will be accessible in the next few decades gives a resource estimate similar to that made in the 2008 provisional assessment. These estimated EGS heat extraction rates far exceed ( greater than 100x) the natural geothermal heat production rate, thus geothermal electricity production at these rates might not be sustainable unless heat is also recharged from other sources (e.g., excess solar energy when supply exceeds demand). In addition to assessment maps and cumulative estimates, the new models of fractured reservoirs developed herein can be used to estimate steady power production given a set of fractures and well spacing, and estimates can be made for setback distances to ensure no thermal interference with nearby powerplants.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS1993

Hydraulic properties of the Madison aquifer system in the western Rapid City area, South Dakota

Available information on hydrogeology, data from borehole geophysical logs, and aquifer tests were used to determine the hydraulic properties of the Madison aquifer. From aquifer-test analysis, transmissivity and storage coefficient were determined for the Minnelusa and Madison aquifers, and vertical hydraulic conductivity (K v ') along with specific storage (Ss') for the Minnelusa confining bed. Borehole geophysical well logs were used to determine the thickness and location of the Minnelusa aquifer, the lower Minnelusa confining bed, and the Madison aquifer within the Madison Limestone. Porosity values determined from quantitative analysis of borehole geophysical well logs were used in analyzing the aquifer-test data. The average porosity at the two aquifer-test sites is about 10 percent in the Minnelusa aquifer, 5 percent in the lower Minnelusa confining bed, and 35 percent in the Madison aquifer. The first aquifer test, which was conducted at Rapid City production well #6, produced measured drawdown in the Minnelusa and Madison aquifers. Neuman and Witherspoon's method of determining the hydraulic properties of leaky two-aquifer systems was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. Analysis of the aquifer test for the Minnelusa aquifer yielded a transmissivity value of 12,000 feet squared per day and a storage coefficient of 3 x 10 -3 . The specific storage of the Minnelusa confining bed was 2 x 10 -7 per foot, and its vertical hydraulic conductivity was 0.3 foot per day. The transmissivity of the Madison aquifer at this site was 17,000 feet squared per day, and the storage coefficient was 2 x 10 -3 . The second aquifer test, which was conducted at Rapid City production well #5 (RC-5) produced measured drawdown only in the Madison aquifer. Hantush and Jacob's method of determining the hydraulic properties of leaky confined aquifers with no storage in the confining bed was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. The analysis of data from the RC-5 aquifer test showed that transmissivity was not equal in all directions. Hantush's method was used to determine the direction of radial anisotropy and magnitude of the major and minor axes of transmissivity. The major axis of transmissivity is at an angle of 42° east of north, and the transmissivity along this axis is about 56,000 feet squared per day. The minor axis of transmissivity is at an angle of 48° west of north, and the transmissivity along this axis is about 1,300 feet squared per day. The major axis of transmissivity intersects Cleghorn Springs, a large resurgent spring on the west edge of Rapid City. The shape of the potentiometric contours of the Madison aquifer near RC-5 agree with the orientation of the transmissivity ellipse. The average value of the storage coefficient from the isotropic analysis of the aquifer-test data was 3.5 x 10 -4 , and the average vertical hydraulic conductivity of the lower Minnelusa confining bed was 9.6 x 10 -3 foot per day.

South Dakota↗