Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Frontiers in Earth Science”

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,612 records · Page 67Linked to original sources
USGS2025

Mineral commodity summaries 2025

Introduction Each mineral commodity chapter of the 2025 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2024 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the MCS 2025 are welcomed.

Mineral Commodity Summaries↗
USGS2026

Mineral commodity summaries 2026

Introduction Each mineral commodity chapter of the 2026 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2025 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the 2026 MCS are welcomed.

Mineral Commodity Summaries↗
USGS2026

Morphodynamics of snow barchans at Concordia Station, Antarctica

Barchans are crescent-shaped dunes that form and migrate under unimodal wind regimes in areas of low sediment availability. While sand barchan morphodynamics has been extensively studied, snow barchans remain poorly documented. Here, using hourly time-series images taken near Concordia Station, Antarctica, we measure the sizes (lengths 0.3 m to 13 m, widths 0.4 m to 7.5 m) and migration rates (0.13 m h −1 to 3.5 m h −1 ) of 44 populations of snow barchans during 2017–2022. We compare this dataset with established height–velocity and length-width relationships for sand barchans in both air and water. Our results show that snow barchans follow the same scaling as sand barchans, allowing us to estimate a characteristic sintering (bonding of snow particles) time of about one day. This sintering process limits the size of snow barchans when active, reduces their aspect ratios, and, once wind falls below the transport threshold, explains their preservation over a few years. During wind events, migrating barchans are too small to maintain their shape and orientation, making them highly responsive to fluctuations in wind direction. Thus, studying snow barchans improves understanding of an elementary bedform encountered in multiple environments and also of the surface properties of snow on different timescales, offering insights into snow accumulation and wind-blown snow in ice cores.

Geomorphica↗
USGS2025

Rupture process of the Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake

The Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake ruptured a km long portion of the east-west trending Mendocino fault zone (MFZ). In order to clarify the rupture process, we assemble three-component seismograms from regional seismic stations, horizontal coseismic displacement vectors derived from Global Navigation Satellite System (GNSS) time series, and a Sentinel-1 ascending interferogram. These data are interpreted with a model of slip distributed on two vertical fault planes representative of the eastern MFZ and spanning the ~70 km length of the aftershock zone. Assuming right-lateral strike slip, we find that the rupture initiates in the oceanic mantle at 20-30 km depth and proceeds unilaterally updip and toward the east. Early aftershocks locate adjacent to the peak slip areas, tracking the coseismic rupture propagation from oceanic mantle to shallower depth and implying a significant role of static stress transfer in driving aftershocks in an ocean plate environment.

California↗
USGS2025

Geochemistry and Soils of the Big Smoky Valley Fens, Nevada

Fens are groundwater-fed wetlands that can provide habitat for plants and animals. Due to anthropogenic activities and climate change, many fens around the world are at risk. This paper presents the results of a study of the hydrology and geochemistry of fens in Big Smoky Valley, central Nevada to support the Bureau of Land Management’s activities in the area. A water sample from the largest fen in the study area was analyzed for its water chemistry and compared to a nearby alluvial aquifer and hot spring. The high SiO 2 concentration of the fen sample implies that the fen water may originate from geothermal water. A soil core was taken to analyze radiocarbon age and soil type. A majority of the core was composed of silt and clay interlayered with water-filled voids. Changes in the character of the clay with depth suggest that there may have been changes in the depositional environment over time. Radiocarbon dating of Ruppia seeds showed longevity of the fen, with the minimum 14 C age of the core as 4,375±40 years. This paper provides reconnaissance-level information on the Big Smoky Valley fens, but further information would be needed to better understand the source of water to the fens or how the fen environment has changed over time with climate.

Nevada↗
USGS2026

Spatially referenced watershed models for the binational Red–Assiniboine River Basin: Bayesian vs frequentist comparison

Excess nutrient loading remains a leading cause of declining water quality in lakes, estuaries, and coastal waters worldwide, with global economic costs of US$200 billion – US$2 trillion annually from impacts on fisheries, tourism, freshwater resources, and water treatment. Our study focuses on total phosphorus (TP) in Lake Winnipeg and its binational Red-Assiniboine River Basin, where nutrient inputs have degraded water quality and increased cyanobacterial blooms. These changes pose ecological, public health, and economic risks. We applied a spatially referenced watershed model with a hybrid statistical-mechanistic structure partitioning annual nutrient loads into land-use export, land-to-water delivery, and in-reservoir decay. Bayesian and traditional frequentist model calibrations were compared. In the frequentist model, coefficients for agricultural inputs, forests /wetlands, stream channels, precipitation, and reservoir losses were statistically significant, whereas coefficient for wastewater was not. In contrast, all variables were successfully calibrated using the Bayesian approach. Model results delineate TP-export hotspots across the basin, showing that 54–62% of TP originates from the U.S., with agricultural sources ranging 62–72%—highlighting the importance of agriculture-focused Best Management Practices. Given the global relevance of nutrient-driven water-quality challenges, our results highlight Bayesian calibration for robust risk assessment and adaptive nutrient management.

Red–Assiniboine River Basin↗
USGS2026

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia↗
USGS2026

Northern bobwhites select for native grasses on working grazing land

Northern bobwhite ( Colinus virginianus ; bobwhite) populations have experienced an 85% decline across most of its range since the 1960s. The most drastic decreases have been in the southeastern United States where biologists attribute the decline to widespread habitat loss, including conversion of native grasslands to tall fescue ( Lolium arundinaceum ) and other exotic cool-season grass (CSG) pasture and hayfields. Including agricultural lands in conservation programs could improve habitat conditions on a regional scale. One working-lands conservation strategy involves the use of native warm-season grasses (NWSGs) rather than cool-season exotic grasses as economically viable cattle forage. To evaluate this management practice for creating bobwhite habitat, we conducted a field study on a bobwhite population in eastern Kentucky, USA, where grazed NWSG and burned NWSG fields were adjacent to exotic grazed and exotic hayed CSG fields. Between April 2019 and October 2022, we used radio telemetry to track bobwhites to evaluate resource selection at second- and third-order scales. Bobwhites used traditionally managed CSGs less than would be expected by chance at the second-order scale in both breeding and non-breeding seasons. Bobwhite use of grazed NWSG and burned NWSG was greater than would be expected by chance at the second-order scale in both breeding and non-breeding seasons. Similarly, at the third-order scale, bobwhites used CSGs less than would be expected by chance, whereas NWSGs and woody stems were used more than by chance alone. Our study suggests that bobwhites select NWSGs at multiple spatial scales despite broad-scale availability of managed exotic CSGs. Under a working-lands framework, the integration of NWSGs into working lands could create bobwhite habitat in the southeastern United States, especially if woody stems are present.

Kentucky↗
USGS2024

Extending the Boore and Abrahamson (2023) modified square-root-impedance method for the development of site amplifications consistent with the full-resonance approach to a range of VS30 values

The square-root-impedance (SRI) method is commonly used to approximate the seismic site amplifications computed using the full-resonance (FR) method for gradient shear-wave velocity ( V S ) profiles that are smoothly varying with depth. The SRI site amplifications have been observed to systematically underpredict the FR site amplifications by a ratio of FR/SRI amplifications around 1.05 to 1.3 across a wide frequency range (Boore, 2013). Recently, Boore and Abrahamson (2023; hereafter, BA23) related this difference in the SRI and FR methods to differences in the exponent η of the ratio of seismic impedances between the two methods. They proposed the implementation of a modified frequency-dependent η in the SRI method to improve its match to the FR site amplifications. This modified η was derived using only five V S profiles. We investigate the performance of the BA23 η for a wide range of realistic gradient V S profiles with V S30 ranging from 180 to 1500 m/s. These gradient V S profiles are constructed using two power-law functions of depth and are constrained by the assigned VS30 value, the depth and velocity of the half-space, and depths to shear-wave velocity horizons of 1.0 and 2.5 km/s ( Z 1.0 and Z 2.5 ) based on western United States sites. Despite observing a V S30 dependence of η, we find that the BA23 η generally works reasonably well for the range of V S profiles analyzed. Using the VS30 -dependent η derived in this study results in improvements in matching the FR site amplification compared to using the BA23 η. These improvements are more pronounced for the soft-site conditions and become modest to negligible for the stiff site conditions

Bulletin of the Seismological Society of America↗
USGS2023

Connecting flood-related fluvial erosion and deposition with vulnerable downstream road-stream crossings

Fluvial erosion is increasingly responsible for infrastructure and building damages associated with floods as the intensity of extreme rainfalls hit rural and urban rivers in a variety of climate settings across the United States. Extreme floods in 2016 and 2018 caused widespread culvert blockages and road failures, including extensive damage along steep tributaries and ravines in the Marengo River, Wisconsin, watershed during 2016 and 2018. A study conducted by the U.S. Geological Survey (USGS), Wisconsin Wetlands Association (WWA), Ashland County, and the Northwest Wisconsin Regional Planning Commission (NWRPC) investigated the special concern of fluvial erosion hazards (FEHs) associated with gullying, streamside landslides, and the loss of wetland storage in headwaters. In 2019, a pilot study was begun to map and classify ephemeral and perennial streams and wetlands in terms of their sensitivity to FEHs. This study combined data from field-based rapid geomorphic assessments (RGAs) coupled with a stream network-wide geographic information system (GIS) approach for mapping stream segments, referred to as fluvial process zones (FPZ), sensitive to erosion, deposition, and channel change. The GIS approach used nationally available 10-meter (m) resolution topology and an extended stream network to map FPZs based on Strahler stream order, stream power, channel slope, presence of adjacent steep valley sides and headwater flats, and adjacent landform setting. Bankfull channel widths derived from RGA-based hydraulic geometry curves combined with drainage areas, an estimate of bankfull flow, and channel slope were used to calculate specific stream power for the FPZs. Lastly, the FPZs were characterized by their location within three major landform settings that affect erosion potential. The resulting vulnerability maps provided a screening framework to identify FPZs that are sensitive to incision, gullying and mass wasting along steep headwater ephemeral channels, as well as downstream perennial channels that have the potential for valley-side landslides, coarse sediment deposition, and channel change. Lastly, each FPZ was characterized in terms of hydrologic alteration associated with ditching. The vulnerability mapping products and rankings of sensitivity of FPZs will ultimately be used by Ashland County and their collaborators to prioritize natural flood management projects that mitigate FEHs, restore hydrology, and reconnect channels with adjacent wetlands and floodplains.

Wisconsin↗
USGS2025

Water-resources inventory and assessment at Katahdin Woods and Waters National Monument

The U.S. Geological Survey, in cooperation with the National Park Service, prepared a water-resources inventory and assessment for Katahdin Woods and Waters National Monument (KAWW). This compilation includes published and publicly accessible hydrologic data and resource assessments of streams, rivers, ponds, lakes, wetlands, vernal pools, and groundwater in and near KAWW. It also includes reports and datasets summarizing attributes of KAWW’s hydrologic infrastructure, such as stream crossings, dams, wastewater discharge plants, groundwater monitoring wells, and U.S. Geological Survey streamflow-gaging stations. Descriptions of data and details of current limitations in available datasets are included. Wetland, groundwater, streamflow, and water-quality information are all limited. Hydrography data are available; however, there are limited ground-truth data. Accurate streamlines within KAWW were developed from light detection and ranging (lidar) as a part of this work. Hydrologic infrastructure information is available from multiple sources; however, differences exist among the datasets. Datasets are summarized in appendix 1.

Maine↗
USGS2024

Simulated ground-motion records for the seismic assessment of monumental masonry structures

Earthquakes are natural disasters that can cause widespread devastation and loss of life. Simulated ground-motion records can be useful in regions with limited seismic stations or a history of damaging but infrequent earthquakes. This is especially true in areas with a high concentration of heritage masonry structures, which are especially susceptible to damage, as simulated records can be crucial in predicting their seismic response and protecting these buildings from seismic damage. Despite the importance of simulated earthquakes, few studies have investigated how effective they are compared to real earthquakes when assessing the structural response of heritage buildings. To address this knowledge gap, we employ two different simulation methods of the Mw 6.2 Faial earthquake, which occurred on July 9th, 1998, in the Azores, to replicate the recorded time-series at four available stations within an epicentral distance of 150 km. The study has two objectives: first, to validate the simulated records of the 1998 Faial earthquake using alternative stochastic ground-motion simulation approaches, and second, to determine how these approaches impact the seismic assessment of historic masonry structures. To accomplish these objectives, this study uses real and simulated ground-motion datasets to conduct non-linear response history analyses of the São Francisco Church, a monumental structure in Horta that sustained damage during the Faial earthquake. The results show that both simulation approaches yield structural responses similar to the observed records.

Conference Paper↗
USGS2026

Rewilding the upper Klamath River Basin: Rapid recolonization of Chinook salmon following the world’s largest dam removal

Habitat fragmentation is a global driver of decline in biodiversity, yet evaluations of ecological response following re-connection remain rare. In 2024, four hydroelectric dams on the Klamath River were removed, reconnecting over 640 km of anadromous habitat and marking the largest dam removal project completed to date. We used a stationary high-resolution multi-beam imaging SONAR to quantify adult Chinook salmon ( Oncorhynchus tshawytscha ) passage into newly reopened reaches during the first two migratory seasons following dam removal. An estimated 7,742 (95% CI: 7,702–7,778) Chinook salmon migrated upstream in 2024 and 13,310 (95% CI: 12,876–13,733) in 2025, representing 18–19% of fish returning to the Klamath River Basin. Within this period, Chinook salmon recolonized 88% of their documented historical distribution, expanding into 345 km of reconnected habitats and reaching elevations of up to 1,250 m. Recolonization required fish to swim through dam removal sites, complex hydraulic corridors and engineered fishways, demonstrating the propensity for rapid recolonization at basin scale despite these challenges. These findings provide the first quantitative assessment of salmon response to the world’s largest dam removal and demonstrate how reconnection can catalyze large-scale recolonization in historically fragmented river networks.

California, Oregon↗
USGS2025

Arctic fold-and-thrust belts

The modern Arctic has been formed through a series of continent–continent collisions, accretion of terranes and phases of crustal extension. The Neoproterozoic Timanian, Paleozoic Caledonian and Uralian, and late Mesozoic Verkhoyansk–Kolyma, Chukotkan and Brookian orogenies formed several large fold-and-thrust belts (FTBs). The FTBs are exposed across vast areas of continents and continue offshore to form a complex tectonic basement for thick sedimentary basins, playing an important role in the history of accumulation and deformation of younger unmetamorphosed sedimentary successions that are the subject of this volume. Recognition of the importance of FTBs in the Arctic geological history and their role as a controlling factor of development of Arctic sedimentary basins resulted in this chapter, in which we review the current state of knowledge about Arctic FTBs and highlight questions that remain to be addressed. Enclosure D , a map showing boundaries of the FTB and their internal first-order structural fabric, is a part of the overview.

Book chapter↗
USGS2026

Incorporating location uncertainty improves inference with stop-level North American Breeding Bird Survey data

Ecological models should account for uncertainty to be most effective and useful. Yet, uncertainty from model covariates—unlike that from other sources, such as sampling error or process variability—is seldom explicitly incorporated. This can cause underestimates of uncertainty to cascade through model parameter estimates, predictions, and downstream uses. Burner et al. proposed a method for quantifying uncertainty in covariates and incorporating it into models using informative Bayesian priors. This method was applied to stop-level Breeding Bird Survey (BBS) analyses, where land cover uncertainty at each stop arises from substantial stop location uncertainty. A limited validation of model-estimated land cover, using stops with known locations, indicated the method’s potential effectiveness, but it was not rigorously evaluated. We conduct a robust simulation-based test, generating stop locations, extracting land cover, and simulating bird communities across 210 BBS routes in the upper Midwest. We compare 3 models: a “known” model with true land cover, a “naive” model assuming consistent 800-m stop spacing, and a “full” model using informative priors to estimate land cover. Species parameter estimates and predicted prevalence patterns across gradients in land cover from the full model approached those of the known model and were substantially closer to the true values used in simulations relative to those from the naive model. Naive model parameters were more biased relative to the other models, and credible intervals of predicted species prevalence rarely included the true simulated values. The full model also produced land cover covariate estimates closer to true simulation values relative to the mean informative priors. Our results show that, for the BBS, informative priors enable more accurate stop-level analyses despite location uncertainty. In contrast, naive models that ignore this uncertainty yield poor inferences. More broadly, we demonstrate empirically the utility of informative priors to account for covariate uncertainty in ecological models.

Michigan, Minnesota, Wisconson↗
USGS2024

Hydroacoustic observations reveal drivers of mixing and salinization of a karst subterranean estuary during intense precipitation

Karst subterranean estuaries within globally ubiquitous carbonate aquifers are coastal groundwater ecosystems that provide an essential water resource for human populations. To understand the drivers of salinization within a coastal aquifer in the Yucatan Peninsula (Mexico), we employed hydroacoustics in flooded caves to observe how oceanic and atmospheric events facilitate mixing between the meteoric lens (fresh-brackish groundwater) and the saline groundwater on tidal and episodic timescales. Precipitation during Tropical Storm Carlotta increased the flow and salinity of the meteoric lens without evidence for vertical mixing across the halocline. We postulate that vertical migration of haloclines in the conduit relative to those within the rock matrix during precipitation creates lateral density gradients that drive mixing, and ultimately creates a brackish layer within the meteoric lens. These results provide a mechanistic explanation for vertical and lateral exchange in a coastal carbonate aquifer, which has implications for groundwater response to future climatic change.

Ox Bel Ha cave system, Yucatan peninsula↗
USGS2025

Analysis of aquifer framework and properties, Alvahs Lane well field, Cutchogue, New York

The U.S. Geological Survey, in cooperation with the Suffolk County Water Authority, evaluated the aquifer transmissivity and storage properties at the Alvahs Lane well field north of the village of Cutchogue, New York. This analysis of aquifer properties provides the Suffolk County Water Authority with hydrogeologic information needed to develop water supplies to meet the increasing water demands of the residents of Suffolk County, New York. An aquifer test was conducted at the Alvahs Lane well field from October 18 through October 21, 2022, when a production well was pumped at 550 gallons per minute for about 24 hours, and groundwater-level drawdown and recovery were measured in two monitoring wells. The three wells are screened in a glaciofluvial aquifer under unconfined (water table) conditions. Drawdown and recovery data were analyzed with an analytical solution for partial penetration and delayed yield in an unconfined aquifer to provide estimates of the glaciofluvial aquifer properties. Inclusion of lateral aquifer boundaries was not necessary for the analysis to result in satisfactory matches with the observed water-level responses. Aquifer transmissivity was estimated at 32,000 feet squared per day. Assuming a saturated aquifer thickness of 120 feet, this result is equivalent to a horizontal hydraulic conductivity value of 270 feet per day. Specific yield was estimated at 0.15 (dimensionless). The estimated properties are consistent with those of a highly transmissive unconfined aquifer.

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