SEARCH · USGS Science

Results for “The Mountain Geologist”

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

SourceResearch connectionFrom yearThrough year

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

1,662 records · Page 27Linked to original sources

USGS2026

Latest Pleistocene to 19th-century earthquakes on bending-moment reverse faults of the Seattle fault zone, Washington

Fault-related folds and their associated secondary faults play a critical yet often underrecognized role in accommodating strain and generating earthquakes in active fold-and-thrust belts. In the Seattle fault zone (SFZ), Washington, USA, we present new paleoseismic, geomorphic, and geophysical evidence for late Pleistocene and Holocene earthquakes on shallow, south-dipping secondary faults—the Lytle Beach and Vasa Park faults—that lie within the hanging wall of the greater SFZ and are on trend with the primary, blind Blakely Harbor fault. Our data show that these structures have ruptured independently, producing localized uplift and deformation at the surface, with the most recent event (RH2) likely occurring in the early nineteenth century. While a temporal overlap between the late Pleistocene RH1 and VP1 earthquakes raises the possibility of a ≥35 km rupture along the Blakely Harbor fault, structural and temporal evidence instead supports independent rupture on individual faults related to folding. We interpret these faults as bending-moment reverse faults that formed within a synclinal hinge zone of the main fault, reflecting mechanical and kinematic influences of the broader fault system. Combined with prior studies, our findings indicate that faulting related to folding dominates the mode of strain release within the SFZ since the late Pleistocene with more frequent earthquake recurrence (∼350 yr) over the past ∼2500 yr.

Washington↗

USGS2025

Near-fault amplification and ground motion variability during the 2019 Ridgecrest, California sequence

We estimate ground-motion variability near the 2019 M 7.1 Ridgecrest earthquake sequence. Accurate seismic hazard estimation requires understanding ground-motion spatial correlations, yet many studies lack the dense station coverage needed to resolve small-scale variability. The 2019 M 7.1 Ridgecrest earthquake sequence presents a unique opportunity to examine ground motions and their spatial correlations at a range of interstation distances. The permanent seismic network was augmented with hundreds of temporary stations including several fault-crossing nodal arrays. We compute the event ( δE i ) and within-event ( δW ij ) residuals from the observed peak ground velocity and peak ground acceleration data to isolate potential sources of ground-motion variability. We then compare δW ij between station pairs that record an event to understand the semivariance of the ground motion versus interstation distance. By fitting an exponential model to the semivariances, we determine a correlation range of 25 km for the Ridgecrest region. Although the exponential model fits the broad-scale increase of semivariance with interstation distance, we also observe smaller-scale trends. We find that ground motions are less correlated for station pairs that are near or across faults that ruptured during the 2019 Ridgecrest sequence. We also find large, positive median δW ij with relative values 2–3 times larger than nearby stations for individual stations’ near-fault traces. Near-fault amplification and greater ground-motion variability can delineate fault zones and may locally increase the seismic hazard.

California↗

USGS2025

Comparative assessment of a restored and natural wetland using 13C-DNA SIP reveals a higher potential for methane production in the restored wetland

Wetlands are the largest natural source of methane (CH 4 ), a potent greenhouse gas produced by methanogens. Methanogenesis rates are controlled by environmental factors such as redox potential, temperature, and carbon and electron acceptor availability and are presumably dependent on the composition of the active methanogen community. We collected intact soil cores from a restored and natural freshwater depressional wetland on Maryland’s Delmarva Peninsula (USA) to assess the effects of wetland restoration and redox shifts on microbial processes. Intact soil cores were incubated under either saturated (anoxic) or unsaturated (oxic) conditions and amended with 13 C-acetate for quantitative stable isotope probing (qSIP) of the 16S rRNA gene. Restored wetland cores supported a distinct community of methanogens compared to natural cores, and acetoclastic methanogens putatively identified in the genus Methanosarcina were among the most abundant taxa in restored anoxic and oxic cores. The active microbial communities in the restored wetland cores were also distinguished by the unique presence of facultatively anaerobic bacteria belonging to the orders Firmicutes and Bacteroidetes . In natural wetland incubations, methanogen populations were not among the most abundant taxa, and these communities were instead distinguished by the unique presence of aerobic bacteria in the phyla Acidobacteria , Actinobacteria , and class Alphaproteobacteria . Iron-reducing bacteria, in the genus Geobacter , were active across all redox conditions in both the restored and the natural cores, except the natural oxic–anoxic condition. These findings suggest an overall higher potential for methanogenesis in the restored wetland site compared to the natural wetland site, even when there is evidence of Fe reduction.

Maryland↗

USGS2025

Extreme plate boundary localization promotes shallow earthquake slip at the Japan Trench

The 2011 Mw9.1 Tohoku-oki earthquake is exceptional among great earthquakes for having peak slip of ~50-70 m on the shallowest portion of the plate boundary megathrust. International Ocean Discovery Program Expedition 405 drilled multiple holes through the megathrust in the large slip region and at a Pacific Plate input site. The megathrust preferentially develops at the top or base of the pelagic clay in the input section where pronounced contrasts in physical properties are present. This results in a narrow, weak fault located at a major mechanical contact between frontal prism mud and subducted clay. Localization imposed by the input section enhances the tendency for shallow seismic slip, showing the Japan Trench may be more susceptible to ruptures with large shallow slip than margins without weak clays.

Japan Trench↗

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↗

USGS2026

Decoding the oxidative digestion mechanism for polystyrene nanoplastic detection in the Great Lakes using a customizable Raman spectral processing algorithm

Despite the concerns surging around nanoplastics (NPs) regarding their prevalence and bioavailability in freshwater systems, robust detection of NPs in complex environmental matrices is hindered by the lack of standardized sample pretreatment and a mechanistic understanding of oxidative digestion. Here, we systematically investigate the interaction between hydrogen peroxide (H 2 O 2 ) and polystyrene (PS) NPs during digestion in deionized (DI) water and four environmental matrices from in and around the Great Lakes Basin. To facilitate high-throughput analysis, we develop Pre_peak, a customizable Raman spectral processing algorithm that achieves >99% accuracy for both NP identification and interference rejection, allowing reliable NP quantification via pixel counting and systematic decoding of the oxidative digestion mechanisms. In DI water, varying H 2 O 2 doses from 0 to 30% has negligible effects on the recovery and Raman signal intensity of PS NPs over 24 hours of digestion. However, morphological changes and aggregation of PS NPs are observed when the H 2 O 2 dose exceeds 20%. Prolonged digestion further leads to progressive NP loss. In natural waters, the optimal dosage and digestion duration depend on matrix characteristics, including dissolved organic matter (DOM) and ion composition. This study provides mechanistic insights into NP–oxidant interactions and underscores the need for matrix-tailored digestion protocols to advance standardized NP detection in freshwater environments.

Great Lakes↗

USGS2025

How do hydrological variability and human activities control the spatiotemporal changes of riverine nitrogen export in the Upper Mississippi River Basin?

Excessive nitrogen export from agricultural watersheds remains a critical water quality challenge, with the Upper Mississippi River Basin (UMRB) significantly contributing to downstream eutrophication and hypoxia in the Gulf. This study investigates the spatiotemporal dynamics of riverine nitrate plus nitrite (NO 3 – + NO 2 – -N) export across the UMRB at high spatial resolution (12-digit Hydrologic Unit Codes or HUC12 subwatershed scale) during 2001–2020 and quantifies the effects of anthropogenic activities and hydrological variability on riverine NO 3 – + NO 2 – -N export changes in the region between 2001–2005 and 2016–2020. Our results revealed hotspots of substantial increases in NO 3 – + NO 2 – -N yields across the UMRB, with distinct regional patterns in driving factors. Over the entire UMRB, NO 3 – + NO 2 – -N yields increased by 9.7 kg/ha/yr on average from 2001–2005 to 2016–2020, with anthropogenic activities contributing 4.8 kg/ha/yr and hydrological variability contributing 4.9 kg/ha/yr. The northern and western UMRB had combined influences from both anthropogenic activities and hydrological variability, while the east-central regions had predominantly hydrologically driven changes. Agricultural sources, including fertilizer, manure, and biological nitrogen fixation, collectively contributed over 80% of NO 3 – + NO 2 – -N loading throughout the basin. This framework for disentangling human and hydrological impacts provides critical insights for developing effective and targeted watershed management strategies to reduce nutrient losses and improve water quality.

Upper Mississippi River Basin↗

USGS2024

What do we know without the catalog? Eliciting prior beliefs from experts for aftershock models

Fitting parametric seismological models to earthquake catalogs often comes with numerical challenges, especially when catalogs are small. An alternative way to quantify parameter values for a seismic region is by eliciting expert opinions on the seismological characteristics that each parameter corresponds to. For instance, expert beliefs on aftershock patterns can be formulated into prior distributions for aftershock parameters, for example, for the epidemic‐type aftershock sequence (ETAS) model. We illustrate such a method by not only eliciting priors for ETAS parameters for the Pacific Northwest (PNW), a subduction zone with a complex tectonic environment, but also a relatively small catalog. We compare these priors with those suggested by the ETAS literature for global subduction zones, discussing implications for aftershock forecasting for the PNW.

The Seismic Record↗

USGS2025

The influence of nodule versus crust morphology on the composition of seamount-hosted ferromanganese minerals

Ferromanganese mineral precipitation in the global oceans is ubiquitous, occurring in the form of both crusts and nodules at a broad range of depths and seafloor terrains. Although ferromanganese crusts and nodules are both composed of ferromanganese minerals, mineralogy and mean element concentrations compiled for regional crust versus nodule occurrences differ. Notably, most published compilations compare nodules from abyssal plain sediments to crusts forming on seamounts, and do not address the question as to whether mineral morphology may affect the composition of crusts and nodules in addition to environmental factors. To address this, we have compiled a dataset of co-located hydrogenetic (seawater sourced) open-ocean ferromanganese crusts and seamount-hosted nodules. The preliminary data presented here support our hypothesis that co-located crusts and nodules exhibit similar compositions, both in terms of their mineralogy and mean element concentrations, as well as overall variance in concentrations for any given element. These data indicate that local processes influence sample-to-sample variation in both crusts and nodules, in contrast to broader ferromanganese variation reported between open-ocean seamounts- hosted crusts and abyssal plain-hosted nodules.

Conference Paper↗

USGS2024

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington↗

USGS2024

Coastal Science Navigator companion guide—Discover the U.S. Geological Survey coastal science products you need

The Coastal Science Navigator is an online gateway to a wide variety of U.S. Geological Survey (USGS) coastal change hazards-related information, data, and tools relevant to stakeholders’ scientific and decision-making needs. The products within the Coastal Science Navigator provide data related to past, present, and future threats to our coastlines. The filter search allows users to see all available products and identify relevant options by time scale, geographic scope, coastal hazard theme, and other filters. The guided search suggests products based on users’ answers to a short series of questions. A comprehensive summary is available for each product. The idea for the Coastal Science Navigator arose in 2020 in response to stakeholder feedback identifying the need for a central source for USGS coastal science information. It was published in July 2023 and initially included 55 products. Regular updates are planned to integrate other existing and new products. This guide introduces some of the many coastal change hazards-related products available through the USGS. In it, we showcase the products included in the Coastal Science Navigator’s initial publication in July 2023. While it is not representative of all the information, tools, and data available, we hope it serves as a compelling snapshot of what the USGS has to offer and encourages you to explore the Coastal Science Navigator to discover more of the products you need. To navigate this guide, the products have been organized by the time scale they are best suited for—past, present, or future—although many products cover multiple time scales. An additional section features software, one of the many product types available as filters within the Coastal Science Navigator. Other products include downloadable data, websites, and geonarratives (web pages that combine text, images, and interactive maps into narratives you can scroll through). Featured geographic scopes are also highlighted within this guide, detailing some of the many regions in which the USGS conducts research and illustrating another way to filter products within the Coastal Science Navigator.

Circular↗

USGS2026

Bottom-up characterization of geologic methane emissions in the San Juan Basin in the southwestern USA

Methane is a potent greenhouse gas that plays an important role in atmospheric chemistry and global warming. The current global methane budget has large uncertainties, and a better understanding of the budget would help to guide strategies for reducing anthropogenic emissions to fight climate change. Natural geologic methane emissions are a particularly poorly constrained source, with top-down estimates from 14 C in ice cores suggesting much lower geologic emissions than bottom-up scaling of direct flux measurements. Our study aims to contribute to resolving this discrepancy through improved bottom-up characterization of geologic methane seepage in the San Juan Basin in southwestern Colorado and northwestern New Mexico, USA. We performed 983 new flux chamber measurements in this basin during summer 2022 and winter 2023 field campaigns. Our results, in combination with prior measurements, suggest that natural seepage in the San Juan Basin only occurs on or near the Fruitland coal outcrop. Specifically, our new measurements confirm previous measurements of seepage along the northwestern exposure of the Fruitland outcrop in Colorado (a known hydrodynamic overpressure region) and for the first time, identified seepage locations along the southernmost Fruitland outcrop exposure in New Mexico, in association with a coal cleat and a fault. Overall, seepage along the Fruitland coal outcrop is heterogeneously distributed, with both positive and negative (interpreted as microbial soil sink) methane fluxes. Features that are hypothesized to be predictive of seepage (e.g., faults) were not associated with positive methane fluxes in areas outside of the Fruitland outcrop. Our best estimate for total geologic methane seepage in the San Juan Basin from spatial interpolation and statistical upscaling is approximately 0.14 Tg CH 4 /yr, with a range from 0.029 to 0.48 Tg CH 4 /yr. This best-estimate value is lower than a previous bottom-up estimate from a gridded seepage inventory, but higher than a previous top-down estimate.

Colorado, New Mexico↗

USGS1976

Diamond drilling at the Ma'milah Mine, Kingdom of Saudi Arabia

The Ma'millah gold mine is about 90 km southeast of At Ta'if at lat 21°03'N., long 41°18'E., in southwestern Saudi Arabia. The deposit, which was worked extensively by ancient miners, consists of several veins and massive lenses of quartz along a regional fault. The extensive ancient mine workings, large size of some quartz outcrops, and significant gold content of several samples indicated that subsurface exploration was warranted. Accordingly, two holes were drilled into the deposit during the interval December 1972 to February 1973. The drill findings were disappointing; the grade of vein material was found to be too low to be of economic interest and the quartz bodies were found to pinch out both along strike and downdip. The drilling program therefore was terminated. The deposit is not worth further exploration

Ma'milah Mine↗

USGS2026

Importance of fish in the diet of Adélie penguins across multiple life stages

Over recent decades, the Adélie penguin ( Pygoscelis adeliae ) population has grown across most of its range, the exception being the northern coast of the western Antarctic Peninsula (WAP). In the Ross Sea, the very large Cape Crozier colony grew to reach ~9% of the global population. Demographic factors driving the Cape Crozier trend have yet to be identified. However, low chick fledging mass - a negative influence in the northern WAP where the diet is mostly less-energy-dense krill - is not showing an effect in the Cape Crozier colony. Previously, we hypothesized that Cape Crozier fledglings, and post-breeding adults, must be finding sufficient higher-quality prey once free of the prey-depleted colony foraging area, with subsequent higher survival. Here we explore penguin diet within and outside the colony foraging area using stable isotope analysis (SIA) of feathers: those grown by post-breeding adults outside the foraging area, prior to moult, compared to those of near-fledged chicks raised on meals obtained within the foraging area. Second, we explore whether a major difference occurs between breeding and post-breeding diet of adults from a small, nearby colony (Cape Royds) that exploits a smaller, little-depleted preyscape. Finally, we compare the results to SIA analyses published by others to explore whether the pre-moult diet of Ross Island adults, mostly foraging within the north-eastern Ross Sea, compares with that of adults nesting in the southern WAP (Bellingshausen Sea coast). At the latter, populations are not decreasing, with diet containing an appreciable contribution of fish. The results confirmed that fish contributed greatly to chick diets at both Ross Island colonies (δ 15 N 10.6–11.4). Once adults were released from central-place foraging, SIA levels (δ 15 N 9.7–12.9) were similar to those of the southern WAP. Foraging on energy-dense fish in coastal Antarctic waters could be one reason Adélie penguins are among the most abundant and therefore most ecologically successful penguin species.

Antarctic Science↗

USGS2024

Isotopic evaluation of the National Water Model reveals missing agricultural irrigation contributions to streamflow across the western United States

The National Water Model (NWM) provides critical analyses and projections of streamflow that support water management decisions. However, the NWM performs poorly in lower-elevation rivers of the western United States (US). The accuracy of the NWM depends on the fidelity of the model inputs and the representation and calibration of model processes and water sources. To evaluate the NWM performance in the western US, we compared observations of river water isotope ratios ( 18 O / "> 16 O and 2 H / "> 1 H expressed in δ notation) to NWM-flux-estimated (model) river reach isotope ratios. The modeled estimates were calculated from long-term (2000–2019) mean summer (June, July, and August) NWM hydrologic fluxes and gridded isotope ratios using a mass balance approach. The observational dataset comprised 4503 in-stream water isotope observations in 877 reaches across 5 basins. A simple regression between observed and modeled isotope ratios explained 57.9 % ( δ 18 O) and 67.1 % ( δ 2 H) of variance, although observations were 0.5 ‰ ( δ 18 O) and 4.8 ‰ ( δ 2 H) higher, on average, than mass balance estimates. The unexplained variance suggest that the NWM does not include all relevant water fluxes to rivers. To infer possible missing water fluxes, we evaluated patterns in observation–model differences using δ 18 O diff ( δ 18 O obs − δ 18 O mod ) and d diff ( δ 2 H diff - 8 ⋅ δ 18 O diff "> ). We detected evidence of evaporation in observations but not model estimates (negative d diff and positive δ 18 O diff ) at lower-elevation, higher-stream-order, arid sites. The catchment actual-evaporation-to-precipitation ratio, the fraction of streamflow estimated to be derived from agricultural irrigation, and whether a site was reservoir-affected were all significant predictors of d diff in a linear mixed-effects model, with up to 15.2 % of variance explained by fixed effects. This finding is supported by seasonal patterns, groundwater levels, and isotope ratios, and it suggests the importance of including irrigation return flows to rivers, especially in lower-elevation, higher-stream-order, arid rivers of the western US.

Hydrology and Earth Systems Science↗

USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗

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↗

USGS2024

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas↗