Lateral Velocity Variations Research Articles

Pre-stack wave-equation time migration

In conventional time migration, the Green’s function is commonly calculated using ray tracing which may produce inaccurate images for complex structures with strongly lateral velocity variations. To mitigate this issue, we develop an accurate prestack wave-equation time migration method, borrowing from reverse time migration (RTM). First, a two-way acoustic wave equation is derived in the image-ray coordinate system. Then, we directly solve this equation using the finite-difference method to extrapolate source- and receiver-side wavefields, which does not require any approximation as in conventional ray-based time migration. Finally, a zero-lag crosscorrelation imaging condition is applied to generate the time migration image. In addition, we derive the dispersion relation and stability condition for the finite-difference method in the image-ray coordinate system, which provides guidance for temporal and spatial grid discretization. The perfectly matched layer (PML) boundary condition is utilized in the image-ray coordinate system to attenuate boundary reflections. Numerical examples of synthetic and field data verify the feasibility and adaptability of the proposed approach for low signal-to-noise land data. In addition, after mapping the images from the time domain to the depth domain, the image results are even better than RTM results.

A probabilistic full waveform inversion of surface waves

AbstractOver the past decades, surface wave methods have been routinely employed to retrieve the physical characteristics of the first tens of meters of the subsurface, particularly the shear wave velocity profiles. Traditional methods rely on the application of the multichannel analysis of surface waves to invert the fundamental and higher modes of Rayleigh waves. However, the limitations affecting this approach, such as the 1D model assumption and the high degree of subjectivity when extracting the dispersion curve, motivate us to apply the elastic full‐waveform inversion, which, despite its higher computational cost, enables leveraging the complete information embedded in the recorded seismograms. Standard approaches solve the full‐waveform inversion using gradient‐based algorithms minimizing an error function, commonly measuring the misfit between observed and predicted waveforms. However, these deterministic approaches lack proper uncertainty quantification and are susceptible to get trapped in some local minima of the error function. An alternative lies in a probabilistic framework, but, in this case, we need to deal with the huge computational effort characterizing the Bayesian approach when applied to non‐linear problems associated with expensive forward modelling and large model spaces. In this work, we present a gradient‐based Markov chain Monte Carlo full‐waveform inversion where we accelerate the sampling of the posterior distribution by compressing data and model spaces through the discrete cosine transform. Additionally, a proposal is defined as a local, Gaussian approximation of the target density, constructed using the local Hessian and gradient information of the log posterior. We first validate our method through a synthetic test where the velocity model features lateral and vertical velocity variations. Then we invert a real dataset from the InterPACIFIC project. The obtained results prove the efficiency of our proposed algorithm, which demonstrates to be robust against cycle‐skipping issues and able to provide reasonable uncertainty evaluations with an affordable computational cost.

Lateral Transport Controls the Tidally Averaged Gravitationally Driven Estuarine Circulation: Tidal Mixing Effects

Abstract In classic models of the tidally averaged gravitationally driven estuarine circulation, denser salty oceanic water moves up the estuary near the bottom, while less dense riverine water flows toward the ocean near the surface. Traditionally, it is assumed that the associated pressure gradient forces and salt advection are balanced by vertical mixing. This study, however, demonstrates that lateral (across the estuary width) transport processes are essential for maintaining the estuarine circulation. This is because for realistic estuarine bathymetry, the depth-integrated salt transport up the estuary is enhanced in the deeper estuary channel. A closed salt budget then requires the lateral transport of this excess salt in the deeper channel toward the estuarine flanks. To understand how such lateral transport affects the estuarine salt and momentum balances, we devise an idealized model with explicit lateral transport focusing on tidally averaged lateral mixing effects. Solutions for the along-estuary velocity and salinity are nondimensionalized to depend only on one single nondimensional parameter, referred to as the Fischer number, which describes the relative importance of lateral to vertical tidal mixing. For relatively strong lateral tidal mixing (greater Fischer number), salinity and velocity variations are predominantly vertical. For relatively weak lateral tidal mixing (smaller Fischer number), salinity and velocity variations are predominantly lateral. Overall, lateral transport greatly affects the estuarine circulation and controls the estuarine salinity intrusion length, which is demonstrated to scale inversely with the Fischer number.

Tomographic imaging of South Rewa basin, Central India: Implications of Gondwana rifting and Late Cretaceous volcanism

We have imaged thick (3.0–5.0 km) sub-trappean Gondwana sediments deposited in a typical rift-environment of the south Rewa basin in Central India masked by the Deccan basalts. The basalts intrude through the basin forming several dykes, sills, and other intrusions along with large-scale undulations of the basement showing complex geology and tectonic settings due to Late Cretaceous volcanism (∼65 Ma). To delineate complex subsurface geological structures and understand the tectonic settings, we use long-offset seismic reflection data along the N-S trending 127 km Kuvari-Shahdol profile of the basin. The conventional stack and pre-stack time migration techniques fail to image the complex subsurface structures dominated by steeply dipping faults. The steeply dipping faults with small-scale structures are influenced by strong lateral and vertical velocity variations that cause problems in imaging the basin. Hence, a robust tomographic inversion with pre-stack depth migration technique is adopted to image fine-scale subtle subsurface geological structures with smooth velocity variations. The tomographic velocity model and pre-stack depth migration seismic image could decipher basaltic trap thickness, numerous dyke intrusions, alternate horsts and grabens with deposition of thick hydrocarbon-bearing sub-trappean Gondwana rocks along with the basement configuration. The deep basinal faults with highly distorted basement structures represent the intense tectonic activity of this Gondwana rift-basin. The deep-seated hydrocarbon reservoir is discovered with an excellent trapping mechanism forming both pre-rift and syn-rift settings of the rift-basin. The seismic image is further constrained and corroborated by the inversion of residual Bouguer gravity anomaly data to obtain corresponding density model derived from tomographic velocity model. The large basement faults, horsts and grabens, several intrusives like dyke, sill, laccolith, and thick sub-trappean Gondwana formations with Deccan trap cover indicate complex geology and tectonic settings of south Rewa basin forming the Late Cretaceous Volcanic Province of India.

Two‐dimensional SH‐wave and acoustic P‐wave full waveform inversion: A Midland Basin case study

AbstractBuilding a reliable S‐wave velocity model remains challenging from P‐wave land seismic data as multi‐parameter elastic full waveform inversion which updates P‐ and S‐wave velocities simultaneously is not yet widely adopted due to its computational cost, highly nonlinear nature and noise contamination from land seismic data. To overcome this challenge, we propose implementing the SH‐wave full waveform inversion to obtain the S‐wave velocity model. The value of this approach has been examined by applying two‐dimensional time‐domain SH‐ and P‐wave acoustic full waveform inversion to a nine‐component three‐dimensional seismic survey acquired in the Midland Basin. The nine‐component seismic survey was acquired by using both vertical and horizontal component vibrators and receivers, which enables us to apply both SH‐ and P‐wave full waveform inversion to the raw and preprocessed shot gathers along a two‐dimensional line. The SH‐ and P‐wave full waveform inversion applied to the raw shot gathers provide more detailed P‐ and S‐wave velocities compared to the velocities from the stacking velocity analysis. Additionally, there are no problematic artefacts in the inverted S‐wave velocity model, which shows the stability of the SH‐wave inversion. Through the SH‐wave full waveform inversion from the preprocessed shot gathers, we reveal the lateral velocity variations in the Grayburg–San Andres interval, which coincides with the depositional environment changes suggested by the previous regional study mainly based on well log and cores. These variations demonstrate that SH‐wave full waveform inversion can identify additional geological features that are not observed in the inverted P‐wave velocity model. The inverted P‐ and SH‐wave velocities are validated by the flattening of the events observed in the common image gathers after Kirchhoff depth migration. The SH‐wave migration image further reveals the irregular geometries of carbonates in the Spraberry formation. Vp/Vs ratios calculated from the independently inverted P‐ and S‐wave velocity models show strong lateral variations in the Wolfcamp interval (key‐producing interval), which could be caused by the organic content variations between the reservoir and non‐reservoir rocks. The inversion results demonstrate that the SH‐wave full waveform inversion can be used to provide an S‐wave velocity model that is comparable to the P‐wave velocity model derived from acoustic full waveform inversion, which is widely used for P‐wave velocity model building from P‐wave land seismic data.

Improving prestack time migration by introducing a new velocity-related parameter: Parameter picking and 3D real data application

Prestack time migration (PSTM), a commonly used tool for seismic imaging, has been widely applied in 3D seismic data processing. However, the conventional PSTM algorithms use only one effective velocity parameter (i.e., root-mean-square [rms] velocity) for each imaging point, which may not be accurate when stronger lateral variations occur in seismic velocities. In this paper, we introduce a new parameter called the velocity variation factor that considers velocity variations in inhomogeneous media to improve PSTM. This new parameter, together with the rms velocity, describes the propagation Green’s function at an imaging point with two effective parameters rather than one effective parameter as in conventional PSTMs. This provides a more accurate traveltime calculation for the wave propagating through media with moderate lateral velocity variation. Unlike the conventional bending-ray PSTM, the additional effective parameter is fully independent of the rms velocities. We estimate the two effective parameters at each imaging point by flattening the neighboring image gathers with a global optimization algorithm. The objective function is built at each imaging point using a selective crosscorrelation-based time shift, which can quantitatively describe the slight bending of events in the local migrated gathers regardless of the quality of the gathers. We estimate the two effective parameters using the very fast simulated annealing algorithm and multiscale approach, thus avoiding the local minimum caused by the noises in the migrated gathers. We apply our two-parameter PSTM to a real 3D land data set to demonstrate its industrial applicability. A comparison of the new imaging result with the conventional prestack depth migration also is presented.

Deep crustal structure and tectonic implications of the Nansha Trough from wide-angle reflection and refraction seismic data

Deep crustal structure of the Nansha Trough (NT) plays a crucial role in understanding its crustal nature and magmatic activity. This study used the data from a wide-angle seismic profile (NDL1), which passes through the Jinghong seamount (JHSM) and Yinqing seamount (YQSM) along the trough. By using 2D forward ray-tracing and joint refraction and reflection travel-time inversion methods, we obtained the velocity structure of the NT. The model results show a thinned continental crust with a thickness of ∼9–18 km in the east of the trough, which is variable drastically below the seamounts. The upper crust slightly varies at a thickness of ∼5 km except at the seamounts, while the lower crust shows lateral variations in thickness and velocity, accompanied by high velocities below the JHSM due to the magmatic addition. Velocity models of the NDL1 also reveal crustal differences between the JHSM and YQSM. Combined with the magnetic anomaly and multichannel seismic data, we infer that the seamounts formed from multiple stages of magmatism after the spreading of the South China Sea and the NT are still active. High velocities in the lower crust result from the mantle upwelling and decompression melting of flexural responses of the NT, which may provide a genesis clue for the formation of the submarine seamounts.

Lateral variations of shear wave velocity (VS) profile and VS30 over gentle terrain

Lateral variations of shear wave velocity (VS) profile and time-averaged shear wave velocity in the upper 30 m (VS30) are important indicators of site condition and significant contributors to uncertainty in earthquake ground motion and seismic hazards. Few previous studies involved in quantifying these lateral variations, especially at short distances, due to limited data at closely spaced sites. This study quantitatively investigated the lateral variation of VS profile and VS30 over plain and piedmont terrains, focusing on short distances ranging from hundreds of meters to several kilometers. 2510 plain and 749 piedmont borehole profiles from the mainland China borehole profile database were employed to examine these variations using the variogram as the geo-statistics tool. The main findings include: 1) Over plain terrain, the variation in site conditions between sites does not significantly increase with separation distance within a certain range (typically from 1 km to 3–5 km). Over piedmont terrain, site condition variation steadily increases with distance. 2) The variation of site condition over piedmont terrain is higher than that over plain terrain. Nonetheless, for both terrains, the dispersion of site condition between sites at distance of <1 km is engineering negligible, and the correlation of site conditions between sites within 5 km remains strong; 3) For both terrains, the lateral variation of VS in deeper layers is higher than that in shallower layers. The study outcomes emphasize the importance of incorporating the information of nearby measurements in site condition estimations. The study outcomes can also be utilized in the optimization of borehole sampling design for microzonation projects and the evaluation of site condition dispersion at inferred sites.

Continental Fragments in the South China Block: Constraints From Crustal Radial Anisotropy

AbstractThe lithospheric architecture of the South China Block (SCB) is crucial to understanding the formation and evolution of this distinctive and highly reworked continental lithosphere with over 3 billion years of tectonic history. However, due to a lack of high‐resolution geophysical datasets, a detailed picture of the SCB lithosphere is absent, and fundamental questions regarding its formation, assembly, and subsequent reworking processes are actively debated. Assuming that unique deformation patterns due to such tectonic processes can be mapped by seismic anisotropy, we present a new crustal radially anisotropic shear‐wave velocity model along a 1500‐km seismic transect that spans the major tectonic domains of the SCB to characterize the past deformation processes. The new seismic models show significant lateral variations in seismic anisotropy and velocity, suggesting that the SCB consists of several separated (micro)continental blocks or terranes that likely have different origins and have survived the prolonged deformation history since the early formation of these continental fragments. Combining available geophysical datasets, we link individual crustal domains of distinct anisotropy to constrain the multiphase deformation processes of the SCB, including the early formation of the Proto‐Yangtze and Cathaysia Blocks, the assembly of the SCB, and the subsequent reactivation of the interior and extensive deformation that have formed the Basin‐and‐Range style tectonics in the Cathaysia Block. We suggest that relict continental fragments have played critical roles in the formation and reactivation of the SCB lithosphere.

Upper-plate conduits linked to plate boundary that hosts slow earthquakes

In shallow subduction zones, fluid behavior impacts various geodynamic processes capable of regulating slip behaviors and forming mud volcanoes. However, evidence of structures that control the fluid transfer within an overriding plate is limited and the physical properties at the source faults of slow earthquakes are poorly understood. Here we present high-resolution seismic velocity models and reflection images of the Hyuga-nada area, Japan, where the Kyushu-Palau ridge subducts. We image distinct kilometer-wide columns in the upper plate with reduced velocities that extend vertically from the seafloor down to 10–13 km depth. We interpret the low-velocity columns as damaged zones caused by seamount subduction and suggest that they serve as conduits, facilitating vertical fluid migration from the plate boundary. The lateral variation in upper-plate velocity and seismic reflectivity along the plate boundary correlates with the distribution of slow earthquakes, indicating that the upper-plate drainage system controls the complex pattern of seismic slip at subduction faults.

Non-equilibrium CFD simulation of the wet-to-dry expansion of the siloxane MM in a converging–diverging nozzle

Wet-to-dry organic Rankine cycles could generate 30% higher power outputs in the temperature range of 150 to 250 °C compared to existing single-phase cycles. Since the expansion is only partially wet, turboexpanders could potentially be applied provided that the wet portion of the expansion is confined to the stator to avoid erosion in the rotor. To assess the feasibility of achieving complete evaporation in the stator, two-dimensional non-equilibrium numerical simulations of the wet-to-dry expansion of siloxane MM in a covering–diverging nozzle are performed for the first time. The simulation setup is first validated against published experimental data, and a sensitivity study is conducted concerning the selected interphase models. The model is then applied to simulate expansions from inlet pressures ranging from 478 to 1250 kPa and vapour qualities from 0.1 to 0.5. Moreover, the droplet number density was varied between 1010 and 1014. The results show that the evaporation rate, the extent of non-equilibrium effects and the flow’s spatial uniformity are predominantly dependent on the droplet size. Expansions beginning with droplets smaller than 20 μm resulted in complete mixture evaporation and negligible non-equilibrium effects in almost all investigated cases. For larger droplets, ranging from 40 to 100 μm, full evaporation could only be achieved for inlet pressures above 750 kPa and inlet qualities above 0.3, whereas for lower pressures, the outlet vapour quality varied between 80 and 90%. For droplets larger than 200 μm, there is a significant delay in evaporation resulting in outlet quality typically between 40 and 70%. Larger droplet flows are characterised by substantial velocity slip, temperature difference, phase separation and lateral velocity variations. Having said this, droplet breakup analysis indicates that droplets larger than 100 μm are likely to undergo breakup, which could enhance the evaporation rate; however, this requires further investigation. In conclusion, high inlet pressures and high inlet qualities are preferred from the perspective of ensuring dry-vapour conditions at the nozzle outlet.

High-resolution P- and S-wave reflection studies of an intraplate structure: The Azambuja fault, Portugal

The Azambuja fault is a NNE trending structure located 50 km north of Lisbon, the capital and most populous city of Portugal. The fault has been considered as a possible source for the historical, large earthquakes. Understanding this fault is a priority in seismic hazard evaluation of this region. The fault has a clear morphological signature. Miocene and Pliocene sediments are tilted eastward and cut by steeply dipping mesoscale fault segments, presenting reverse and normal offsets with a net downthrow to the east. Neotectonic studies indicate a Quaternary slip on the fault of 0.05–0.06 mm/year. However, no direct evidence of the Azambuja fault affecting the Pleistocene or Holocene sediments was found so far. Here, we present the findings from high-resolution seismic reflection studies using both P- and S-waves over the Holocene deposits. The detection of small-throw faulting in ductile sediments is a challenging task. We show that multiple signatures, like perturbations in the reflection hyperbolae visible in shot and CMP gathers, interruptions of reflectors in stacked sections, lateral seismic velocity variations obtained by horizon velocity analysis, all at coincident locations, strongly suggest that the activity of the Azambuja fault has affected the Holocene sediments in the study area. The lateral velocity variations are corroborated by wavepath eikonal traveltime tomography and velocity analysis supported by seismic modeling. By means of 2D viscoelastic modeling, we explain the absence of fault-related diffractions and negligible back-scattered energy from the fault. Using data from nearby boreholes, we find that the 15 ka old alluvium cover has indeed been disturbed by the presence of shallow fault strands. Considering the estimated vertical throws and the empirical relationships between fault length, co-seismic rupture and magnitude, a slip rate of 0.07 mm/y, slightly larger than previously thought, is expected for this fault.

Lateral Variations of P-Wave Velocity from Deep Borehole Data in the Southern Apennines, Italy

We have selected 28 deep wells in the Southern Apennine area, most of which are located along and around the Val d’Agri Basin. The Southern Apennines, one of the most seismically active regions of the Italian peninsula, is a NE-verging fold-and-thrust belt characterised by the Meso–Cenozoic Apulia carbonate duplex system overlain by a thick column of Apennine carbonate platform and Lagonegro basin units. These units are unconformably covered by Neogene siliciclastic successions. Among the many Quaternary tectonic basins in the area, the Val d’Agri Basin is the most important intramontane depression, and is bordered by a ~ NW–SE-trending active fault system that represents one of the main seismogenic structures of the region. Moreover, the Val d’Agri Basin is the largest onshore oil field basin in Europe. In this context, we have analysed sonic log records from 28 deep wells and compared them with the corresponding stratigraphy and the other geophysical logs. We have obtained detailed measurements of the P-wave velocity (Vp) for each well from 0 to ~ 6 km depth, and found important lateral variations of Vp over very small distances. From these values, we have retrieved the densities of the main units crossed by the wells and the range of the overburden gradient in this area.

An effective scheme of joint migration inversion in the pseudo‐time domain

AbstractTraditional full‐waveform inversion is a non‐linear and ill‐posed inversion problem. To reduce the non‐linearity of it, joint migration inversion (joint migration inversion) was proposed as an alternative. Joint migration inversion tries to minimize the mismatch between measured and modelled reflection data. One key feature of joint migration inversion is its parameterization: two separate parameters, reflectivity (for the amplitudes of reflected events) and propagation velocity (for the phase effects). This separation helps to reduce the non‐linearity of the inversion. During joint migration inversion, with the velocity being updated, the reflectors in the updated image are also shifting in depth accordingly, this phenomenon is called depth–velocity ambiguity. This interaction between the two parameters during inversion is desired to keep the image time consistent with the measured data but may lead to non‐robustness of joint migration inversion due to the presence of local minima. Therefore, we propose a more robust joint migration inversion scheme, which parameterizes the models with vertical time, termed pseudo‐time joint migration inversion. In pseudo‐time, the updates of velocity will not result in the associated vertical location changes of reflectors in the estimated image. Instead, the reflectors are mainly getting more focused. One limitation is that the depth‐pseudo‐time conversion process assumes a simple linear relationship between depth and pseudo‐time, which might cause some artefacts in the converted models when there exist strong lateral velocity variations. One subsequent round of depth joint migration inversion is recommended to resolve this issue. We demonstrate the effectiveness of our proposed method with a two‐dimensional synthetic example in an extreme scenario, where the initial velocity model is homogeneous, a realistic offshore two‐dimensional synthetic example, a two‐dimensional field example from the Vøring basin in Norway and a simple three‐dimensional synthetic example. In all examples, pseudo‐time joint migration inversion manages to recover more reasonable updates in the inverted velocity and invert more focused reflectors in the inverted image, compared to depth joint migration inversion.

Seismic Survey in Lesser Himalayan Thrust Belt, Western Nepal

Two hundred km of 2D seismic survey was carried out at the Lesser Himalayan Thrust Belts in Dailekh district, western Nepal. The main motivation is to elucidate the geologic relationship between the known oil and gas seeps, subsurface structure, and stratigraphy in the area. This is a challenging task which is from its extreme structural and geological complexity such as thrust faulting, tight folding, steep dip layers, and strong lateral variations in seismic velocity. Seismic data were acquired with SERCEL 428XL system and processed by GEOEAST computer software. In order to increase signal-to-noise ratio (SNR), suppress interference, and search for optimum acquisition parameters, a series of comparative tests on the different charge depth and size, group interval, CDP fold, geophone array, and single high-sensitivity geophone were conducted. We also tested 2S3L (two lines shooting and three lines receiving) wide line profiling. The results indicate that single hole with charge depth of 12 m, 4-16 kg charge size (less charge size for the densely populated areas), single high-sensitivity geophone, and 1S2L wide line profiling with 132 folds are the optimum acquisition parameters. On the basis of comparative process experiment, data processing workflow consisting of data preparation, prestack denoising, amplitude compensation, deconvolution, tomography static correction, velocity analysis, residual static correction, CRS stack, poststack migration, prestack time migration (PSTM), and prestack depth migration (PSDM) was selected. Maybe affected by problem of conflicting dip in complex media, CRS stack section does not show satisfactory geological characteristics. PSTM profile has moderate signal-to-noise (S/N) ratio; the shallow, medium, and deep continuous reflections can be observed in section. More details of the geological structures can be observed in PSDM section, especially in medium and shallow layers (less than 3000 ms or 4000 m), but PSDM method is more expensive and highly time consuming than that of CRS stack and PSTM. So, the PSTM section can be reasonably used for geological interpretation. By reference to field mapping, thrust characteristics, and MT data, the final interpretation to the PSTM section identified the interfaces of 6 geological units (Paleoproterozoic Nabhisthan Fm., Paleoproterozoic Dubidanda Fm., Neogene to Late Cretaceous Surkhet group, Late Carboneferous to Early Cretaeous Gondwana group, Mesoproterozoic Upper Lakharpata group, and Lower Lakharpata group) and delineated Main Boundary Thrust (MBT), Ramgarh Thrust (RMT), Padukasthan Thrust (PT), and Dailekh Thrust (DT). The bottom of Surkhet group which is our top target zone is about 4250 meters deep.

An efficient 3D wave-equation pre-stack time migration for high-density and wide-azimuth data

Abstract The sharp increase in data size of the high-density and wide-azimuth seismic acquisition makes seismic processing routines very challenging. To reduce the turn-around time of pre-stack time migration (PSTM) for large-scale field data (terra to petabyte size), we present an efficient wave-equation PSTM that is an approximate version of the double-square-root (DSR) migration for moderate lateral velocity variation. By converting the DSR dispersion relation into the vertical-time domain and ignoring an offset-dip term, a novel common-offset ray parameter wave-equation PSTM is developed. Meanwhile, the common-offset ray parameter data (i.e. offset plane waves) can be formed by applying the linear Radon transform to the common middle point CMP) gathers. To enhance the signal-to-noise ratio of the conventional Radon transform, we integrate the Radon transform result over azimuth, yielding azimuth-independent offset plane waves that are well suited to irregular sampling and uneven offset distribution. Furthermore, the split-step Fourier method is introduced to extrapolate the wavefield efficiently. Since the pre-stack migration of common-offset ray parameter data is migrated independently, the common image gathers (CIGs) can be generated naturally by sorting the migrated sections by ray parameters. A 3D high-density field data test demonstrates the efficiency and effectiveness of the proposed method, making it very promising for industrial-sized data processing.

Surface-wave tomography using SeisLib: a Python package for multiscale seismic imaging

Summary To improve our understanding of the Earth’s interior, seismologists often have to deal with enormous amounts of data, requiring automatic tools for their analyses. It is the purpose of this study to present SeisLib, an open-source Python package for multiscale seismic imaging. At present, SeisLib includes routines for carrying out surface-wave tomography tasks based on seismic ambient noise and teleseismic earthquakes. We illustrate here these functionalities, both from the theoretical and algorithmic point of view and by application of our library to seismic data from North America. We first show how SeisLib retrieves surface-wave phase velocities from the ambient noise recorded at pairs of receivers, based on the zero crossings of their normalized cross-spectrum. We then present our implementation of the two-station method, to measure phase velocities from pairs of receivers approximately lying on the same great-circle path as the epicentre of distant earthquakes. We apply these methods to calculate dispersion curves across the conterminous United States, using continuous seismograms from the transportable component of USArray and earthquake recordings from the permanent networks. Overall, we measure 144 272 ambient-noise and 2055 earthquake-based dispersion curves, that we invert for Rayleigh-wave phase-velocity maps. To map the lateral variations in surface-wave velocity, SeisLib exploits a least-squares inversion algorithm based on ray theory. Our implementation supports both equal-area and adaptive parametrizations, with the latter allowing for a finer resolution in the areas characterized by high density of measurements. In the broad period range 4–100 s, the retrieved velocity maps of North America are highly correlated (on average, 96 per cent) and present very small average differences (0.14 ± 0.1 per cent) with those reported in the literature. This points to the robustness of our algorithms. We also produce a global phase-velocity map at the period of 40 s, combining our dispersion measurements with those collected at global scale in previous studies. This allows us to demonstrate the reliability and optimized computational speed of SeisLib, even in presence of very large seismic inverse problems and strong variability in the data coverage. The last part of the manuscript deals with the attenuation of Rayleigh waves, which can be estimated through SeisLib based on the seismic ambient noise recorded at dense arrays of receivers. We apply our algorithm to produce an attenuation map of the United States at the period of 4 s, which we find consistent with the relevant literature.

Analytical model for predicting the lateral profiles of velocities through a partially vegetated channel

In an open channel, an emergent canopy of vegetation enhances local resistance, producing a complicated velocity distribution within and around the canopy. Laboratory experiments were performed to understand the flow velocity variations in a partially vegetated channel, revealing both streamwise and lateral velocity variations as the flow entered the canopy. An analytical model was proposed to predict the lateral velocity profiles across the vegetated region and bare channel in both the flow adjustment region and the fully developed flow region. The flow velocity variations in the streamwise and lateral directions were parameterized and included in the new proposed model, allowing the model to capture both the streamwise and the lateral velocity variations. Twenty sets of experimental data from our experiments and data from published studies were used to verify the new model. The predicted velocity profiles exhibited good agreement with the measurements under a wide range of vegetation and flow conditions, suggesting that the proposed model is capable of accurately predicting the velocity profiles in a partially vegetated channel. Finally, the proposed model was applied to evaluate the critical frontal area per canopy volume at which Kelvin–Helmholtz vortices form along the emergent canopy side edge. With the same spatially averaged velocity, the critical frontal area exhibited no dependence on stem diameter.

Moment Tensor Solutions for Earthquakes in the Southern Korean Peninsula Using Three-Dimensional Seismic Waveform Simulations

Precise estimates of earthquake source properties are crucial for understanding earthquake processes and assessing seismic hazards. Seismic waveforms can be affected not only by individual event properties, but from the Earth’s interior heterogeneity. Therefore, for accurate constraints on earthquake source parameters, the effects of three-dimensional (3D) velocity heterogeneity on seismic wave propagation need evaluation. In this study, regional moment tensor solutions for earthquakes around the southern Korean Peninsula were constrained based on the spectral-element moment tensor inversion method using a recently developed high-resolution regional 3D velocity model with accurate high-frequency waveform simulations. Located at the eastern margin of the Eurasian plate, the Korean Peninsula consists of complex geological units surrounded by thick sedimentary basins in oceanic areas. It exhibits large lateral variations in crustal thickness (&amp;gt; 10 km) and seismic velocity (&amp;gt;10% dlnVs) at its margins in the 3D model. Seismic waveforms were analyzed from regional earthquakes with local magnitudes &amp;gt; 3.4 that occurred within and around the peninsula recorded by local broadband arrays. Moment tensor components were inverted together with event locations using the numerically calculated Fréchet derivatives of each parameter at periods ≥ 6 s. The newly determined solutions were compared with the results calculated from the one-dimensional (1D) regional velocity model, revealing a significant increase in a double-couple component of &amp;gt; 20% for earthquakes off of the coastal margins. Further, compared to initial solutions, ≤ 5 km change in depth was observed for earthquakes near the continental margin and sedimentary basins. The combination of a detailed 3D crustal model and accurate waveform simulations led to an improved fit between data and synthetic seismograms. Accordingly, the present results provide the first confirmation of the effectiveness of using 3D velocity structures for accurately constraining earthquake source parameters and the resulting seismic wave propagation in this region. We suggest that accurate 3D wave simulations, together with improved source mechanisms, can contribute a reliable assessment of seismic hazards in regions with complex continental margin structures and sedimentary basins from offshore earthquakes whose seismic waveforms can be largely affected by 3D velocity structures.

Crustal structure and tectonic implications of the southernmost Merida Andes, Venezuela, from wide-angle seismic data analysis

The crustal structure of the Merida Andes in southwestern Venezuela is crucial to our understanding of both, the complex regional geodynamics of the northern part of South America, as well as the development of the hydrocarbon basins located along its flanks. Symmetric and asymmetric models of the Merida Andes have been proposed based on geological and geophysical (principally gravimetric) surveys, without counting detailed structures from wide-angle seismic. This work examines the crustal structure of the southwestern part of the Merida Andes in Venezuela using controlled seismic sources along a 500 km wide-angle seismic profile (named ANDES SUR). Forward modeling is carried out to obtain a layer-based, P-wave, crustal velocity model. An alternative velocity model is obtained by applying a first arrival tomography for the upper 20 km depth, and a joint PmP and Pn inversion for the deeper crust. Important lateral P-wave velocity variations and low-velocity zones (LVZ) are observed within the upper crust. To the northwest, these lateral velocity variations are probably related to the deformation zones of Icotea and Tarra Faults, and the southeast to the Mantecal and Apure grabens. The integration of seismological data to the profile provides some clues about a possible slow-shear movement between the upper and the lower crust (brittle-ductile transition zone) at the NW end of the profile, as well as the possible northwestward incipient subduction of the South American Plate beneath the Maracaibo Triangle Block. We interpret the Boconó Fault's main trace as a take-off surface detaching from the interface between the upper and lower crust and dipping to the SE. A slight NW displaced crustal root of the Merida Andes is observed in all models, being the root ∼150 km wide with an approximately 9 km downward deflection (from ∼46 to 55 km depth) of the Moho. The crustal structure along the profile ANDES SUR, together with the observations along two profiles located further north, provide information about the tectonic implication along the strike of the Merida Andes, favoring an asymmetric model with possible NW directed South American lower crust and upper mantle underthrusting or incipient subduction beneath the Maracaibo Triangle Block beneath the center of the Merida Andes.