USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Planetary Science Journal”

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,665 records · Page 32Linked to original sources
USGS2000

Estimating the probability of elevated nitrate (NO2+NO3-N) concentrations in ground water in the Columbia Basin Ground Water Management Area, Washington

Logistic regression was used to relate anthropogenic (man-made) and natural factors to the occurrence of elevated concentrations of nitrite plus nitrate as nitrogen in ground water in the Columbia Basin Ground Water Management Area, eastern Washington. Variables that were analyzed included well depth, depth of well casing, ground-water recharge rates, presence of canals, fertilizer application amounts, soils, surficial geology, and land-use types. The variables that best explain the occurrence of nitrate concentrations above 3 milligrams per liter in wells were the amount of fertilizer applied annually within a 2-kilometer radius of a well and the depth of the well casing; the variables that best explain the occurrence of nitrate above 10 milligrams per liter included the amount of fertilizer applied annually within a 3-kilometer radius of a well, the depth of the well casing, and the mean soil hydrologic group, which is a measure of soil infiltration rate. Based on the relations between these variables and elevated nitrate concentrations, models were developed using logistic regression that predict the probability that ground water will exceed a nitrate concentration of either 3 milligrams per liter or 10 milligrams per liter. Maps were produced that illustrate the predicted probability that ground-water nitrate concentrations will exceed 3 milligrams per liter or 10 milligrams per liter for wells cased to 78 feet below land surface (median casing depth) and the predicted depth to which wells would need to be cased in order to have an 80-percent probability of drawing water with a nitrate concentration below either 3 milligrams per liter or 10 milligrams per liter. Maps showing the predicted probability for the occurrence of elevated nitrate concentrations indicate that the irrigated agricultural regions are most at risk. The predicted depths to which wells need to be cased in order to have an 80-percent chance of obtaining low nitrate ground water exceed 600 feet in the irrigated agricultural regions, whereas wells in dryland agricultural areas generally need a casing in excess of 400 feet. The predicted depth to which wells need to be cased to have at least an 80-percent chance to draw water with a nitrate concentration less than 10 milligrams per liter generally did not exceed 800 feet, with a 200-foot casing depth typical of the majority of the area.

Washington↗
USGS2026

Methodology for construction of a three-layer geologic model of the conterminous United States using land surface, top of bedrock, and top of basement

This report describes the methodology used for the construction of a digital three-layer geologic model of the conterminous United States by mapping the altitude of three surfaces: land surface, the top of bedrock, and the top of basement. These surfaces are mapped through the compilation and synthesis of published stratigraphic horizons from numerous topical studies. The mapped surfaces create a three-layer geologic model with three geomaterial-based subdivisions: unconsolidated to weakly consolidated sediment; layered consolidated rock strata that constitute bedrock; and crystalline rocks that are described as “basement,” consisting of either igneous, metamorphic, or highly deformed rocks. The data compilation and synthesis are highly dependent on the definition of the informal terms “bedrock” and “basement,” which may describe different ages or types of rock in different parts of the conterminous United States. This report presents the conceptualization of the three mapped layers, describes the datasets used, and summarizes the decisions made while compiling the three-layer model from the various sources. This digital dataset was created as part of efforts by the U.S. Geological Survey to develop subsurface geologic data in geospatial form as part of a broad directive to develop two-dimensional and three-dimensional geologic information at detailed, national, and continental scales. This digital dataset partly fulfills the goal of the U.S. Geological Survey’s National Cooperative Geologic Mapping Program to construct a national-scale three-dimensional geologic model.

conterminous United States↗
USGS2024

Mapping karst groundwater flow paths and delineating recharge areas for springs in the Little Sequatchie and Pryor Cove watersheds, Tennessee

The Little Sequatchie River and Pryor Cove Branch, in southern Tennessee, drain the eastern escarpment of the Cumberland Plateau to the Sequatchie River near the southern end of the Sequatchie Valley. The Little Sequatchie River is the largest tributary to the Sequatchie River by drainage area, covering over 120 square miles. The hydrology of the two drainage areas has been largely altered by karst processes, which has caused the majority of the streams to sink at the contact between the Mississippian Pennington Formation and the underlying Mississippian Bangor Limestone. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service and Tennessee Department of Environment and Conservation, initiated a study in 2021 to map the karst groundwater pathways in both watersheds in order to delineate recharge areas for several springs. One of these springs, Sequatchie Cave, represents a significant habitat for two Species of Greatest Conservation Need, the Glyphopsyche sequatchie (Sequatchie caddisfly) and the federally endangered Marstonia ogmorhaphe (royal marstonia). Springs and springflow-dominated streams in the Little Sequatchie River valley and Pryor Cove also provide water for agricultural practices and serve as a drinking water source for nearby communities. During the study, a total of 25 dye injections were conducted over eight rounds from January 2022 through March 2023. Dye traces from these injections helped to delineate recharge areas for six major springs, ranging from 7.3 to 65.2 square miles in area. The majority of the dye traces remained subsurface (from sinkpoint to recovery site) for long distances, with karst groundwater travelling nearly 8 miles before resurfacing. The dye traces also had rapid traveltimes, often travelling hundreds to thousands of feet per hour. The goal of this project was to provide scientific data related to karst groundwater pathways and spring recharge areas to aid State and Federal agencies in making informed decisions to protect and preserve this unique and vulnerable karst system.

Tennessee↗
USGS2026

Standardized method for logging drill core at the Idaho National Laboratory, Idaho

The U.S. Geological Survey’s (USGS) Lithologic Core Storage Library (CSL) at the Idaho National Laboratory stores more than 120,000 feet of drill core that is accessible to the public for research and sampling. To effectively convey the physical and descriptive properties of the drill core, USGS staff at the Idaho National Laboratory Project Office log the drill core and publish the lithologic logs as data releases. The logs provide essential data on the lithology, texture, mineralogy, alteration, and other physical properties of the core, which serve as valuable information for researchers to guide their research and sampling efforts. To ensure consistent, quality, and dependable lithologic logs, this document outlines the procedures and expectations for logging drill core at the CSL. This document describes the processes for storing, photographing, and logging core, and includes a variety of resources, reference materials, and appendixes designed to standardize and aid the logging process. Following the procedures outlined in this document will produce consistent, detailed logs that facilitate dependable observations and serve as an easy reference for researchers and other interested parties.

Idaho↗
USGS2025

The Hardscrabble Creek complex: A newly discovered, mostly buried, Mesoproterozoic mafic-ultramafic pluton in the Wet Mountains, Colorado, USA

The origin of prolific ca. 1.4 Ga ferroan magmatism between the southwestern USA and eastern Canada is enigmatic and has been explained by various models, including extensional, mantle plume, and convergent plate-margin models. Rare mafic plutons are associated with the ferroan plutons, which may help constrain their mantle source and tectonic setting. In the southwestern USA, only two such mafic plutons are known to exist. We present the first evidence for a third, mostly buried, potentially layered, mafic-ultramafic Mesoproterozoic pluton, informally referred to as the Hardscrabble Creek complex, in the central Wet Mountains of Colorado, USA. Recent geophysical data show an elliptical magnetic and gravity high spatially coincident with local gabbroic outcrops. New field and petrographic analyses of these exposed rocks reveal that they consist of ultramafic to mafic cumulates, including orthopyroxenite, olivine norite, norite, and anorthosite. High-precision U-Pb dating of zircon from orthopyroxenite and norite yield weighted mean 206 Pb/ 238 U dates of 1352.36 ± 1.60 Ma and 1352.37 ± 1.71 Ma, respectively. These dates indicate that the complex formed over a narrow timeframe, after the adjacent 1362 ± 7 Ma ferroan San Isabel Granite, and during the waning stages of the regional ca. 1.4 Ga ferroan magmatism. Whole-rock geochemistry and Nd-Sr-Pb isotope compositions of samples from the Hardscrabble Creek complex are similar to those of the San Isabel Granite, suggesting that they were derived from the same or a similar mantle source. The mineral chemistry of the samples is comparable to Proterozoic massif-type anorthosites and related mafic intrusions, indicating that the Hardscrabble Creek complex and San Isabel Granite together represent a rare anorthosite-mangerite-charnockite-granite (AMCG) suite in the southwestern USA. The Hardscrabble Creek complex is unique because it formed ~80 m.y. after the other few mafic plutons in the southwestern USA, and it contains an ultramafic section that is absent from these plutons and rare to the AMCG suite in general. A combination of arc-like whole-rock geochemistry, chondrite uniform reservoir-like Nd-Sr-Pb isotopes, and ocean island basalt (OIB)-like zircon trace element chemistry suggests that the complex was derived from a partial melt of OIB-like mantle and interacted with metasomatically enriched lithospheric mantle. The enriched lithospheric mantle signature, combined with the long ~160 m.y. duration of magmatism in the region, is consistent with a period of protracted convergent tectonism.

Colorado↗
USGS2024

Summary of annual repeat magnetotelluric surveys of the Geysers geothermal field

As part of a multi-year project, funded by the California Energy Commission, annual magnetotelluric (MT) surveys have been collected at The Geysers geothermal field in northern California with the goal of measuring temporal changes within the steam field. The repeat surveys started in 2021 and repeated a 2017 survey (Peacock et al., 2020) with further extension to the southern part of the geothermal field. Temporal variations in the MT transfer functions are observed to be spatially coherent and compartmentalized. Mapping residual phase tensor ellipses demonstrates the direction of maximum change is often aligned with existing fracture orientations. Three dimensional inversion of the MT data, using the inversion results from previous years as the starting model, indicates that the steam reservoir has generally become more resistive over time (~10%), suggesting more steam in the field. A few pockets within the steam field have become more electrically conductive over time and are collocated with injection wells, suggesting either more fluid content in those zones, less steam, or more saline fluids.

California↗
USGS2024

Oblique contraction along the fastest ocean-continent transform plate boundary focuses rock uplift west of the Fairweather fault, southeast Alaska

Contraction along the Yakutat–North America plate boundary drives 4.6–9.0 mm/year Holocene rock uplift rates along Earth's fastest slipping (≥49 mm/year) ocean–continent transform fault, the Fairweather Fault. Between Icy Point and Lituya Bay, the near-vertical Fairweather fault focuses rock uplift and rapid right-lateral slip by accommodating both vertical and fault-parallel strain during oblique-slip and separate, predominantly strike-slip ruptures. Unusually high uplift rates, indicated by radiocarbon and luminescence dating, result from a 10-km-wide, asymmetric, positive flower structure along a 20°, ∼30-km-long restraining double bend in the Fairweather fault. The principal reverse fault in the flower structure, the offshore, blind Icy Point–Lituya Bay fault, ruptures no more than every 460–1040 years evidenced by uplifted Holocene shorelines. Maximum 3–5 m coseismic uplifts imply 3.1–10 m dip slip per event and earthquake magnitudes of M w 7.0–7.5. The Yakutat block collides obliquely into North America, and our model entails oblique slip on the Fairweather fault with and without corupture on the reverse fault. Oblique slip is evident by vertically offset (>25 m) fluvial and marine terraces and by the primary Fairweather fault strand that strikes >20° to the west of plate-boundary motion.

Alaska↗
USGS1988

Geology of an area near Brentwood, Williamson County, Tennessee

The lithology and geologic structure of an area near Brentwood, Williamson County, Tennessee, were studied to define the potential aquifers and confining units of the ground-water flow system. Four formations were identified. They are, in descending order, the Bigby-Cannon Limestone, the Hermitage Formation, the Carters Limestone, and the Lebanon Limestone. The Bigby-Cannon Limestone and the Carters Limestone are potential aquifers. The Hermitage Formation and the Lebanon Limestone are confining units. The Bigby-Cannon Limestone and the Hermitage Formation have been affected by recent erosion. Any variation of the Carters Limestone is controlled by pre-Carters erosion of the top of the Lebanon Limestone. The thickness of the Carters Limestone ranges from 65 to 79 feet. Structurally, the area reflects the regional northwest dip. A small scale anti-cline-syncline pair also is evident. This feature is not a result of erosion and also occurs in the T-3 bentonite bed in the Carters Limestone.

Tennessee↗
USGS2025

Phase 1 technical implementation plan for the expansion of the ShakeAlert earthquake early warning system to Alaska

Executive Summary The conference report accompanying the fiscal year (FY) 2022 Consolidated Appropriations Act (Public Law 117–103) for the U.S. Department of the Interior and related agencies directed the U.S. Geological Survey (USGS) to “work with the State of Alaska to develop an implementation plan to be completed within two years in order to put ShakeAlert/Earthquake Early Warning in Alaska” (p. 29). Congress included $1 million in the FY 2022 appropriation to conduct this effort. The USGS Earthquake Hazards Program, along with partner organizations, has developed the ShakeAlert earthquake early warning (EEW) system for the West Coast, which currently operates in California, Oregon, and Washington. The purpose of the system and its alert delivery partners is to reduce the impact of earthquakes and save lives and property by delivering ShakeAlert-powered alerts that are transmitted to the public via mass notification technologies, and by providing more detailed data streams to institutional users and commercial service providers to trigger automated, user-specific, protective actions. ShakeAlert was designed in such a way that it could be expanded to other U.S. regions with high earthquake risk, after the build-out of seismic and geodetic networks to support ShakeAlert in a specified region is completed and the necessary funding is secured for long-term operation and maintenance. When an earthquake occurs, seismic waves radiate from the rupturing fault like waves on a pond. It is these waves that people feel as earthquake shaking and that can cause damage to structures. Using networks of ground-motion sensors and sophisticated computer algorithms, ShakeAlert can detect an earthquake seconds after it begins, calculate its location and magnitude, and estimate the resulting intensity of shaking. Early warnings of impending shaking are then sent to people and systems that may experience damaging shaking, allowing them to take appropriate protective actions. Depending on the user’s distance from the earthquake, alerts may be delivered before, during, or after the arrival of strong shaking. There will almost always be a region near the earthquake epicenter where alerts arrive after damaging shaking has begun. The ShakeAlert system updates its ground-motion estimates as an earthquake grows larger. In response to the FY 2022 congressional direction, the USGS worked with the State of Alaska to devise this implementation plan for ShakeAlert expansion to Alaska. The USGS engaged with the Alaska Division of Homeland Security and Emergency Management (DHS&EM) and the Alaska Division of Geological and Geophysical Surveys (DGGS). A cooperative agreement was awarded to the Alaska Earthquake Center (AEC) at the University of Alaska Fairbanks (UAF) for their contributions to the plan and their work coordinating with other networks in Alaska. The USGS engaged with the Alaska Seismic Hazards Safety Commission (ASHSC) throughout the process. The USGS also held a series of Alaska stakeholder engagements. The process of developing the implementation plan was facilitated by contracted staff from Corner Alliance, which is a government consulting firm. This implementation plan describes the details and estimates the costs for a Phase 1 expansion of the ShakeAlert system to Alaska. A geographically limited Phase 1 goal was chosen that covers the highest risk and most populated areas of Alaska. The areas proposed encompass the State’s main population centers and 90 percent of the State’s population. This Phase 1 design is considered very challenging and ambitious from the viewpoint of network operators. The lessons learned if this plan is implemented could be used to consider subsequent phases to expand EEW beyond Phase 1 in Alaska in the future. ShakeAlert is built on the foundation of the sensor networks and data processing infrastructure of the USGS-led Advanced National Seismic System (ANSS). This implementation plan calls for a total of 450 high-quality, real-time EEW-capable ANSS seismic stations in Alaska: 270 new stations, 160 upgraded stations, and 20 existing stations. These seismic station numbers are based on a station spacing of 10 kilometers (km) in urban areas, 20 km in seismic source areas that endanger population centers, and 40 km in other areas. The associated costs also include support for some EEW-capable global navigation satellite system (GNSS) stations, with a focus on improving warnings for large subduction zone earthquakes. For effective EEW, ShakeAlert requires low-latency, high-availability, robust telemetry links to deliver continuous, real-time data from field stations to the data centers. The Alaska data processing hardware infrastructure would follow the general design for fail-safe operation that is used for the ShakeAlert system on the West Coast. The ShakeAlert architecture uses two independent layers: the production layer for earthquake processing and the alert layer to make alerting decisions and serve alerts to users. This implementation plan includes two geographically separated data centers in Alaska, each with two fully independent production and alert layers using the same system design developed for the West Coast. As of March 2024, the ShakeAlert system is at version 3.0.1, with more advanced versions in the development and testing pipeline. ShakeAlert originally used two algorithms to determine the location and magnitude of earthquakes using seismic data. A third algorithm that can calculate very large magnitudes of very large earthquakes with geodetic data was added in March 2024. ShakeAlert publishes several data and alert products to meet the needs of different users. All messages include the location of the earthquake, either as a point or a line, and its magnitude. Ground-shaking estimates are published in two forms, as ground-motion contours and a map grid. Providing adequate warning time for strong shaking (the “target threshold”) requires sending alerts at a threshold lower than that strong shaking level (the “alert threshold”). The thresholds for public alerting in Alaska would be a joint USGS and State decision. To have the greatest benefit, ShakeAlert-powered alerts would be delivered to institutional users and individuals by all practical pathways. The USGS alert layer can support thousands of institutional users and alert redistributors, but the USGS does not have the mission nor the infrastructure and expertise to perform mass notifications to the public or implement automatic actions for end users of the alerts. To meet this need, ShakeAlert recruits private sector “technology enablers” that have the necessary expertise to develop end-user implementations using EEW alerts with the goal of stimulating an EEW industry. Earthquake early warning alerts are useless if people do not know how to respond to them. Although the alert messages include instructions about what to do (drop, cover, and hold on), alerts are more effective if people have been trained in advance. Messages about ShakeAlert’s capabilities, limitations, and benefits could be integrated with existing earthquake education programs, including State-run programs. Therefore, ShakeAlert would coordinate with both public and private partners and stakeholders through various partnerships and agreements to accomplish consistent and ongoing public earthquake hazard education. The estimated capital cost of completing the computing infrastructure and sensor networks for the Phase 1 ShakeAlert expansion to Alaska is approximately $66 million in 2024 dollars. The annual operation and maintenance cost of the completed system is estimated to be $12 million per year in 2024 dollars when fully built out.

Alaska↗
USGS2025

Migration of seismicity from the mantle to the upper crust beneath Harrat Lunayyir volcanic field, Saudi Arabia

Harrat Lunayyir is a volcanic field in Saudi Arabia that experienced a Mw~5.4 earthquake driven by an upper-crustal dike intrusion in May 2009. This volcanic field has exhibited numerous forms of volcanic seismicity both prior to and since the 2009 dike intrusion. Significantly, earthquakes within the lithospheric mantle and, rarely, the lower crust are present in the two-decade long seismicity catalog of Harrat Lunayyir. Here we analyze 24 years of volcanic seismicity at Harrat Lunayyir from 1998 to 2022. We find that: 1) precursory seismicity began at least eight years prior to the 2009 event, with a particularly notable seismic episode one year prior; 2) lithospheric mantle seismicity is highly localized in space and in time, largely occurring in discrete sequences lasting on the order of a few hours to a few days; 3) one seismic sequence clearly migrates upward from the lithospheric mantle to the upper crust, including seismicity within the nominally ductile lower crust; 4) crustal seismicity has been slowly declining over time; and 5) lithospheric-mantle seismicity does not show any apparent decline with time. From these observations we infer that the seismicity is driven by magmatic fluids or volatiles, and seismic monitoring of this volcanic field should continue into the future.

Harrat Lunayyir volcanic field↗
USGS2026

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona↗
USGS1997

Hydrologic and water-quality conditions in the Horse Creek basin, west-central Florida, October 1992-February 1995

A baseline study of the 241-square-mile Horse Creek basin was undertaken from October 1992 to February 1995 to assess the hydrologic and water-quality conditions of one of the last remaining undeveloped basins in west-central Florida. During the period of the study, much of the basin remained in a natural state, except for limited areas of cattle and citrus production and phosphate mining. Rainfall in 1993 and 1994 in the Horse Creek basin was 8 and 31 percent, respectively, above the 30-year long-term average. The lowest and highest maximum instantaneous peak discharge of the six daily discharge stations occurred at the Buzzard Roost Branch and the Horse Creek near Arcadia stations with 185 to 4,180 cubic feet per second, respectively. The Horse Creek near Arcadia station had the lowest number of no-flow days with zero days and the Brushy Creek station had the highest number with 113 days. During the study, the West Fork Horse Creek subbasin had the highest daily mean discharge per square mile with 30.6 cubic feet per second per square mile, and the largest runoff coefficient of 43.7 percent. The Buzzard Roost Branch subbasin had the lowest daily mean discharge per square mile with 5.05 cubic feet per second per square mile, and Brushy Creek and Brandy Branch shared the lowest runoff coefficient of 0.6 percent. Brandy Branch had the highest monthly mean runoff in both 1993 and 1994 with 11.48 and 19.28 inches, respectively. During the high-baseflow seepage run, seepage gains were 8.87 cubic feet per second along the 43-mile Horse Creek channel. However, during the low-baseflow seepage run, seepage losses were 0.88 cubic foot per second. Three methods were used to estimate average annual ground-water recharge in the Horse Creek basin: (1) well hydrograph, (2) chloride mass balance, and (3) streamflow hydrograph. Estimated average annual recharge using these three methods ranged from 3.6 to 8.7 inches. The high percentage of carbonate plus bicarbonate analyzed at the Carlton surficial aquifer well could indicate an upward ground-water flow from the underlying intermediate aquifer system. Based on constituent concentrations in water samples from the six daily discharge stations, concentrations generally are lower in the upper three subbasins, West Fork Horse Creek, Upper Horse Creek, and Brushy Creek than in the lower three subbasins. Typically, concentrations were highest for major ions at Buzzard Roost Branch and nutrients at Brushy Creek.

Florida↗
USGS2025

Ordovician stratigraphy, structure, and karst of the Falling Spring Valley, Alleghany County, Virginia, USA

This one-day trip highlights new findings on a preliminary bedrock geologic map that shows results from ongoing geologic mapping in the Falling Spring Valley of Alleghany County, Virginia, USA, which is the southern end of the larger Warm Springs Valley, an elongated anticlinal valley rimmed by Ordovician and Silurian siliciclastic rocks, and which is famous for its thermal springs. This mapping includes stratigraphic, structural, and karst field and lab research focused on the Ordovician strata exposed in the area, the oldest of which is the dolomitic upper part of the Beekmantown Formation (Lower Ordovician, Darriwilian), and the youngest of which is the Juniata Formation (Upper Ordovician, Katian), a sequence of siliciclastic redbeds. Warm Springs Valley is the location of the only known caves in the eastern United States—three at present—with thermal waters flowing in some of their passages. Stops on the trip will highlight key details from mapping efforts, primarily within the structurally deformed Ordovician carbonate sequence that is exposed in the core and limbs of the anticline, as well as the associated karst features that are developed in those carbonate rocks, including results of recent dye traces and water temperature monitoring that have improved our understanding of the karst hydrogeologic systems developed in these strata.

Virginia↗
USGS2024

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California↗
USGS1998

Assessment of the hydraulic connection between ground water and the Peace River, west-central Florida

The hydraulic connection between the Peace River and the underlying aquifers along the length of the Peace River from Bartow to Arcadia was assessed to evaluate flow exchanges between these hydrologic systems. Methods included an evaluation of hydrologic and geologic records and seismic-reflection profiles, seepage investigations, and thermal infrared imagery interpretation. Along the upper Peace River, a progressive long-term decline in streamflow has occurred since 1931 due to a lowering of the potentiometric surface of the Upper Floridan aquifer by as much as 60 feet because of intensive ground-water withdrawals for phosphate mining and agriculture. Another effect from lowering the potentiometric surface has been the cessation of flow at several springs located near and within the Peace River channel, including Kissengen Spring, that once averaged a flow of about 19 million gallons a day. The lowering of ground-water head resulted in flow reversals at locations where streamflow enters sinkholes along the streambed and floodplain. Hydrogeologic conditions along the Peace River vary from Bartow to Arcadia. Three distinctive hydrogeologic areas along the Peace River were delineated: (1) the upper Peace River near Bartow, where ground-water recharge occurs; (2) the middle Peace River near Bowling Green, where reversals of hydraulic gradients occur; and (3) the lower Peace River near Arcadia, where ground-water discharge occurs. Seismic-reflection data were used to identify geologic features that could serve as potential conduits for surface-water and ground-water exchange. Depending on the hydrologic regime, this exchange could be recharge of surface water into the aquifer system or discharge of ground water into the stream channel. Geologic features that would provide pathways for water movement were identified in the seismic record; they varied from buried irregular surfaces to large-scale subsidence flexures and vertical fractures or enlarged solution conduits. Generally, the upper Peace River is characterized by a shallow, buried irregular top of rock, numerous observed sinkholes, and subsidence depressions. The downward head gradient provides potential for the Peace River to lose water to the ground-water system. Along the middle Peace River area, head gradients alternate between downward and upward, creating both recharging and discharging ground-water conditions. Seismic records show that buried, laterally continuous reflectors in the lower Peace River pinch out in the middle Peace River streambed. Small springs have been observed along the streambed where these units pinch out. This area corresponds to the region where highest ground-water seepage volumes were measured during this study. Further south, along the lower Peace River, upward head gradients provide conditions for ground-water discharge into the Peace River. Generally, confinement between the surficial aquifer and the confined ground-water systems in this area is better than to the north. However, localized avenues for surface-water and ground-water interactions may exist along discontinuities observed in seismic reflectors associated with large-scale flexures or subsidence features. Ground-water seepage gains or losses along the Peace River were quantified by making three seepage runs during periods of: (1) low base flow, (2) high base flow, and (3) high flow. Low and high base-flow seepage runs were performed along a 74-mile length of the Peace River, between Bartow and Nocatee. Maximum losses of 17.3 cubic feet per second (11.2 million gallons per day) were measured along a 3.2-mile reach of the upper Peace River. The high-flow seepage run was conducted to quantify losses in the Peace River channel and floodplain between Bartow and Fort Meade. Seepage losses calculated during high-flow along a 7.2-mile reach of the Peace River, from the Clear Springs Mine bridge to the Mobil Mine bridge, were approximately 10 percent of the river flow, or 118 cubic feet per second. Calculated seepages along the Peace River in Hardee and De Soto Counties were inconclusive, because most seepages were within the range of discharge measurement error. Two continuous aerial thermal infrared imagery surveys were conducted to locate sites of ground-water discharge along the Peace River. Although temperature and hydrologic conditions were ideal to observe spring flow using thermal infrared imaging techniques, no sources of ground-water discharge were identified using this method. Diffuse ground-water seepage may, however, provide significant ground-water discharge.

Florida↗
USGS2025

Predictive regressive models of recent marsh sediment thickness improve the quantification of coastal marsh sediment budgets

Coastal marsh wetlands experience variations in vertical gains and losses through time, which have allowed them to infill relict topography and record variations in drivers. The stratigraphic unit associated with the development of the marsh also reflects the long-term importance of key ecosystem services supplied by the marsh environment, including carbon storage and storm mitigation. Mapping these coastal wetland sediments and the marsh unit thickness is challenging as traditional coastal geophysical tools are not easily deployable (acoustic methods) or are unreliable in saline-soil environments (e.g., ground-penetrating radar), leaving core-based methods the most viable mapping method. In the present study, we utilized prior information on the geologic architecture of the region to select spatial and physical metrics that likely persisted throughout evolution of the marsh during the late Holocene. We then assessed the individual and collective power of these metrics to predict marsh thickness observed from cores. Employing regressive predictive models powered by these data, we improve the quantification of marsh thickness for a coastal fringing marsh within the Grand Bay estuary in Mississippi and Alabama (USA). The information gained from this approach yields improved estimates of the carbon stocks in this environment. Additionally, the stored sediment masses reflect the past, and potential future, persistence of the Grand Bay marsh under historical and present marsh-estuarine sediment exchange fluxes. Such improvements to both the sediment budget of recent marsh stratigraphic units and the spatial extent provide new resources for comparison with large-scale landscape models, the latter of which may be used, when validated, to predict future change and ecosystem transformations.

Alabama. Mississippi↗
USGS2024

Characteristics of the fault damage zone From high-resolution seismic imaging along the Palos Verdes Fault, California

The distribution and intensity of fault damage zones provides insight into fault activity and its relationship to fluid flow in the crust. Presently, measures of the in-situ distribution of fault damage remain limited and along-strike studies are rare. This study focuses on an offshore section Palos Verdes Fault damage zone that spans 28 km, near Los Angeles, California. To investigate the previously unresolved shallow (∼400 m below the seafloor) fault damage zone we use densely spaced (∼500 m line separation) newly collected sparker multichannel seismic lines and sub-bottom profiles. The combination of high-resolution acquisition methods and specialized seismic processing workflows provide improved imaging of shallow faulting. We apply a multi-trace similarity technique to identify discontinuities in the seismic data that may be attributed to faults and fractures. This fault detection approach reveals diverse fault damage patterns on adjacent seismic profiles. However, a discernible damage zone pattern emerges by stacking multiple damage detection profiles along strike. We find that peak damage identified in this way corresponds to the active main fault strand, confirmed in this study, and thus the technique may be useful for identifying active fault strands elsewhere. Additionally, we observe that the variable width of the damage zone along strike is controlled by fault obliquity. Furthermore, our observations reveal a correlation between fault damage and seafloor fluid seeps visible in the water column, suggesting that damage plays a role in controlling fluid flow around the fault.

California↗
USGS2025

Using gridded seismicity to forecast the long-term spatial distribution of earthquakes for the 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

Gridded (or background) seismicity models are a critical component of probabilistic seismic hazard assessments, accounting for off‐fault and smaller‐magnitude earthquakes. They are typically developed by declustering and spatially smoothing an earthquake catalog to estimate a long‐term seismicity rate that can be used to forecast future earthquakes. Here, we present new gridded seismicity models for use in the 2025 National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The previous PRVI NSHM was released in 2003, and our new models incorporate updates to both data and methodology. We utilize an updated earthquake catalog based on improved Puerto Rico Seismic Network data with newly characterized completeness epochs. The catalog is divided into crustal, subduction interface, and intraslab seismicity using new methods and Slab2 subduction zone geometries. To forecast the long‐term spatial distribution of earthquakes, we use an updated methodology developed for the 2023 U.S. 50‐state NSHM, considering three declustering methods and two spatial smoothing methods based on 2D Gaussian kernels. To adapt it for the complex seismotectonics of the region, we also adopt probabilistic methods to account for events with unknown depths and uncertainties in tectonic classification, and develop a new method for spatial scaling to counteract the effects of spatial variability in network coverage while maintaining the use of smaller events. Finally, we test the performance of these spatial models in forecasting the location of M w ≥ 5earthquakes in the region. Our updated methodology improves the representation of epistemic uncertainty relative to the 2003 model, and our results demonstrate the effectiveness of the new measures we have introduced to address heterogeneities in network detection and systematically evaluate forecast performance.

Puerto Rico, U.S. Virgin Islands↗