Geophysical Journal International

SummaryWe refined the Attenuated ProPagation of Local Earthquake Shaking (APPLES) ground-motion-based earthquake early warning (EEW) approach, and directly compare APPLES performance with that of the source-characterization-based U.S. ShakeAlert EEW system for a suite of historical earthquakes in the U.S. West Coast and Japan. APPLES is an extension of the Propagation of Local Undamped Motion (PLUM) algorithm in which observed shaking intensity at seismic stations is used to forward-predict intensity distributions to surrounding areas using an attenuation model derived from an intensity prediction equation. We test new configuration options within APPLES, such as using the second highest estimated ground motion rather than the maximum, to better match median ground-motion observations and reduce alerts for small magnitude earthquakes, both of which are key alerting priorities within ShakeAlert. We evaluate these configurations alongside ShakeAlert by comparing the ground-motion estimation accuracy and available warning times relative to station observations and ShakeMap distributions. Our preferred APPLES configuration produces accurate ground-motion estimates and corresponds better with median observations compared to ShakeAlert’s estimates. This preferred configuration substantially reduces alert issuance for M < 5.0 earthquakes compared to the previous APPLES configuration, and alert-release criteria can further restrict alerts to primarily M ≥ 5.5 earthquakes without requiring magnitude estimation. Prioritizing matching median-observed ground motions may reduce APPLES warning times compared to configurations that were tuned to avoid missed alerts (such as those that use the maximum estimated ground motions), which can lead to shorter warning times compared to ShakeAlert for the same alert threshold. However, station-based warning time assessments demonstrate that APPLES can outperform ShakeAlert for high target thresholds. APPLES is a simple, independent EEW approach that may improve the robustness of EEW for the West Coast of the U.S.

SummaryMetallic infrastructure, such as steel sheet situated within landfills, poses significant challenges to accurate tracking of leachate using induced polarization (IP) methods. The application of IP method is efficient to delineate leakage; however, the presence of metallic structures can cause an interference on the survey and generate high-chargeability anomalies as observed in field survey. To comprehensively validate the interference caused by steel sheets, both numerical and empirical field tests were conducted. As expected, both results demonstrate that interference diminishes as the distance between survey line and metallic structure increases. Additionally, at consistent intervals, the chargeability values inverted using integral chargeability (IC) exhibit a monotonic increase with depth. Moreover, the interference induced by metallic structures is also affected by the controlling factors (i.e. depth, width and thickness) of the structure alongside the intrinsic resistivity and chargeability. Strategic utilization of the size, chargeability, and spatial positioning of metallic structures relative to survey lines can significantly enhance background polarization. This approach offers a promising framework for improving the spatial resolution of subsurface targets exhibiting low polarization effects. The optimization of survey line placement, which must consider the dimensions and electrical properties of metallic structures such as steel sheets, is essential for accurately characterizing landfill leachate using the IP method.

SummaryNorthern Europe experiences vertical land motion and sea level changes that deviate from the average as a consequence of past changes in ice sheet cover in Fennoscandia and the British Isles. The process, called Glacial Isostatic Adjustment (GIA), is controlled by the subsurface structure. Numerical models of GIA can be compared to observations of uplift or past sea level changes to constrain the subsurface structure, and such models can also be used to correct present-day sea level observations to reveal sea level changes due to climate change. GIA models for northern Europe usually adopt a homogeneous upper mantle viscosity even though seismic studies indicate contrasting elastic lithosphere thickness and upper mantle structure between Northwestern Europe and Eastern Europe. This raises the question whether the effect of lateral variations in structure (3D viscosity) can be detected in observations of GIA and whether including such variations can improve GIA model predictions. In this study we compare model output from a finite element GIA model with 3D viscosity to observations of paleo sea level and current vertical land motion. We use two different methods to derive 3D viscosities, based on seismic velocity anomalies and upper mantle temperature estimates. We use three different reconstructions of the Eurasian ice sheet, one based on an inversion using a 1D model, and two others based on glacial geology and modelling. When we use these two reconstructions, we find that the data are fit better using 3D viscosity models. Models with two separate 1D viscosities for Fennoscandia and for the British Isles cannot replicate a 3D model because a 3D model redistributes GIA-induced stresses differently from a combination of models with separate 2D viscosities. The fit to data across Fennoscandia is improved when, as indicated by seismic models, the upper mantle viscosity is higher than for the rest of Northern Europe. The best fit is obtained with a model with dry olivine rheology, in agreement with other evidence from Fennoscandia.

SummaryIn both onshore and offshore seismic exploration, seismic source localization plays a crucial role in ensuring operational safety and environmental protection. With the continuous advancement of the Marchenko method in the fields of seismic migration and internal multiple elimination, this paper investigates a seismic source localization method based on the Marchenko method, aiming to further extend application domain of this method. The key to this method lies in the data reconstruction based on convolution operations. The conventional Marchenko method is then applied to obtain a seismic profile, which includes the location of the seismic source. In the experiments, this study first uses an anticline model to simulate seismic source localization in onshore seismic exploration. The results show that the proposed method can accurately estimate both the distance to the seismic source and its depth. Furthermore, in large-scale marine model experiments, the method is also able to reliably determine the distance between the seismic source and the observation stations.

SummaryFor the interpretation of Spectral Induced Polarization spectra, the determination of the Relaxation Time Distributions (RTD) can be useful, for instance to extract the grain size distribution. However, this is an ill-posed problem, and retrieving the RTD often requires regularization during the inversion process. In this note, we use Bézier curves and simulated annealing to determine the RTD. The procedure that does not require any regularization nor smoothing, by reducing the number of parameters thanks to Bézier curves which are intrinsically continuous and infinitely derivable. We successfully applied our methodology to three examples (Cole-Cole model, Davidson-Cole model, and an experimental spectrum), demonstrating its interest and efficiency.

SummaryThe extension of direct current resistivity methods to induced polarization methods has enriched the tools available for subsurface exploration. This enrichment involves an increase in the number of parameters used in the models, as well as addressing different physical phenomena than those observed with direct current. Accounting for non-linearities, if they exist, can further enhance the sophistication of our models. Non-linearities are often observed, particularly in laboratory experiments. However, we question their origin, and the experiment described here suggests that the non-linearities observed under typical experimental conditions may be artifacts related to the electrodes, rather than reflecting the actual response of the subsurface. Indeed, we first replaced the polarizable injection electrodes with non-polarizable electrodes. The non-linearities observed due to the presence of harmonics were significantly reduced. Then, we replaced the voltage control with a current control, which completely eliminated the non-linearities still present.We know that it is impossible to prove the non-existence of a phenomenon that does not exist. This fundamental epistemological principle (as pointed out by Russell and Popper) means that we are not claiming that nonlinearity does not exist. We are simply describing an experiment that can raise doubts about its existence.

SummaryWe present and validate an efficient GPU-accelerated solver for seismic wave propagation in three-dimensional elastic media. The solver achieves up to a 372× speedup relative to a CPU implementation and supports forward simulations on grids ranging from 100 million to 1 billion cells. It is based on a velocity-stress, first-order formulation of the elastodynamic wave equation and supports kilometer-scale models with layered isotropic and anisotropic structure. We validate the solver by comparing synthetic seismograms to analytical solutions from a propagator matrix method in axisymmetric media. Simulations include moment-tensor sources for a 2017 nuclear explosion and collapse in North Korea, and a magnitude ∼4 earthquake near Linthal, Switzerland (6 March 2017). Anisotropic effects for the Swiss event are modeled using rotated orthorhombic stiffness tensors derived from laboratory measurements of gneiss. Projection onto orthorhombic symmetry enables solver compatibility. We find that anisotropy changes waveform polarity, amplitude, and phase at near-source stations. Unscaled laboratory values produce polarity reversals, while velocity-rescaled tensors correct them. These results demonstrate the impact of anisotropy on waveform modeling and indicate that simplified 1D isotropic models may be insufficient for complex crustal settings. We review how structural effects, including anisotropy and 3D heterogeneity, contribute to transverse-component energy in the 2017 DPRK explosion and discuss implications for seismic source classification.