USGS Science⌕ Search

SEARCH · USGS Science

Results for “Deep Sea Research Part II: Topical Studies in Oceanography”

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

Timing and geometry of the Chemehuevi Formation reveal a late Pleistocene sediment pulse into the Lower Colorado River

The Chemehuevi Formation is a distinctive 50−150-m-thick wedge-shaped Pleistocene sedimentary unit deposited by the Colorado River. It lines the perimeters of the river’s floodplains and bedrock canyons for more than 600 km between the mouth of the Grand Canyon and the delta region in the Gulf of California. The formation is composed of a basal tan to light-yellowish-brown and pale-orange mud-dominated facies overlain and interbedded by a light-yellow-brown sand-dominated facies. The unit is one of two extensively exposed aggradational packages in the Lower Colorado River corridor, in addition to a series of other smaller alluvial terrace deposits. The Chemehuevi Formation appears to represent the response of a fully integrated Colorado River system to a significant perturbation, in contrast to the Bullhead Alluvium, which is likely a unique result of Pliocene river integration. The aggradation of the Chemehuevi Formation in the Lower Colorado River corridor may be similarly due to a unique event in the Colorado River system, or it may instead be a well-preserved sedimentary sequence recording typical behavior of the Colorado River below the Grand Canyon in the late Pleistocene. As such, multiple causal mechanisms have been proposed, but no study to date has conclusively explained the Chemehuevi Formation. To help resolve its timing, duration, and origin, we applied post-infrared infrared stimulated luminescence, carbonate U-Th series, and zircon sensitive high-resolution ion microprobe U-Th series geochronology to determine the ages of key exposures of the unit over a wide spatial area. These new data demonstrate that the Chemehuevi Formation was deposited ca. 110−90 ka. The depositional ages collectively overlap, suggesting that deposition occurred rapidly relative to the resolution of the geochronometers. The new depositional timing coincides with a shift from glacial to interglacial conditions after the marine isotope stage 5-6 transition. This observation is consistent with a climate-induced sediment pulse as a causal mechanism, yet correlations with similar deposits in the Colorado River headwaters or in neighboring catchments appear elusive. Potentially, climate transitions between glacial and interglacial periods induced a sediment pulse from hillslopes of the Colorado River system that resulted in the Chemehuevi Formation. An alternative or additional explanation is that the Chemehuevi Formation represents release of lava dam−impounded sediment in the Grand Canyon. The surface geometry of the Chemehuevi Formation projects upstream to the approximate location of lava dams, and the largest possible lava dam impoundment (the Upper Prospect dam) is comparable in volume to the formation. The lava dam hypothesis appears to be a possible explanation for the Chemehuevi Formation. However, tying deposition to a specific lava dam or series of lava dams remains challenging due to discrepancies in timing and volume. The combined effects of a series of lava dams may have led to the Chemehuevi Formation, as the last Pleistocene lava dam eruption coincides with the onset of deposition. Alternatively, the formation may result from the combined effects of both regional climate transitions and the lava dams that created a transient reservoir to compound a climate transition−driven sediment pulse. The geochronologic data presented here do not allow us to distinguish between the lava dam or climate transition hypotheses but will need to be reconciled with any future proposed depositional model.

Arizona, California, Nevada↗
USGS2024

Spatial distribution of API gravity and gas/oil ratios for petroleum accumulations in Upper Cretaceous strata of the San Miguel, Olmos, and Escondido Formations of the south Texas Maverick Basin—Implications for petroleum migration and charge history

The Maverick Basin of south Texas is currently undergoing active exploration and production of gas and oil from tight sandstone reservoirs. The most productive tight sandstones in the basin are in the Upper Cretaceous San Miguel, Olmos, and Escondido Formations. These units are second only to the Eagle Ford Shale in terms of cumulative production volumes. The structural history of the Maverick Basin, from rifting to subsidence to exhumation, has had a profound effect on the characteristics of these reservoirs and the petroleum resources contained therein. This U.S. Geological Survey review of the production history of these strata reflects a recent shift from conventional production to horizontal drilling (unconventional) that exploits low permeability reservoirs in previously overlooked areas of existing oil and gas fields in southern Texas, typically outside of established field boundaries. To investigate the physical properties of the Maverick Basin hydrocarbon accumulations, this case study compiled American Petroleum Institute (API) gravity measurements and calculated cumulative gas/oil ratios (GOR) for thousands of producing wells from the San Miguel, Olmos, and Escondido Formations. Maps were generated from the compiled well production data to show the spatial heterogeneity of API gravity and GOR values for the three formations within the Maverick Basin and immediately outside the basin to the northeast. Within the Maverick Basin, the spatial patterns of API gravity values indicate lighter oils downdip towards the southern basin edge. GOR values indicative of wet and dry gases within the basin are seen interspersed, with values that correspond to black and heavy oils. Differences in the spatial patterns of the petroleum properties within the Maverick Basin are interpreted as effects of Eocene basin inversion caused by Laramide orogenic deformation, and the resulting reservoir exhumation of basin strata. East of the Maverick Basin, spatial distributions of API gravity and GOR values show progressively heavier oils updip to the northwest, grading to dry gases downdip to the southeast, which correlates to the oil and gas windows of the underlying Eagle Ford Shale. Correlation of API gravity and GOR values from the San Miguel, Olmos, and Escondido Formations with thermal maturity data from the Eagle Ford Shale suggests that the Eagle Ford Shale is the petroleum source, and that petroleum migration was approximately vertical for areas to the east of the Maverick Basin. The discontinuity of API gravity and GOR properties within the Maverick Basin implies a complex petroleum charge history, possibly involving the remigration of petroleum and the addition of petroleum from other source intervals in Mexico, to the southwest. Depressurization of exhumed, overpressured reservoirs of the San Miguel, Olmos, and Escondido Formations can explain the intermittent occurrence of gas production throughout the southern Maverick Basin by exsolution of gas from formation brines and the resulting dry gas flushing of hydrocarbon-charged reservoirs. The introduction of dry gas through flushing can, in turn, explain why the patterns of API gravity and GOR values are so dissimilar in the Maverick Basin. This process has implications for possible future production of unconventional resources from undiscovered tight-gas reservoirs in strata of the San Miguel, Olmos, and Escondido Formations, and a different approach to petroleum exploration may be needed in the Maverick Basin relative to exploration techniques applied in other basins within the northern Gulf of Mexico.

Texas↗
USGS2025

Origin and evolution of mafic volcanism associated with 3 m.y. of andesite production at the Goat Rocks volcanic cluster, southern Washington Cascade Range

More than 3 m.y. of mafic volcanism near the Goat Rocks volcanic cluster in the southern Washington Cascade Range, USA, lends insight into the evolution of basalts and the subarc mantle at a long-lived, major arc volcanic locus. We contribute field observations, 40 Ar/ 39 Ar dates, paleomagnetic directions, and bulk rock and mineral compositions to characterize nine mafic units that erupted in association with the Goat Rocks volcanic cluster. The time frame of mafic volcanism, ca. 3.6 Ma to 60 ka, encompasses the lifespan of the central volcanic cluster (3.1 Ma to 115 ka), with a lull from ca. 2.7 Ma to 1.4 Ma. A climactic period of voluminous mafic activity and far-traveled lava flows, including construction of the Hogback Mountain shield volcano, coincided with voluminous andesite eruptions from the central volcanic cluster. The basaltic rocks in the Goat Rocks area are calc-alkaline to barely tholeiitic and have high field strength element depletion relative to large-ion lithophile elements characteristic of calc-alkaline basalts (CAB) of the Cascade volcanic arc. Unlike at neighboring andesitic volcanic centers (Mounts Adams, St. Helens, and Rainier), no other mafic end members such as high-aluminum olivine tholeiite (HAOT) or intraplate-type basalt (IPB) are present at or near the Goat Rocks volcanic cluster, although some of the calc-alkaline basalts in this study have IPB-like affinities. The Goat Rocks mafic units exhibit two main temporal trends in composition: (1) the most primitive basalts erupted earlier, compared to less primitive and more evolved compositions later, and (2) high field strength element concentrations are higher in the younger basalt units relative to the oldest two. In contrast to these temporal trends, the mafic units define two compositional groups that recur through time, a low-Sr and a high-Sr group, each with distinct trace element and Sr and Nd isotope ratios. Although radiogenic isotope ratios are generally aligned with High Cascades CAB and HAOT, some extend toward IPB of Mount Adams and Simcoe Mountains volcanic field. Olivine-dominated crystal fractionation at shallow pressure from a small range of parent magma compositions accounts for much of the variation among the basalts and basaltic andesites. A high-pressure fractionation model is plausible for only one of the youngest basalt units (basalt of Walupt Lake volcano). Mafic recharge and crustal assimilation accounts for the incompatible-element enriched composition of basaltic andesites erupted during construction of the largest andesitic centers, further supporting sustained basalt mass flux and thermal energy driving andesite genesis. We model the most primitive members of the Goat Rocks mafic units as partial melts of successively less depleted mantle in time. Variable degrees of fluxing with fluids and melts from subduction explain the distinction between high-Sr and low-Sr groups. We propose that mantle metasomatism by ancestral subduction and fluid-flux melting is heterogeneously distributed through the local subarc mantle and played a greater role in the genesis of the high-Sr basalt group. The limited range of primitive basalt types around the Goat Rocks volcanic cluster contrasts with the much greater diversity of basalts throughout the southern Washington to northern Oregon Cascade arc. On the other hand, the central volcanic cluster encompasses nearly the entire diversity observed at neighboring composite volcanoes. In the case of the Goat Rocks area at least, and perhaps attributable to the entire region, this means that the genesis of diverse intermediate magmas is independent from and does not require vastly different parental basalt compositions.

California, Oregon, Washington↗
USGS2025

Timing of hydrothermal alteration and Au-Sb-W mineralization, Stibnite-Yellow Pine district, Idaho

The Stibnite-Yellow Pine district of central Idaho was mined from the early 1900s until the 1990s, extracting gold, antimony, tungsten, and mercury from veins and disseminated and replacement ores in mountainous terrain along the headwaters of the Salmon River. Mining during the two World Wars supplied critical antimony and tungsten to the war efforts. Recent exploration has delineated mineral resources of over 187 metric tons Au, 274 metric tons Ag, and 93,000 metric tons Sb. Mineralization is hosted in Cretaceous Idaho batholith granitic rocks and a sequence of Neoproterozoic to Paleozoic metasedimentary strata of carbonate and siliciclastic compositions. Historical studies outlined some of the complex paragenesis but debated the absolute age of mineralization. New petrographic and geochronologic work documents a sequence of five hydrothermal events in the Stibnite-Yellow Pine district. Event 1 is related to Cretaceous magmatic and hydrothermal activity and includes events ranging in age from 86 to 75 Ma, including sparse quartz-molybdenite veins dated at 86 Ma. Disseminated gold mineralization of event 2 is associated with sericitic alteration and sulfidation of igneous biotite and replacement of plagioclase by potassium feldspar, largely in granodiorite. Gold is present in zoned arsenian pyrite in both disseminated ores and in crosscutting carbonate-quartz veins containing pyrite and arsenopyrite. The large Yellow Pine deposit, localized at a dilatant bend in the Meadow Creek fault, hosts such disseminated and vein gold. Event 2 is interpreted as the major gold-forming event; 40 Ar/ 39 Ar ages of sericite and potassium feldspar alteration range, respectively, from 70 to 59 and 66 to 56 Ma. The long span is interpreted to reflect the age of gold mineralization and local overprinting by event 3. A narrower range from 66 to 61 Ma is interpreted to date the peak of gold mineralization and alteration. Event 3, tungsten mineralization with scheelite, is texturally later than event 2 gold and localized along the Meadow Creek structure. Event 3 scheelite has been dated by isotope dilution-thermal ionization mass spectrometry (ID-TIMS) and laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U-Pb methods at 57 Ma. Event 4, best developed in the West End area, includes gold-silver–bearing quartz-carbonate-pyrite veins and breccias with epithermal textures and potassium feldspar alteration envelopes. Adularia from event 4 yields 40 Ar/ 39 Ar plateau ages of 52 to 51 Ma. Event 5 antimony and mercury mineralization consists of stibnite veins and breccia cements at the Yellow Pine and Hangar Flat deposits as well as cinnabar veins and replacements at the peripheral Fern and Hermes deposits; it is constrained by an LA-ICP-MS U-Pb date on scheelite (ca. 47 Ma) intergrown with stibnite. Minor propylitic and argillic alteration is evident in 47 Ma igneous dikes, which do not contain economic mineralization. The Au-Sb-W ores in the Stibnite-Yellow Pine mining district formed over an extended time period from about 70 to 45 Ma in multiple pulses that were localized along the Meadow Creek fault zone. Each event corresponds to episodes of magmatism and/or hydrothermal activity in the region. Insignificant event 1 skarn and molybdenum mineralization is similar in age to the Thompson Creek porphyry molybdenum deposit in central Idaho. Event 2 gold mineralization occurred during a magmatic gap in central Idaho but was synchronous with magmatism in the Bitterroot lobe further north; event 2 is similar in age to orogenic gold-arsenic mineralization at the Beartrack mine in eastern Idaho. Event 3 scheelite mineralization coincides with tungsten mineralization at the Quartz Creek deposit, late magmatism in the Bitterroot lobe, and rapid exhumation of the Atlanta lobe of the Idaho batholith. Event 4 gold mineralization is coincident with the onset of regional Challis magmatism and extension. Event 5 antimony and mercury mineralization is time-equivalent to epithermal gold mineralization in the nearby Thunder Mountain volcanic field and the peak of Challis magmatism.

Idaho↗
USGS2024

Volcanic gas monitoring

Introduction As magma rises through the crust, decreasing pressure conditions allow volatiles to exsolve from the magma. These volatiles then migrate upward through the crust, where they can be stored at shallower levels or escape to the atmosphere. Rising magma also heats rock masses beneath volcanic centers, causing water in shallow aquifers and hydrothermal systems to boil and release additional gases and steam (see chapter F , this volume; Ingebritsen and Hurwitz, 2024). The chemistry and quantity of gases that reach the surface during periods of quiescence or volcanic unrest can reveal that gas-rich magma is ascending, crystallizing, or alternatively stalling, with important implications for volcanic hazard (for example, Sutton and others, 1992; Aiuppa and others, 2007, 2021; Werner and others, 2009, 2011, 2012; Moretti and others, 2013; de Moor and others, 2016; Lewicki and others, 2019; Edmonds and others, 2022; Kern and others, 2022; Kunrat and others, 2022). Most volcanoes in Alaska and the western United States are characterized by weak degassing, with one or more low-temperature fumaroles (typically near the local boiling temperature of water) and connect to a deeper and sometimes extensive hydrothermal system (for example, McGee and others, 2001; Symonds and others, 2003a, b). Hydrothermal systems will affect the chemistry of rising gases exsolved from deeper magma (Symonds and others, 2001), including sulfur dioxide (SO 2 ), hydrogen chloride (HCl), and water vapor (for example, Doukas and Gerlach, 1995; Gerlach and others, 1998, 2008; Symonds and others, 2001; Werner and others, 2013). As an example, depending on factors such as temperature, pressure, and oxidation state, rising SO 2 will react with groundwater to form hydrogen sulfide (H 2 S) gas, dissolved sulfate (SO 4 2− ), or elemental sulfur (Christenson, 2000; Symonds and others, 2001; Werner and others, 2008). The reaction and dissolution of SO 2 into shallow groundwater is commonly referred to as scrubbing, and can reduce the likelihood that ascending, degassing magma can be detected. Carbon dioxide, however, in addition to exsolving from magma early in the ascent process, is not easily removed by hydrothermal fluids (Lowenstern, 2001). As scrubbing and other processes take place, the SO 2 /H 2 S, CO 2 /SO 2 , and CO 2 /H 2 S ratios may change. High rates of SO 2 emission indicate that magma has moved to relatively shallow levels in the volcano and that the system has heated up enough to establish dry pathways from depth to the surface. Monitoring multiple gas species and the total output of those species is thereby useful for volcano monitoring during both periods of quiescence, to establish background degassing conditions, and during unrest, when gas geochemistry and emission rates can provide information on changing conditions, such as magma ascent. To provide context for multidisciplinary volcano forecasts, we focus on the following two key required capabilities: (1) characterizing baseline geochemistry and gas discharge from volcanoes and volcanic regions and (2) monitoring changes in gas geochemistry and discharge to inform forecasts of volcanic eruptions and their effects. Sufficient baseline data must be collected to identify and interpret anomalous degassing associated with volcanic unrest (for example, Sorey and others, 1998; Rouwet and others, 2014). Differences in volcano type, baseline degassing rates, local hydrology, and geography (for example, high versus low latitude) will result in a different baseline for each volcano. Volcanoes of any threat level that exhibit one or more degassing phenomena would ideally be monitored by techniques needed to establish baseline degassing data, with the sampling frequency of baseline data dictated by the threat level (table E1). Additional monitoring techniques become necessary during periods of unrest. In general, three of the most important techniques for gas monitoring are (1) direct sampling of fumarole, spring, and soil gases for laboratory geochemical measurements, (2) measurements of the chemical composition of the volcanic plume and emission rates of major gas species (for example, H 2 O, CO 2 , SO 2 , and H 2 S) by satellite, airborne, or ground-based techniques, and (3) measurements of diffuse emissions of CO 2 and other gases through soils. Various methods and instruments may be useful both for baseline studies and during unrest.

Scientific Investigations Report↗
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↗
USGS2024

Dissolved arsenic concentrations in surface waters within the upper portions of the Klamath River Basin, Oregon and California

Arsenic toxicity is an environmental health problem. Levels of arsenic in surface waters at some locations in the Klamath River Basin in southern Oregon and northern California can exceed the U.S. Environmental Protection Agency (EPA) standard for drinking water. There are both anthropogenic and natural sources of arsenic. The Klamath River Basin consists primarily of volcanic deposits and contains an underground geothermal system with hot springs and warm water wells, all known natural sources of arsenic. Anthropogenic sources of arsenic are related to the agricultural use of herbicides, fungicides, and insecticides. Surface water arsenic levels can also be affected by fertilizer amendments, evaporative concentration, oxygen-level depletion, and various geochemical transformations that can increase arsenic mobilization. In this study by the U.S. Geological Survey and the Bureau of Reclamation, dissolved concentrations of arsenic, copper, and lead were measured in surface waters at 39 unique sites within the upper portions of the Klamath River Basin between 2018 and 2022. In every year, except 2022, sites were sampled four times between April and November. Surface-water arsenic concentrations varied up to four-orders of magnitude among sites. Median arsenic concentration was lowest at Cherry Creek (0.03 micrograms per liter [μg/L]) and highest at Wood Kimball Spring (36.7 μg/L), two sites located north of Upper Klamath Lake. The highest arsenic concentrations (17.4±4.9 μg/L, n =3) were found in drain sites (defined here as a waterbody returning used irrigation water) while the lowest arsenic concentrations were found in an artesian well (0.8 μg/L, n =1). The elevated arsenic concentrations of the drain sites suggest that arsenic might be concentrated or mobilized by agricultural activities, water re-use practices, and (or) by geochemical processes occurring around water stored in drains (that is, in the water column and across sediment water boundaries). A source of arsenic in drain water in the Klamath Strait Drain area includes water used for irrigation originating from Ady Canal. Other potential sources include groundwater, geothermal water, and local soils and sediments. Seasonal differences in surface-water arsenic concentrations were detected at 13 sites, 10 of which had higher arsenic concentrations in summer than in either spring or fall. The sites sampled around Upper Klamath Lake, the impounded rivers, one of the two canal sites, and 5 of the 14 river sites had higher surface-water arsenic concentrations in the summer than in either spring or fall. Surface-water arsenic concentrations from groundwater sources (that is, springs and in the artesian well) did not vary significantly among seasons (p-values greater than 0.1). Median surface-water concentrations of copper and lead ranged from 0.03 to 3.7 μg/L, and from 0.013 to 0.175 μg/L ( n =2–18), respectively. Dissolved concentrations of both metals were below acute toxicity endpoints reported by the EPA for freshwater animals. Surface-water arsenic concentrations varied independently from corresponding changes in surface-water lead or copper concentrations. However, arsenic concentrations measured in bed-sediment samples collected from a subset of sites located north of Upper Klamath Lake correlated strongly and significantly with the corresponding sedimentary lead concentrations ( p =0.015). Aqueous arsenic speciation measured in a subset of sites in 2019 and 2022 showed that all the arsenic existed as arsenic (V), the most oxidized arsenic species, and presumably, the least toxic. The highest proportions of arsenite (As(III)), the presumably most toxic arsenic species, relative to total arsenic concentrations were found at drain sites. Our assessment of dissolved arsenic concentrations in various surface-water bodies in the Upper Klamath River Basin reveals geographical areas of consistently low (below 2.1 μg/L), moderate (below 10 μg/L) and high (above 10 μg/L) surface-water arsenic concentrations. South of Upper Klamath Lake, surface-water arsenic concentrations were consistently higher than 20 μg/L at two drain sites located in an area of predominant agricultural land use with extensive water re-use practices. North of Upper Klamath Lake, surface-water arsenic concentrations greater than 20 μg/L were consistently measured at sites with limited nearby agricultural activities, suggesting a geogenic source. The consistently high arsenic levels from the Wood River at Jackson F. Kimball State Park, Fort Creek, and Crooked Creek, which are sites located at or near headwater spring sources, suggest a natural background source of arsenic. Water flowing downstream from this area could be a potential source of arsenic to Upper Klamath Lake and the Upper Klamath River.

California, Oregon↗
USGS1991

Reconnaissance of the occurrence of agricultural chemicals in ground water in Haywood, Lake, Obion and Shelby Counties, Tennessee

Data on the occurrence of agricultural chemicals in ground wafer in Tennessee are sparse. The surficial alluvial aquifer is an important source of domestic water supply in West Tennessee, and potentially is subject fo contamination from the application of agricultural chemicals in the area. Nineteen shallow wells completed in the alluvial aquifers in areas of high density agricultural use were sampled in the winter and again in the summer of 1988 to ascertain the occurrence of agricultural chemical in ground water. Although no triazine herbicides or organophosphorus insecticides were detected in any of the wells sampled, elevated nitrite plus nitrate (as nitrogen) concentrations were detected. Results from the winter sampling period indicate a range of nitrite plus nitrate (as nitrogen) concentrations of less than 0.1 to 7.8 milligrams per liter with a median concentration of 2.6 milligrams per liter. Results from the summer sampling period indicate a range of nitrite plus nitrate (as nitrogen) concentrations of less than 0.1 to 8.9 milligrams per liter, median, 2.5 milligrams per liter. The highest concentrations occurred in the shallowest wells, and, in one instance, in a shallow well near a heavily irrigated field.

Tennessee↗
USGS2025

Long-term trends in microseismicity during operational shut-ins at the Coso Geothermal Field, California

Pausing injection and production can lead to induced seismicity in a variety of settings, with some of the largest events occurring during these so-called shut-ins. In geothermal fields, shut-ins are periodically conducted for maintenance on wells and surface infrastructure, thereby offering recurring means of estimating stress changes in the subsurface that lead to increased seismicity rates. Here, we examine daily production and microseismicity data from the Coso Geothermal Field (CGF) in California between 1996 and 2010 to investigate the repetitive trends of operational shut-in microseismicity. Using 24 local seismic stations, we first analyze spatial and temporal trends of over 60,000 earthquakes with magnitudes between -0.4 to 3.8. We find that the northern region exhibits no significant seismicity changes during shut-ins, whereas the rest of the field experiences induced seismicity during almost every shut-in with an increasing intensity towards the southern and eastern portions of the field, highlighting local differences in stress within the CGF. Additionally, we cluster the seismicity using waveform cross-correlation, revealing several earthquake clusters primarily occurring during shut-in periods. These observations suggest that certain fracture and fault sections respond quicker to changes in pore pressure and poroelastic stresses within the geothermal system, possibly highlighting main fluid pathways.

California↗
USGS2026

Monitoring changes in Landsat thermal features in urban and non-urban interfaces from 1986 to 2023 in two international urban centers: Implications for climate and global issues

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.

Remote Sensing↗
USGS2025

Estimating the hypothetical endowment of critical minerals and other commodities in porphyry copper mine waste in the Four Corners states, USA

Society is fundamentally dependent upon commodities that are used in end-use products for the aerospace, defense, energy, telecommunication, and transportation sectors, resulting in centuries of mining to supply these commodities and materials. Waste from these mining operations can remain on the landscape indefinitely, but there is a lack of national understanding of the distribution and scale of such waste features. The renewable energy transition will continue to increase demand for critical minerals and will result in increasing volumes of mine waste on the Earth’s surface. Reprocessing mine waste can reduce environmental risks and recover needed commodities to match growing demand for societal growth. Therefore, understanding the approximate abundance of commodities that may be available for recovery within mine waste features can be an important piece of domestic critical mineral supply.

Arizona, Colorado, New Mexico, Utah↗
USGS2025

Earthquake stress-drop values delineate spatial variations in maximum shear stress in the Japanese forearc lithosphere

Earthquake stress drop (Δσ) may increase with depth and stress in the brittle lithosphere. However, the range of uncertainty in Δσ and the lack of constraints on absolute stress make it difficult to establish whether they are correlated. Here, we investigate Δσ dependence on depth and maximum shear stress ( τ max ) based on ~11 years of seismicity in the northeastern Japanese forearc following the 2011 Tohoku-Oki megathrust earthquake. We interpret Δσ estimates computed using both individual spectra and spectral-ratio methods and find that Δσ exhibits a clear depth dependence within the seismically active upper ~60 km of the forearc lithosphere ( ~ 0.8 MPa per 10 km). We further compare Δσ values with quantitative τ max estimates from finite-element models of force balance. We find that median Δσ values increase with τ max in the brittle forearc lithosphere and that earthquake stress release is proportional to τ max . The dependence of Δσ on τ max explains the apparent depth dependence of Δσ and suggests that average Δσ values provide a relative measure of the stress at failure. In the northeastern Japanese forearc, Δσ values remained roughly constant in the decade following the Tohoku-Oki earthquake, suggesting negligible changes in failure stress in the forearc since the mainshock.

Communications Earth and Environment↗
USGS2024

Deformation by pressure solution and grain boundary sliding in a retrograde shear zone in southern New England, USA

Alleghanian phyllonites in a shear zone in southern New England were formed by the retrogression and hydration of a high-grade Acadian pelitic schist. The retrogression was locally incomplete, resulting in both heterogeneous mineralogy and mineral compositions, and in many arrested reaction textures. These features, documented by backscattered electron and element map imaging, and by 40 Ar/ 39 Ar dating, provide a record of the processes active as the schist was converted to a phyllonite. On the margins of the shear zone relic biotite and garnet grains survive but are partially replaced by chlorite. Relic high-grade muscovite flakes containing up to 35 % paragonite (pg) also persist but are entrained in swarms of pg-poor muscovite flakes producing folia with an anastomosing network of muscovite grains of variable composition. Single crystals of pg-rich muscovite are truncated by muscovite in successively younger cross-cutting folia that contain decreasing pg contents, falling to <5% pg in the youngest folia. In the core of the shear zone all high-grade minerals are destroyed, but recrystallization continues with earlier greenschist facies folia replaced by lower grade and pg-poorer muscovite in younger folia. Muscovite grains in truncated folia are commonly kinked and selective recrystallization of muscovite in kink bands to pg-poor-compositions demonstrates that strain energy helped drive recrystallization. The boundaries of these new pg-poor muscovite folia are sharp, and truncations indicate that grain boundary sliding (GBS) involved dissolution. GBS also enabled folding by flexural slip along muscovite grain boundaries. Fold amplitudes increased as chlorite dissolved from fold limbs precipitated in fold hinges producing crescent-shaped saddle-reef-like structures. Together these observations of truncation and replacement by dissolution and precipitation demonstrate that pressure solution and GBS facilitated both the formation of the phyllonitic shear zone and slip along it.

Connecticut↗
USGS2026

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines↗
USGS1999

Effects of snowmobile use on snowpack chemistry in Yellowstone National Park, 1998

Snowmobile use in Yellowstone National Park has increased substantially in the past three decades. In areas of greatest snowmobile use, elevated levels of by-products of gasoline combustion such as ammonium and benzene have been detected in snowpack samples. Annual snowpacks and snow-covered roadways trap deposition from local and regional atmospheric emissions. Snowpack samples representing most of the winter precipitation were collected at about the time of maximum annual snow accumulation at a variety of locations in the park to observe effects of a range of snowmobile traffic levels. Concentrations of organic and inorganic compounds in snow samples from pairs of sites located directly in and off snow-packed roadways used by snowmobiles were compared to concentrations in samples collected at nearby off-road sites. Concentrations of ammonium were 2 to 5 times higher for the in-road snow compared to off-road snow for each pair of sites. Thus, concentrations decreased rapidly with distance from roadways. In addition, concentrations of ammonium, nitrate, sulfate, benzene, and toluene in snow were positively correlated with snowmobile use.

Wyoming↗
USGS2025

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗