Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Computational Physics”

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,659 records · Page 13Linked to original sources
USGS1989

Techniques for simulating flood hydrographs and estimating flood volumes for ungaged basins in east and west Tennessee

A dimensionless hydrograph developed for a variety of basin conditions in Georgia was tested for its applicability to streams in East and West Tennessee by comparing it to a similar dimensionless hydrograph developed for streams in East and West Tennessee. Hydrographs of observed discharge at 83 streams in East Tennessee and 38 in West Tennessee were used in the study. Statistical analyses were performed by comparing simulated (or computed) hydrographs, derived by application of the Georgia dimensionless hydrograph, and dimensionless hydrographs developed from Tennessee data, with the observed hydrographs at 50 and 75% of their peak-flow widths. Results of the tests indicate that the Georgia dimensionless hydrography is virtually the same as the one developed for streams in East Tennessee, but that it is different from the dimensionless hydrograph developed for streams in West Tennessee. Because of the extensive testing of the Georgia dimensionless hydrograph, it was determined to be applicable for East Tennessee, whereas the dimensionless hydrograph developed from data on streams in West Tennessee was determined to be applicable in West Tennessee. As part of the dimensionless hydrograph development, an average lagtime in hours for each study basin, and the volume in inches of flood runoff for each flood event were computed. By use of multiple-regression analysis, equations were developed that relate basin lagtime to drainage area size, basin length, and percent impervious area. Similarly, flood volumes were related to drainage area size, peak discharge, and basin lagtime. These equations, along with the appropriate dimensionless hydrograph, can be used to estimate a typical (average) flood hydrograph and volume for recurrence-intervals up to 100 years at any ungaged site draining less than 50 sq mi in East and West Tennessee.

Tennessee↗
USGS2026

Effects of a low-lipid diet on the gut microbiome and head kidney transcriptome of juvenile Chinook Salmon

Objective Pacific salmon Oncorhynchus spp. reared in production hatcheries are typically fed high-lipid, energy-dense diets to achieve large size and high body condition prior to release. In contrast, juveniles in natural environments tend to consume low-lipid, high-protein diets, and fish reared for research or conservation purposes are sometimes fed diets that are formulated to mimic natural diets and promote wild-like phenotypes. Understanding how these alternative diets affect fish health beyond growth and body condition could ultimately contribute to improving hatchery fish fitness. Methods In this work, we evaluated changes in the fecal microbiome and gene expression of juvenile Chinook Salmon O. tshawytscha on a standard high-lipid hatchery diet versus a low-lipid diet formulated to mimic the nutrition profile of natural-origin fish. To evaluate the time scale at which diet alters the fecal microbiome, we collected longitudinal samples over a 12-week period and switched the diets of a subset of fish twice during the experiment. We used 16S ribosomal RNA gene amplicon sequencing to characterize fecal microbiome differences between fish on the two diets as well as hatchery-reared fish at a production hatchery, hatchery fish that had been captured after release into a stream, and natural-origin, stream-reared fish of similar ages. Additionally, we conducted RNA sequencing on head kidney samples from laboratory-reared fish to evaluate changes in gene expression in this important immune organ. Results We found that the low-lipid diet and the hatchery diet resulted in microbiomes that differed from the microbiome of natural-origin fish and from each other and that diet-driven changes to the microbiome could occur in under 14 d. The low-lipid diet did not result in a microbiome that resembled the microbiome of naturally produced fish. Instead, the low-lipid diet resulted in a microbiome community that was distinct from those of fish reared on the hatchery diet and fish sampled from the wild. The RNA sequencing results indicated differential enrichment of pathways related to immunity, metabolism, and hormone synthesis between fish that were fed the two experimental diets. Conclusions The results suggest that additional environmental factors influence the microbiome more strongly than diet formulation or that the low-lipid diet has a smaller effect on the microbiome than a natural, ­invertebrate-based diet. Given that the gut microbiome and systemic immune function contribute significantly to disease resistance, our findings highlight the importance of understanding how diets fed to fish in captivity may affect fish health beyond growth and body condition metrics.

Journal of Aquatic Animal Health↗
USGS2025

Flood-inundation maps for 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, 2023

Digital flood-inundation map libraries for two reaches that constitute 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, were created by the U.S. Geological Survey (USGS) in cooperation with the Papio-Missouri River Natural Resource District. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://www.usgs.gov/mission-areas/water-resources/science/flood-inundation-mapping-fim-program , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at Little Papillion Creek at Irvington, Nebr. (USGS station 06610750), Little Papillion Creek at Ak-Sar-Ben at Omaha, Nebr. (USGS station 06610765), and Big Papillion Creek at Q Street at Omaha, Nebr. (USGS station 06610770) streamgages. Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System database at https://doi.org/10.5066/F7P55KJN or from the National Weather Service Advanced Hydrologic Prediction Service at https://water.weather.gov/ahps/ . Flood profiles were computed for two different reaches that constitute 14.8 miles of stream length in the study area by using hydraulic models. The models were calibrated by adjusting roughness coefficients to best represent the current (2022) stage-streamflow relation at the streamgages within the study reach. The hydraulic models were then used to compute water-surface profiles at 1-foot stage intervals for selected stage ranges to represent various flooding scenarios at the streamgages in each reach. The simulated water-surface profiles then were combined with a digital elevation model using a geographic information system, which had a 10-foot grid spacing to delineate the flooding extents and water depths for each stage. The availability of these flood-inundation maps, along with information regarding current stage from the USGS streamgages, can provide emergency management personnel and residents with information that is critical for flood response activities and post flood recovery efforts.

Nebraska↗
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

Using the D-Claw software package to model lahars in the Middle Fork Nooksack River drainage and beyond, Mount Baker, Washington

Lahars, or volcanic mudflows, are the most hazardous eruption-related phenomena that will affect communities living along rivers that originate on Mount Baker. In the past 15,000 years, the largest lahars from Mount Baker have affected the Middle Fork Nooksack River drainage and beyond. Here we use the physics-based D-Claw software package to model nine lahar scenarios that are initiated as water-saturated landslides between Sherman Crater and the Roman Wall on the Mount Baker edifice and flow down the Middle Fork Nooksack River. The scenarios range in volume from 1 to 260 million cubic meters and have an initial hydraulic permeability from 10 −12 to 10 −10 meters squared. Model output includes data such as flow depth, velocity, runout distance, area inundated, arrival time, and sediment concentration as well as information that allows scientists to calculate other important hydrologic characteristics such as lahar discharge. These data are important to officials who have the responsibility to plan for, or take mitigation measures against, future Mount Baker lahars. To check the validity of the D-Claw results, we compare the scenarios to known geologic information. We also compare D-Claw results with empirical models that have been used in the past to determine potential inundation areas, runout distances, and arrival times. These comparisons highlight similarities and differences between empirical and physics-based models. We also present D-Claw scenario-based animations to help scientists, officials, and lay people alike to visualize how future lahars could affect communities.

Washington↗
USGS2025

Lessons in business recovery following the 2023 Kahramanmaraş earthquake sequence, Türkiye informed by women entrepreneurs

On 6 February 2023, Southern Türkiye was hit by devastating earthquakes, directly affecting over 14 million people in 11 cities, causing more than 50,000 deaths and the destruction of more than 800,000 buildings. This article goes beyond the physical damage imposed by the catastrophe to discuss the effects of the earthquakes on the operations of women-owned businesses. The mixed-method study with entrepreneurs belonging to a women’s business association operating in a moderately disrupted part of the region explores their struggles and recovery expectations. Thirty-five questionnaires were analyzed to identify the reasons for business closure, challenges, and needs faced in the post-disaster period and their recovery strategies. In addition, 23 entrepreneurs participated in roundtable discussions to provide a broader context to their responses to survey topics as well as lessons learned. Across both the survey and roundtables, while many respondents reported minor physical damage to their building, they also experienced financial and personal challenges from disruption to equipment, infrastructure, services, supply chains, institutional decisions, employee well-being, and customer base. Many used their business resources and personal savings to assist employees and others in the community. The women entrepreneurs often felt their recovery needs were ignored by government and private relief organizations and encountered barriers to receiving assistance from public and private institutions. Organizing together as women in business, even informally, provided mutual support during the crisis and recovery periods and catalyzed their role in support of their communities. The results illuminate functional community recovery as a balance of recovery of built infrastructure functionality and recovery of the broader social and economic fabric of the community.

southern Turkey↗
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↗
USGS2024

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗
USGS2026

Metallurgical coal—Deposits, production, resources, market dynamics, and supply chain risks

Plain Language Summary Metallurgical coal (met coal; consumed to produce coke for steelmaking) must meet specific chemical and physical specifications. In 2023, the conterminous United States produced 66 million short tons (mst) of met coal, consumed 15.85 mst domestically, exported 51.1 mst, and imported 0.7 mst. Most met coal was produced in the Appalachian Basin, but there are also resources that meet the specifications for met coal in the Western United States.

conterminous United States↗
USGS2026

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

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

Puerto Rico, U.S. Virgin Islands↗
USGS2025

Seismic moment and local magnitude scales in Ridgecrest, CA from the SCEC/USGS Community Stress Drop Validation Study

We illustrate the systematic difference between moment magnitude and local magnitude caused by underlying earthquake source physics, using seismic moments submitted to the Statewide California Earthquake Center/United States Geological Survey Community Stress Drop Validation Study 2019 Ridgecrest data set. While the relationship between seismic moment and moment magnitude ( M or M w ) of log 10 ( M 0 ) ~ 1.5* M is uniformly valid for all earthquake sizes by definition (Hanks and Kanamori, 1979), the relationship between local magnitude M L and moment is itself magnitude dependent. For moderate events, ~3< M < ~6, M and M L are coincident; for earthquakes smaller than ~3, M L ~ 1.0 log 10 M 0 (Hanks and Boore, 1984). This is a physical consequence of the corner frequency fc becoming larger than the upper frequency of observation and implies that M L and M differ systematically by a factor of 1.5 for these small events. While this idea is not new, we propose a new, continuous relationship between local magnitude and moment, for magnitudes 2 to 6 which extrapolates to smaller and larger magnitudes, applicable to southern California specific to the Ridgecrest region. We make use of the plethora of seismic moments as submitted by many participants of the Community Stress Drop study, compared to the Southern California Seismic Network (SCSN) catalog magnitudes. Overall, the seismic moments in the Community Study recover moment magnitude well, so we use our new M L - M 0 to convert M L to M , refining the SCSN operational M Lr scale. This systematic difference of 50% in slope between local and moment magnitude at small magnitudes has implications for spectral stress drop estimates, earthquake ground motion modeling, as well as other magnitude scales and earthquake occurrence statistics.

California↗
USGS1986

Water-quality appraisal of NASQAN stations below impoundments, eastern Tennessee

The National Stream Quality Accounting Network (NASQAN) is a network of stations at which systematic and continuing water-quality data are collected. Major objectives of this U.S. Geological Survey program are (1) to depict areal variability of streamflow and water-quality conditions nationwide on a year-by-year basis and (2) to detect long-term changes in streamflow and stream quality. Several NASQAN stations in East Tennessee are downstream from impoundments which have a significant effect on water quality. NASQAN data obtained from the Tennessee River below Watts Bar Dam and the Clinch River below Melton Hill Dam were compared to water-quality data from the basins upstream. The comparison indicates that NASQAN data obtained below impoundment may not be adequate to describe a composite picture of water quality in the accounting unit. Detention time of storage in the impoundments is believed to moderate the range of constituent values observed at the NASQAN stations. Data obtained upstream and downstream from Watts Bar Dam indicate that the water sampled at the NASQAN station comes from stratified layers of the impoundment and is not representative of an integrated sample of water from the impoundment. Values of total recoverable iron suggest that, because of adsorption to sediments in impoundments, some constituents are not accurately described by sampling below impoundments. Relations between water-quality constituents and flow at stations on the Clinch River and Tennessee River are not well defined due to regulation. Direct load computations for many constituents were therefore not possible, which diminished the utility of data from these NASQAN stations to account for quantity versus quality of the water. Load computations were only possible for ionic constituents through use of a continuous specific-conductance record as an intermediary. Compensation for the effects of discharge prior to application of the Seasonal Kendall test for trends could not be done and identification of trends in water-quality constituents caused by some process (source) change was not possible. Some water-quality trends indicated by data from the Clinch and Tennessee Rivers might reflect the decreasing trend in discharge during the 1972-82 water years. Thus the stations below Watts Bar Dam and below Melton Hill Dam do not adequately meet the NASQAN objective to detect and assess long-term changes in stream quality.

Tennessee↗
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↗
USGS2026

Estimating aftershock risk for entry into earthquake-damaged buildings

We present a simple method to estimate the risk of experiencing strong shaking from aftershocks during entry into earthquake-damaged buildings. We compute wait times until the probability of strong ground shaking from aftershocks reaches a predefined risk threshold; for example, a 0.4 percent probability of experiencing Modified Mercalli Intensity 7 or greater shaking during the planned building entry. We also develop a relation between aftershock probability and the probability of strong shaking, so that users can reference the U.S. Geological Survey aftershock forecast during an ongoing aftershock sequence to determine if the risk threshold has been met. We apply our method to active continental regions (for example, the Western United States), stable continental regions (for example, the Central and Eastern United States), and subduction zones (for example, Cascadia or Alaska).

Open-File Report↗
USGS2024

StreamStats—A quarter century of delivering web-based geospatial and hydrologic information to the public, and lessons learned

StreamStats is a U.S. Geological Survey (USGS) web application that provides streamflow statistics, such as the 1-percent annual exceedance probability peak flow, the mean flow, and the 7-day, 10-year low flow, to the public through a map-based user interface. These statistics are used in many ways, such as in the design of roads, bridges, and other structures; in delineation of floodplains for land-use zoning and setting of insurance rates; for regulatory purposes, such as the permitting of wastewater discharges; and for hydrologic and climate change studies. StreamStats was first developed for Massachusetts and released in 2001. The application provided users with the ability to obtain streamflow statistics computed from data collected at USGS streamgages and to obtain estimates of streamflow statistics for user-selected ungaged sites. Massachusetts StreamStats used geographic information system software and digital mapping to compute drainage-basin characteristics, which were then used in statistical models to estimate streamflow statistics for the user-selected sites. The statistical models were in the form of equations that were developed through a process known as regression analysis. StreamStats was the first known web application with the ability to do interactive geoprocessing. The utility of Massachusetts StreamStats was instantly apparent, leading the USGS to develop a version of StreamStats that could be implemented nationally. USGS State offices normally were required to develop custom regression equations and prepare local digital mapping data needed for implementing StreamStats for their States. Funding needed to complete this work usually was provided through cooperative agreements between the USGS and State agencies. In 2004, Idaho became the first to be released in the national version of StreamStats. By 2023, 44 States were fully implemented and six were undergoing implementation. StreamStats has undergone many modifications over the years to keep up with changes to the underlying software and to add functionality. Customized functionality and separate linked StreamStats applications were developed for several States. Meeting the high demand for additions and improvements to StreamStats while also adhering to budgetary constraints has, at times, been challenging. The StreamStats development team has identified numerous additional improvements that could be made to provide better performance and more functionality. The lessons learned from the experience of building and operating StreamStats for nearly 25 years could be relevant to others interested in pursuing efforts of a similar scale.

Circular↗
USGS2026

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗
USGS2024

Insights on gas hydrate formation and growth within an interbedded sand reservoir from well logging at the Qiongdongnan Basin, South China Sea

Although variable well log resolution and its control on saturation estimation has been studied, it has not been directly applied to a specific location to explore the nature of gas hydrate within a sand reservoir. We applied in-situ measurements of resistivities, neutron porosity, and gamma ray at two sites in the Qiongdongnan Basin, South China Sea (QDN-W05–2021 and QDN-W08–2021) to investigate the reservoir parameters of a hydrate-bearing sand reservoir. Our results show that gas hydrate is distributed in 5 zones with a total thickness of 10.7 m and an average saturation of 69% at the QDN-W05–2021 site, while they are distributed in 2 zones with a total thickness of 4.3 m and an average saturation of 49% at the QDN-W08–2021 site. We found that variances in saturations estimated from lateral-extra deep button (RX), phase shift (P40H-P40L), and attenuation (A40H-A40L) resistivities within the laterally mapped continuous sand body were affected by the nature of gas hydrate occurrences. Results indicate gas hydrate forms and accumulates at the center of the sand layer and tends to be less or not present toward the top and base. Integrated with seismic data, the in-situ measurements provide insights in the evolution of a mushroom-shaped, hydrate-gas reservoir system. In the system, free gas is likely horizontally transported from the top-center of the gas chimney to the surrounding areas in the early stage dominated by a warm-gas environment, whereas hydrate forms in the opposite pathway starting from the surrounding areas in the following stage with temperature reducing. Our study suggests that high-resolution in-situ measurements not only are a tool to identify the physical properties, but also can be used to help explain the physical process of hydrate growth and accumulation.

Marine Geology↗
USGS2024

Assessment of fresh groundwater discharge and saline surface-water intrusion at Operable Unit 2, North Chevalier Field Disposal Area (Site 11), Naval Air Station Pensacola, Florida, 2018–22

Site 11 is a former landfill at North Chevalier Field Disposal Area in Operable Unit 2 at Naval Air Station Pensacola, in northwest Florida. Site 11 is adjacent to Bayou Grande, a shallow, tidally influenced, saline estuary of the Pensacola Bay watershed. Federal and Florida regulators have expressed concern that contaminants detected in groundwater beneath the inland parts of Site 11 may discharge to Bayou Grande. In 2017, the Department of Defense, U.S. Navy, Naval Facilities Engineering Systems Command Southeast asked the U.S. Geological Survey to assess the occurrence of fresh groundwater discharge to Bayou Grande at Site 11 and to delineate to the extent practicable the location of groundwater discharge. Between 2018 and 2022, the U.S. Geological Survey used a multiple-lines-of-evidence approach that included a visual method and three physical methods based on the temperature difference between groundwater and surface water to assess groundwater discharge. One of the physically based methods also used the difference in specific conductance between fresh groundwater and brackish to saline surface water. Combined, the data indicate that fresh groundwater from across Site 11 discharges primarily along the shoreline of the northern and northeastern part of Site 11. The data also indicate that saline surface water from Bayou Grande intrudes tens of feet into the shallow aquifer beneath Site 11. The combined data indicate that the interface between fresh groundwater and saline surface water changes over space and time. Any new monitoring wells proposed for installation near the shoreline of Site 11 should include approaches to monitor the changes in the location of the freshwater/saltwater interface. Care would need to be taken to collect any groundwater samples at the correct season and tidal period to provide the highest probability of collecting a representative sample of Site 11 groundwater unaffected by saltwater.

Florida↗