Geophysical Journal International

SummaryApplications of spectral induced polarization (SIP) require electrodes that maintain hydrologic contact with surrounding soils to capture small electrical responses, often observed as phase shifts in milliradians. For unsaturated soils, electrodes must overcome the increased electrical contact impedance due to reduced pore fluid. Traditional designs use a ceramic membrane electrode (CME) with a water reservoir and metal conductor, requiring periodic maintenance to retain electrolytic solution. For field applications where maintenance is impractical, alternative designs are needed. This study evaluated a new electrode design (silica flour electrode, SFE) alongside a CME design. SFEs use packed silica flour to store water via capillary forces against a metal conductor. The study examined both designs in three variably saturated soils at soil suctions up to 700 mbar and soil water contents below 1 per cent, with SIP measurements across 0.01 to 10,000 Hz frequencies. SFEs match CMEs at high frequencies and perform better at lower frequencies, without requiring ongoing maintenance, making them ideal for field use. In water-only experiments, CMEs produced errors and high noise below 1.5 Hz, whereas SFEs were more accurate. However, CMEs performed better above 300 Hz. In fine sand, SFEs performed better due to the relatively lower contact impedance as compared to CMEs. Both electrode types performed comparably in silty sand and silt loam soils, although CMEs required ongoing maintenance, suggesting potential for long-term reliability issues.

SummaryThe finite-difference method (FDM), limited by uniform grids, often encounters severe oversampling in high-velocity regions when applied to multi-scale subsurface structures, leading to reduced computational efficiency. A feasible solution to this issue is the use of non-uniform grids. However, previous discontinuous grid approaches required careful consideration of interpolation operations in transition regions, while single-block continuous grids lacked flexibility. This paper proposes a novel approach using multi-block stretched grids with positive and negative singularities to achieve non-uniform grids, the numerical simulation of seismic waves is realized by combining it with the curvilinear grid finite-difference method (CGFDM). Our method facilitates seamless information exchange between coarse and fine grids without additional interpolation or data processing and allows for flexible grid configurations by adjusting singularity pairs.The effectiveness of our approach is verified through comparisons with the generalized reflection/transmission method (GRTM) and the finite-element method (FEM). Numerical experiments demonstrate the method's reliable accuracy and significant reduction in grid points compared to uniform grids. Although the stability of our method has not been rigorously mathematically proven, we demonstrate that the algorithm remains applicable for sufficiently long simulations to address realistic scenarios.

SummaryFor a weakly anisotropic medium, Rayleigh and Love wave phase speeds at angular frequency ω and propagation azimuth ψ are given approximately by V(ω, ψ) = A0 + A2ccos 2ψ + A2ssin 2ψ + A4ccos 4ψ + A4ssin 4ψ. Earlier theories of the propagation of surface waves in anisotropic media based on non-degenerate perturbation theory predict that the dominant components are expected to be 2ψ for Rayleigh waves and 4ψ for Love waves. This paper is motivated by recent observations of the the 2ψ component for Love waves and 4ψ for Rayleigh waves, referred to here as “unexpected anisotropy”. To explain these observations, we present a quasi-degenerate theory of Rayleigh-Love coupling in a weakly anisotropic medium based on Hamilton’s Principle in Cartesian coordinates, benchmarking this theory with numerical results based on SPECFEM3D. We show that unexpected anisotropy is expected to be present when Rayleigh-Love coupling is strong and recent observations of Rayleigh and Love wave 2ψ and 4ψ anisotropy can be fit successfully with physically plausible models of a depth-dependent tilted transversely isotropic (TTI) medium. In addition, when observations of the 2ψ and 4ψ components of Rayleigh and Love anisotropy are used in the inversion, the ellipticity parameter ηX, introduced here, is better constrained, we can constrain the absolute dip direction based on polarization measurements, and we provide evidence that the mantle should be modeled as a tilted orthorhombic medium rather than a TTI medium. Ignoring observations of unexpected anisotropy may bias the estimated seismic model significantly. We also provide information about the polarization of the quasi-Love waves and coupling between fundamental mode Love and overtone Rayleigh waves in both continental and oceanic settings. The theory of SV-SH coupling for horizontally propagating body waves is presented for comparison with the surface wave theory, with emphasis on results for a TTI medium.

SummaryExtensional faults in Southern Calabria (Italy) have been widely studied for their capability of generating high magnitude earthquakes (Mw 7-7.2). An example is the historical seismic sequence occurred in 1783, which caused numerous fatalities near the villages located along the longest faults of this region: the Cittanova and the Serre faults. In this work, we estimated the seismic potential of these two faults by a kinematic block modeling approach using GNSS data of both campaign points and permanent stations. Our results indicate that both faults are accommodating the recognized extensional velocity gradient (∼ 1 mm/yr) by long-term slip rates (∼ 2 mm/yr). To estimate the back slip distribution and the interseismic coupling degree of the Cittanova and Serre faults, we discretised these by a triangular dislocation elements (TDEs) mesh. This approach has allowed us to distinguish the fault areas where elastic seismic rupture is more likely to happen from those affected by aseismic creeping behaviour. The obtained results show that the highest values of coupling are located near the shallow portion of the fault planes and near the southern tip of the Cittanova fault. We therefore estimated a set of possible rupture scenarios finding that the Southern Calabria domain is accumulating an interseimic moment rate at most equal to 2.16 ×1016 Nm/yr, the equivalent of an earthquake of Mw 4.86 for each year.

SummaryChemical potentials are defined as the partial derivatives of the Helmholtz energy with respect to moles of chemical components under conditions of zero domain strain and fixed temperature. Under hydrostatic conditions, chemical potentials are dependent only on state properties. Under nonhydrostatic conditions, they also depend on a “chemical expansivity tensor” - a second-order tensor with unit trace that characterises how the elastic network is compressed to accommodate new material within the local domain element. The five degrees of freedom of this tensor generate a class of chemical potentials. An important group within this class are the “uniaxial chemical potentials”, which quantify the Helmholtz energy change when new material is incorporated via compression along a single axis. Chemical and mechanical equilibrium is achieved when all uniaxial chemical potentials remain constant along their respective axes. The derived expressions apply to both crystalline and amorphous materials. Their utility is demonstrated through solutions to classic phase-equilibrium problems.