Geophysical Journal International

SummaryThis study introduces a new method for calculating acoustic-gravity waves in a spherically layered atmosphere. The method introduces a model assumption and divides the atmosphere into finely stratified layers to solve the PDE with respect to the radial coordinate. The time-domain synthetic signal is obtained by summing over the orders of the associated Legendre functions and then applying the FFT. The method is applied to numerically simulate wave behaviour, including Earth curvature effects, and compares with the horizontally layered model (HLM). Results show that at near-field distances, our method aligns closely with HLM, but significant differences emerge in the far field, particularly beyond an epicentral distance of 50°, where Earth curvature becomes critical. Our method successfully simulates head waves of seismic phases, and Rayleigh waves, even for waves travelling multiple times around the Earth, which HLM cannot achieve. Simulations using a homogeneous Earth model reveal head wave characteristics consistent with previous studies, with the strongest energy observed in Rayleigh head waves. The application of the AK135 Earth model highlights the visibility of seismic phases through the Earth’s core. We validate our method by comparing synthetic records with actual data from the 1999 Chi-Chi earthquake. The synthetic records show good agreement with observed seismic signals and ionospheric perturbations in terms of arrival time and wave envelope. These results demonstrate the accuracy of our method in simulating acoustic-gravity waves at large epicentral distances.

SummaryThe widespread, multi-year crustal deformation induced by megathrust earthquakes (Mw8+) is primarily controlled by the combined effects of continuous aseismic slip on the fault plane (afterslip) and viscoelastic relaxation driven by coseismic stress perturbations in the upper mantle. However, till today it remains a considerable challenge to separate these two mechanisms in geodetic observations. We derived the first 3-year GNSS observations following the 2021 Chignik Mw8.2 earthquake to investigate the mechanisms of postseismic deformation. We established a model capable of simultaneously simulating afterslip and viscoelastic relaxation, and constrained the upper mantle rheology beneath the Alaska Peninsula. The best-fit model effectively reproduces the GNSS observations and reveals a notable viscosity difference between the mantle wedge and the oceanic asthenosphere, with steady-state viscosities of $3 \times {{10}^{18}}$ Pa s and $4 \times {{10}^{19}}$ Pa s, respectively. The inferred mantle wedge viscosity beneath the Alaska Peninsula is lower than the values reported for south-central and southeastern Alaska, suggesting an eastward increase in viscosity along the subduction zone. Two main patches of afterslip are identified during the first 3 years. The patch of up-dip afterslip overlaps with the 1938 Chignik Mw8.3 earthquake rupture zone, and demonstrates a close spatial correlation with the slow slip event in 2018. The above new results enhance our insights into the spatial variability of regional rheology and slip behavior along the Alaska-Aleutian subduction zone.

SummaryNew Q (1/attenuation) models of the Central Alpine Fault provide unprecedented resolution to 20 km depth by incorporating new t* measurements from dense temporary seismograph deployments in the area. The models reveal significant heterogeneity in the crust, with the main Q features broadly similar along-strike the Alpine Fault but varying at length scales of 10-30 km. Accounting for heterogeneity is an important step towards understanding the seismic cycle of M7+ Alpine Fault earthquakes. Our models show the Alpine Fault as a southeast-dipping zone of very (<300) to moderately (600-900) low Q, contrasting sharply with high Q values (Qp>600, Qs>1000) within the Western Province bedrock and high Q values (Qp∼900, Qs∼1200) associated with uplifted Alpine schists to the east. The wealth of previous geologic and geophysical studies along this section of the Alpine Fault support a detailed interpretation of the observed Q values. We interpret the low Q values along the Alpine Fault as resulting from enhanced fracturing within the brittle crust with a proportion of these fractureslikely filled with fluids, which further enhance seismic attenuation through viscous dissipation. In the ductile crust (below ∼8 km depth), low Q values (<400) are likely predominantly caused by grain-size reduction from very high total shear strain and by small amounts of metamorphic fluids. Low Q values of 200-400 at 20-40 km depth downdip of the Alpine Fault and the generally low Q (<600) within the crustal root farther from the Alpine Fault, suggest increasing role of metamorphic fluids relative to that of grain-size reduction with depth and distance from the fault. The updated model also reveals a newly identified zone of low Q east of the Main Divide, approximately 40 km southeast of the Alpine Fault trace. This zone of low Q indicates significant strain accumulation on faults striking along the eastern flank of the Southern Alps, some of which have produced M6+ earthquakes in recent history. These faults represent a considerable seismic hazard for the South Island. The improved dataset and recent velocity models from temporary deployments also allow us to investigate the influence of the initial velocity model on the resulting t* measurements and Q models.

SummaryThe bulk component of the electrical conductivity of a porous material is related to the (connected) porosity and saturation by power-law functions defining the first and second Archie's laws. Recently, it was shown that for porous materials with fractal characteristics, the power-law exponent of Archie's law could be related to the fractal dimension of such materials. Similarly, the real and imaginary parts of the complex-valued surface conductivity are not just proportional to the specific surface area and saturation of the material but to power law functions of these properties defining two additional “interfacial” Archie's laws, which are called the third (saturated case) and fourth (unsaturated case) Archie's laws in this paper. These new laws have been poorly recognized and studied so far. A number of porous materials and especially clay-rich media are multiscale materials characterized by broad distributions of particle and pore sizes. We extend Archie's laws concept to describe the complex conductivity of such materials. We use both numerical simulations in fractal porous materials as well as published experimental datasets to propose a unified physical interpretation of the exponents entering the four Archie's type power-law relationships, which offer an updated complex conductivity model for natural porous media.

SummeryThe seismic hazard due to higher magnitude Himalayan earthquakes largely depends on the geometry of the underthrusting Indian Plate beneath the Himalayas, i.e., the Main Himalayan Thrust (MHT). For an objective assessment of seismic hazard in the central Himalayan seismic gap, we determine the geometry of the Main Himalayan Thrust (MHT) along 4 ∼SW-NE oriented arc-normal seismic profiles covering the central Himalayan seismic gap. We use teleseismic P- wave coda autocorrelation on waveforms recorded at 117 broadband seismic stations spread along these profiles, with an interstation spacing of 3-5 km. The results show that along these seismic profiles, the MHT is mostly of flat-ramp-flat geometry. However, the mid-crustal ramp of the MHT shows variations in its location, dip angle, and width. We also observe variations in the MHT near the Main Frontal Thrust (MFT) and Main Boundary Thrust (MBT). The observed variations in the MHT geometry within the central Himalayan seismic gap thus suggest the possibility of along-strike segmentation of the Himalayan arc, and different seismic hazard scenarios may be present during any possible higher magnitude earthquake in the central Himalayan seismic gap.

SummaryRayleigh wave is widely used for characterizing shallow subsurface structures. The conventional Rayleigh-wave methods rely on the manual picking of dispersion curves, and the dispersion curves of multi-component data are usually merged manually. The manual processing of multi-component Rayleigh waves reduces the efficiency of the method, especially when the data size and the number of modes are large. To overcome these limitations, we develop an energy-based clustering method, namely the Energy-Density-Based Spatial Clustering of Applications with Noise (E-DBSCAN) algorithm. The E-DBSCAN algorithm extracts energy clusters and dispersion curves from a single dispersion image. It considers the dispersion-energy values of the surface wave and is able to pick the dispersion curve more reliably compared with the conventional DBSCAN algorithm. We propose a two-step clustering approach for the automatic picking of multi-mode dispersion curves from multi-component data: we first extract the energy clusters in the dispersion spectra of horizontal- and vertical-component data using E-DBSCAN, respectively, and combine them in the frequency-velocity domain. Then we extract multi-modal dispersion curves from the combined multi-component energy clusters with E-DBSCAN or DBSCAN. Numerical results show that our proposed method has fairly high accuracy and estimates more abundant multi-modal dispersion curves than the single-component method. Two field examples, including an active-source and an ambient-noise dataset, prove the validity of our method and the outperformance of multi-component results compared with the single-component results. Our proposed method has a relatively low dependence on parameter selection and is also applicable to multi-offset data, which is valuable for picking multi-modal dispersion curves.

SummaryRecently, there has been an increasing interest in employing rotational motion measurements for seismic source inversion, structural imaging and ambient noise analysis. We derive reciprocity and representation theorems for rotational motion. The representations express the rotational motion inside an inhomogeneous anisotropic earth in terms of translational and rotational motion at the surface. The theorems contribute to the theoretical basis for rotational seismology methodology, such as determining the moment tensor of earthquake sources.

SummaryDetermining earthquake focal mechanisms is essential for understanding fault geometry and the stress field in the Earth's crust. When focal mechanisms are estimated using P-wave first-motion polarities, accurate polarity determination is critical. In recent years, deep-learning-based polarity-determination models have been developed. However, the estimation of focal mechanisms using P-wave polarities is often not robust. When automating this process using deep learning models, it is crucial to identify and utilize only those polarity predictions that the model infers with high accuracy and low uncertainty. In this study, we propose a novel deep learning model, PoViT-UQ, that combines a Vision Transformer (ViT) with Monte Carlo Dropout (MCD) to estimate high-precision initial P-wave polarity classification and arrival time detection with uncertainty quantification. Using seismic waveform data sampled at 100 and 250 Hz, the model classifies polarities into three classes (Up, Down, and Noise) and simultaneously estimates P-wave arrival times. The results showed a classification accuracy exceeding 98% and a standard deviation of 0.027 s in the arrival time estimation using the 250 Hz model. By integrating MCD, we evaluate prediction uncertainty and apply an interquartile range threshold of ≤0.15 to improve the accuracy of focal mechanism estimates. Validation using aftershock data from the 2016 Central Tottori Earthquake confirmed that our approach contributes to efficient and high-precision focal mechanism estimates. Our model advances automated initial P-wave polarity determination and enables reliable data selection based on uncertainty quantification.

SummaryThe flows within Earth’s fluid outer core push and pull on the core-mantle boundary (CMB) through dynamic pressure variations, potentially leading to deformation of the CMB. It is therefore crucial to obtain a realistic estimate of the pressure associated with flows within the fluid core. In many studies, it is commonly assumed that the flow tangent to the CMB is in balance between Coriolis and pressure gradient forces, known as a tangentially geostrophic (TG) flow. A static pressure field is thereby associated kinematically to the flow field at the core’s surface. We run direct numeric simulations of the magnetohydrodynamic equations in the Boussinesq approximation that can solve for the pressure field and allow for a comparison between a fully dynamic solution and the TG pressure estimate. An excellent agreement between the two pressure fields is found for a steady image of the core surface dynamics. However, the performance of the TG pressure estimate is not without limitations. Although it effectively captures most of the temporal dynamics associated with the fluid flow, discrepancies arise, particularly near the equator and for rapid changes in flow dynamics.

AbstractSeismology has entered the petabyte era, driven by decades of continuous recordings of broadband networks, the increase in nodal seismic experiments, and the recent emergence of Distributed Acoustic Sensing (DAS). This review explains how cloud platforms, by providing object storage, elastic compute, and managed databases, enable researchers to “bring the code to the data,” thereby providing a scalable option to overcome traditional HPC solutions’ bandwidth and capacity limitations. After literature reviews of cloud concepts and their research applications in seismology, we illustrate the capacities of cloud-native workflows using two canonical end-to-end demonstrations: 1) ambient noise seismology that calculates cross-correlation functions at scale, and 2) earthquake detection and phase picking. Both workflows utilize Amazon Web Services, a commercial cloud platform for streaming I/O and provenance, demonstrating that cloud throughput can rival on-premises HPC at comparable costs, scanning 100 TBs to 1.3 PBs of seismic data in a few hours or days of processing. The review also discusses research and education initiatives, the reproducibility benefits of containers, and cost pitfalls (e.g., egress, I/O fees) of energy-intensive seismological research computing. While designing cloud pipelines remains non-trivial, partnerships with research software engineers enable converting domain code into scalable, automated, and environmentally conscious solutions for next-generation seismology. We also outline where cloud resources fall short of specialised HPC-most notably for tightly coupled petascale simulations and long-term, PB-scale archives-so that practitioners can make informed, cost-effective choices.

SummaryThe 1985 Mw 6.9 Wuqia earthquake, one of the strongest instrumentally recorded seismic events in the Pamir foreland thrust system, caused significant surface ruptures. The Pre-earthquake KH-9 and post-earthquake WorldView-3 and SPOT-6 satellite images are used to investigate the fault rupture and slip behavior of this earthquake. We revealed a more detailed ∼22 km long displacement belt beyond the previously documented ∼15 km rupture, using optical image correlation with sophisticated error post-processing. Several new fractures in western segment are identified which are confirmed in the displacement map. A comprehensive analysis of the strike change, near-surface dip and cross-fault offsets shows a ∼1.6 km dextral strike-slip tearing fault resulted from the heterogeneous strain release. Based on the empirical scaling relationship, a down-dip rupture width of 10.55 km is estimated using the observed rupture length and inverted slip. Combined with the previously published 3D fault geometry based on seismic imaging, we suggest that the 1985 Wuqia earthquake ruptured only the upper ramp. This study provides precise constraints on surface rupture characteristics, and new insights into the complex rupture pattern of a thrust-type earthquake within the tectonically active Pamir foreland region.

SummaryUnderstanding the geophysical drivers of seasonal geocenter motion (GCM) variations remains challenging due to the complexity of Earth system interactions, limited data on individual mass redistribution components, and model uncertainties. This study presents a comprehensive investigation of seasonal GCM signals from April 2002 to January 2024 using the Fingerprint Approach (FPA), which enables direct quantification of contributions from distinct Earth system components. Additionally, Multichannel Singular Spectrum Analysis (MSSA) is applied to quantify the influence of terrestrial water storage (TWS), atmosphere (ATM), and ocean (OCN) variability on seasonal GCM fluctuations. Correlation and lag analyses are employed to explore their temporal relationships and underlying geophysical linkages. The results reveal that TWS, ATM, and OCN jointly explain 97.9 per cent, 98.1 per cent, and 90.8 per cent of the seasonal variance in the X, Y, and Z components of GCM, respectively. TWS exerts as the dominant contributor in the Y (66.4 per cent) and Z (67.9 per cent) components, while ATM and OCN each contribute less than 49 per cent in all components. Further analysis indicates that ATM, OCN, and TWS exhibit varying lag relationships with GCM in the X and Z components, while TWS demonstrates a notably stronger correlation with GCM in the Y component. Importantly, an approximately 120-day periodic signal identified in GCM is, for the first time, linked to global precipitation variability, providing a novel geophysical interpretation. These findings enhance our understanding of climate-driven geophysical mass redistribution and offer new insights into the processes governing seasonal GCM variations.

SummaryCrashing ocean waves, or surf, have previously been identified as persistent generators of coherent infrasound signals from 0.5 to 20 Hz. Here, we demonstrate that infrasonic and seismic (seismo-acoustic) signals from surf are composed of repetitive transient events which can be detected and characterized using template matching. Using data collected from a series of field experiments designed to study seismo-acoustic surf signals in Santa Barbara, California, we show that source regions of these events can be constrained primarily to just offshore of a local coastal headland using a reverse-time-migration implementation on a small spatial scale (<5 km2). Our data include one continuously running infrasound sensor (September 2022–July 2023) to examine temporal signal evolution, complemented by several short-duration campaigns involving various infrasound arrays, co-located seismometers, and video recordings. Throughout varied oceanographic and atmospheric conditions, we detect up to tens of thousands of independent surf repeaters per day over the course of a year. The amplitudes of detected infrasound signals are correlated with offshore significant wave height and local wind speed. We identify coincident arrivals of seismic and infrasound signals with similar spectral characteristics, suggesting a linked source mechanism locally producing both the seismic and acoustic transient signals. Source regions estimated from array- and network-based methods correspond to the surf zone as seen in video footage, and the directions of selected transient signals align with the location of a rocky reef shelf nearshore. This work showcases the ability to extract near-real-time information about the coastal sea state from seismic and acoustic signal features.

SummaryCurrent geodetic velocities show that over half (up to 10 mm/yr) of Arabia-Eurasia shortening in the west is accommodated within a relatively narrow zone across the Kura basin of Azerbaijan, in which the most prominent active structure is the Kura fold-and-thrust belt, bordering the southern margin of the Greater Caucasus. The GNSS velocities furthermore suggest equivalent amounts of north-south right-lateral shear across the eastern Kura basin along the West Caspian fault zone that is accommodating relative motion between the Kura basin and the South Caspian basin. Although destructive historical earthquakes are known to have occurred, their spread is restricted geographically and their moment release accounts for only half of the accumulated deformation. These observations can be explained by incompleteness of the historical record, that the faults fail in rare larger earthquakes, or that they slip aseismically. To distinguish between these hypotheses we produce an InSAR velocity field using Sentinel-1 SAR data to image active tectonic deformation within the Kura basin of Azerbaijan. Tectonic signals are superimposed on those relating to non-tectonic processes, including widespread mud volcano inflation that highlights the important role of fluid flow within the basin sediments. We show aseismic creep occurs on two parallel faults of the West Caspian fault zone, and infer this also on the Kura Fold and Thrust Belt from sharp gradients in velocity indicating active fold growth. Recent paleoseismic studies of the faults imaged here indicate discrete slip events, and we speculate that the creep may be episodic, perhaps triggered by deeper earthquake events or by periods of enhanced fluid mobilisation. Together, the right-lateral and left-lateral faults appear to accommodate a large-scale expulsion of the Absheron region towards the South Caspian basin, perhaps driven by gravitational potential energy contrasts.

SummaryWhile global sea level rise is a major concern for ocean-connected coastal regions, inland seas such as the Black Sea exhibit water level changes primarily governed by regional hydrological and climatic factors. Understanding the drivers of water level variability in the Black Sea is essential due to its sensitivity to river inflow, evaporation, and limited connection to the global ocean. This study, for the first time, integrates satellite altimetry and satellite gravity data from 2003 to 2023 to analyze the long-term and seasonal variations in the Black Sea water levels, as well as local sea level variability and its driving factors. The results indicate that the sea level in the Black Sea experiences a positive trend of 1.04 ± 0.39 mm/yr. This comes, however, with a negative trend on a seasonal scale during autumn (–1.14 ± 0.27 mm/yr), which contrasts with the rising trends observed in other seasons. We found that the loading deformation induced by global mass redistribution contributes to 39 per cent of the Black Sea level trend, leading to an overestimation of actual climate-induced sea level change by 0.41 mm/yr. We further found that the reduction of precipitation and river inflow from surrounding basins leads to an increase in the salinity of the Black Sea, driving the decline in steric sea level. On the other hand, it has also increased the water mass of the Black Sea, compensating for the steric sea level drop.

SummaryDetecting marsquakes is crucial for understanding the interior structure of Mars. However, the detection and analysis of such seismic events is challenging due to the strong background noise and low signal-to-noise ratios (SNR) of marsquake seismograms. Moreover, the volume of the seismic data needed to create an efficient method to detect marsquakes, and the required amplitudes of such events is limited by the low energy of the seismic. Here, we develop a new deep learning-based detection algorithm, MarsConvNet (MANet), which can fully exploit the latent information of seismic sequences by upsampling and downsampling for detecting marsquakes in the continuous seismic records. Our approach is to construct a deep encoder structure and nonlinear mapping relationship between seismic signals and arrival times of P and S waves by training on the STanford EArthquake Dataset (STEAD) from Earth. Application to the continuous waveform data from the InSight seismic data from 2019 to 2022 shows that the algorithm can detect marsquakes with weak amplitude in addition to reproducing all events detected manually by the Marsquake Service (MQS). Using our method, we detected 67 previously undetected marsquakes.

SummaryA new approach for the parallel forward modelling of transient electromagnetic (TEM) fields is presented. It is based on a family of uniform-in-time rational approximants to the matrix exponential that share a common denominator independent of the evaluation time points. The partial fraction decomposition of this family is exploited to devise a fast solver with high parallel efficiency. The number of shifted linear systems that need to be solved in parallel does not depend on the number of required time channels nor the spatial discretization. We also argue that similar parallel efficiency gains can be expected when solving the inverse TEM problem.

AbstractDesert strata consist of dune and interdune deposits, and they have different magnetic properties. Contrasting magnetic properties of these two types of strata is a foundation for extracting paleoclimate information from dune-interdune strata in these eco-fragile regions, but few have done so. Here we compared magnetic properties between dune and interdune strata from the late Miocene desert deposits in the northwestern Tarim Basin. We found that interdune strata displayed enhanced magnetic properties mainly due to content increase of stable single-domain (SSD) and small pseudo single-domain magnetite (PSD) within chlorite layers. But magnetic enhancement in interdune strata does not correspond with content increase of superparamagnetic (SP) ferrimagnetic minerals. Our high-resolution transmission electron microscope results suggest pedogenic SSD and small PSD magnetite form within chlorite layers, and we propose weathering of chlorites resulted in high concentration of Fe (II) and Fe (III) ions within layers, facilitating the formation of SSD and small PSD magnetite. By contrast, dry climate in the Tarim Basin results in micropore spaces between mineral particles having limited concentration of Fe (II) and Fe (III) ions, inhibiting SP magnetite formation. This pattern is different from Chinese loess-paleosols sequence, where magnetic enhancement is accompanied with content increase of both SP and SSD or small PSD ferrimagnetic mineral content. Therefore, magnetic enhancement pattern differences (whether magnetic enhancement is accompanied with SP magnetite formation) may inform climate conditions. To infer environmental variations in arid regions using magnetic parameters, reliance on proxies that indicate the relative content of ferrimagnetic grains or content of magnetic grains with stable remanence is essential.

SUMMARYSeismology has advanced significantly with ambient noise interferometry, enhancing the extraction of surface waves for detailed Vs imaging. However, direct Vp measurements remain critical for geological and engineering applications. Guided P waves offer potential for Vp insights, yet their detection from ambient noise is challenging, especially in urban environments. This study explores the feasibility of using underground urban seismic noise, specifically tunnel traffic, for guided P wave extraction and subsurface imaging. We introduce an adaptive workflow for preprocessing and enhancing guided P wave signals, applied to data from the Zijingang Tunnel beneath the campus of the Zhejiang University. Our refined data processing workflow significantly improves the quality of retrieved empirical Green’s functions of guided P waves. The joint inversion of surface waves and guided P waves provides high-resolution Vs and Vp profiles, revealing detailed subsurface structures. We also examine the impacts of source depth and type on guided P wave retrieval. We observed distinct differences in urban ambient noise characteristics between daytime and nighttime, suggesting that daytime traffic noise, with higher frequency content, is more suitable for high-frequency guided P wave retrieval. Our results demonstrate the effectiveness of using underground urban seismic noise for guided P wave extraction, highlighting the potential applications of urban ambient noise imaging in geotechnical engineering, urban planning, and seismic hazard assessment.

SUMMARYAirborne natural source electromagnetics (EM) provide spatially densely gridded, frequency-dependent inter-site transfer function estimates referring airborne vertical and horizontal magnetic field data to the horizontal components of a ground site. Common processing schemes for airborne natural source EM data are conceptually based on approaches developed for ground-based magnetotelluric sounding curve estimation and yield statistically independent transfer function estimates for selected times windows and separately for each component of the airborne recording. However, the magnetic field is spatially smooth and, moreover, its components are mutually dependent in space. Here, we propose to exploit spatial potential field relations directly in the processing workflow to refine transfer functions. The self-consistent response function estimator is applied after a first statistical estimation of transfer functions and supports the determination of a spatially smooth and self-consistent set of vertical and horizontal magnetic transfer functions that reflects physical constraints. Self-consistency is enforced by fitting sources distributed within an equivalent sheet to a preliminary set of transfer function estimates determined on a statistical basis. Smoothness is controlled by the depth to the equivalent sources. Once estimated, the distribution of vertical and horizontal magnetic transfer functions can be predicted from the equivalent sheet at any point in the air half-space. Additionally, the scheme allows for the identification of anomalous parts of the horizontal fields at the ground site. The procedure is demonstrated on a large-scale field data set from Gobabis (Namibia). We observe that transfer functions derived by using the equivalent source layer interpolator improve processing results when compared to purely statistically determined transfer function averages.