USGS Science⌕ Search

SEARCH · USGS Science

Results for “Contributions to Mineralogy and Petrology”

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,674 records · Page 83Linked to original sources
USGS2025

A probabilistic assessment methodology for the evaluation of geologic energy storage capacity—Natural gas storage in depleted hydrocarbon reservoirs

The need for energy storage, particularly underground, where capacity and duration may far exceed battery storage technologies, is especially relevant given the increasing demands for reliable power alongside the development of intermittent renewable electricity sources. Geologic energy storage facilities already exist, and expanded use would enable storing gases such as methane and hydrogen. In 2018, a National Academies of Sciences, Engineering, and Medicine report, “Future Directions for the U.S. Geological Survey's Energy Resources Program,” recommended that the U.S. Geological Survey (USGS) prioritize assessing underground energy storage in geologic formations in the United States. The U.S. Geological Survey has since developed a methodology for assessing natural gas storage capacities in depleted hydrocarbon reservoirs on a national scale. The methodology introduced in this report prescribes three approaches for calculating gas storage capacity. This methodology relies on the availability of input data, including cumulative hydrocarbon production records, reservoir petrophysical properties, and reservoir pressure data. Assessment inputs can be obtained from public, State-level databases and propriety national-scale databases, although the use of analogs could be warranted for estimating input parameters. Probabilistic assessment results are aggregated to play, petroleum province, regional, and national scales. The steps defined in this report are demonstrated on the Michigan Basin Province, which includes the Mississippian Sandstone Gas Play and the Clinton Structural Play. This methodology could be used to systematically and consistently assess hydrocarbon plays and provinces for natural gas storage capacity across the United States.

Michigan↗
USGS2025

The GorDAS Distributed Acoustic Sensing experiment above the Cascadia locked zone and subducted Gorda Slab

The southernmost portion of the Cascadia Subduction zone in Northern California produces high rates of moderate and large earthquakes owing to subduction of the Gorda slab and deformation associated with the Mendocino Triple Junction. Distributed Acoustic Sensing (DAS) is rapidly advancing as a method for detecting earthquakes and imaging crustal structure. We have begun a long-term DAS monitoring experiment on buried telecom fiber in Arcata, California, with the goal of increasing the available recordings of moderate to large earthquakes as well as imaging seismogenic structures. We have recorded over a year's worth of data, including most aftershocks of the 2022 M w 6.4 Ferndale earthquake, though not the mainshock itself. The dataset includes numerous magnitude 3.5 and larger earthquakes including the 2023/01/01 M w 5.4 Rio Dell earthquake. Here we present initial results comparing an earthquake detection algorithm, run in real-time on the processing unit of the interrogator system, with both the ShakeAlert earthquake early warning system as well as a post-processed earthquake catalog developed with deep-learning phase-picker algorithms. The rapid onboard processing of the detector demonstrates the potential utility of DAS-based edge computing for earthquake early warning. We also verify the quality of the strain waveforms both in terms of peak amplitudes and waveform similarity using about five months of nodal seismometer data. These instruments were deployed roughly every 300 m along the ~15km long cable and validate large variations in peak strain over short distances that are seen in the DAS data. All data from time windows surrounding both the local and teleseismic earthquakes are publicly available, which will improve our understanding of both the performance of DAS systems in moderate earthquakes and earthquake hazards associated with the Gorda subduction zone.

California↗
USGS1976

Trace elements in some rocks and slag from the Yemen Arab Republic and their bearing on the iron ore near Sa'dah

Thirty-seven samples of rocks, saprolite, laterite, gossan, and slag were collected in the Yemen Arab Republic during June-July 1975 when field work was conducted to refine an interpretation of the regional geology made from LANDSAT-1 images by M. J. Grolier. These samples we reanalyzed in the laboratories of the U. S. Geological Survey, and the results were interpreted for use in programs of economic development as a guide for the evaluation of mineral resources in the Yemen Arab Republic. The work was undertaken as part of the program of the U. S. Agency for International Development with the cooperation of the Mineral and Petroleum Authority, Ministry of Economy, Yemen Arab Republic.

Open-File Report↗
USGS2024

Rock sample photogrammetry

This step-by-step protocol describes the photogrammetry process used by the U.S. Geological Survey Spokane Imaging Lab (SPIMG) lab to create 3D models of geologic samples. Steps related to photographing small objects are applicable to photogrammetry in general, however, SPIMG-specific steps involving lab hardware and software may not be.

Protocols.io↗
USGS2024

Updating the crustal fault model for the 2023 National Seismic Hazard Model for Alaska

We present the crustal fault model for Alaska, based on geologic observations, as a primary input for the 2023 revision of the U.S. Geological Survey National Seismic Hazard Model. We update the 2013 Alaska Quaternary fault and fold database to produce a simplified model of 105 fault sections and four fault zone polygons with basic geologic parameters including slip sense and rate. Significant updates include the following: (1) a slip rate of ∼53 mm/year on the Queen Charlotte Fault indicating it accommodates all of the plate boundary motion; (2) quantified slip rates on megathrust splay faults in the southern Prince William Sound region and near Kodiak Island; (3) improved details of structures in the Chugach-St. Elias orogen; (4) revision of the Castle Mountain Fault from right-lateral slip to a predominantly reverse fault; (5) improved Interior Alaska tectonic models that clarify relationships between the Denali, Totschunda, and thrust faults on both sides of the Alaska Range; (6) identified large earthquake sources in the eastern Brooks Range; and (7) omission of the Chatham Strait section of the Denali Fault. The fault model underscores that the collision of the Yakutat microplate is the dominant driver of active crustal faulting in most of Alaska.

Alaska, Yukon↗
USGS2025

Strong shaking from past Cascadia Subduction Zone earthquakes encoded in coastal landforms

Strong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large-magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty-four deep-seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

California, Oregon, Washington↗
USGS2026

Groundwater quality near an oil field in a stream-dominated recharge setting, California, USA

Alluvial valley aquifers are important sources of water supply in many areas but effects of co-located oil and gas development on these resources have not been widely reported, especially in settings where recharge is dominated by stream infiltration. Interpreting the presence of geochemical indicators in the context of hydrology, geology, and other factors provides a more complete understanding of the relations between groundwater and sources of oil-field fluids and aids in identifying risks associated with oil and gas development. Groundwater and Salinas River water samples were collected in an alluvial valley near the San Ardo Oil Field in Monterey County, California and analyzed for a wide range of dissolved chemical, gas, and isotopic constituents to determine if oil-field fluids (water and gas from oil-producing and non-producing zones) have mixed with fresh groundwater used for supply. Hydraulic gradients, age-dating tracers, and other geochemical indicators show that recharge from the Salinas River has the potential to dilute oil-field fluids that might migrate or seep into the aquifer. Groundwater and Salinas River water collected downgradient of the San Ardo Oil Field showed little or no evidence of mixing with oil-field fluids. Some samples within the oil field contained trace amounts of hydrocarbons or elevated temperatures, indicating that any potential effects from oil-field activities are minor or have been diluted by recharge from the Salinas River. The two samples with the most geochemical evidence of potential mixing with oil-field fluids (SP-18 and GW-17) were collected west of or along the Los Lobos fault, where naturally occurring hydrocarbons are near the land surface. Those samples are also near active or inactive oil-field wells, and so anthropogenic activities and pathways cannot be ruled out as a cause of trace detections of hydrocarbons and elevated temperatures in the aquifer.

California↗
USGS1983

Use of digital land-cover data from the Landsat satellite in estimating streamflow characteristics in the Cumberland Plateau of Tennessee

Characteristics of Streamflow at ungaged sites in coal-mining areas of the Cumberland Plateau may be estimated with regression equations that relate streamflow characteristics to the physiographic and climatic characteristics of the corresponding drainage basins. An experiment was performed to compare the accuracy of equations using basin characteristics derived from maps and climatological records (control group equations) with the accuracy of equations using basin characteristics derived from digital processing of Landsat spectral data as well as maps and climatological records (experimental group equations). Results of this experiment show that (with the exception of low flows and four of six annual peak logarithms) drainage area can explain more than 90 percent of the variance in all streamflow characteristics in both groups of equations. Seventeen of 39 experimental group equations that have two basin characteristics each are different from the corresponding control group equations. Five of the 17 differing experimental group equations have no counterpart in the control group because of lack of significance in control group basin characteristics, seven are slightly more accurate, and five are not measurably different in accuracy. When the equations in both groups are arranged into six flow categories, there is no substantial difference in accuracy between equations using basin characteristics derived from maps and climatological records (control group) and equations using basin characteristics derived from Landsat tapes as well as maps and climatological records (experimental group) for this particular study area, the Cumberland Plateau of Tennessee.

Tennessee↗
USGS2026

Assessing environmental drivers of denitrification in restored riverine floodplains

Restoration of impaired floodplains is an increasingly prevalent strategy for alleviating water quality concerns and reducing downstream flooding at watershed scales. Floodplains temporarily store water and slow flow velocity to promote sedimentation during overbank flooding and remove inorganic nitrogen from floodwater and groundwater via denitrification. Evaluating the impacts of different restoration strategies on denitrification can inform more strategic investments into floodplain modifications that improve water quality outcomes. Our research investigates how denitrification rates in floodplains respond to environmental factors that are actionable from an engineering perspective through design and water resources management. We seasonally measured soil denitrification enzyme activity and various environmental characteristics in 4 floodplains with different restoration design and management approaches at the confluence of the Wabash and Tippecanoe Rivers in Indiana, United States. Our results showed that denitrification rates in an agricultural floodplain were significantly lower than in restored floodplains with native vegetation. Certain soil conditions characteristic of floodplain wetlands were associated with higher denitrification, particularly elevated total nitrogen, moisture, silt, and organic matter contents. Vegetation species composition was correlated with denitrification rates. This link may reflect the direct effects of vegetation on soil conditions, such as supplying labile organic carbon, or indirect effects, such as vegetation acting as an indicator of hydrologic regime and land use. Denitrification seasonally varied, peaking in winter when nitrate supply from rivers draining agricultural watersheds in the region is also high. Substrate limitation of soil denitrification enzyme activity was most significant during the summer when overbank flooding, which replenishes soil nitrogen stocks, rarely occurs. Our findings indicate that denitrification capacity will likely be maximized in riverine floodplains that are restored as wetlands with diverse native vegetation and enhanced hydrologic connectivity. Such restoration activities promote higher denitrification rates via elevated moisture, fine sediment deposition, and soil organic matter.

Tippecanoe River, Wabash River↗
USGS2026

Biome-scale spatial patterns of avian abundance reveal proactive conservation opportunities in North American grasslands

North American grassland birds have experienced steeper population declines than any other avian guild, yet conservation efforts remain largely reactive and fragmented. We used nearly four decades of North American Breeding Bird Survey data to identify biome-scale spatial patterns (clustering) of grassland bird abundance for the Great Plains. Our results reveal an ecological core in the north-central Plains where community-level abundance is either increasing by >100% or remains high and stable, providing a strategic roadmap for a “Defend the Core” conservation approach. This approach flips the script from reactive triage centered on isolated population fragments to a proactive strategy of maintaining large-scale ecosystem integrity. Conversely, we found that population losses are more spatially clustered than wins, reflecting the relentless, one-way movement of woody encroachment and agricultural conversion. This asymmetry supports prioritizing intact landscapes, as current restoration rates are often outpaced by the scale of habitat loss. Notably, we found that community-level spatial clustering is a more robust indicator of biome condition than trends of individual flagship species, suggesting that managing for ecosystem integrity provides a more effective multi-species umbrella. Given our results, there is an opportunity for operationalizing a Great Plains Conservation Design that is ecosystem-centric and rooted in the sustainability of the private-land cattle production that maintains these open spaces. By leveraging avian abundance as a biological sensor, managers and producers can deploy a shared vision that matches the spatial scale of the threats, moving from reactive triage to proactive defense of core working grasslands in North America.

Great Plains biome↗
USGS2025

Petrogenesis and mineralization potential of the Bradley Peak komatiitic basalts, Wyoming Province

The ca. 2.72 Ga Bradley Peak greenstone terrane in the Wyoming Province contains spinifex- and cumulate-textured komatiitic volcanic rocks that may host Ni-Cu-PGE sulfide mineralization. New whole rock geochemistry classifies these rocks as Al-undepleted komatiitic basalts derived from a parental melt with ~17 wt. % MgO. Isotopic data (Sm-Nd, Re-Os) and REE profiles suggest a plume source and possible assimilation of juvenile footwall volcanic rocks. Calculated sulfur concentrations at sulfide saturation and low-detection PGE geochemistry indicate these magmas did not saturate a sulfide liquid resulting in sulfide mineralization. Although the Bradley Peak komatiitic basalt flows are unlikely to host magmatic sulfide mineralization, they are likely the source of Au for regional orogenic gold deposits.

Wyoming↗
USGS2025

Potentiometric surface maps and groundwater-level hydrographs for confined aquifers of the New Jersey Coastal Plain, 2018

The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection (NJDEP), prepared potentiometric surface maps for 10 confined aquifers of the New Jersey Coastal Plain physiographic province based on water-level measurements collected during late 2018 and early 2019 from 951 wells in New Jersey and parts of Pennsylvania and Delaware. Maps were prepared for the confined Cohansey aquifer, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system. Potentiometric surface maps indicate regional cones of depression in the following aquifers and the counties in which they are centered: Atlantic City 800-foot sand in Atlantic County, the Piney Point aquifer in Cumberland County, the Wenonah-Mount Laurel aquifer and Englishtown aquifer system in Monmouth and Ocean Counties, the Wenonah-Mount Laurel aquifer in Camden and Gloucester County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Ocean County, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system in Camden County. Cones of depression with smaller areal extents were in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand centered in Cape May County, the Piney Point aquifer centered in Ocean County, the Wenonah-Mount Laurel aquifer in Salem and Burlington Counties, the Englishtown aquifer system in Camden County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth County, and the middle aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth, Ocean, and Salem Counties. No cone of depression was interpreted in the Vincentown aquifer. Long-term hydrographs are presented for 75 wells spanning each of the 10 confined aquifers, and contain a mix of discrete water-level measurements and daily mean water levels based on continuously recorded 15-minute data. Changes of water levels during 2014–19, as indicated by the hydrographs, were compared with those of previous periods to assess any departures from historical data. During 2014–19, water levels were stable and fluctuated within similar ranges as previous periods in the following aquifers and locations: all wells in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Vincentown aquifer, the Englishtown aquifer system, the Piney Point aquifer wells in Burlington and Ocean Counties, six of eight wells in the Wenonah-Mount Laurel aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system except those within NJDEP Critical Area II. Increasing water levels in 2014–19, ongoing since historical periods, were indicated in the following aquifers and locations: Atlantic County wells in the Piney Point aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system within NJDEP Critical Area II. Water levels in the Atlantic City 800-foot sand also increased during 2014–19 in wells in Atlantic County and northern Cape May County closer to the center of the cone of depression in that aquifer, which is a response unique to this period and absent from previous periods. During 2013–19, continued decreasing water levels, ongoing since previous periods, were indicated by hydrographs of Atlantic City 800-foot sand wells in southern Cape May County, Piney Point aquifer wells in Cumberland County where the regional cone of depression is located, and two wells in the Wenonah-Mount Laurel aquifer—070478, which in 2014–19 departed from previous periods, and 330020, which continued a gradual decrease throughout its period of record.

New Jersey↗
USGS2024

A history of cryohydrogeology modeling and recent advancements through the integration of solute transport

Groundwater flow systems and permafrost are interrelated because permafrost thaw enhances permeability, while groundwater flow can advect heat and accelerate permafrost thaw (McKenzie et al. 2021). Given amplified climate change in cold regions, there is renewed interest in ‘cryohydrogeology’, the study of groundwater in cold regions. Many data-driven studies have shown that permafrost thaw is leading to activated aquifers and increased baseflow across the pan-Arctic region (e.g. Walvoord and Striegl 2007, Evans et al. 2020). Empirical evidence of a subsurface ‘replumbing’ (Walvoord and Kurylyk 2016) in permafrost regions raises questions about the fate of sequestered contaminants in the North (Langer et al. 2023). We will discuss the history of and emerging opportunities in cryohydrogeological modeling, with a focus on recent contaminant transport modeling.

Conference Paper↗
USGS2026

Special issue “The 2024 M7.6 Noto Peninsula earthquake and seismic swarm”

The 2024 Mw7.5 (M JMA 7.6) Noto Peninsula earthquake struck the Noto Peninsula, central Japan on January 1, 2024, causing widespread damage. This Mw7.5 mainshock, together with the prolonged seismic swarm that persisted in the northeastern Noto Peninsula since the end of 2020, of which the largest event was the 2023 M JMA 6.5 earthquake, constitutes one of the most significant seismic sequences in Japan in recent decades. The sequence provides a rare opportunity to investigate the complex interplay among swarm activity, fluid migration, fault geometry and distribution, and the preparatory conditions required for a large inland crustal earthquake. Motivated by this scientific context, this special issue brings together studies examining the sequence from diverse disciplinary perspectives, incorporating seismic, geodetic, and tsunami-related observations and modelling.

Noto Peninsula↗
USGS2025

The effects of line simplification on planform geometry

Data on maps should retain accuracy regardless of scale. Yet, as cartographic lines are generalized, there can be impacts on properties such as topology, density, and planform geometry. Here, we investigate the use of the Scale Specific Sinuosity (S3) metric (Stanislawski et al., 2023) to evaluate the effects of line simplification on planform geometry, which is the bends of streams in map view. We employ an open-source Python S3 workflow to characterize the geometry of five diverse stream channels in the United States. The original data are extracted from the U.S. Geological Survey National Hydrography Dataset 1:24,000-scale vector data (U.S. Geological Survey, 2000) (Table 1), and the simplification is done using the Visvalingam and Whyatt method (2017) with a simplification tolerance of 0.5, 1.0, 1.5, and 2.0 km. The S3 analysis is calculated at each level of simplification and S3 derivatives are generated. Derivatives include measures of sinuosity, fractal dimension, and the dominant bend wavelength. The findings show that the change in planform geometry is scale-dependent, though simplification will have little effect on straighter lines. The change becomes more apparent in complex lines as the degree of simplification aligns with the scale of the dominant bend geometries. These logical conclusions are evidence that the S3 is a useful metric for automated characterization of bend geometry regardless of line complexity.

Conference Paper↗
USGS2025

Analysis of the potential effects of Uzbekistan’s mineral endowment on the critical mineral supply of tungsten

Tungsten appears on the 2018 and 2022 U.S. Geological Survey critical mineral lists in part because of a very high global production concentration in China, which produces almost 83 percent of the world’s mined tungsten. Using known parameters and values from other tungsten mining operations, we created hypothetical scenarios in which three tungsten deposits in Uzbekistan are considered for development. Our results show that all three deposits are likely to be economically viable to develop under 2024 market conditions. If the three studied tungsten deposits were put into production, Uzbekistan could become the third-leading tungsten-producing country in the world and increase world output of tungsten by 2.7 percent. Putting these tungsten deposits in Uzbekistan into production could slightly reduce the tungsten global market concentration, therefore reducing the supply disruption potential for tungsten.

Open-File Report↗
USGS2025

AAPG Energy and Minerals Division Tight Oil and Gas Committee Activities and Commodity Report for 2021-2022: Tuscaloosa Marine Shale, Gulf Coast basin, Louisiana and Mississippi

The Upper Cretaceous Tuscaloosa marine shale (TMS) potential production area encompasses 20.4 million acres across central Louisiana (LA), southern Mississippi (MS), southwestern Alabama (AL), and a small southwestern section of the Florida panhandle (Hackley et al., 2018). It remains a minor and largely undeveloped unconventional shale oil play with production from the TMS confined along the east-west LA-MS State boundary (Fig. 1). The mean undiscovered, technically recoverable resources in the TMS are estimated at 1.5 billion barrels of oil and 4.6 trillion cubic feet of gas (Hackley et al., 2018). Geochemical analyses of source rock solvent extracts and oil samples indicate that, in the play area, the TMS is the primary source of shale oil produced from the TMS as well as of accumulated oils in the underlying conventional reservoirs of the lower Tuscaloosa (Hackley et al., 2020). The API gravity of TMS oils ranges from approximately 34 to 46 degrees (Hackley et al., 2020; Croke et al., 2020).

AAPG Wiki↗
USGS2025

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington↗