Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Applied Geophysics”

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.

At least 55 records · Page 3Linked to original sources
USGS1999

The color of the Martian sky and its influence on the illumination of the Martian surface

The dust in the atmosphere above the Mars Pathfinder landing site produced a bright, red sky that increases in redness toward the horizon at midday. There is also evidence for an absorption band in the scattered light from the sky at 860 nm. A model of the sky brightness has been developed [ Markiewicz et al. , this issue] and tested against Imager for Mars Pathfinder (IMP) observations of calibration targets on the lander. The resulting model has been used to quantify the total diffuse flux onto a surface parallel to the local level for several solar elevation angles and optical depths. The model shows that the diffuse illumination in shadowed areas is strongly reddened while areas illuminated directly by the Sun (and the blue forward scattering peak) see a more solar-type spectrum, in agreement with Viking and IMP observations. Quantitative corrections for the reddening in shadowed areas are demonstrated. It is shown quantitatively that the unusual appearance of the rock Yogi (the east face of which appeared relatively blue in images taken during the morning but relatively red during the afternoon) can be explained purely by the changing illumination geometry. We conclude that any spectrophotometric analysis of surfaces on Mars must take into account the diffuse flux. Specifically, the reflectances of surfaces viewed under different illumination geometries cannot be investigated for spectral diversity unless a correction has been applied which removes the influence of the reddened diffuse flux.

Journal of Geophysical Research E: Planets↗
USGS2010

Introduction to the JEEG Agricultural Geophysics Special Issue

Near-surface geophysical methods have become increasingly important tools in applied agricultural practices and studies. The great advantage of geophysical methods is their potential rapidity, low cost, and spatial continuity when compared to more traditional methods of assessing agricultural land, such as sample collection and laboratory analysis. Agricultural geophysics investigations commonly focus on obtaining information within the soil profile, which generally does not extend much beyond 2 meters beneath the ground surface. Although the depth of interest oftentimes is rather shallow, the area covered by an agricultural geophysics survey can vary widely in scale, from experimental plots (10 s to 100 s of square meters), to farm fields (10 s to 100 s of hectares), up to the size of watersheds (10 s to 100 s of square kilometers). To date, three predominant methods—resistivity, electromagnetic induction (EMI), and ground-penetrating radar (GPR)—have been used to obtain surface-based geophysical measurements within agricultural settings. However, a recent conference on agricultural geophysics (Bouyoucos Conference on Agricultural Geophysics, September 8–10, 2009, Albuquerque, New Mexico; www.ag-geophysics.org) illustrated that other geophysical methods are being applied or developed. These include airborne electromagnetic induction, magnetometry, seismic, and self-potential methods. Agricultural geophysical studies are also being linked to ground water studies that utilize deeper penetrating geophysical methods than normally used.

Journal of Environmental & Engineering Geophysics↗
USGS2017

On extracting sediment transport information from measurements of luminescence in river sediment

Accurately quantifying sediment transport rates in rivers remains an important goal for geomorphologists, hydraulic engineers, and environmental scientists. However, current techniques for measuring long-time scale (102–106 years) transport rates are laborious, and formulae to predict transport are notoriously inaccurate. Here we attempt to estimate sediment transport rates by using luminescence, a property of common sedimentary minerals that is used by the geoscience community for geochronology. This method is advantageous because of the ease of measurement on ubiquitous quartz and feldspar sand. We develop a model from first principles by using conservation of energy and sediment mass to explain the downstream pattern of luminescence in river channel sediment. We show that the model can accurately reproduce the luminescence observed in previously published field measurements from two rivers with very different sediment transport styles. The model demonstrates that the downstream pattern of river sand luminescence should show exponential-like decay in the headwaters which asymptotes to a constant value with further downstream distance. The parameters from the model can then be used to estimate the time-averaged virtual velocity, characteristic transport lengthscale, storage time scale, and floodplain exchange rate of fine sand-sized sediment in a fluvial system. The sediment transport values predicted from the luminescence method show a broader range than those reported in the literature, but the results are nonetheless encouraging and suggest that luminescence demonstrates potential as a sediment transport indicator. However, caution is warranted when applying the model as the complex nature of sediment transport can sometimes invalidate underlying simplifications.

Journal of Geophysical Research F: Earth Surface↗
USGS2016

A new strategy for earthquake focal mechanisms using waveform-correlation-derived relative polarities and cluster analysis: Application to the 2014 Long Valley Caldera earthquake swarm

In microseismicity analyses, reliable focal mechanisms can typically be obtained for only a small subset of located events. We address this limitation here, presenting a framework for determining robust focal mechanisms for entire populations of very small events. To achieve this, we resolve relative P and S wave polarities between pairs of waveforms by using their signed correlation coefficients—a by-product of previously performed precise earthquake relocation. We then use cluster analysis to group events with similar patterns of polarities across the network. Finally, we apply a standard mechanism inversion to the grouped data, using either catalog or correlation-derived P wave polarity data sets. This approach has great potential for enhancing analyses of spatially concentrated microseismicity such as earthquake swarms, mainshock-aftershock sequences, and industrial reservoir stimulation or injection-induced seismic sequences. To demonstrate its utility, we apply this technique to the 2014 Long Valley Caldera earthquake swarm. In our analysis, 85% of the events (7212 out of 8494 located by Shelly et al . [2016]) fall within five well-constrained mechanism clusters, more than 12 times the number with network-determined mechanisms. Of the earthquakes we characterize, 3023 (42%) have magnitudes smaller than 0.0. We find that mechanism variations are strongly associated with corresponding hypocentral structure, yet mechanism heterogeneity also occurs where it cannot be resolved by hypocentral patterns, often confined to small-magnitude events. Small (5–20°) rotations between mechanism orientations and earthquake location trends persist when we apply 3-D velocity models and might reflect a geometry of en echelon, interlinked shear, and dilational faulting.

California↗
USGS2012

The detection of coastal-trapped waves

We outline a simple method for estimating the cross-spectral matrix of coastal-trapped wave amplitudes,A, from a set of oceanographic observations. Specifically, we propose that A may be estimated by (M'M) -1 M'UM(M'M) -1 where a prime denote conjugate transpose, U is the sample cross-spectral matrix of observations and M is a matrix which has the spatial form of the waves for columns. In general, M will be complex and frequency-dependent. We discuss the bias of this estimator and show how to estimate the variance of the power and cross spectra of wave amplitudes. We also outline an ad hoc scheme for assessing the predictive skill of the coastal trapped wave representation and finally give some advice on how to interpret A. Although the method is presented in the context of shelf circulation and coastal trapped waves, it may be applied to any linear system where the spatial forms of the waves are known and the cross-spectral matrix of their amplitudes is required.

Journal of Geophysical Research C: Oceans↗
USGS2021

Assessment of a claimed ultra-low frequency electromagnetic (ULFEM) earthquake precursor

The rate of occurrence of anomalous ultra-low frequency electromagnetic (ULFEM) pulses has been claimed to have increased days to weeks prior to the M5.4 2007 and M4.0 2010 Alum Rock earthquakes. We re-examine the previously reported ultra-low frequency (ULF: 0.01–10 Hz) magnetic data recorded at a QuakeFinder site located 9 km from the earthquake hypocentre, and compare to data from a nearby Stanford-USGS site located 42 km from the hypocentre, to analyse the characteristics of the pulses and assess their origin. Using pulse definitions and pulse-counting algorithms analogous to those previously reported, we corroborate the increase in pulse counts before the 2007 Alum Rock earthquake at the QuakeFinder station, but we note that the number of pulses depends on chosen temporal and amplitude detection thresholds. These thresholds are arbitrary because we lack a clear physical model or basis for their selection. We do not see the same increase in pulse counts before the 2010 Alum Rock earthquake at the QuakeFinder or Stanford-USGS stations. In addition, the majority of pulses in the QuakeFinder data and Stanford-USGS data do not match temporally, indicating the pulses lack a common origin and are not from lightning or solar-driven ionospheric/magnetospheric disturbances. Our assessment of the temporal distribution of pulse counts shows pulse counts increase during peak human activity hours, suggesting these pulses result from local cultural noise and are not tectonic in origin. The many unknowns about the character and even existence of precursory earthquake pulses means that standard numerical and statistical tests cannot easily be applied. Yet here we show that exhaustive investigation of many different aspects of ULFEM signals can be used to properly characterize their origin.

California↗
USGS2024

Hyperspectral imaging predicts differences in carbon and nitrogen status among representative biocrust functional groups of the Colorado Plateau

Biological soil crusts (biocrusts) are widespread soil photosynthetic communities covering about 12% of Earth's land surface and play crucial roles in terrestrial carbon (C) and nitrogen (N) cycles, yet scalable quantifications of biocrusts and their biogeochemical contributions are notably lacking. While remote sensing has enormous potential to assess, scale, and contextualize biocrusts and their functions, the applicability of hyperspectral data in predicting C- and N-related biocrust traits remains largely unexplored. We address this issue by evaluating the potential of in situ hyperspectral data to predict C and N across a range of biocrust species and different environmental conditions. We found that in situ hyperspectral reflectance measurements can be used to predict biocrust tissue C/N ratios and N concentrations with relatively high accuracy but to a lesser extent for potential biocrust N2 fixation rates. Critical wavelength domains included the visible region of the spectrum from roughly 490–600 nm, which most effectively captured variations in biocrust tissue C, and the shortwave infrared region from 1,150 to 1,350 nm and 1,550–1,650 nm, which most effectively captured biocrust tissue N and N2 fixation potential. Finally, we provide evidence that multi- and hyperspectral missions with targeted band placement, such as the proposed 26-band Landsat Next, could be effective in predicting biocrust traits. This work provides a critical step in understanding how to apply data from new and upcoming satellite missions to the monitoring of biocrusts.

Utah↗
USGS2010

The perfect debris flow? Aggregated results from 28 large-scale experiments

Aggregation of data collected in 28 controlled experiments reveals reproducible debris-flow behavior that provides a clear target for model tests. In each experiment ∼10 m 3 of unsorted, water-saturated sediment composed mostly of sand and gravel discharged from behind a gate, descended a steep, 95-m flume, and formed a deposit on a nearly horizontal runout surface. Experiment subsets were distinguished by differing basal boundary conditions (1 versus 16 mm roughness heights) and sediment mud contents (1 versus 7 percent dry weight). Sensor measurements of evolving flow thicknesses, basal normal stresses, and basal pore fluid pressures demonstrate that debris flows in all subsets developed dilated, coarse-grained, high-friction snouts, followed by bodies of nearly liquefied, finer-grained debris. Mud enhanced flow mobility by maintaining high pore pressures in flow bodies, and bed roughness reduced flow speeds but not distances of flow runout. Roughness had these effects because it promoted debris agitation and grain-size segregation, and thereby aided growth of lateral levees that channelized flow. Grain-size segregation also contributed to development of ubiquitous roll waves, which had diverse amplitudes exhibiting fractal number-size distributions. Despite the influence of these waves and other sources of dispersion, the aggregated data have well-defined patterns that help constrain individual terms in a depth-averaged debris-flow model. The patterns imply that local flow resistance evolved together with global flow dynamics, contradicting the hypothesis that any consistent rheology applied. We infer that new evolution equations, not new rheologies, are needed to explain how characteristic debris-flow behavior emerges from the interactions of debris constituents.

Journal of Geophysical Research F: Earth Surface↗
USGS2023

Flow development and entrainment in turbulent particle-laden jets

Explosive eruptions expel volcanic gases and particles at high pressures and velocities. Within this multiphase fluid, small ash particles affect the flow dynamics, impacting mixing, entrainment, turbulence, and aggregation. To examine the role of turbulent particle behavior, we conducted an analogue experiment using a particle-laden jet. We used compressed air as the carrier fluid, considering turbulent conditions at Reynolds numbers from approximately 5,000 to 20,000. Two different particles were examined: 14-μm diameter solid nickel spheres and 13-μm diameter hollow glass spheres. These resulted in Stokes numbers between 1 and 35 based on the convective scale. The particle mass percentage in the mixture is varied from 0.3% to more than 20%. Based on a 1-D volcanic plume model, these Stokes numbers and mass loadings corresponded to millimeter-scale particle diameters at heights of 4–8 km above the vent during large, sustained eruptions. Through particle image velocimetry, we measured the mean flow behavior and the turbulence statistics in the near-exit region, primarily focusing on the dispersed phase. We show that the flow behavior is dominated by the particle inertia, with high Stokes numbers reducing the entrainment by more than 40%. When applied to volcanic plumes, these results suggest that high-density particles can greatly increase the probability of column collapse.

Journal of Geophysical Research: Atmospheres↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Enhanced visualization for the interpretation of Magellan radar data: Supplement to the Magellan special issue

The differences of radar data from more familiar photographic images, coupled with the unique geologic characteristics of Venus's surface, can make interpretation of the standard Magellan data products difficult for both planetary scientists and nonspecialists. We describe a set of digital processing techniques for transforming individual remote‐sensing datasets in order to make the information they contain more apparent, and for combining multiple datasets of diverse resolution and information content into synthetic images that allow the viewer to explore correlations among the data. The creation of derived data products by techniques like those we describe is becoming increasingly widespread and important, both for scientific research and for use in educational and media presentations. We apply our processing techniques to Magellan synthetic aperture radar images, altimetry, and microwave emissivity measurements of four areas that exemplify some of the most important styles of geologic activity on Venus.

Journal of Geophysical Research E: Planets↗
USGS1999

Stress sensitivity of fault seismicity: A comparison between limited-offset oblique and major strike-slip faults

We present a new three-dimensional inventory of the southern San Francisco Bay area faults and use it to calculate stress applied principally by the 1989 M = 7.1 Loma Prieta earthquake and to compare fault seismicity rates before and after 1989. The major high-angle right-lateral faults exhibit a different response to the stress change than do minor oblique (right-lateral/thrust) faults. Seismicity on oblique-slip faults in the southern Santa Clara Valley thrust belt increased where the faults were undamped. The strong dependence of seismicity change on normal stress change implies a high coefficient of static friction. In contrast, we observe that faults with significant offset (>50–100 km) behave differently; microseismicity on the Hayward fault diminished where right-lateral shear stress was reduced and where it was undamped by the Loma Prieta earthquake. We observe a similar response on the San Andreas fault zone in southern California after the Landers earthquake sequence. Additionally, the offshore San Gregorio fault shows a seismicity rate increase where right-lateral/oblique shear stress was increased by the Loma Prieta earthquake despite also being clamped by it. These responses are consistent with either a low coefficient of static friction or high pore fluid pressures within the fault zones. We can explain the different behavior of the two styles of faults if those with large cumulative offset become impermeable through gouge buildup; coseismically pressurized pore fluids could be trapped and negate imposed normal stress changes, whereas in more limited offset faults, fluids could rapidly escape. The difference in behavior between minor and major faults may explain why frictional failure criteria that apply intermediate coefficients of static friction can be effective in describing the broad distributions of aftershocks that follow large earthquakes, since many of these events occur both inside and outside major fault zones.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Locating non-volcanic tremor along the San Andreas Fault using a multiple array source imaging technique

Non-volcanic tremor (NVT) has been observed at several subduction zones and at the San Andreas Fault (SAF). Tremor locations are commonly derived by cross-correlating envelope-transformed seismic traces in combination with source-scanning techniques. Recently, they have also been located by using relative relocations with master events, that is low-frequency earthquakes that are part of the tremor; locations are derived by conventional traveltime-based methods. Here we present a method to locate the sources of NVT using an imaging approach for multiple array data. The performance of the method is checked with synthetic tests and the relocation of earthquakes. We also applied the method to tremor occurring near Cholame, California. A set of small-aperture arrays (i.e. an array consisting of arrays) installed around Cholame provided the data set for this study. We observed several tremor episodes and located tremor sources in the vicinity of SAF. During individual tremor episodes, we observed a systematic change of source location, indicating rapid migration of the tremor source along SAF. ?? 2010 The Authors Geophysical Journal International ?? 2010 RAS.

Geophysical Journal International↗
USGS2016

Probabilistic assessment of erosion and flooding risk in the northern Gulf of Mexico

We assess erosion and flooding risk in the northern Gulf of Mexico by identifying interdependencies among oceanographic drivers and probabilistically modeling the resulting potential for coastal change. Wave and water level observations are used to determine relationships between six hydrodynamic parameters that influence total water level and therefore erosion and flooding, through consideration of a wide range of univariate distribution functions and multivariate elliptical copulas. Using these relationships, we explore how different our interpretation of the present-day erosion/flooding risk could be if we had seen more or fewer extreme realizations of individual and combinations of parameters in the past by simulating 10,000 physically and statistically consistent sea-storm time series. We find that seasonal total water levels associated with the 100 year return period could be up to 3 m higher in summer and 0.6 m higher in winter relative to our best estimate based on the observational records. Impact hours of collision and overwash—where total water levels exceed the dune toe or dune crest elevations—could be on average 70% (collision) and 100% (overwash) larger than inferred from the observations. Our model accounts for non-stationarity in a straightforward, non-parametric way that can be applied (with little adjustments) to many other coastlines. The probabilistic model presented here, which accounts for observational uncertainty, can be applied to other coastlines where short record lengths limit the ability to identify the full range of possible wave and water level conditions that coastal mangers and planners must consider to develop sustainable management strategies.

Journal of Geophysical Research C: Oceans↗
USGS1985

Inversion of seismic refraction data in planar dipping structure

A new method is presented for the direct inversion of seismic refraction data in dipping planar structure. Three recording geometries, each consisting of two common-shot profiles, are considered: reversed, split, and roll-along profiles. Inversion is achieved via slant stacking the common-shot wavefield to obtain a delay time—slowness (tau— p ) wavefield. The tau— p curves from two shotpoints describing the critical raypath of refracted and post-critically reflected arrivals are automatically picked using coherency measurements and the two curves are jointly used to calculate velocity and dip of isovelocity lines iteratively, thereby obtaining the final two-dimensional velocity model. This procedure has been successfully applied to synthetic seismograms calculated for a dipping structure and to field data from central California. The results indicate that direct inversion of closely-spaced refraction/wide-aperture reflection data can practically be achieved in laterally inhomogeneous structures.

Geophysical Journal International↗
USGS2013

Seismotectonic framework of the 2010 February 27 M w 8.8 Maule, Chile earthquake sequence

After the 2010 M w 8.8 Maule earthquake, an international collaboration involving teams and instruments from Chile, the US, the UK, France and Germany established the International Maule Aftershock Deployment temporary network over the source region of the event to facilitate detailed, open-access studies of the aftershock sequence. Using data from the first 9-months of this deployment, we have analyzed the detailed spatial distribution of over 2500 well-recorded aftershocks. All earthquakes have been relocated using a hypocentral decomposition algorithm to study the details of and uncertainties in both their relative and absolute locations. We have computed regional moment tensor solutions for the largest of these events to produce a catalogue of 465 mechanisms, and have used all of these data to study the spatial distribution of the aftershock sequence with respect to the Chilean megathrust. We refine models of co-seismic slip distribution of the Maule earthquake, and show how small changes in fault geometries assumed in teleseismic finite fault modelling significantly improve fits to regional GPS data, implying that the accuracy of rapid teleseismic fault models can be substantially improved by consideration of existing fault geometry model databases. We interpret all of these data in an integrated seismotectonic framework for the Maule earthquake rupture and its aftershock sequence, and discuss the relationships between co-seismic rupture and aftershock distributions. While the majority of aftershocks are interplate thrust events located away from regions of maximum co-seismic slip, interesting clusters of aftershocks are identified in the lower plate at both ends of the main shock rupture, implying internal deformation of the slab in response to large slip on the plate boundary interface. We also perform Coulomb stress transfer calculations to compare aftershock locations and mechanisms to static stress changes following the Maule rupture. Without the incorporation of uncertainties in earthquake locations, just 55 per cent of aftershock nodal planes align with faults promoted towards failure by co-seismic slip. When epicentral uncertainties are considered (on the order of just ±2–3 km), 90 per cent of aftershocks are consistent with occurring along faults demonstrating positive stress transfer. These results imply large sensitivities of Coulomb stress transfer calculations to uncertainties in both earthquake locations and models of slip distributions, particularly when applied to aftershocks close to a heterogeneous fault rupture; such uncertainties should therefore be considered in similar studies used to argue for or against models of static stress triggering.

Maule↗
USGS2014

Predictions of barrier island berm evolution in a time-varying storm climatology

Low-lying barrier islands are ubiquitous features of the world's coastlines, and the processes responsible for their formation, maintenance, and destruction are related to the evolution of smaller, superimposed features including sand dunes, beach berms, and sandbars. The barrier island and its superimposed features interact with oceanographic forces (e.g., overwash) and exchange sediment with each other and other parts of the barrier island system. These interactions are modulated by changes in storminess. An opportunity to study these interactions resulted from the placement and subsequent evolution of a 2 m high sand berm constructed along the northern Chandeleur Islands, LA. We show that observed berm length evolution is well predicted by a model that was fit to the observations by estimating two parameters describing the rate of berm length change. The model evaluates the probability and duration of berm overwash to predict episodic berm erosion. A constant berm length change rate is also predicted that persists even when there is no overwash. The analysis is extended to a 16 year time series that includes both intraannual and interannual variability of overwash events. This analysis predicts that as many as 10 or as few as 1 day of overwash conditions would be expected each year. And an increase in berm elevation from 2 m to 3.5 m above mean sea level would reduce the expected frequency of overwash events from 4 to just 0.5 event-days per year. This approach can be applied to understanding barrier island and berm evolution at other locations using past and future storm climatologies.

Louisiana↗
USGS2024

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↗