USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Waterway, Port, Coastal, and Ocean Engineering”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,722 records · Page 20Linked to original sources
USGS2026

Uncertainty quantification of geologic energy storage in depleted gas reservoirs using material balance equations embedded in a hierarchical errors-in-variables model

The storage potential for gas in geologic settings, such as depleted hydrocarbon reservoirs and solution-mined salt caverns, is becoming salient to future energy infrastructure planning. Technologies such as carbon capture, utilization, and storage, carbon dioxide-enhanced oil recovery, and natural gas and hydrogen storage help meet growing energy demands, reduce carbon emissions to meet climate goals, and provide energy security amid geopolitical uncertainties. Therefore, estimates of underground gas storage capacity could be useful for efficiently navigating the energy transitions. Material balance is a fundamental method in reservoir engineering for estimating original gas in place and potential storage capacity at the scale necessary for national assessments of subsurface pore space resources. However, the deterministic method cannot accommodate multiple data sources or quantify uncertainty in predictions. In this study, a method that embeds material balance equations within a hierarchical errors-in-variables model is proposed which allows the estimation of the distributions of reservoir properties needed for assessments. Uncertainties associated with these reservoir properties have traditionally been expert-elicited, whereas the uncertainty estimates from the proposed models are data-driven. Capacity and uncertainty estimates can be used in a probabilistic resource assessment, supplementing information traditionally used by assessors or even replacing this expert elicitation step when data are unavailable. Various regression models are compared in a case study of the Michigan Basin, a large contributor to the United States’ current natural gas storage capacity. In particular, errors-in-variables models help ameliorate regression dilution and can quantify uncertainty in predictions of pressure in addition to storage capacity. Overfitting is addressed by quantifying generalization error and model averaging in simple and stratified cross-validation against reported working gas capacity, representing the varying quality and quantity of available data. Incorporating a statistical framework into existing numerical methods in reservoir engineering can improve the quality of estimation, and in particular, this method brings rigor to uncertainty quantification as part of a larger effort by the U.S. Geological Survey to assess domestic energy gas storage resources in depleted hydrocarbon reservoirs.

Mathematical Geosciences↗
USGS2024

Plan to coordinate post-earthquake investigations supported by the National Earthquake Hazards Reduction Program (NEHRP)

Introduction This report presents a plan supported by the National Earthquake Hazards Reduction Program (NEHRP) to coordinate domestic and international post-earthquake investigations (herein called “the Plan”). Post-earthquake scientific and engineering investigations are undertaken to capture critical information to understand the causes and impacts of the event, lessons from which can substantially improve the Nation’s resilience after future earthquakes. NEHRP is the Federal Government’s coordinated nationwide program to reduce risks to life and property from earthquakes. The Plan describes the activation and coordination of the four designated NEHRP Agencies in the Federal Government: Federal Emergency Management Agency, National Institute of Standards and Technology, U.S. National Science Foundation, and U.S. Geological Survey (USGS). The Plan also describes coordination between NEHRP Agencies and other organizations that may participate in pre-event and post-earthquake investigations, including non-NEHRP Federal agencies; State, regional, local, Tribal, and territorial agencies; domestic nongovernmental organizations; academic institutions and affiliated organizations; private companies; foreign governmental agencies and nongovernmental organizations; and international organizations. The Plan delineates the coordination of NEHRP post-earthquake scientific and engineering investigations to document the direct, indirect, and cascading physical and societal impacts from fault rupture and ground shaking hazards and from secondary hazards such as landslides, liquefaction, and tsunamis. In addition, the Plan identifies pre-event activities necessary to ensure that post-earthquake investigations are executed effectively. The USGS is the lead NEHRP Agency for activating and coordinating NEHRP post-earthquake investigations and for implementing this Plan. The USGS also leads coordination of the NEHRP Agencies in completing the pre-event activities identified in the Plan. The Plan has new information and supersedes USGS Circular 1242, “The Plan to Coordinate NEHRP Post-Earthquake Investigations,” which was published in 2003. This second-generation Plan was developed with the assistance of the Applied Technology Council (ATC) of Redwood City, California, under USGS contract 140G0121P0309, ATC-155. A 17-member Project Review Panel provided guidance on plan development, and input was solicited from subject matter experts representing key stakeholder groups and from participants at a public workshop.

Circular↗
USGS1976

Availability of ground water in the area surrounding the Trident submarine construction facility, Kitsap County, Washington

This report provides general information on the water resources--with emphasis on ground-water occurrence and availability--in that part of Kitsap County (referred to as Trident Impact Area) that would be most affected by the development of the Trident submarine-construction facility at Bangor. Ground water occurs principally in a lower, confined aquifer near and below sea level and an upper, unconfined aquifer above sea level; both aquifers are extensive beneath the study area and are composed of sand and gravel layers separated by poorly permeable silt, clay, and till layers. Locally throughout the area efficiently constructed wells tapping each of the two aquifers are capable of yielding 100 to 1,300 or more gallons per minute. The estimated 1970 water use in the study area averaged about 13 million gal/day (gallons per day); of this amount about 9 million gal/day came from surface-water sources--from a large reservoir outside the study area--and about 4 million gal/day came from ground water pumped from the two aquifers. Anticipated water use soon will be about 18 to 21 million gal/day; virtually all the additional quantity required (about 5 to 8 million gal/day) above present use must come from ground-water sources. Preliminary evaluation of the aquifers suggests that an additional 1.5 million gal/day can be developed from the upper aquifer and 7 million gal/day from the lower aquifer. Existing wells tapping the lower aquifer might yield additional water and increase the total yield in the area by 3.5 million gal/day, and new wells drilled in selected areas could produce an additional 3.5 million gal/day from this aquifer. However, additional, large-scale ground-water withdrawal from the lower aquifer could induce saltwater intrusion into wells situated in coastal areas. Recommended future studies include additional data collection, monitoring of chloride concentrations in water from large-yield coastal wells, confirmation of well yields and capacities, and test-well drilling and analysis.

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

Bathymetric maps, surface areas, and storage capacities of Council Grove Lake and Marion Reservoir, Kansas, and Pine Creek Lake, Oklahoma, 2024

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, completed high-resolution multibeam bathymetric surveys to compute new elevation-area-capacity tables for Council Grove Lake and Marion Reservoir, Kansas, and Pine Creek Lake, Oklahoma. Elevation-area-capacity tables identify the relation between the water-surface elevation, surface area, and storage capacity of the lake. The surface areas and storage capacities of each lake were computed from bathymetric surfaces combining multibeam echo sounder data collected in 2024 and light detection and ranging point-cloud data collected in 2016 and 2018.

Kansas, Oklahoma↗
USGS2024

Geophysical characterization of an alkaline‑carbonatite complex using gravity and magnetic methods at Magnet Cove, Arkansas, USA

The Magnet Cove alkaline‑carbonatite complex (MCC), located in the Ouachita Mountains of south-central Arkansas in the United States, hosts an extensive variety of rare rock types and critical mineral resources with physical properties (density and magnetization) that contrast significantly with the sedimentary rocks into which they have intruded. Newly acquired ground-based gravity and magnetic data were used to develop two-dimensional and three-dimensional geophysical models of the Cretaceous-aged Magnet Cove intrusive complex. The models reveal that the MCC: (1) widens out at middle crustal depths to as much 22 km across, and may reach a depth of 20 km; (2) has a total volume (exposed and subsurface) that may be over 800 km 3 ; (3) is likely connected at depth to other intrusions in the Arkansas alkaline province; and (4) has a geometry that is aligned with pre-existing structures such as the Reelfoot rift and the Ouachita orogenic belt, some of which were likely structurally controlled by the Precambrian crystalline basement and the continent-ocean transition zone buried beneath the Ouachita orogen. For the first time, the magnetic models of the MCC account for the presence of strong remanent magnetization. This results in a geophysical workflow necessary to accurately interpret magnetic anomalies over the much larger Arkansas alkaline province, its geologic and structural framework, and critical mineral potential.

Arkansas↗
USGS2024

Integrating Sr isotopes, microchemistry, and genetics to reconstruct Salmonidae species and life history

Recent approaches to fisheries research emphasize the importance of the coproduction of knowledge in building resilient and culturally mindful fisheries management frameworks. Despite widespread recognition of the need for Indigenous knowledge and historical reference points as baseline data, archaeological data are rarely included in conservation biology research designs. Here we propose a novel multiproxy method to learn from former fisheries stewards by generating archaeological data on past salmonid population parameters. We used a newly developed, high throughput qPCR (HT-qPCR) chip, originally designed for environmental DNA (eDNA), for species identification of archaeological salmonid vertebrae. We combine this with the laser ablation split-stream (LASS) approach to identify ocean-migration versus freshwater residency. We test this multidisciplinary approach using both contemporary and archaeological salmonid samples and new radiocarbon dates from the Tronsdal Site on the Skagit River, Washington State, USA. This is a useful approach for extracting information about Salmonidae species and life history diversity from archaeological remains to reconstruct historic baselines for several population parameters in anadromous species with long periods of freshwater residency. The approach outlined in this paper may be particularly useful for research investigating past fisheries dynamics, offering hundreds to thousands of years of temporal depth for modern fisheries management, harvest policies, restoration ecology, and conservation biology.

Idaho, Oregon, Washington↗
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↗
USGS2024

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map↗
USGS2025

A geospatial analysis of water-quality threats from orphan wells in principal and secondary aquifers of the United States

Throughout the history of oil and gas production in the United States, millions of wells have been drilled for exploration and energy production. Hundreds of thousands of unplugged wells are no longer actively producing and are currently under orphan status, with no responsible party obligated for plugging. Orphan wells can pose threats to water resources by providing pathways for contaminants such as hydrocarbons and brines to migrate into water-supply aquifers. In this study, we investigate the potential threats to groundwater resources posed by orphan wells at the national scale. Water-quality data is extremely sparse in relation to orphan wells nationally and may not be suitable for identifying contamination from oil and gas development. We used geospatial and statistical methods to evaluate which principal and secondary aquifer systems may be most susceptible to contamination from orphan wells. Analysis involved three sets of susceptibility factors including: 1) factors related to the number and density of orphan wells; 2) factors that can threaten well integrity and contribute to transport of contaminants; and 3) factors related to groundwater withdrawal rates and the affected populations/communities in the event of water quality disturbances. From a dataset of 117,672 documented orphan wells, 64,203 fall within a principal aquifer system, while the remainder fall within a secondary aquifer system. By assessing the combination of well integrity and hydrogeologic factors within these aquifer systems, five groupings of principal aquifers were identified, where groups ranged from aquifer systems with high numbers of orphan wells, multiple well integrity threats and high withdrawals, to aquifers with a relatively low number of orphan wells, limited well integrity threats and minimal water use. Three regions of the country emerge containing aquifers with higher susceptibility to contamination from orphan oil and gas wells. These regions include 1) The Appalachian Basin (including the Pennsylvanian Aquifer System), 2) The Gulf Coast Aquifers (including the Coastal Lowlands Aquifer system) and 3) The California Aquifers (including the California Coastal Basin Aquifer system). This work is the first multivariate geospatial investigation of orphan wells and groundwater resources on a national scale, and sheds light on which aquifers are most susceptible to groundwater contamination from orphan wells.

Science of the Total Environment↗
USGS2023

Thickness of the Saudi Arabian crust

As part of a joint Saudi Geological Survey (SGS) and U.S. Geological Survey (USGS) project, we analyzed P-wave receiver functions from seismic stations covering most of the Kingdom of Saudi Arabia to map the thickness of the crust across the Arabia Plate. We present an update of crustal-thickness estimates and fill in gaps for the western Arabian Shield and the rifted margin at the Red Sea (the coastal plain), as well as the eastern Arabian Platform. We applied a conventional H-k stacking algorithm and included careful attention to stacking weights, two forms of sedimentary corrections for stations located on the Arabian Platform, and additional processing for noisy stations. We obtained useful results at 154 stations from 898 teleseismic events over a 2-year period from 1995–1997 (for non-SGS stations) and a 6-year period from 2008–2014 (for SGS stations). Average crustal thickness (that is, depth to the Mohorovičić discontinuity [Moho] below the surface) beneath the Red Sea coastal plain (the rift margin) is 29 kilometers (km), beneath the volcanic fields (known in Arabic as harra [plural] or harrat [singular]) is 35 km, beneath the Arabian Shield (excluding harrats) is 37 km, and beneath the Arabian Platform is 38 km. Crustal thinning appears not to extend east of the rift escarpment, suggesting uniform extension that is no broader at depth than at the surface. In contrast to some previous claims that the Arabian Platform crust is thicker than that of the Arabian Shield, we find no statistically significant difference between their whole crustal thicknesses. However, the average subsedimentary crustal thickness (that is, the crystalline crust) for stations on the Arabian Platform is 34 km, 3 km thinner than the crust of the Arabian Shield. Individual station P-wave (pressure) velocity and S-wave (shear) velocity ratios ( V P / V S ) are highly variable for the Arabia Plate, ranging from 1.60 to 1.97 and averaging 1.75, with a standard deviation of 0.07. There are no statistically significant differences between V P / V S ratios of the different geologic regions of Saudi Arabia. Similar V P / V S ratios, coupled with similar crustal thicknesses for harrats and the Arabian Shield, indicate that Cenozoic magmatism has contributed negligibly to crustal growth.

Professional Paper↗
USGS2026

Bascom Laser Diffraction Sedimentology Laboratory, Reston, Virginia

Introduction At the Bascom Laser Diffraction Sedimentology Laboratory, which is located in the Florence Bascom Geoscience Center at U.S. Geological Survey (USGS) headquarters in Reston, Virginia, scientists use physical sedimentology and particle characterization techniques to conduct detailed sediment characterization. Scientists address research problems in collaboration with other USGS science centers, State geological surveys, commercial industry, universities, and other partners. Laboratory capabilities include laser diffractometry for quantitative particle-size analysis, portable x-ray fluorescence (XRF) analysis for determining elemental abundances in rock or sediment samples, petrographic analysis of geologic media, and mechanical sieve analysis. These methods are used to analyze soil and sediment core material from terrestrial, marine, and lacustrine environments, surface sediments from coastal regions, and calcareous materials. Work done by the laboratory supports geologic mapping, resource assessments, land change studies, and geohazard analyses.

Fact Sheet↗
USGS2024

Land-use interactions, Oil-Field infrastructure, and natural processes control hydrocarbon and arsenic concentrations in groundwater, Poso Creek Oil Field, California, USA

Like many hydrocarbon production areas in the U.S., the Poso Creek Oil Field in California includes and is adjacent to other land uses (agricultural and other developed lands) that affect the hydrology and geochemistry of the aquifer overlying and adjacent to oil development. We hypothesize that the distributions of hydrocarbons and arsenic in groundwater in such areas will be controlled by complex interactions between mixed land uses, oil-field infrastructure, and natural processes. In 2020–2021, samples of groundwater and surface water were collected and analyzed for a large suite of inorganic and organic chemicals and isotope and gas tracers to test this hypothesis. Those data are supplemented with ancillary data on historical geochemistry, hydrology, geology, and oil-field infrastructure. Hydrocarbons in groundwater (e.g., methane through pentane gases and benzene) are associated with natural processes (e.g., fault offsets or transition in sediment depositional environment) and oil-field infrastructure (e.g., fluid-migration pathways associated with uncemented annulus in oil wells or unlined pits). Arsenic concentrations >10 μg per liter (μg/L; maximum concentration 12.9 μg/L) are associated with natural processes in old, high-pH groundwater, and more recent recharge of water from natural and/or engineered recharge processes. Along the southwest margin of the oil field, pumping for drinking-water and irrigation supplies in combination with engineered groundwater recharge produce a depression in groundwater elevations where groundwater with elevated sulfate concentrations from agricultural areas and groundwater with hydrocarbons from the oil field mix to produce a zone of sulfate reduction that removes hydrocarbons and arsenic from groundwater but produces elevated sulfide (S 2- ) concentrations (maximum concentration 29 mg per liter, mg/L). In this study, multiple approaches were required to resolve the overlapping effects of land uses, oil-field infrastructure, and natural processes on the distributions of hydrocarbons and arsenic in groundwater. The combined use of geographic, historical, physical, chemical, isotopic, and other information to constrain processes could be a useful approach for studies in other hydrocarbon-production areas. This is particularly important where land uses affect aquifer hydrology to an extent that causes mixing of groundwaters with different chemical compositions.

California↗
USGS2026

Regional variations in sea ice and primary productivity in the Bering Sea during Marine Isotope Stage 11

Marine Isotope Stage (MIS) 11 (424-374 ka) has long been an analog for Holocene climate, because it is the most recent interglacial period with similar orbital conditions. However, there is significant global and regional variability in the climate response to MIS 11 warmth. Here, we review sediment core records from across the Bering Sea to investigate changes in paleoceanographic conditions during Marine Isotope Stages 12-10. Sea ice was present over much of the Bering Sea during MIS 11, but today, none of the sites investigated are ever ice-covered. This suggests that sea ice regimes in the Bering Sea during MIS 11 were different to those of the Holocene. There are also regional differences in the response of sea ice to MIS 11 warming. At the Umnak Plateau, Southeastern Bering Sea, sea ice concentrations decline during deglaciation, but they remain high at the slope sites until Peak MIS 11. Sea ice re-advances over the Umnak Plateau during peak interglacial warmth, at the same time that it declines over the slope sites. Late MIS 11 is characterized by high concentrations of seasonal sea ice at the Umnak Plateau, whilst sea ice at the slope sites fluctuates between consolidated and unconsolidated ice cover. This east-west dichotomy may be explained by changes in the behavior of the Aleutian Low. Productivity increases dramatically during deglaciation due to increased upwelling and sea level rise bringing fresh nutrients into the oceans. This is characterized by increased diatom productivity, increased terrestrial carbon deposition, and laminations at all sites.

Alaska↗
USGS2024

Deep structure of Siletzia in the Puget Lowland: Imaging an obducted plateau and accretionary thrust belt with potential fields

Detailed understanding of crustal components and tectonic history of forearcs is important due to their geological complexity and high seismic hazard. The principal component of the Cascadia forearc is Siletzia, a composite basaltic terrane of oceanic origin. Much is known about the lithology and age of the province. However, glacial sediments blanketing the Puget Lowland obscure its lateral extent and internal structure, hindering our ability to fully understand its tectonic history and its influence on modern deformation. In this study, we apply map-view interpretation and two-dimensional modeling of aeromagnetic and gravity data to the magnetically stratified Siletzia terrane revealing its internal structure and characterizing its eastern boundary. These analyses suggest the contact between Siletzia (Crescent Formation) and the Eocene accretionary prism trends northward under Lake Washington. North of Seattle, this boundary dips east where it crosses the Kingston arch, whereas south of Seattle the contact dips west where it crosses the Seattle uplift (SU). This westward dip is opposite the dip of the Eocene subduction interface, implying obduction of Siletzia upper crust at this southern location. Elongate pairs of high and low magnetic anomalies over the SU suggest imbrication of steeply-dipping, deeply rooted slices of Crescent Formation within Siletzia. We hypothesize these features result from duplication of Crescent Formation in an accretionary fold-thrust belt during the Eocene. The active Seattle fault divides this Eocene fold-thrust belt into two zones with different structural trends and opposite frontal ramp dips, suggesting the Seattle fault may have originated as a tear fault during accretion.

Washington↗
USGS2025

Application of Hydrologic Simulation Program—FORTRAN (HSPF) as part of an integrated hydrologic model for the Salinas Valley, California

The U.S. Geological Survey (USGS), in cooperation with the Monterey County Water Resources Agency, conducted studies to help evaluate the surface-water and groundwater resources of the Salinas Valley study area, consisting of the entire Salinas River watershed and several smaller, adjacent coastal watersheds draining into Monterey Bay. The Salinas Valley study area is a highly productive agricultural region that depends on the coordinated use of surface water and groundwater to meet demand for irrigation and public water supply. To continue to meet these demands, a better understanding of the historical water balance and the effects of water-resource development on the long-term sustainability of water resources in the Salinas Valley study area is needed.

California↗
USGS2026

Recycling radiogenic osmium by crustal foundering in subduction zones: Evidence from pyroxenite xenoliths in the north Andean arc

Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts

Earth and Planetary Science Letters↗
USGS2026

Miocene evolution of the Humboldt Current

Diatom records from the East Pisco Basin (EPB) of southern Peru and offshore Ocean Drilling Program (ODP) Hole 682 A reveal stepwise increases in the primary productivity of the Humboldt Current during the middle and late Miocene. Although diatoms are present back through the late middle Eocene, successively enhanced diatom production occurs during the Miocene in four steps. The first step between 14.2 and 13.8 Ma marks the onset of diatom deposition in the Pisco-0 sequence. A second step at ∼12.9 Ma coincides with a major drop in global sea level during which diatom deposition ceased in the EPB but continued in ODP 682 A as well as in Chile, Ecuador, and Colombia. Beginning at 10.4 to 10.0 Ma, a major increase in diatom sedimentation rates in both the EPB and in ODP 682 A signals an abrupt intensification of the Humboldt Current productivity, which may be related to the closure of the Central American Seaway to deeper water circulation. A fourth step coincides with the transgressive deposition of the diatom-rich Pisco-2 sequence at 8.4 Ma and is enhanced by the onset of the global Late Miocene Biogenic Bloom (8.0–4.5 Ma). Between 7 and 6 Ma, common subtropical diatoms in ODP 682 A contrast with high diatom deposition rates in the EPB, suggesting variable El Niño-La Niña conditions. During this time, fossil vegetation evidence from southern Peru supports the presence of seasonal periods of enhanced rainfall typical of El Niño conditions. The Miocene trend of increasing sea mammal diversity in the EPB parallels these steps, further supporting stepwise enhancement of primary productivity in the Humboldt Current during the Miocene.

Humboldt Current↗