USGS Science⌕ Search

SEARCH · USGS Science

Results for “Philosophical Transactions of the Royal Society B: Biological Sciences”

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

The projected exposure and response of a natural barrier island system to climate-driven coastal hazards

Accelerating sea level rise (SLR) and changing storm patterns will increasingly expose barrier islands to coastal hazards, including flooding, erosion, and rising groundwater tables. We assess the exposure of Cape Lookout National Seashore, a barrier island system in North Carolina (USA), to projected SLR and storm hazards over the twenty-first century. We estimate that with 0.5 m of SLR, 47% of current subaerial barrier island area would be flooded daily, and the 1-year return period storm would flood 74%. For 20-year return period storms, over 85% is projected to be flooded for any SLR. The modelled groundwater table is already shallow (< 2 m deep), and while projected to shoal to the land surface with SLR, marine flooding is projected to overtake areas with emergent groundwater. Projected shoreline retreat reaches an average of 178 m with 1 m of SLR and no interventions, which is over 60% of the current island width at narrower locations. Compounding these hazards is subsidence, with one-third of the study area currently lowering at > 2 mm/yr. Our results demonstrate the difficulty of managing natural barrier systems such as those managed by federal park systems tasked with maintaining natural ecosystems and protecting cultural resources.

North Carolina↗
USGS2024

Correlation analysis of groundwater and hydrologic data, Kaloko-Honokōhau National Historical Park, Hawai‘i

Designated in 1978, Kaloko-Honokōhau National Historical Park is located on the west coast of the Island of Hawaiʻi. The Kaloko-Honokōhau National Historical Park encompasses about 1,200 acres of coastal land and nearshore ecosystems, which include wetlands, anchialine pools (landlocked bodies of brackish water with hydrologic connections to the ocean), fishponds, a fishtrap, and coral reefs. These nearshore ecosystems are dependent on groundwater discharge with a freshwater component and provide habitat for threatened and endangered, endemic species, such as the orangeblack Hawaiian damselfly ( Megalagrion xanthomelas ) and the Hawaiian coot (ʻAlae keʻokeʻo, Fulica alai ). The populations of these native species, however, are threatened because of habitat loss related to urban development and environmental changes. Kaloko-Honokōhau National Historical Park is within the Keauhou aquifer system and the North Kona District, which experienced a 52 percent resident-population increase between 2000 and 2020 and a 41 percent visitor increase between 2008 and 2019. To support the current water demand associated with this growing population, groundwater is the primary source of freshwater used in the North Kona District, with about 15 million gallons of groundwater withdrawn from the Keauhou aquifer system per day since 2009. With anticipated development, future (2015–35) groundwater withdrawal from the Keauhou aquifer system is projected to be about 55 percent greater than recent (2012–14) withdrawal. Because Kaloko-Honokōhau National Historical Park is located within a coastal aquifer, natural and human-induced changes can affect the quality and quantity of groundwater, which can threaten groundwater-dependent ecosystems. To improve understanding of recent groundwater conditions, the U.S. Geological Survey, in cooperation with the National Park Service, undertook this study to document correlations between hydrologic time-series datasets from sites in and near Kaloko-Honokōhau National Historical Park using the nonparametric (distribution-free) Kendall’s tau statistical test. For the statistical analyses, dependent variables representing the groundwater system include groundwater level, the groundwater-level difference between pairs of sites, and specific conductance, and independent variables include datasets of sea level, rainfall, and groundwater withdrawal. About 34 percent of the 140 non-time-lagged Kendall’s tau statistical tests evaluated in this report are statistically significant ( p -value ≤ 0.050) with generally weak (0.1 ≤ tau ≤ 0.2) to moderate (0.2 ≤ tau ≤ 0.3) correlations. Groundwater levels measured at monitoring sites have the strongest correlation with the multivariate El Niño–Southern Oscillation index and withdrawal from production wells at the nearby Kohanaiki Private Club Community. Specific conductance is not consistently and significantly correlated with the independent hydrologic variables investigated in this report. Because the relations between hydrologic variables are commonly not instantaneous, a second set of correlations was evaluated after applying a range of time lags to the independent variable datasets. Relative to the non-time-lagged case (the set of correlations that did not use time-lagged independent variables), some of the time-lagged independent variables improved correlations with some of the dependent variables. For a particular independent variable, similar time lags were expected between the independent variable and dependent variable at all four monitoring sites. However, different time lags among the four sites sometimes produced the strongest correlations. This study identified several correlations that are statistically significant and hydrologically plausible, but the correlations could indicate that multiple concurrent factors are controlling the observed groundwater-system response, which might be better addressed using multivariate analyses. This study only investigates bivariate correlations, which may not explain all the variance in the data. The correlations analyzed in this report are limited by the quantity of available hydrologic data in the area near Kaloko-Honokōhau National Historical Park and are based on 14 years of time-series data, which were aggregated to a relatively coarse monthly temporal resolution that represents the minimum resolution common to all datasets.

Hawaii↗
USGS2025

Characterization of stream water quality and groundwater levels in the Central Pine Barrens region, Suffolk County, New York, 2017–23

The area locally known as the “Central Pine Barrens” region, located in Suffolk County, New York, contains most of Long Island’s preserved and undeveloped land. This region overlays an aquifer system that provides potable groundwater for residents of Suffolk County. Between 2017 and 2023, the U.S. Geological Survey, in cooperation with the Central Pine Barrens Joint Planning & Policy Commission and the Town of Brookhaven, monitored groundwater levels and stream water quality in this region. Groundwater levels were measured monthly at five wells and continuously (15-minute intervals) at a sixth well. Water quality was monitored at five locations in the Carmans River and at two locations in the Peconic River, and samples were analyzed for major ions, trace elements, nutrients, pharmaceuticals, and pesticides. The major ion compositions at the sites were mainly sodium-chloride type waters, and compositions varied the most at the furthest upstream sites in both streams. Concentrations above aquatic-life criteria thresholds also occurred most frequently at the furthest upstream sites. The seasonal patterns of nutrient loads and concentrations varied between the Carmans and Peconic Rivers. Several organic compounds including pharmaceuticals, domestic use products, and pesticides were detected at low concentrations in both streams. Metformin was the most frequently detected pharmaceutical compound, and herbicides were the most frequently detected pesticide class. Water-quality conditions influenced by anthropogenic contributions are a result of current and historical land use, and these contributions include onsite wastewater disposal systems, commercial or domestic fertilizers and pesticides, and urban or industrial contaminants in road runoff. This study characterizes and improves understanding of the current hydrologic conditions in the Central Pine Barrens region and the study findings can help inform the development of plans to manage, protect, and restore water resources.

New York↗
USGS2024

Using geologic mapping to understand temporal and spatial relations of closely clustered to concurrent latest Holocene surface ruptures on two intersecting faults, south-central Mojave Desert, California

The Pinto Mountain Fault Zone (PMFZ) marks a major structural boundary between east-oriented sinistral faults of the eastern Transverse Ranges (to the south) and northwest-oriented dextral faults of the south-central Mojave Desert (to the north). These structural fault systems comprise sinistral and dextral deformational domains of the Eastern California Shear Zone (ECSZ) that intersect one another in the Copper Mountain and Twentynine Palms areas. The U.S. Geological Survey (USGS) is conducting detailed geologic mapping and geochronologic investigations designed to clarify geometric, kinematic, and temporal relations among the two domains, that are focused on the central portion of the left-lateral PMFZ near its intersection with the right-lateral Copper Mountain Fault (CMF) and Mesquite Lake Fault Zone (MLFZ).

California↗
USGS2026

Deaf, deafblind, and hard of hearing university student experiences with earthquake early warning in the United States: Evaluating language planning and technology access

The growing literature on deaf and hard of hearing (DHH+) populations and disasters demonstrates that emergency communication (including alerts) is not reaching global DHH + individuals with dangerous impacts for morbidity and mortality. This is the first research study in the U.S. to qualitatively explore the experiences of DHH + persons with earthquake early warning (EEW) through group-based dialogue sessions. The study investigates eight DHH + university students'past earthquake experiences, access to EEW alerts, and perceptions of ShakeAlert Ⓡ , an EEW system for detecting earthquakes and alerting residents of California, Oregon, and Washington. Findings highlight key gaps in disaster alert usability within four thematic areas: lack of messaging in participants' language(s), unclear alert messaging, deficient message delivery mechanisms for deafblind persons, and insufficient access to earthquake information and training that leads to dependence on informal information networks. Weaknesses identified in these four themes reduce DHH + trust in EEW systems and compromise the capacity of alert recipients to take swift protective action or to mentally prepare before shaking starts. The study also underscores structural factors such as insufficient linguistic representation in disaster language planning and technology design, which ignores the linguistic and sensory access needs of DHH + individuals. Building on disaster language planning frameworks, we recommend involving DHH + populations to co-develop EEW alerts. By centering DHH + perspectives, this research contributes to ongoing efforts to ensure that EEW systems reach everyone.

California↗
USGS2024

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗
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↗
USGS1977

Leachate plumes in ground water from Babylon and Islip landfills, Long Island, New York

Landfills operated by the Towns of Babylon and Islip in southwest and central Suffolk County contain urban refuse, incinerated garbage, and scavenger (cesspool) waste; some industrial refuse is deposited at the Babylon site. The Islip landfill was started in 1933, the Babylon landfill in 1947. The landfills are in contact with and discharge leachate into the highly permeable upper glacial aquifer (hydraulic conducitvity 190 to 500 ft/d [feet per day]). The aquifer is 74 feet thick at the Babylon landfill and 170 feet thick at the Islip landfill. The leachate-enriched water occupies the entire thickness of the aquifer beneath both landfills, but hydrologic boundaries retard downward migration of the plumes to deeper aquifers. The Babylon plume is 1,900 feet wide at the landfill and narrows to about 700 feet near its terminus 10,000 feet from the landfill. The Islip plume is 1,400 feet wide at the landfill and narrows to 500 feet near its terminus 5,000 feet from the landfill. Hydrochemical maps and sections show the distribution of the major chemical constituents of the plumes. The most highly leachate-enriched ground water obtained was from the Babylon site; it contained 860 mg/L (milligrams per liter) sodium, 110 mg/L potassium, 565 mg/L calcium, 100 mg/L magnesium, 2,700 mg/L bicarbonate, and 1,300 mg/L chloride. Sulfate was notably absent or in low concentration in most parts of both plumes. Nitrogen in plume water was mostly in the form NH4+, and concentrations as high as 90 mg/L were found; concentrations of nitrogen as N in the plume were less than 10 mg/L. As much as 440 mg/L iron and 190 mg/L manganese were found in the leachate-enriched water. Samples were also tested for arsenic, boron, cadmium, cobalt, chromium, copper, mercury, nickel, lead, selenium, strontium, and zinc. Boron was more or less ubiquitous and was found in concentrations as high as 2 mg/L. Organic carbon was found in concentrations as high as 2,250 mg/L in the most highly leachate-enriched water but attenuated rapidly to less than 20 mg/L. Dissolved-solids concentrations near the landfills were between 400 and 3,000 mg/L at Babylon and between 500 and 1,500 mg/L at Islip. Ground-water temperatures near the landfills exceed those in ambient water by as much as 7°C at Babylon and 16°C at Islip. Heat contributed by the landfills was mostly dissipated within 0.4 mi of the landfill, butat Islip, the warm leachate-enriched water extended 0.5 mi downgradient. The entrance of leachate into the less dense ground water as pulsations after rainfall may explain the presence of high leachate enrichment near the bottom of the aquifer. A comparison of the physical characteristics of leachate-enriched ground water with those of ambient water suggests that the downward movement of leachate results from its greater density. Simulation of the movement and dispersion of the Babylon plume with a mathematical dispersion model indicated the coefficient of longitudinal dispersion to be about 60 ft 2 /d (feet squared per day) and the ground-water velocity to be 1 ft/d. However, the velocity determined from the hydraulic gradient and public-supply wells in the area was 4 ft/d, which would cause a plume four times as long as that predicted by the mathmatical dispersion model. At the Islip site, the plume was one-third the length calculated on the basis of the age of the landfill. The shortness of the plumes has not been explained; it may be a result of the leachate's having been too dilute to form a plume during the early years of the landfills.

New York↗
USGS2024

Post-glacial stratigraphy and late Holocene record of great Cascadia earthquakes in Ozette Lake, Washington, USA

Ozette Lake is an ~100-m-deep coastal lake located along the outer coast of the Olympic Peninsula (Washington, USA); it is situated above the locked portion of the northern Cascadia megathrust but also relatively isolated from active crustal faults and intraslab earthquakes. Here we present a suite of geophysical and geological evidence for earthquake-triggered mass transport deposits (MTDs) and related turbidite deposition in Ozette Lake since ca. 14 ka. Comprehensive high-resolution bathymetry data, seismic reflection profiles, and sediment cores are used to characterize the post-glacial stratigraphic framework and examine paleoseismic evidence in the lacustrine sediments. Stacked sequences of MTDs along the steep eastern flanks of the lake appear to grade basin-ward from thick, chaotic, blocky masses to thin, parallel-bedded turbidite beds. The discrete turbidite event layers are separated by fine-grained (silt and clay) lake sedimentation. The event layers are observed throughout the lake, but the physical characteristics of the deposits vary considerably depending on proximity to primary depocenters, steep slopes, and subaqueous deltas. A total of 30–34 event deposits are observed in the post-glacial record. Radiometric dating was used to reconstruct a detailed sedimentation history over the last ~5.5 k.y., develop an age model, and estimate the recurrence (365–405 yr) for the most recent 12 event layers. Based on sedimentological characteristics, temporal overlap with other regional paleoseismic chronologies, and recurrence estimates, at least 10 of the dated event layers appear to be sourced from slope failures triggered by intense shaking during megathrust ruptures; the recurrence interval for these 10 events is 440–560 yr. Thus, Ozette Lake contains one of the longest and most robust geological records of repeated shaking along the northern Cascadia subduction zone.

Washington↗
USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗
USGS2025

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗
USGS2026

Computing flow-field distortion coefficients from well-construction and formation properties

Direct measurements of groundwater velocity made with borehole flowmeters in screened wells must be compensated for the effects of flow-field distortion (also known as borehole acceleration). A theoretical equation developed by Drost et al. (1968) and simple inputs describing hydraulic properties of well construction and geologic formation were programmed into an Excel workbook to facilitate computation by groundwater-flowmeter users. Tables describing the physical and hydraulic properties for well constructions and gravel pack media are provided with an example to facilitate use of the workbook. Groundwater flowlines converge or diverge as they pass from a geologic formation, through a gravel pack and well screen. The extent of flowline convergence or divergence and the value of the flow-field distortion coefficient is related to the relative changes in hydraulic conductivity of the well screen, gravel pack, and geologic formation. Convergence or divergence is accompanied by acceleration or deceleration of groundwater. Direct measurements of groundwater velocity at the center of the monitoring well can be adjusted to provide a more accurate estimate of velocity in the formation by applying a correction for flow-field distortion. Variables required to compute the flow-field distortion coefficient include the hydraulic conductivity of the gravel pack, well screen, and the geologic formation surrounding the well screen; the borehole radius, and the inside radius and outside radius of the well screen.

Groundwater↗
USGS2026

Simulation of groundwater flow in the Silurian aquifer, eastern Iowa, 2020–45

The Silurian aquifer is an important water source for municipalities, industry, and rural households and communities in eastern Iowa, including Johnson County. Increasing demand for groundwater from the aquifer indicated that analytical tools may be beneficial to quantify groundwater resources and inform water-management decisions for the aquifer. The U.S. Geological Survey, in cooperation with the Johnson County Board of Supervisors, developed conceptual and numerical groundwater models to simulate water levels in the Silurian aquifer, determine groundwater budgets, and forecast changes in groundwater levels through 2045. A MODFLOW numerical model was constructed to match water levels for 2020 through 2022. Model performance was assessed by comparing simulated and observed water-level hydrographs at several pumping and monitoring wells. Hydrograph comparison indicates that the model simulates water levels accurately at monitoring wells that are not near areas of pumping and generally overpredicts water levels near pumping wells, while accurately simulating the drawdown trend over time at pumping wells. This transient model framework was then used to predict water levels through 2045 in a scenario of moderate drought and increased groundwater withdrawals. Results from the predictive model indicate as much as 13 meters of additional drawdown in 2045 from 2020 water levels in the area of pumping wells.

Iowa↗
USGS2025

The impact of burial diagenesis on soil-formed minerals in paleosols using stable isotopes of phyllosilicates and carbonate clumped isotopes

To understand the effects of burial diagenesis on the stable isotope geochemistry of soil-formed clay and carbonate minerals in paleosols, samples were collected from seven cores, spanning middle- to upper-Pennsylvanian strata of the Illinois Basin, with varied maximum burial depths of 1–3 km. Mixed-layer illite-smectite and kaolinite mixtures give δ 2 H and δ 18 O values of −83 ‰ to −36 ‰ and 11.9 ‰ to 21.1 ‰ (VSMOW), respectively. After carbonates were screened petrographically for diagenetic textures using transmitted light and cathodoluminescence, measured clumped isotope Δ 47 values range from 0.504 to 0.563 ‰ (I-CDES). Resulting mineral formation temperatures for phyllosilicate mineral mixtures are 28 to 66 °C (mean = 47 °C), whereas T(Δ 47 ) estimates for calcites are 36 to 61 °C (mean = 45 °C). Calculated δ 18 O water values from which phyllosilicate minerals and calcites precipitated under isotopic equilibrium ranges from −7.1 to −1.2 ‰ and − 1.4 to +4.9 ‰, respectively. Closed and open-system phyllosilicate-fluid exchange modeling indicates that phyllosilicate alteration occurred in the presence of a low temperature brine or meteoric water and is interpreted to occur in a layer-by-layer illitization transformation. Due to the lack of diagenetic textures and positively correlated T(Δ 47 ) and δ 18 O water , calcites are interpreted to have undergone solid-state bond reordering. Despite low to moderate temperatures (<125 °C) and varying depths of shallow burial (1–3 km), solid-state transformation of phyllosilicates and calcites indicates paleosols had prolonged exposure to burial conditions which has implications for the use of paleosol minerals for paleoenvironmental reconstructions.

Illinois, Indiana, Kentucky↗
USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the northern Appalachian orogen, USA

Lithium demand is projected to increase more than 48 times by 2040 due to electric vehicle production and other energy storage needs. Most lithium production is outside of the USA, thereby increasing supply chain vulnerability. The combined end use importance and heightened supply risk of lithium make this lightest metallic element a critical commodity to the USA. To mitigate this supply risk, the US Geological Survey is actively assessing lithium deposits in the USA. Herein, we detail an assessment for lithium-mineralized pegmatites in the US northern Appalachian Mountains. Permissive tracts were generated by cross-referencing tectonic and geologic maps and mineral occurrence data with mappable criteria derived from generalized and region-specific lithium pegmatite ore deposit models; tracts were then ranked as having high, medium, or low permissibility. Available geophysical and geochemical data were found to be of minimal utility for this deposit type at the scale of the assessment. The number of undiscovered deposits were estimated and integrated into probabilistic simulations, which included an expanded and updated global grade and tonnage model of pegmatite-hosted lithium ore. The estimated total amount of undiscovered resources for the northern Appalachian Orogen has a median value of 1,410,000 metric tons of Li 2 O when considering moderate correlation across sub-regions. At a confidence level of 90%, a resource of at least 90,000 metric tons of Li 2 O remains undiscovered, and at a 10% confidence level, a resource of as much as 7,380,000 metric tons Li 2 O remains undiscovered. After applying an up-to-date economic filter to convert median contained lithium to recoverable material, a correlated total of 900,000 metric tons of Li 2 O may be economically extractable, equating to enough Li 2 O to provide the current annual US lithium supply deficit (presently obtained through net imports) for 127 years at 2025 rates of apparent consumption. This period of provision will inevitably shorten with projected increasing consumption rates, emphasizing that further research could be completed to better delineate regions of high lithium resource potential and support exploration and domestic production.

Connecticut, Delaware, Maine, Massachusetts, New H↗
USGS1998

Transmissivity and water quality of water-producing zones in the intermediate aquifer system, Sarasota County, Florida

The intermediate aquifer system is an important water source in Sarasota County, Florida, because the quality of water in it is usually better than that in the underlying Upper Floridan aquifer. The intermediate aquifer system consists of a group of up to three water-producing zones separated by less-permeable units that restrict the vertical movement of ground water between zones. The diverse lithology, that makes up the intermediate aquifer system, reflects the variety of depositional environments that occurred during the late Oligocene and Miocene epochs. Slight changes in the depositional environment resulted in aquifer heterogeneity, creating both localized connection between water-producing zones and abrupt culmination of water-producing zones that are not well documented. Aquifer heterogeneity results in vertical and areal variability in hydraulic and water-quality properties. The uppermost water-producing zone is designated producing zone 1 but is not extensively used because of its limited production capability and limited areal extent. The second water-producing zone is designated producing zone 2, and most of the domestic- and irrigation-supply wells in the area are open to this zone. Additionally, producing zone 2 is utilized for public supply in southern coastal areas of Sarasota County. Producing zone 3 is the lowermost and most productive water-producing zone in the intermediate aquifer system. Public-supply well fields serving the cities of Sarasota and Venice, as well as the Plantation and Mabry Carlton Reserve well fields, utilize producing zone 3. Heads within the intermediate aquifer system generally increase with aquifer depth. However, localized head-gradient reversals occur in the study area, coinciding with sites of intense ground-water withdrawals. Heads in producing zones 1, 2, and 3 range from 1 to 23, 0.2 to 34, and 7 to 42 feet above sea level, respectively. Generally, an upward head gradient exists between producing zones 3 and 2. However, an upward head gradient between producing zones 2 and 1 does not consistently occur throughout Sarasota County, probably the result of greater ground-water withdrawals from producing zone 2 than from producing zone 1. The transmissivity of the intermediate aquifer system is spatially variable. Specific-capacity data from selected wells penetrating producing zones 2 and 3, were used to estimate transmissivity. Estimated transmissivity values for producing zones 2 and 3 range from about 100 to 26,000 feet squared per day and from about 1,300 to 6,200 feet squared per day, respectively. Because the capacity of specific water-producing zones is highly variable from site to site, estimating the performance of a specific water-producing zone as a water resource is difficult. Water samples collected during the study were analyzed for major-ion concentrations. Generally, bicarbonate type water from rock interaction occurs in northern Sarasota County; enriched calcium-magnesium-sulfate type water from deeper aquifers occurs in central Sarasota County; and sodium-chloride type water from saltwater mixing occurs in southern Sarasota County. In some areas of northern Sarasota County, the major-ion concentrations in water are lower in producing zone 2 than in producing zone 1. Major-ion concentrations in water are higher in producing zone 3 throughout the study area. A major objective of the study was to evaluate hydraulic and water-quality data to determine distinctions that could be used to characterize a particular producing zone. However, data indicate that both hydraulic and water-quality properties are highly variable within and between zones, and are more related to the degree of connection between and areal extent of water-producing zones than to aquifer depth and distance from the coast.

Florida↗
USGS2025

Evidence for fluid pressurization of fault zones and persistent sensitivity to injection rate beneath the Raton Basin

Subsurface wastewater injection has increased the seismicity rate within the Raton Basin over more than two decades, with the basin-wide injection rate peaked between 2009-2015. To understand the evolution of injection-induced earthquakes, we systematically analyzed 2016-2024 broadband recordings with a machine-learning-based phase picker and constructed a catalog with 95,993 earthquakes (-1≤ M L ≤4.3). We then inverted for full centroid moment tensors (CMT) for 90 M L ≥ 2 events, with a special interest in constraining the non-double-couple components via probabilistic metrics. Both relocations and CMT solutions support basement-rooted normal faults, including graben and half-graben structures. Furthermore, we observe the non-double-couple components that imply elevated pore pressure in the fault zones. An earthquake cluster emerged in the north-central basin in 2023, preceded by ~1-yr of increased injection volume from wells within 15km. Despite a basin-wide decrease in the injection volume, we highlights the persistence of seismicity that remains to sensitive to injection rates within the Raton Basin.

Colorado, New Mexico↗
USGS2025

Quantitative mineral resource assessment of lithium pegmatite deposits in the Appalachian Orogen, USA

Lithium is classified as a U.S. critical mineral commodity, and its demand is projected to drastically increase through 2040, driven by electric vehicle production and energy storage applications (IEA 2021).Most global lithium production is not in the United States increasing vulnerability to a supply disruption. The U.S. Geological Survey is actively assessing domestic lithium deposits including lithium-bearing pegmatites in the Appalachian orogen. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical, and mineral occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. Estimates were then integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was used to estimate the amount of potentially recoverable undiscovered resources. Preliminary computations for the northern Appalachians, including application of the economic filter to the median recoverable contained resource, yields 900,000 metric tons of Li 2 O that correspond to enough Li 2 O to replace 127 years of import reliance at the current rate (7,100 t Li 2 O/yr; USGS, 2025). For the southern Appalachians, preliminary computations yielded 1,430,000 metric tons of Li 2 O, which corresponds to 201 years of import reliance.

Alabama, Connecticut, Delaware, Georgia, Maine, Ma↗