USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Geochemical Exploration”

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

Taking heat (downstream): Simulating groundwater and thermal equilibrium controls on annual paired air–water temperature signal transport in headwater streams

Headwater stream temperature often exhibits spatial variation at the kilometer-scale, but the relative importance of the underlying hydrogeological processes and riverine perturbations remains poorly understood. In this study, we investigated the relative importance of groundwater (GW) and other processes on downstream annual stream temperature signal characteristics using deterministic heat budget model (HFLUX) scenarios within an idealized stream reach representative of mountainous forested conditions. We summarized annual stream thermal regimes from the relationship of paired sinusoidal air and water temperature signals (amplitude ratio, phase lag, and mean ratio). Results showed that downstream changes in annual temperature depended on the thermal gradient between water and the hypothetical equilibrium temperature (where all heat fluxes sum to zero). GW inflow, riparian shading, and the boundary input signal were the most significant factors affecting downstream annual water temperature signals, while flow volume and channel dimensions impacted how quickly annual temperature signals changed. Effects of GW were dominated by advective rather than conductive heat exchange processes, but conduction played a larger role when GW input was more spatially diffuse. Our results indicated several mechanisms by which local processes may affect stream thermal resilience to disturbances and can help guide management of wildfire and climate change.

Journal of Hydrology

Post-fire soil hydrologic response and recovery in northern California (USA)

Background Wildfires abruptly change landscapes by altering soil properties and vegetation cover. These changes are thought to reduce soil infiltration capacity, making landscapes susceptible to runoff and erosion. However, post-fire soil response is complex and likely varies across locations and time. Aims Here, we aim to understand regional post-fire soil response and recovery by tracking changes across different northern California (USA) lithology and vegetation types. Methods We conducted repeat in situ soil infiltration tests for 3 years post-fire at 31 burned and 10 unburned sites spanning the 2021 Dixie, 2020 LNU Lightning Complex, 2020 Walbridge and 2020 Glass fires. Key results Our two main findings are: (1) burned chaparral soils have increased hydraulic conductivity compared with unburned sites, and (2) infiltration rates return to pre-fire conditions within 3 years across most lithologies and vegetations. Conclusions Recovery might be generalizable by vegetation and lithology but differ regionally, making it important to identify meaningful hydrologic response units (HRUs). Multi-year studies with paired burned and unburned measurements can constrain the recovery timeline and provide information missed by observations solely of burned soils. Implications Understanding where, and for how long, soil remains susceptible to runoff and erosion can help prioritize areas and time periods most in need of mitigation.

California

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington

Temporal, spatial, and chemical evolution of Quaternary high-silica rhyolites in the Mineral Mountains, Utah

The Mineral Mountains in southwestern Utah are a structurally controlled core complex at the confluence of the Colorado Plateau and the Basin and Range physiographic provinces. Aside from hosting Utah’s largest batholith, the Mineral Mountains host some of the State’s youngest high-silica rhyolites, which have been linked to a magma source that is presently being utilized as an enhanced geothermal system. The high-silica rhyolites take the form of effusive lavas and domes, and explosive products are rare. Previous K-Ar dating of these Pleistocene rhyolites placed eruptions between about 790 and 500 kilo-annum (ka) with contemporaneous basalts erupting in the valley to the east of the Mineral Mountains. Large uncertainties on these ages obscured the tempo of eruptions and thus hindered attempts to constrain the timescales of the petrogenetic processes that produced the rhyolites. In this study, we build on previous studies conducted in the 1970s and 1980s by using new geochronologic and geochemical data to investigate the temporal and spatial evolution of the youngest phase of volcanism in the Mineral Mountains. We identify two major eruptive periods, from approximately 850 to 750 ka and from approximately 590 to 480 ka. The older phase is characterized by the eruption of several basaltic lavas, two obsidian flows, and a series of coalescing porphyritic rhyolite domes. The younger phase included the eruption of six evolved high-silica rhyolite domes and one pyroclastic deposit, followed by the eruption of trachyandesite in the adjacent valley to the east. Whole-rock geochemical data indicate that the rhyolites can be divided into three chemical groups, with more evolved compositions erupting through time. The youngest rhyolites along the range crest have the lowest total iron and TiO 2 concentrations and the highest incompatible element concentrations, indicative of increasing differentiation with time and elevation. Improved precision on the eruption ages indicates a recurrence interval of approximately 20 thousand years. The eruptive flux for both periods of rhyolitic volcanism is about 0.01 cubic kilometers per thousand years, which is less than the magma resurgence flux rates for syn-caldera and post-caldera eruptions of the Valles Caldera and Yellowstone Caldera volcanic systems. Collectively, these geochemical, geochronological, and volumetric data may facilitate a better understanding of heat flux and the longevity of magmatic sources related to geothermal resources in similar small-volume, silicic systems.

Utah

Storm impact scale for barrier islands

A new scale is proposed that categorizes impacts to natural barrier islands resulting from tropical and extra-tropical storms. The proposed scale is fundamentally different than existing storm-related scales in that the coupling between forcing processes and the geometry of the coast is explicitly included. Four regimes, representing different levels of impact, are defined. Within each regime, patterns and relative magnitudes of net erosion and accretion are argued to be unique. The borders between regimes represent thresholds defining where processes and magnitudes of impacts change dramatically. Impact level 1 is the 'swash' regime describing a storm where runup is confined to the foreshore. The foreshore typically erodes during the storm and recovers following the storm; hence, there is no net change. Impact level 2 is the 'collision' regime describing a storm where the wave runup exceeds the threshold of the base of the foredune ridge. Swash impacts the dune forcing net erosion. Impact level 3 is the 'overwash' regime describing a storm where wave runup overtops the berm or, if present, the foredune ridge. The associated net landward sand transport contributes to net migration of the barrier landward. Impact level 4 is the 'inundation' regime describing a storm where the storm surge is sufficient to completely and continuously submerge the barrier island. Sand undergoes net landward transport over the barrier island; limited evidence suggests the quantities and distance of transport are much greater than what occurs during the 'overwash' regime.

Journal of Coastal Research

Integrating marine historical ecology into management of Alaska’s Pacific cod fishery for climate readiness

The Pacific cod ( Gadus macrocephalus ) fishery was closed in 2020 after a rapid decline in biomass caused by the marine heat waves of 2014–2019. Pacific cod are exceptionally thermally sensitive and management of this fishery is now challenged by increasingly unpredictable climate conditions. Fisheries monitoring is critical for climate readiness, but short-term monitoring data may be inadequate for recognizing and anticipating change under rapid climate changes. We propose an interdisciplinary, marine historical ecology framework that looks to long-term records (local and traditional knowledge, history, archaeology, and paleoclimatology) to capture a long range of ecological variability and provide historical context for management. In order to connect to contemporary fisheries management, this framework must be built on a common vocabulary and an understanding of the key metrics used in fisheries stock assessments. Here, we propose metrics derived from Pacific cod stock assessment and synthesize information relevant to understanding the effects of past warming periods on cod populations across the Gulf of Alaska and Bering Sea. This case study provides a framework for thinking about how to use these historical records in the context of fisheries management under rapidly changing climate conditions.

Alaska

Borehole geophysical time-series logging to monitor passive ISCO treatment of residual chlorinated-ethenes in a confining bed, NAS Pensacola, Florida

In-situ chemical oxidation (ISCO) is a common method to remediate chlorinated ethene contaminants in groundwater. Monitoring the effectiveness of ISCO can be hindered because of insufficient observations to assess oxidant delivery. Advantageously, potassium permanganate, one type of oxidant, provides the opportunity to use its strong electrical signal as a surrogate to track oxidant delivery using time-series borehole geophysical methods, like electromagnetic (EM) induction logging. Here we report a passive ISCO (P-ISCO) experiment, using potassium permanganate cylinders emplaced in boreholes, at a chlorinated ethene contamination site, Naval Air Station Pensacola, Florida. The contaminants are found primarily at the base of a shallow sandy aquifer in contact with an underlying silty-clay confining bed. We used results of the time-series borehole logging collected between 2017 and 2022 in 4 monitoring wells to track oxidant delivery. The EM-induction logs from the monitoring wells showed an increase in EM response primarily along the contact, likely from pooling of the oxidant, during P-ISCO treatment in 2021. Interestingly, concurrent natural gamma-ray (NGR) logging showed a decrease in NGR response at 3 of the 4 wells possibly from the formation of manganese precipitates coating sediments. The coupling of time-series logging and well-chemistry data allowed for an improved assessment of passive ISCO treatment effectiveness.

Florida

A process-based model for forecasting wave runup along the coast of Georgia

Wave runup is an important nearshore process that impacts total water level, sediment transport, and coastal design. Current methods for forecasting wave runup implement an empirical model that considers offshore wave height, wave period, and generalized beach slope. In this study, the authors generated wave runup forecasts from offshore wave conditions and a system of polynomial equations derived from numerical simulations at three different still water datums for each beach profile. They developed a process-based methodology that incorporated site-specific cross-shore topobathy into the phase-resolving numerical model. A comparison between the system of equations, deterministic hydrodynamic simulations, and observed high-water marks was made using Hurricanes Matthew (2016) and Irma (2017) for 12 cases, and it showed that the polynomials were capable of being consistent with the results from full simulation runs, while not requiring hours of runtime when a forecast was needed—the differences between the polynomial and the observed high water marks ranged from 3 to 32 cm for the Irma hindcast and 9–70 cm for Matthew. Then, using forcings from Hurricanes Ian and Nicole (2022), the model predicted the occurrence of dune collision, overwash, and inundation for the coast of Georgia and suggested that wave runup was impacted by the still water level and local topobathy.

Georgia

Long-term monotonic trends in water budget components in the contiguous United States: Insights from two hydrologic models

Characterizing changes to water availability for domestic, industrial, agricultural, and other uses is essential to support water management. To better quantify these changes, the U.S. Geological Survey and National Science Foundation National Center for Atmospheric Research produced two hydrologic models simulating water budget components from 1980 to 2021 over the contiguous United States (CONUS). Both hydrologic models were driven by a common atmospheric forcing dataset and aggregated to common spatial and temporal scales, which enables a novel evaluation of congruency between the models. We present annual and seasonal trends in six water budget components (precipitation, evapotranspiration, streamflow, groundwater recharge, soil saturation, and snow water equivalent) based on the Mann–Kendall test for monotonic trend and Theil-Sen slope estimate for the water year 1983–2021 period for ~86,000 catchments in CONUS. Additional components and metrics from our analysis pipeline are available in an associated published dataset, which contains more than 46 million trend results. The water budget trends showed broad agreement with prior observational and modeling studies that indicate increasing trends in the northeast and decreasing trends in southwestern CONUS. We found the seasonal variability in water budget trends was greatest in the southern, central, and northwest CONUS. These findings support integrated trend assessments when coupled with trends in water quality and use.

Contiguous United States

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

Evaluation of submersible pressure transducers for streamflow monitoring in small streams

Compact streamgages requiring minimal infrastructure and equipped with submersible pressure transducers (PTs) are increasingly used to monitor small streams, yet disparate implementations obscure their accuracy under real-world conditions. This study isolated instrumentation-derived uncertainty in stage monitoring by co-locating various combinations of commercial vented and unvented PTs with seven U.S. Geological Survey (USGS) reference streamgages on small streams. Multi-year PT stage records from compact streamgages, collected and corrected following USGS protocols, were compared to concurrent reference observations. Vented PTs demonstrated an average measurement uncertainty of ±0.005 m and mean absolute percent error (MAPE) of ±0.2%. Unvented PTs exhibited higher uncertainty, averaging ±0.009 m and ±0.4% MAPE. Although both sensor types had stage errors up to ±0.3 m, 95% of vented PT and unvented PT observations were within 0.01 and 0.02 m of reference stage, respectively. Analysis of additional unvented PTs revealed stage errors of up to ±1.2% when using barometric sensors within 15 km of the in-water sensor. Propagation of stage error to discharge using reference rating models resulted in cumulative discharge MAPEs of ±4.5% for vented and ±5.5% for unvented PTs. These findings highlight PTs as practical alternatives to reference instrumentation when deployed with standardized procedures, potentially expanding access to reliable streamflow data.

conterminous United States

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research

A new water temperature modeling approach to predict thermal habitat suitability for nonnative cichlids in Florida rivers

As global temperatures increase, the spatiotemporal arrangement of thermal habitats in Florida rivers may shift, creating the potential for greater dispersal and establishment of nonnative tropical freshwater fishes. To understand how water temperature changes may affect the spatial distribution of these nonnative species, more effective water temperature prediction models are necessary. Currently, most models employ either a generalized air–water temperature relationship or require expensive and complicated tools to measure hydrometeorological factors (e.g. groundwater input). Thus, we developed a novel modeling approach that is accurate, accessible, and cost-effective in allowing fisheries managers to project water temperatures in rivers across Central and North Florida. To characterize the potential for nonnative fishes to spread northward, we evaluated two hardy and abundant species currently found primarily in South Florida: Mayan Cichlid ( Mayaheros urophthalmus ) and Oscar ( Astronotus ocellatus ). Our results show an increase in thermally suitable winter days for both species in 10 of 11 rivers studied, consistent with predicted water temperature warming under 16 climate-change scenarios spanning different levels of air temperature warming (+1 °C, +2 °C, +3 °C, +4 °C) and precipitation/groundwater thermal sensitivity (0, 0.33, 0.66, 1). Considering resource limitations, fisheries managers can use our water temperature modeling approach to predict effects of climate change on Mayan Cichlid and Oscar survival, growth, and dispersal and take actions to manage potential northward movement of these species.

Florida

PFAS remediation in a bioelectrochemical system inoculated with the west branch consortium (WBC-2)

Groundwater contamination by per - and polyfluoroalkyl substances (PFAS) poses a persistent environmental and public health concern. This study evaluates a two-chambered bioelectrochemical system (BES) inoculated with the West Branch Consortium (WBC-2) for PFAS remediation. Under an applied cathodic potential of −450 mV (versus Ag/AgCl), the BES with active WBC-2 achieved >99.0% perfluorooctanesulfonic acid (PFOS) removal within 21 days in deionized water with culture medium and > 98.9% removal of PFOS, perfluorooctanoic acid (PFOA), perfluorohexanoic acid (PFHxA), and perfluorohexanesulfonic acid (PFHxS) in contaminated groundwater after 102 days. Intermediate formation (e.g., PFOA, 6:2 fluorotelomer sulfonate (6:2 FTS), perfluoropropionic acid (PFPrA), perfluorobutanoic acid (PFBA)) and background-corrected fluoride release were consistent with PFOS transformation under anaerobic reducing conditions potentially involving defluorination. Following repeated PFOS spikes (100 μg/L on Days 0, 50, and 399), PFOA, PFPrA, and PFBA accumulated over 664 days. Despite being the dominant accumulated compound, PFOA accounted for <1.8% of the total spiked PFOS mass. Minimal PFOS transformation occurred in controls without active WBC-2, highlighting the importance of microbial metabolism. Biofilm analysis revealed dense colonization of rod-shaped bacteria on carbon fiber brushes. Enrichment of Bacillus , Agrobacterium , and other low-abundance taxa suggests selective adaptation to BES and PFAS conditions. These findings highlight BES driven by electrochemically stimulated microbial activity as a promising strategy for PFAS remediation.

Journal of Water Process Engineering

Drone-based radiometric surveys provide high-resolution mine waste characterization

Airborne radiometric surveys use passive geophysical techniques to characterize geochemical variations at or near earth’s surface. These methods have been used for a variety of mapping applications, including mineral resource evaluation. However, detailed characterization of smaller geologic targets, including mine waste features, requires flying at lower altitudes and with tighter line spacing than is feasible with traditional aircraft. Here, a small uncrewed aircraft system (sUAS) equipped with a radiometric sensor was used to acquire high-resolution gamma-spectrometry over small mine waste features and a low-grade stockpile in southwestern New Mexico. The sUAS radiometric system mapped local variability within each survey area and revealed ~2–10 m wide zones where radioelements K, Th, and U may be elevated 2–10× the surrounding material. Additionally, the sUAS radiometric data revealed radioelement variability across survey sites, which correlated reasonably well with variability seen in geochemical samples at each survey site, even though samples collected from individual sites showed high local variability. The sUAS data characterized local heterogeneity within mine waste and other small geologic targets at scales of a few meters to tens of meters, which is not possible with traditional crewed aircraft, and with continuity of coverage that is not possible with ground surveys, thus filling a key gap in geophysical survey spatial resolution.

New Mexico

Systematic approach to prioritize wells for effective groundwater monitoring and management in the Arkansas Headwaters Basin, Colorado, USA

Study region The Arkansas Headwaters Basin, an intermountain basin in the Southern Rocky Mountains of North America. Study focus Our specific focus is choosing a set of wells to support a possible future regional groundwater-surface water model that would support water management. We present a three-step process using multiple criteria to score, predict, and choose prioritized wells that capture the full distribution of data including extremes. The three-step process provides accessible visualizations, fiscally efficient well prioritization, and screening useful for subsequent groundwater modeling. The novelty of the proposed methodology is the systematic approach integrating a scoring and a predictive approach to support a selection path. The systematic approach may be broadly adapted for other basins. New hydrological insights for the region Understanding regional hydrology hinges on efficient collection of hydrologic data that captures the relevant dynamics including extremes. The present study, a case study for a particular basin in the Southern Rocky Mountains, is the first use of a scripted (R software) strategy to select an economical and representative set of monitoring wells. Our findings suggest caution when using proximity as a proxy for correlation, because proximal wells in the same geologic formation and similar depths are not always correlated. In the Arkansas Headwaters Basin, subsurface geology may be less influential on groundwater elevations than broader hydrologic influences, such as regional drought.

Colorado

Episodic sediment-discharge events in Cascade Springs, southern Black Hills, South Dakota

Cascade Springs is a group of artesian springs in the southern Black Hills, South Dakota, with collective flow of about 19.6 cubic feet per second. Beginning on February 28, 1992, a large discharge of red suspended sediment was observed from two of the six known discharge points. Similar events during 1906-07 and 1969 were documented by local residents and newspaper accounts. Mineralogic and grain-size analyses were performed to identify probable subsurface sources of the sediment. Geochemical modeling was performed to evaluate the geochemical evolution of water discharged from Cascade Springs. Interpretations of results provide a perspective on the role of artesian springs in the regional geohydrologic framework. X-ray diffraction mineralogic analyses of the clay fraction of the suspended sediment were compared to analyses of clay-fraction samples taken from nine geologic units at and stratigraphically below the spring-discharge points. Ongoing development of a subsurface breccia pipe(s) in the upper Minnelusa Formation and/or Opeche Shale was identified as a likely source of the suspended sediment; thus, exposed breccia pipes in lower Hell Canyon were examined. Upper Minnelusa Formation breccia pipes in lower Hell Canyon occur in clusters similar to the discrete discharge points of Cascade Springs. Grain-size analyses showed that breccia masses lack clay fractions and have coarser distributions than the wall rocks, which indicates that the red, fine-grained fractions have been carried out as suspended sediment. These findings support the hypothesis that many breccia pipes were formed as throats of abandoned artesian springs. Geochemical modeling was used to test whether geochemical evolution of ground water is consistent with this hypothesis. The evolution of water at Cascade Springs could not be suitably simulated using only upgradient water from the Minnelusa aquifer. A suitable model involved dissolution of anhydrite accompanied by dedolomitization in the upper Minnelusa Formation, which is caused by upward leakage of relatively fresh water from the Madison aquifer. The anhydrite dissolution and dedolomitization account for the net removal of minerals that would lead to breccia pipe formation by gravitational collapse. Breccia pipes in the lower Minnelusa Formation are uncommon; however, networks of interconnected breccia layers and breccia dikes are common. These networks, along with vertical fractures and faults, are likely pathways for transmitting upward leakage from the Madison aquifer. It is concluded that suspended sediment discharged at Cascade Springs probably results from episodic collapse brecciation that is caused by subsurface dissolution of anhydrite beds and cements of the upper Minnelusa Formation, accompanied by replacement of dolomite by calcite. It is further concluded that many breccia pipes probably are the throats of artesian springs that have been abandoned and exposed by erosion. The locations of artesian spring-discharge points probably have been shifting outwards from the center of the Black Hills uplift, essentially keeping pace with regional erosion over geologic time. Thus, artesian springflow probably is a factor in controlling water levels in the Madison and Minnelusa aquifers, with hydraulic head declining over geologic time, in response to development of new discharge points. Development of breccia pipes as throats of artesian springs would greatly enhance vertical hydraulic conductivity in the immediate vicinity of spring-discharge points. Horizontal hydraulic conductivity in the Minnelusa Formation also may be enhanced by dissolution processes related to upward leakage from the Madison aquifer. Potential processes could include dissolution resulting from leakage in the vicinity of breccia pipes that are abandoned spring throats, active spring discharge, development of subsurface breccias with no visible surface expression or spring discharge, as well as general areal leakage from the Madison aquifer into the Minnelusa Formation.

South Dakota

Cross-fade sampling: Extremely efficient Bayesian inversion for a variety of geophysical problems

This paper introduces cross-fade sampling, a computationally efficient Markov Chain Monte Carlo simulation method that uses a semi-analytical approach to quickly solve Bayesian inverse problems that do not themselves have an analytical solution. Cross-fading is efficient in two ways. First, it requires fewer samples to obtain the same quality simulation of the target probability density function (PDF). Secondly, it is much faster to evaluate the posterior probability of each sample than conventional sampling methods for simulating Bayesian posterior PDFs. Conventional methods require evaluating the prior probability (which describes your a priori constraints) and data likelihood (which describes the fit between the observations and the predictions of the model) for each sample model. However, cross-fading does not require evaluating the data likelihood, meaning that ‘big data’ can be fit with zero additional computational cost. Further, the cross-fading approach can be used to calculate the marginal likelihood associated with a model design, facilitating model comparison and Bayesian model averaging. Topics covered in this paper include derivation of the cross-fade approach and how it can be used to simulate Bayesian posterior PDFs and compute the marginal likelihood, discussion of the class of problems to which cross-fading can be applied (with examples from earthquake statistics, earthquake ground motion modelling, volcanic eruption forecasting, and finite fault slip modelling), demonstration of efficiency relative to existing sampling methods and discussion of how cross-fading can be used to account for prediction errors (i.e. epistemic errors) as part of the geophysical inverse problem.

Geophysical Journal International