Sharp Velocity Research Articles

Double-Pair Double-Difference Relocation Improves Depth Precision and Highlights Detailed 3D Fault Geometry for Induced Seismicity in Alberta, Canada

Abstract Precise earthquake locations with well-constrained uncertainties can improve our understanding of faulting. Double-difference relocation methods, particularly event-pair double-difference relocations, are well established and have been applied to large earthquake catalogs to provide fault geometries. Previous adaptations of the event-pair double-difference method include data space extensions to use additional information from station pairs, referred to as double-pair double-difference relocation. We apply double-pair double-difference relocation to data from a dense network of borehole geophones for induced seismicity monitoring. This experiment was acquired in an area with strong lithological variation and sharp velocity contrasts, and most previous studies using this dataset are subject to poorly constrained focal depths. We compare the double-pair double-difference to event-pair double-difference relocations and study the effectiveness and uncertainties of both methods. Although double-pair double-difference relocation does not improve absolute locations, substantially improved relative locations and reduced uncertainties are obtained. The method reduces the impact of path effects in the source region, which is essential for applications where reservoir units in the source region can exhibit strong velocity contrasts, anisotropy, and fractures. From the improved relocation, we produce a detailed 3D fault interpretation of the dataset that constrains the geological interpretation. The improved catalog shows excellent depth constraints with seismicity that is restricted to specific geological units. We interpret that seismicity activated pre-existing faults in the reservoir layer and adjacent units. Notably, the results show no evidence of induced seismicity activating basement structures.

Deep plutonic bodies over low-frequency earthquakes revealed from receiver-side Green's functions

Seismological heterogeneity in subduction zones provides insights into slow earthquakes and potential megathrust earthquakes. Studies at the Kii Peninsula in the Nankai subduction zone suggest that there are high-density and high-velocity plutonic bodies in the accretionary prism over the subducting slab, potentially influencing megathrust earthquakes. The lateral variation of heterogeneity and the spatial extent of plutonic bodies remain to be investigated well. Our passive-source imaging of receiver-side Green's functions, from widely distributed campaign seismic observations, reveals a sharp negative S-wave velocity contrast on the top surface of the subducting Philippine Sea plate common to all along-dip profiles and a positive phase tilted upward in the forearc crust. The low permeability of the forearc crust prevents the infiltration of slab-dehydrated fluid further into the upper crust. In the western area, we also found positive phases tilted upward in the forearc crust. The negative phase extends towards the deeper extent of slow-earthquake sources. Meanwhile, the positive phase likely represents the top surface of plutonic rocks of the Kumano and Ohmine plutons that span all the way down to the plate interface. Together with observations of gravity anomaly, intraslab seismicity, and seismic tomography, our interpretation supports the presence of plutonic bodies which extend deep beneath the forearc crust as well as laterally over the subducting PHS slab, rather than a serpentinized mantle wedge. The upper plate is generally low in permeability, but areas with localized high permeability may exist on the updip side of tremor sources. This condition, wherein fluid can infiltrate upwards locally, may maintain the relatively less active slow earthquakes in the western area. The lateral variation of the upper-plate lithology likely influences fluid processes and slow earthquake activities.

ATOMS: ALMA three-millimetre observations of massive star-forming regions–XVII. High-mass star-formation through a large-scale collapse in IRAS 15394-5358

ABSTRACT Hub-filament systems are considered as natural sites for high-mass star formation. Kinematic analysis of the surroundings of hub-filaments is essential to better understand high-mass star formation within such systems. In this work, we present a detailed study of the massive Galactic protocluster IRAS 15394$-$5358, using continuum and molecular line data from the ALMA three-millimetre observations of massive star-forming regions (ATOMS) survey. The 3 mm dust continuum map reveals the fragmentation of the massive ($\rm M=843~{\rm M}_{\odot }$) clump into six cores. The core C-1A is the largest (radius = 0.04 pc), the most massive ($\rm M=157~{\rm M}_{\odot }$), and lies within the dense central region, along with two smaller cores ($\rm M=7~and~3~{\rm M}_{\odot }$). The fragmentation process is consistent with the thermal Jeans fragmentation mechanism and virial analysis shows that all the cores have small virial parameter values ($\rm \alpha _{vir}\lt \lt 2$), suggesting that the cores are gravitationally bound. The mass versus radius relation indicates that three cores can potentially form at least a single massive star. The integrated intensity map of $\rm H^{13}CO^{+}$ shows that the massive clump is associated with a hub-filament system, where the central hub is linked with four filaments. A sharp velocity gradient is observed towards the hub, suggesting a global collapse where the filaments are actively feeding the hub. We discuss the role of global collapse and the possible driving mechanisms for the massive star formation activity in the protocluster.

Segregation‐Induced Flow Transitions in Rock‐Ice Mixtures: Implications for Rock‐Ice Avalanche Dynamics

AbstractGlobal climate change has been intensifying the scale and frequency of rock‐ice avalanches and similar catastrophic mass movements in high‐mountain regions. The difference in the physical characteristics of rock and ice particles leads to mixing and segregation during flow. Although, both particle segregation and the presence of ice fundamentally alter flow behavior, the joint influence and feedback of these two aspects are overlooked in state‐of‐the‐art rock‐ice avalanche models. Using discrete element simulations, we show that by controlling the distribution of inter‐particle frictional interactions within the mixture, segregation patterns resulting from the size, density, concentration, and surface friction differences of rock and ice phases can induce sharp velocity gradients along the flowing thickness. Flowing layers where low friction contacts with ice are abundant tend to flow faster and can induce slow creeping motion in an otherwise static basal layer dominated by more frictional rocks. Based on these observations, we find that the effective friction of rock‐ice flows for various mixture concentrations and size ratios can be obtained as a sum of the single‐phase rheologies of rocks and ice weighted according to their microscopic contact probabilities. This effective friction for rock‐ice mixtures allows us to extend a recent non‐local granular fluidity framework that captures the complex segregation‐flow feedback mechanism in rock‐ice flows. The findings provide a deeper micromechanical understanding of how particle interactions influence rock‐ice avalanche mobility, which ultimately improves flow models needed for hazard assessment and mitigation.

Is There a Carbonated Mid‐Lithosphere Discontinuity in Cratons?

AbstractThe mid‐lithosphere discontinuity (MLD), identified by a sharp velocity drop at ∼70–100 km depths within the cratonic lithosphere is key to comprehending the chemical composition and thermal structure of the cratonic lithosphere. The MLD is widely accepted to be caused by composition anomalies, such as hydrous minerals, which show low velocities and high electrical conductivities. However, noticeable high‐electrical conductivity anomalies have not been detected in the most cratonic lithosphere. Dolomite has an electrical conductivity similar to olivine and can be originated by carbonatitic melts trapped at ∼80–140 km depths. Here we investigated the elasticity of dolomite under mantle conditions using ab initio calculations and found dolomite exhibits significantly lower velocities than the primary minerals in the lithospheric mantle. Therefore, the dolomite enrichment might provide a good explanation for the observed velocity drop of the MLD in cratonic regions where no high‐conductivity anomaly has been detected, such as the northern Slave craton.

Seismic evidence for melt-rich lithosphere-asthenosphere boundary beneath young slab at Cascadia

The Lithosphere-Asthenosphere Boundary (LAB) beneath oceanic plates is generally imaged as a sharp seismic velocity reduction, suggesting the presence of partial melts. However, the fate of a melt-rich LAB is unclear after these plates descend into the mantle at subduction zones. Recent geophysical studies suggest its persistence with down-going old and cold slabs, but whether or not it is commonly present remains unclear, especially for young and warm slabs such as in the Cascadia subduction zone. Here we provide evidence for its presence at Cascadia in the form of a large (9.8±\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\pm$$\\end{document}1.5%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\%$$\\end{document}) decrease in shear-wave velocity over a very small (<3 km) depth interval. Similarly large and sharp seismic velocity reduction at the bottom of both old and young slabs, as well as along the base of oceanic plates before subduction, possibly represents widespread presence of melts. The melt-rich sub-slab LAB may strongly influence subduction dynamics and viscoelastic earthquake cycles.

Hamiltonian Monte Carlo based elastic full-waveform inversion of wide-angle seismic data

SUMMARY Full-waveform inversion (FWI) of seismic data provides quantitative constraints on subsurface structures. Despite its widespread success, FWI of data around the critical angle is challenging because of the abrupt change in amplitude and phase at the critical angle and the complex waveforms, especially in the presence of a sharp velocity contrast, such as at the Moho transition zone (MTZ). Furthermore, the interference of refracted lower crustal (Pg) and upper mantle (Pn) arrivals with the critically reflected Moho (PmP) arrivals in crustal and mantle studies makes the application of conventional FWI based on linearized model updates difficult. To address such a complex relationship between the model and data, one should use an inversion method based on a Bayesian formulation. Here, we propose to use a Hamiltonian Monte Carlo (HMC) method for FWI of wide-angle seismic data. HMC is a non-linear inversion technique where model updates follow the Hamiltonian mechanics while using the gradient information present in the probability distribution, making it similar to iterative gradient techniques like FWI. It also involves procedures for generating distant models for sampling the posterior distribution, making it a Bayesian method. We test the performance and applicability of HMC based elastic FWI by inverting the non-linear part of the synthetic seismic data from a three-layer and a complex velocity model, followed by the inversion of wide-angle seismic data recorded by two ocean bottom seismometers over a 70 Ma old oceanic crustal segment in the equatorial Atlantic Ocean. The inversion results from both synthetic and real data suggest that HMC based FWI is an appropriate method for inverting the non-linear part of seismic data for crustal studies.

The Marmara Sea basin as a regional depression constrained from ambient noise correlation tomography

SUMMARY We computed a 3-D shear wave velocity model of the Marmara Sea region from ambient noise tomography. The correlations of up to 8 yr of vertical-component seismic recordings from 80 broad-band stations provided Rayleigh wave group velocity measurements in the period band 6–21 s at more than 1400 selected virtual source–receiver pairs. Rayleigh wave group velocity maps were used to derive a shear wave velocity model through simulated annealing inversion. The resulting crustal model provides coverage of the Marmara Sea along with its surrounding regional tectonic features. This allows for an investigation of the spatial extents of the Marmara Sea on a scale larger than that of basins. The low-velocity structures of the Marmara Sea and the Thrace Basins are coeval to a depth of approximately 9 km. The crustal velocities beneath the Marmara Sea basins exhibit a low vertical gradient and smooth horizontal variations. The regional tectonic structures, such as Istranca Massif, Istanbul and Sakarya Zones, display sharp velocity contrasts with the lower velocity crust beneath the Marmara Sea. The observed low crustal velocities, along with depth variations of the velocity isosurfaces (i.e. 3.4 km s−1) indicate that the Marmara region is a structural depression much deeper and larger than the three basins of the North Marmara Trough. The North Anatolian Fault Zone is unlikely to be the primary factor contributing to the origin of this significant depression, as the basin's development appears to have occurred before the fault propagated into the region.

Effect of surface roughness on large-scale downburst-like impinging jets

Downbursts are cold descending winds that develop from thunderstorm clouds and, after impingement on the ground, produce an intense low-level horizontal front characterized by an axisymmetric toroidal vortex structure. Surface roughness is a key factor in the characterization of mean and turbulent wind speed features of synoptic-scale stationary atmospheric boundary layer winds. The goal of the present research is to physically assess whether the same can apply to the surface layer produced during thunderstorms, which are non-stationary, highly time-transient, and spatially limited phenomena. Downburst-like flows were produced through the impinging jet technique at the WindEEE Dome, at Western University in Canada. Three different surfaces were tested, and an equivalent full-scale roughness length (z0,eq) was determined. Experimental records are made publicly available. The large geometric and kinematic scales produced high Reynolds numbers, which enabled us to classify the flow as “fully turbulent” and therefore representative of full-scale downbursts. Results indicate a weak dependency on the Reynolds number, which suggests no relevant flaws in extending the results to the natural environment. The overall wind speed maxima weakly depend on z0, whereas a sharp velocity decrease is observed beyond the radial position of the maxima with increasing z0. Surface roughness enhances the boundary layer separation and consequently elevates the height of maximum wind speed above the surface. Vertical profiles of the horizontal velocity return a quite clear nose shape. Turbulence intensity shows a C-like shape with maxima in the near proximity of the ground that increase with z0.

Martian Crustal Model From a Joint Inversion of Receiver Functions and Apparent Shear Wave Velocity

AbstractReceiver function (RF) based studies have revealed the layered Martian crustal structure beneath the InSight, which provides key information about the composition, rock alteration, deep magma process, and thermal history of Mars. However, limited by the frequency and trade‐offs between the thickness and velocity, inverting RFs alone is difficult to provide an accurate image of the upper crust. Here, by achieving stable measurements of the apparent shear wave velocities (Vs,app) at different frequencies for high‐quality marsquakes, we perform a joint inversion with the RFs. Our crustal models differ from previous models by having an extra shallow structure in the upper crust, either a layer with a sharp velocity jump interface at 4.0 ± 0.9 km depth or a gradient structure with shear wave velocity increasing gradually from 1.1 ± 0.1 km/s near the surface to 1.7 ± 0.3 km/s at a depth of 6.0 ± 2.4 km, which well fits both RFs and Vs,app's. Such an upper crustal velocity model is consistent with the depth‐decaying of porosity due to self‐compaction and indicates a highly porous upper crust with a surface porosity larger than 30%. Besides, a much higher porosity decay constant than that scaled from the Moon is derived, suggesting different mechanical properties of the lunar and Martian crust.

Analysis of Mass Wasting Processes in the Slumgullion Landslide Using Multi-Track Time-Series UAVSAR Images

The Slumgullion landslide is a large translational debris slide whose currently active part has likely been sliding for approximately 300 years. Its permanent motion and evolutionary processes have attracted the attention of many researchers. In order to study its mass wasting processes and evolution trend, the spatial–temporal displacement of the Slumgullion landslide was retrieved using an adaptive pixel offset tracking (POT) method with multi-track Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) images. Based on three-dimensional displacement and slope information, we then revealed the spatial–temporal distribution of surface mass depletion or accumulation in the landslide, which provides a new perspective to analyze the evolutionary processes of landslides. The results indicate that the Slumgullion landslide had a spatially variable displacement, with a maximum displacement of 35 m. The novel findings of this study mainly include two parts. First, we found that the surface mass accumulated in the toe of the landslide and depleted in the top and middle area during the interval, which could increase the resisting force and decrease the driving force of the Slumgullion landslide. This result is compelling evidence which indicates the Slumgullion landslide should eventually tend to be stable. Second, we found that the distribution of geological structures can well explain some of the unique mass wasting in the Slumgullion landslide. The larger local mass depletion in the landslide neck area verifies that the sharp velocity increase in this region is not only caused by the reduction in width but is also significantly affected by the local normal faults. In summary, this study provides an insight into the relation between the landslide motion, mass volume change, and geological structure. The results demonstrate the great potential of multi-track airborne SAR for displacement monitoring and evolutionary analysis of landslides.

PREM-like velocity structure in the outermost core from global SKS and ScS waveform modeling

Seismic structure in the Earth's outermost core is important for our understanding of thermochemical stratification in the outer core, and is yet heavily debated. Here we study the compressional velocity structure in the Earth's outermost core based on waveform modeling of a unique untapped SKS and ScS dataset near bifurcation distances, collected from global seismic arrays for earthquakes occurring from 2000 to 2020. Using the SKS-ScS array dataset minimizes the effects of many uncertainties associated with earthquake source parameters and seismic heterogeneities in the mantle, and affords an opportunity to study and assess the seismic structure in the outermost core. We study outer core structure by testing two end-member models: 1) the D″ model that attributes any anomalous seismic observations to the effects of the lowermost mantle structure and is paired with a PREM (the Preliminary Reference Earth's Model) structure in the outermost core and 2) the outer core model that is paired with either a PREM or a tomographic structure in the lowermost mantle and attributes other unexplained seismic signals to an outer core structure. The results of the outer core models present unreasonable large lateral variations of >3.1% in the outermost core, while the inferred D″ models exhibit large-scale seismic anomalies that are consistent with the tomographic models and small-scale anomalies that are confirmed by further analysis of the seismic array data. Our analyses suggest a PREM-like seismic velocity structure and a lack of strong thermochemical anomalies in the topmost ∼200 km of the outer core, placing bounds on possible thermal and compositional conditions in the region of the Earth's outermost core. Our study also identifies the existence of small-scale seismic anomalies and sharp velocity variations in the lowermost mantle beneath the south coast of Alaska, northwestern Atlantic and the middle of Central America.

Migration of Mechanical Perturbations Estimated by Seismic Coda Wave Interferometry During the 2018 Pre‐Eruptive Period at Kīlauea Volcano, Hawaii

AbstractWe use seismic ambient noise correlation and coda wave interferometry to estimate velocity variations at high temporal resolution, during the pre‐eruptive period and the onset of the 2018 eruption of Kīlauea volcano. A progressive velocity increase is observed from March to the end of April. It is followed by rapid decrease starting a few days before the onset of the East Rift Zone (ERZ) eruption and then by sharp velocity drop when the eruption started. The change of trend from velocity increase to decrease is progressively delayed by a few days from the summit caldera toward the ERZ. The location of the velocity perturbations shows a migration of the sources of velocity changes from the summit caldera toward the ERZ before the eruption. Using a model of pressure source, we show that the simultaneous caldera inflation and velocity increase probably result from an anisotropic distribution of fault and crack orientations. The velocity decrease could be due to damaging processes above the shallow reservoir and to plastic deformations around the caldera. We introduce a forward model of rock damage associated with the volcano‐tectonic seismicity to calculate the velocity decrease. The good agreement between the calculated and the observed velocity variations shows that a large part of the velocity decrease results from damage of the medium. The delayed onsets of velocity decrease and the source migration of velocity perturbations are probably related to progressive fault openings in the Southern and Eastern parts of the caldera and to magma transfer toward the ERZ.

Effect of water on seismic attenuation of the upper mantle: The origin of the sharp lithosphere–asthenosphere boundary

Oceanic lithosphere moves over a mechanically weak layer (asthenosphere) characterized by low seismic velocity and high attenuation. Near mid-ocean ridges, partial melting can produce such conditions because of the high-temperature geotherm. However, seismic observations have also shown a large and sharp velocity reduction under oceanic plates at the lithosphere-asthenosphere boundary (LAB) far from mid-ocean ridges. Here, we report the effect of water on the seismic properties of olivine aggregates in water-undersaturated conditions at 3 GPa and 1,223 to 1,373 K via in-situ X-ray observation using cyclic loading. Our results show that water substantially enhances the energy dispersion and reduces the elastic moduli over a wide range of seismic frequencies (0.5 to 1,000 s). An attenuation peak that appears at higher frequencies (1 to 5 s) becomes more pronounced as the water content increases. If water exists only in the asthenosphere, this is consistent with the observation that the attenuation in the asthenosphere is almost constant over a wide frequency range. These sharp seismic changes at the oceanic LAB far from mid-ocean ridges could be explained by the difference in water content between the lithosphere and asthenosphere.

Enhanced heat transfer in a microchannel with pseudo-roughness induced by Onsager-Wien effect

A three-dimensional numerical analysis of flow and conjugate heat transfer in a microchannel in the presence of the electric-field-induced Onsager–Wien effect is performed. A novel design is proposed to induce a pseudo-roughness effect in the microchannel and thereby increasing the heat transfer. A series of thin plate electrode pairs are flushed along the bottom wall of the microchannel. The electric-field-enhanced dissociation of ions induces the Onsager–Wien effect, and generates small flow vortices near the bottom wall of the channel. These flow vortices with sharp local velocity gradients effectively disrupt the viscous and thermal boundary layers, and thus, introduce a pseudo-roughness effect. This disruption of the boundary layers improves the heat transfer between the channel wall and the working fluid. The thermal and hydraulic performances in the microchannel are quantified as a function of the flow Reynolds number Re and electric Reynolds number ReEL. In general, the performance factor PF is higher when the flow and electric Reynolds numbers are higher. However, the associated pressure drop penalty reduces the PFPF<1 in viscous-effect-dominated flows at Re=250. By contrast, the Onsager–Wien effect increases the PF to a maximum value of 1.26 in inertia-dominated flows at Re=1000 and ReEL=15. A significant heat transfer enhancement is realized, even at a low applied voltage of ∼1 kV. The results of this study indicate that the pseudo-roughness produced by the electric-field-induced Onsager–Wien effect can substantially enhance the heat transfer in a microchannel with a trivial amount of additional power consumption. Results presented in this study will serve as a benchmark to design microchannels with enhanced heat transfer by using active vortex generation based on Onsager-Wien effect.

Interferometry Observations of the Gravity Wave Effect on the Sporadic E Layer

On the basis of interferometry measurement made with the Chung-Li VHF radar, we investigated the effects of upward propagating gravity waves on the spatial structures and dynamic behavior of the 3 m field-aligned irregularities (FAIs) of the sporadic E (Es) layer. The results demonstrate that the quasi-periodic gravity waves oscillating at a dominant wave period of about 46.3 min propagating from east-southeast to west-northwest not only modulated the Es layer but also significantly disturbed the Es layer. Interferometry analysis indicates that the plasma structures associated with gravity wave propagation were in clumpy or plume-like structures, while those not disturbed by the gravity waves were in a thin layer structure that descended over time at a rate of about 2.17 km/h. Observation reveals that the height of a thin Es layer with a thickness of about 2–4 km can be severely modulated by the gravity wave with a height as large as 10 km or more. Moreover, sharply inclined plume-like plasma irregularities with a tilted angle of about 55° or more with respect to the zonal direction were observed. In addition, concave and convex shapes of the Es layer caused by the gravity wave modulations were also found. Some of the wave-generated electric fields were so intense that the corresponding E × B drift velocities of the 3 m Es FAIs approximated 90 m s−1. Most interestingly, sharp Doppler velocity shear as large as 68 m/s/km of the Es FAIs at a height of around 108 km, which bore a strong association with the result of the gravity wave propagation, was provided. The plausible mechanisms responsible for this tremendously large Doppler velocity shear are discussed.

3-D imaging of the Balmuccia peridotite body (Ivrea–Verbano zone, NW-Italy) using controlled source seismic data

SUMMARY We provide new results from a controlled-source seismic experiment on the deepest part of the Val Sesia crust–mantle section of the Ivrea–Verbano zone (IVZ) in the Italian Alps. The IVZ is a tilted, almost complete section through the continental crust and exposes gabbros and peridotites in the structurally deepest level, coinciding with high-resolution gravity anomalies imaging the Ivrea geophysical body. The seismic experiment SEIZE (SEismic imaging of the Ivrea ZonE) was conducted along two crossing profiles: an NNE-SSW profile of ∼11 km length and an E-W profile of ∼16 km length. 432 vibration points were recorded with 110 receivers resulting in 24 392 traveltime picks. Inversion methods using Markov chain Monte Carlo techniques have been used to derive an isotropic 3-D P-wave velocity model based on first break traveltimes (refracted phases) from controlled source seismic data. Resulting seismic P-wave velocities (Vp ) range from 4.5 to 7.5 km s−1, with an expected general trend of increasing velocities with depth. A sharp velocity change from low Vp in the West to high Vp in the East marks the Insubric Zone (ISZ), the Europe–Adria plate boundary. The most prominent feature of the 3-D tomography model is a high-velocity body (Vp increases from 6 to 7.5 km s−1) that broadens downwards. Its pointy shape peaks the surface East of Balmuccia at a location coincident with the exposed Balmuccia peridotite. Considering rock physics, high-resolution gravity and other geophysical data, we interpret this high-velocity body as dominantly composed of peridotite. The dimension of this seismically imaged peridotite material is far bigger than interpreted from geological cross-sections and requires a revision of previous models. The interpretation of ultramafic bodies in the IVZ as fragments of mantle peridotites interfingered in the crust during pre-Permian accretion is not supported by the new data. Instead, we revive a model that the contact between the Balmuccia peridotite and the Permian mafic magmas might represent a fossil continental crust–mantle transition zone.

Extracting high‐resolution P‐wave reflectivity of the shallow subsurface by seismic interferometry based on autocorrelation of blast mining signals

AbstractBody‐wave reflections are sensitive to sharp velocity contrasts, making them useful for lithological imaging. We analysed seismic data from natural earthquake, ambient noise and mine blasts to map P‐wave reflection profiles at the Hishikari mine area by autocorrelation analysis. Because fissure‐filled gold veins are dominant in this area, we evaluated the potential of autocorrelation analysis for investigating the shallow subsurface, including the ore deposits. Seismic interferometry is commonly performed based on the autocorrelation of ambient noise or natural earthquake signals; here, we instead used blasting in the mine because blast events include high‐frequency signals that boost the spatial resolution of the imaging. To effectively extract P‐wave reflections from seismic signals including blast events, we applied Gaussian smoothing and spectral whitening to remove source effects and then investigated the optimum frequency band. We successfully obtained auto‐correlograms showing high‐resolution seismic reflectors at shallow formation depths. These reflections are interpreted to be lithological boundaries shallower than 500 m. A comparison with profiles obtained from ambient noise and earthquake data showed that blasting signals yielded highly spatially consistent reflections that would not be achievable with natural or ambient seismic sources. This study highlights the potential of using blast autocorrelation seismic analysis during short survey periods. By using single‐blast shots and dense seismic station spacings, we successfully achieved higher resolution 3D reflection images of lithological interfaces, possibly including ore veins.

Seasonal variability of eddy kinetic energy in the north Indian Ocean

The seasonality of eddy kinetic energy (EKE) is analyzed in the north Indian Ocean by adopting high-resolution ocean reanalysis data. Significant eddy energy can be mainly spotted in six regions, including the Somali Current (SC) region, the Gulf of Aden, the Laccadive Sea, the east of Sri Lanka, the East Indian Coastal Current (EICC) region, and the northwest of Sumatra. As the most energetic region, the EKE averaged above 200 m could exceed 0.15 m2·s-2 in the SC region, whereas the mean EKE above 200 m is less than 0.04 m2·s-2 in the other regions. The barotropic and baroclinic instabilities are vital to eddy energy, and the contribution of each term in the barotropic/baroclinic equations varies with season and region. In the SC region and EICC region, EKE is primarily generated by barotropic conversion due to the sharp velocity shear caused by the strong SC during the summer monsoon and the EICC from March to June. For the other regions, the leading source of EKE is the eddy potential energy (EPE), which is extracted from available potential energy of mean flow via baroclinic conversion, and then the EPE is converted into EKE through vertical density flux. Once generated, EKE will be redistributed by pressure work and advection via eddy energy flux, which varies in sync with the monthly variation of total EKE, transporting EKE to the adjacent region or deeper layer. From the vertical aspect, eddy energy conversions are more prominent above 200 m. The maximal EKE and barotropic conversion mostly occur at the surface, whereas the EPE and baroclinic conversion may have two peaks, which lie at the surface and in the thermocline. Using the satellite altimeter data and wind data, we further investigate the impact of geostrophic eddy wind work, which reveals a slightly dampening effect to EKE in the north Indian Ocean.

A new approach to construct 3-D crustal shear-wave velocity models: method description and application to the Central Alps

We develop a new inversion approach to construct a 3-D structural and shear-wave velocity model of the crust based on teleseismic P-to-S converted waves. The proposed approach does not require local earthquakes such as body wave tomography, nor a large aperture seismic network such as ambient noise tomography, but a three-component station network with spacing similar to the expected crustal thickness. The main features of the new method are: (1) a novel model parametrization with 3-D mesh nodes that are fixed in the horizontal directions but can flexibly vary vertically; (2) the implementation of both sharp velocity changes across discontinuities and smooth gradients; (3) an accurate ray propagator that respects Snell’s law in 3-D at any interface geometry. Model parameters are inverted using a stochastic method composed of simulated annealing followed by a pattern search algorithm. The first application is carried out over the Central Alps, where long-standing permanent and the temporary AlpArray Seismic Network stations provide an ideal coverage. For this study we invert 4 independent parameters, which are the Moho discontinuity depth, the Conrad discontinuity depth, the P-velocity change at the Conrad and the average Vp/Vs of the crust. The 3-D inversion results clearly image the roots of the Alpine orogen, including the Ivrea Geophysical Body. The lower crust's thickness appears fairly constant. Average crustal Vp/Vs ratios are relatively higher beneath the orogen, and a low-Vp/Vs area in the northern foreland seems to correlate with lower crustal earthquakes, which can be related to mechanical differences in rock properties, probably inherited. Our results are in agreement with those found by 3-D ambient noise tomography, though our method inherently performs better at localizing discontinuities. Future developments of this technique can incorporate joint inversions, as well as more efficient parameter space exploration.