Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geoenergy Science and Engineering”

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

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗
USGS2025

Simulation of the impacts of projected climate change on groundwater resources in the urban, semiarid Yucaipa Valley watershed, southern California using an integrated hydrologic model

Managing water resources in semiarid watersheds is challenging due to limited supply and uncertain future climate conditions. This paper examines the impact of future climate changes on an urban watershed in southern California using an integrated hydrologic model. GSFLOW modeling software is used to simulate the nonlinear relationships between climate trends and precipitation partitioning into ET, runoff, and subsurface storage. Four global circulation models (GCMs), each with two greenhouse-gas scenarios, RCP45 and RCP85 are used to project future climate conditions. GCMs include the CanESM2, CNRM-CM5, HadGEM2-ES, and MIROC5 models. The model's simulated hydrologic conditions are compared with historical data to assess changes in water budgets and groundwater supply. Results indicate decreased groundwater storage in most scenarios due to increased natural evapotranspiration, vegetation consumptive use, and streamflow out of the watershed. Only scenarios with substantially increased future precipitation show increased groundwater storage. The study also highlights increased future aridity despite the rise in precipitation and large precipitation events forecast by GCMs, which increase the risk of urban floods and decrease stream leakage and water available to vegetation.

California↗
USGS2025

Gaps in water quality modeling of hydrologic systems

This review assesses gaps in water quality modeling, emphasizing opportunities to improve next-generation models that are essential for managing water quality and are integral to meeting goals of scientific and management agencies. In particular, this paper identifies gaps in water quality modeling capabilities that, if addressed, could support assessments, projections, and evaluations of management alternatives to support ecosystem health and human beneficial use of water resources. It covers surface water and groundwater quality modeling, dealing with a broad suite of physical, biogeochemical, and anthropogenic drivers. Modeling capabilities for six constituents (or constituent categories) are explored: water temperature, salinity, nutrients, sediment, geogenic constituents, and contaminants of emerging concern. Each constituent was followed through the coupled atmospheric-hydrologic-human system, with prominent modeling gaps described for a diverse array of relevant inputs, processes, and human activities. Commonly identified modeling gaps primarily fall under three types: (1) model gaps, (2) data gaps, and (3) process understanding gaps. In addition to potential solutions for addressing specific individual modeling limitations, some broad approaches (e.g., enhanced data collection and compilation, machine learning, reduced-complexity modeling) are discussed as ways forward for tackling multiple gaps. This gap analysis establishes a framework of diverse approaches that may support improved process representation, scale, and accuracy of models for a wide range of water quality issues.

Water↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2026

Site-specific amplifications in Northwestern Turkiye: A generic approach

In this study, we derive generic site amplification functions applicable to Northwestern Türkiye and follow guidance on site classifications established by the United States National Earthquake Hazards Reduction Program (NEHRP). We employ the one-dimensional (1-D) equivalent linear ground response analysis method and focus on recordings from a select suite of free-field strong motion (SM) stations. These particular SM stations are associated with a comprehensive geotechnical site conditions dataset comprising 76 shear-wave velocity profiles. Despite their varying resolutions and depths, this diversity allows for a more representative analysis across different site conditions, enhancing the robustness of our findings. For analyses, we utilize a robust set of global ground motion records, selected from both within and outside of Türkiye, to serve as strong and weak input bedrock motions. Site amplification is determined by propagating these bedrock records through 76 site-specific Afet ve Acil Durum Yönetimi Başkanlığı (AFAD) station profiles and calculating the spectral acceleration ratio between the ground surface and the input motion. Our analyses reveal substantial disparities between the amplification functions derived from 1-D site response analyses using strong and weak input motions. These disparities are central to our analysis and inform the subsequent comparison with site amplification functions reported in the literature. We compare the generic site amplification functions derived in this study with those reported in the literature for regions with similar geological and seismic conditions to Northwestern Türkiye. For example, site amplification functions for NEHRP Site Classes C and D have been documented for Greece and the United States utilizing frequency-dependent methods such as the quarter-wavelength approximation. Our findings reveal that the generic site amplification functions derived in this study exhibit differences in peak amplitudes and predominant frequencies compared to those in the literature, reflecting regional variations in geotechnical and seismic properties. This study represents an initial effort to formulate generic site amplification functions for Türkiye, with a particular emphasis on the Northwestern region. A thorough literature review indicated no prior studies have proposed such functions specifically for this area.

Northwestern Turkiye↗
USGS2026

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California↗
USGS2025

Paleoproterozoic vein graphite mineralization caused by decarbonation in the Ruby Range, Montana, USA

Hydrothermal graphite veins are a possible source for modern battery materials and require better understanding of their carbon source(s) and absolute timing to develop mapable criteria for exploration models. We present new observations of graphite vein and alteration paragenesis and U-Pb LA-ICP-MS titanite age data from the Ruby prospect, Montana, USA, that constrain mineralization timing and source. The graphite veins cut high-temperature metamorphic rocks of the lower Christensen Range suite and are associated with intense diopside (Di0.69Hd0.27Jo0.04) alteration of marble. The oldest titanite ages in fresh marble and partially altered calc-silicate gneiss are ca. 2,500 – 2,450 Ma, show elevated REE values, and Eu/Eu* <1.5, consistent with growth during early regional metamorphism. Titanite in diopside-graphite alteration, interpreted as part of the hydrothermal vein-forming episode, cluster around 1,750 Ma, are characterized by lower REE values, and show Eu/Eu* >2; suggesting a low fO 2 fluid generated from metamorphism during the Big Sky orogeny. Our paragenetic observations and titanite ages indicate graphite vein formation via skarnoid decarbonation reactions in marble late in regional orogenesis. Granulite-facies, carbonate-bearing, supracrustal rock terranes are thus favorable for hydrothermal graphite vein deposits.

Montana↗
USGS2026

Increasing artificial light at night enhances salmonid predator foraging efficiency in an urbanized lake

Artificial light at night (ALAN) poses a threat to ecosystems globally. It includes both direct and indirect light, or skyglow, which occurs when ALAN scatters in the atmosphere, extending beyond its original source. We analyzed ALAN trends in Lake Washington, WA, from 2014 to 2023 using Visible Infrared Imaging Radiometer Suite (VIIRS) nighttime light measurements, evaluated the relationship between in situ and satellite measurements, and modeled predator search volumes for a juvenile salmon predator, cutthroat trout ( Oncorhynchus clarki ), as a function of ambient light and turbidity conditions. Open water regions experienced significant increases in ALAN, while nearshore areas primarily showed no or negative trends, revealing the role of skyglow in shaping open water light environments. Using a visual foraging model for cutthroat trout, we found that juvenile salmon at the shallow southern pelagic site experienced light 28 times brighter, resulting in a 168% greater predation vulnerability than those at the northern site. In the Ship Canal, a narrow corridor for outmigrating salmon, predator search volumes were 249% higher than at the southern site. These contrasts in predation vulnerability demonstrate how local conditions influence predator–prey dynamics and provide critical insight for targeting mitigation of both nearshore and distant light sources.

Washington↗
USGS2026

Latest Pleistocene to 19th-century earthquakes on bending-moment reverse faults of the Seattle fault zone, Washington

Fault-related folds and their associated secondary faults play a critical yet often underrecognized role in accommodating strain and generating earthquakes in active fold-and-thrust belts. In the Seattle fault zone (SFZ), Washington, USA, we present new paleoseismic, geomorphic, and geophysical evidence for late Pleistocene and Holocene earthquakes on shallow, south-dipping secondary faults—the Lytle Beach and Vasa Park faults—that lie within the hanging wall of the greater SFZ and are on trend with the primary, blind Blakely Harbor fault. Our data show that these structures have ruptured independently, producing localized uplift and deformation at the surface, with the most recent event (RH2) likely occurring in the early nineteenth century. While a temporal overlap between the late Pleistocene RH1 and VP1 earthquakes raises the possibility of a ≥35 km rupture along the Blakely Harbor fault, structural and temporal evidence instead supports independent rupture on individual faults related to folding. We interpret these faults as bending-moment reverse faults that formed within a synclinal hinge zone of the main fault, reflecting mechanical and kinematic influences of the broader fault system. Combined with prior studies, our findings indicate that faulting related to folding dominates the mode of strain release within the SFZ since the late Pleistocene with more frequent earthquake recurrence (∼350 yr) over the past ∼2500 yr.

Washington↗
USGS2026

Long-term intermittent connection between the western Snake River Plain and Columbia basin: A two-phased incision history of Hells Canyon

For more than a century, researchers have debated the route of the Snake River across the northern Cordillera and U.S. Pacific Northwest, including the associated incision history of Hells Canyon. Here, we use detrital zircon U-Pb provenance analysis of Miocene strata upstream and downstream of Hells Canyon to constrain the evolution in fluvial pathways over time. Downstream of Hells Canyon, we show that the ca. 10−6 Ma Clarkston Heights gravel was dominantly sourced from the nearby Salmon and Clearwater rivers but also with substantial input from drainages that now flow into the western Snake River Plain. These same western Snake River Plain drainages were contributing to time-equivalent Lake Idaho strata, upstream of Hells Canyon, indicating a Hells Canyon fluvial connection between the western Snake River Plain and Columbia basin before 6 Ma. Results from Lake Idaho strata suggest a two-phased lacustrine history in the western Snake River Plain: Before ca. 6.0 Ma, lake strata were derived from local tributaries, suggesting that the western Snake River Plain was isolated from the modern upper Snake River system. Between ca. 4.3 Ma and 2.2 Ma, the source area expanded to include eastern tributaries in conjunction with migration of high-standing topography of the Yellowstone hotspot. Our results challenge the long-held hypothesis of “capture” of the Snake River through Hells Canyon. Instead, we infer that Hells Canyon was a long-established route for outflow of lakes occupying the western Snake River Plain, possibly intermittently, during a ca. 10−2 Ma phase of slow regional incision, followed by rapid incision in Hells Canyon starting ca. 2 Ma in conjunction with erosion of the topographic barrier impounding Lake Idaho.

Idaho, Oregon↗
USGS2025

A time-space model of graphite mineral systems

Increasing demand for graphite in energy storage systems warrants review of graphite ore genesis in a mineral systems framework. Orogenic graphite encompasses the metamorphic and orogenic mineral systems that produce flake graphite and hydrothermal vein (lump and chip) graphite deposits, respectively. A common feature of orogenic graphite deposits is an association with upper amphibolite- to granulite-facies metasedimentary rocks in continent-continent or continent-island arc collisional orogens. Orogenic flake graphite deposits form primarily through graphitization of organic carbon during regional metamorphism, but strain localization and partial melting of pelitic protoliths are likely important processes for graphite grade and quality enrichment. Orogenic vein graphite deposits precipitate from hydrothermal fluids, possibly derived from metamorphism or anatexis at depth. Decarbonation reactions in mixed calcareous-carbonaceous metasedimentary sequences are the most likely carbon sources for the veins. In contrast, intrusion-related graphite includes magmatic-hydrothermal and metamorphic mineral systems that form primarily in continental arc settings via the interaction of magmas with carbonaceous sedimentary packages. Magmatic-hydrothermal flake graphite deposits are hosted in plutonic and volcanic rocks, and result from the exsolution of CO 2 -CH 4 -rich fluids from melts contaminated by such packages. Contact metamorphism of carbonaceous sedimentary rocks by plutons produces some microcrystalline (amorphous) graphite deposits, including many in China. Compilation of geologic data from known graphite deposits globally suggests that pulses of carbon deposition in the Paleoproterozoic, Mesoproterozoic, and Neoproterozoic provided source carbon material. Subsequent supercontinent orogenesis at ca. 2,100 to 1,700 Ma (Columbia), ca. 1,300 to 1,000 Ma (Rodinia), and ca. 650 to 500 Ma (Gondwana) resulted in the genesis of orogenic flake and vein graphite deposits, where favorable geologic components overlapped with organic ± carbonate carbon-rich strata. Cryogenian deposition of graphite protoliths and Cryogenian – Cambrian metamorphic mineralization account for nearly 75% of all known resources globally and coincide with profound carbon isotope excursions and climate variability, implying a link with the global carbon budget. Comparatively few graphite deposits are associated with Pangea-forming orogens, attributed to less exhumation and/or denudation. High-temperature metasedimentary belts containing organic carbon-rich protoliths are most favorable for hosting orogenic flake graphite deposits, whereas sequences that also contain carbonate protoliths are favorable for orogenic graphite veins. Continent-scale orogenic belts may host both deposit types along with vanadium deposits. Use of a time-space mineral systems framework for graphite deposits can improve exploration models needed to ensure future supply of this critical mineral and provide insights into Earth’s long-term carbon cycle.

Mineralium Deposita↗
USGS2025

Seamless Geologic Map Database for the Intermountain West, United States: A foundational dataset for mineral systems analysis

The Intermountain West has a complex geologic history, resulting in the formation of a diverse array of mineral deposits. Effective mineral exploration requires understanding the spatial and temporal relationships among geologic processes and events, a key focus of the mineral systems approach to exploration. This paper presents the Intermountain West Seamless Geologic Map Database, a unified dataset designed to support mineral exploration. Integrating geologic provinces, structural settings, and hydrothermal alteration, the database leverages the Seamless Integrated Geologic Mapping (SIGMa) extension to the USGS Geologic Map Schema (GeMS) to standardize geologic data from varied sources. SIGMa's hierarchical stratigraphic organization and feature-level metadata enhance data interoperability and reusability, enabling seamless query, analysis, and visualization of lithology, structural features, mineral deposits, geochronology, hydrothermal alteration. volcanic activity, and By providing a regionally consistent and dynamically evolving geologic map, this database provides a foundational framework for mineral exploration and geologic research. It also allows for an efficient workflow that expedites the publication of integrated geologic map databases.

Conference Paper↗
USGS2025

Mineralization processes at Escanaba Trough

Escanaba Trough off northern California is the southernmost and only sediment covered segment of the Gorda Ridge. Hydrothermal mineralization was first discovered in Escanaba Trough in 1985, and subsequent investigations included the delineation of the NESCA (northern Escanaba) and SESCA (southern Escanaba) sections of Escanaba Trough, as well as drilling by Ocean Drilling Program (ODP) Leg 169 in 1996. Preliminary research conducted in the 1980s distinguished two classes of sulfide mineral assemblages within the Escanaba Trough: one dominated by pyrrhotite (Fe 1-x S) and one containing ‘mixed sulfides’ including chalcopyrite (CuFeS 2 ) and sphalerite (ZnS). In 2022, we conducted a research expedition to the Escanaba Trough to more thoroughly characterize mineralization processes and consider preservation of precipitated minerals throughout hydrothermally active and inactive sections of Escanaba Trough, with the specific objective of characterizing the host minerals, and processes of enrichment for specific elements which are currently listed as critical elements and can also be toxic under certain conditions (Sb, As, Ba, Co, Ga, Ge, Te). The following will include initial results from that work.

California↗
USGS2025

Fingerprinting magmatic REE deposit sources with zircon petrochronology

Carbonatites and associated alkaline silicate rocks are of considerable economic interest due to their enrichments in rare earth elements. The petrogenesis and source(s) of these complexes, however, are poorly understood. Models propose either mantle plume-derived carbon-rich melts or a mantle source enriched by subduction-related metasomatism. We use zircon trace elements to fingerprint carbonatite-alkaline complex sources, focusing on the economically significant 1.4 Ga Mountain Pass intrusive suite (MPIS). Autocrystic zircon from MPIS alkaline silicate rocks are enriched in Th and U relative to Nb suggesting a subduction influenced, oxidized source region. Lower Sc/Yb, higher Ti concentrations, and the absence of Eu anomalies in autocrystic zircon suggest derivation from a less hydrous, hotter, and deeper mantle source relative to arc-related 1.8–1.6 Ga inherited zircon. These data are best explained by syn- to post-collisional Mesoproterozoic reactivation of a subduction-metasomatized mantle source. The source and tectonic setting of the MPIS contrasts with Mesoproterozoic Gifford Creek and Bayan Obo carbonatites which exhibit plume-like compositions and are associated with rifts, suggesting different sources and tectonic settings for these economically significant deposits. We demonstrate that zircon petrochronology is a robust method for distinguishing carbonatite sources and can inform more targeted exploration strategies for critical mineral resources.

Conference Paper↗
USGS2024

Potential corrosivity of untreated groundwater in Louisiana

Corrosive groundwater can cause lead, copper, and other metals to leach from pipes and plumbing fixtures in water distribution systems. Metals, if ingested, could lead to serious health implications to the nearly 2.9 million people in Louisiana who obtain their drinking water from groundwater sources. Four indices—the Langelier Saturation Index (LSI), Ryznar Stability Index (RSI), Puckorius Scaling Index (PSI), and the Potential to Promote Galvanic Corrosion (PPGC)—in addition to an analysis which normalized the results from the existing indices, the Combined Index (CI), were used to assess the corrosivity of groundwater in Louisiana and identify areas within eight major aquifers and aquifer systems with moderate to high corrosivity potential. The purpose of this study is to provide State and local governments, public water system managers, and the nearly 500,000 private well owners in Louisiana with information needed to manage drinking-water supplies and mitigate potential health risks related to leaching of metals from water pipes and fixtures. The average scores of untreated groundwater samples from approximately 375 wells by index are as follows: LSI, −1.28; RSI, 9.78; PSI, 9.34; and CI, 4.14. The PPGC does not produce a numerical score, but the total percentage of class counts can be used to assign a classification; overall, samples in Louisiana were classified as significant concern. The percentages of groundwater samples from wells classified as potentially corrosive, by index, are as follows: LSI, 53 percent; RSI, 94 percent; PSI, 81 percent; PPGC, 98 percent; and CI, 81 percent. The percentages of samples classified as indeterminate, by index, are as follows: LSI, 46 percent; RSI, 5 percent; PSI, 12 percent; PPGC, 0 percent; and CI, 18 percent.

Louisiana↗
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↗
USGS2024

Recent applications of the USGS National Crustal Model for Seismic Hazard Studies

The U.S. Geological Survey is developing the National Crustal Model (NCM) for seismic hazard studies to facilitate modeling site, path, and source components of seismic hazard across the conterminous United States. The NCM is composed of a 1km grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a threedimensional (3D) geologic framework and geophysical rules that use (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters (VS30), the depths to 1.0 and 2.5 km/s shear-wave velocity (Z1.0 and Z2.5), and sediment thickness, can be computed from the NCM. As GMMs continue to improve in the future, other metrics could also be extracted or derived from the NCM, such as fundamental period, site attenuation (ko), a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis, including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

conterminous United States↗
USGS2025

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report↗