Velocity Model Research Articles

Extracting Multimodal Surface-Wave Dispersion Curves from Ambient Seismic Noise Using High-Resolution Linear Radon Transform

Abstract In the past two decades or so, ambient noise tomography (ANT) has emerged as an established method for imaging subsurface seismic velocity structures. One of the key steps in ANT is to extract surface-wave dispersion curves. The predominant approach for subsurface shear-wave velocity structure inversion involves utilizing fundamental-mode surface waves in ANT. Nevertheless, a notable challenge encountered is the issue of nonuniqueness when employing the dispersion information of fundamental-mode surface waves to invert for shear-wave velocity models. The inclusion of higher-mode dispersion curves in the inversion offers several benefits, including the reduction of nonuniqueness, enhancement of inversion stability, and decreased dependence on the initial model. In this study, we illustrate the applicability of the high-resolution linear radon transform method (HRLRT) for extracting multimodal surface-wave dispersion energy from ambient seismic noise data. We apply HRLRT to both the synthetic noise data and real data recorded by USArray. Our results of applications show that the HRLRT method can extract multimodal surface-wave dispersion information. Compared with established methods such as the frequency–Bessel transform and multicomponent frequency–Bessel transform, the HRLRT exhibits an advantage in suppressing “crossed” artifacts and the second-/third-type artifacts caused by sparse spatial sampling, and the resulting dispersion energy from HRLRT has narrower peaks, meaning high resolution of dispersion curves based on the HRLRT method.

Investigation of mesoscopic boundary conditions for lattice Boltzmann method in laminar flow problems

For use with the lattice Boltzmann method, the macroscopic boundary conditions need to be transformed into their mesoscopic counterparts. Commonly used mesoscopic boundary conditions use the equilibrium density function, which introduces undesirable artifacts into the numerical solution, especially near interfaces with other types of boundary conditions. In this work, several variants of the mesoscopic boundary conditions are summarized and numerically investigated by means of benchmark problems for the incompressible Navier-Stokes equations with known analytical solutions. To improve the numerical approximation of the velocity and pressure fields, moment-based boundary conditions are extended for the D3Q27 velocity model. Furthermore, the interpolated boundary conditions are improved. These newly developed boundary conditions are shown to produce results with a substantially smaller numerical error.

Machine-learning based location of the 2021 MW 7.4 Maduo, Qinghai, China earthquake sequence: Insight into intraplate seismogenesis

On 22 May 2021, an MW 7.4 earthquake occurred in Maduo County, Qinghai Province, China, which is located on the Kunlun Mountain Pass-Jiangcuo fault inside the Bayan Har block, providing a good opportunity to investigate seismogenesis of large intraplate earthquakes. We analyze two years of continuous seismic data from June 2021 to June 2023, which were recorded at 34 portable seismic stations of the MaduoArray deployed in the source zone by our group. The LOC-FLOW workflow of automatic detection and location is applied to construct a complete and high-precision seismic catalog for the region, which includes machine-learning phase picking (PhaseNet), earthquake phase association (REAL), velocity model updating and station correction (VELEST), absolute earthquake location (HypoInverse), and relative location (HypoDD). As a result, 78,832 Maduo aftershocks and other local earthquakes are detected and relocated precisely. Our results show that the length of the Maduo aftershock zone is ∼170 km, which is mainly distributed along the NWW-SEE oriented Kunlun Mountain Pass-Jiangcuo fault, and there is a horsetail bifurcation feature at the eastern end of the aftershock sequence. The seismogenic fault is nearly vertical, and local seismicity occurs on both sides of the fault. Our results also show that there is no seismic gap or aftershock sparse area in the region. Previous studies have revealed a low-velocity and high-conductivity anomaly below the source zone, reflecting fluids ascending from the lower crustal flow. These results provide new insights into the cause of the 2021 Maduo earthquake.

Acoustic wave simulation in strongly heterogeneous models using a discontinuous Galerkin method

In recent years, the discontinuous Galerkin method (DGM) has been rapidly developed for the numerical simulation of seismic waves. For wavefield propagation between two adjacent elements, it is common practice to apply a numerical flux to the boundary of each element to propagate waves between adjacent elements. Several fluxes, including the center, penalty, local Lax-Friedrich (LLF), upwind, and Rankine-Hugoniot jump condition-based (RH condition) fluxes, are widely used in numerical seismic wave simulation. However, some fluxes do not account for media differences between adjacent elements. Although different fluxes have been successfully used in DGM for many velocity models, it is unclear whether they can produce sufficiently accurate or stable results for strongly heterogeneous models, such as checkerboard models commonly used in tomographic studies. We test different fluxes using the acoustic wave equation. We analyze the accuracy of the penalty, LLF, upwind, and RH-condition fluxes based on the results of the numerical simulations of the homogeneous and two-layer models. We conduct simulations using checkerboard models, and the results indicate that the LLF, penalty, and upwind fluxes may have instability problems in heterogeneous models with long-time simulations. We observe instability issues in the LLF, penalty, and upwind fluxes when the wave-impedance contrast is high at the media interface. However, the results of the RH-condition flux remain consistently stable. The series of numerical examples presented in this work provide insights into the characteristics and application of fluxes for seismic wave modeling.

Self‐organizing cooperative hunting for unmanned surface vehicles with constrained kinematics

SummaryThe article aims at solving a cooperative hunting problem for multiple unmanned surface vehicles (USVs) subject to constrained kinematics. In order to cooperatively trap the evader into the hunting domain, a velocity model with control variable for the pursuers is firstly proposed according to the Apollonius circle. Then, a flexible self‐organizing control strategy is developed, which enables the pursuers to approach the evader while forming an encirclement. The pursuers can dynamically adapt their strategies in real‐time by choosing the optimal control variable. Additionally, take into account the limitation imposed on the vessel's motion, the optimal control variable with constraint can be obtained by using the particle swarm optimization with log‐barrier method. The simulation results ultimately demonstrate the validity and superiority of the proposed cooperative hunting algorithm.

Analysis of Bubble Flow in an Inclined Tube and Modeling of Flow Prediction

The lubricating oil system is a significant component of aviation engine lubrication and cooling, and the scavenge pipe is an essential component of the lubricating oil system. Accurately identifying and understanding the flow state of the scavenge pipe is very important. This article establishes a visualization test bench for a 45-degree inclined scavenge pipe, with upward and downward flow directions, respectively. The test temperature is 370 K, and a high-speed camera captures the changes in the two-phase flow inside the pipeline. Based on high-speed photography photos, we develop software for analyzing the flow characteristics of bubbles inside the tube and explore the influence of gas phase conversion velocity and liquid phase conversion velocity on the apparent velocity of bubbles inside the tube. Multiple algorithms were used to develop the model by combining machine learning with speed and accuracy to establish a data regression prediction model for the apparent velocity of bubbles inside the tube. Through calculation and analysis, it was found that the root mean square error of the prediction model using the BP neural network algorithm was the lowest, and the decision coefficient of the prediction model using the support vector machine algorithm was the highest.

An Improved Longitudinal Driving Car-Following System Considering the Safe Time Domain Strategy.

Car-following models are crucial in adaptive cruise control systems, making them essential for developing intelligent transportation systems. This study investigates the characteristics of high-speed traffic flow by analyzing the relationship between headway distance and dynamic desired distance. Building upon the optimal velocity model theory, this paper proposes a novel traffic car-following computing system in the time domain by incorporating an absolutely safe time headway strategy and a relatively safe time headway strategy to adapt to the dynamic changes in high-speed traffic flow. The interpretable physical law of motion is used to compute and analyze the car-following behavior of the vehicle. Three different types of car-following behaviors are modeled, and the calculation relationship is optimized to reduce the number of parameters required in the model's adjustment. Furthermore, we improved the calculation of dynamic expected distance in the Intelligent Driver Model (IDM) to better suit actual road traffic conditions. The improved model was then calibrated through simulations that replicated changes in traffic flow. The calibration results demonstrate significant advantages of our new model in improving average traffic flow speed and vehicle speed stability. Compared to the classic car-following model IDM, our proposed model increases road capacity by 8.9%. These findings highlight its potential for widespread application within future intelligent transportation systems. This study optimizes the theoretical framework of car-following models and provides robust technical support for enhancing efficiency within high-speed transportation systems.

On accounting for the effects of crust and uppermost mantle structure in global scale full-waveform inversion

Summary Fundamental mode surface wave data have often been used to construct global shear velocity models of the upper mantle under the so-called “path average approximation”, an efficient approach from the computational point of view. With the advent of full-waveform inversion and numerical wavefield computations, such as afforded by the spectral element method, accounting for the effects of the crust accurately becomes challenging. Here, we assess the merits of accounting for crustal and uppermost mantle effects on surface and body waveforms using fundamental mode dispersion data and a smooth representation of the shallow structure. For this we take as reference a model obtained by full waveform inversion and wavefield computations using the spectral element method, model SEMUCB-WM1 (French and Romanowicz, 2014) and compare the waveform fits of synthetics to different parts of three component observed teleseismic records, in the period band 32-300 s for body waves and 40-300 s for surface waves and their overtones for three different models. The latter are: a dispersion-only based model (model Disp_20s_iter5), and two models modified from SEMUCB-WM1 by successively replacing the top 200 km (model Merged _200km) and top 80 km (model Merged _80km), respectively, by a model constrained solely by fundamental mode surface wave dispersion data between periods of 20 and 150 s. The crustal part of these 3 models (resp. SEMUCB-WM1) is constrained from dispersion data in the period range 20-60 s (resp. 25-60 s), using the concept of homogenization (e.g., Backus 1962, Capdeville & Marigo 2007) which is tailored to simplify complex geological features, enhancing the computational efficiency of our seismic modeling. We evaluate the fits to observed waveforms provided by these 3 models compared to those of SEMUCB-WM1 by computing three component synthetics using the spectral element method for 5 globally distributed events recorded at 200+ stations, using several measures of misfit. While fits to waveforms for model 3 are similar to those for SEMUCB-WM1, the other two models provide increasingly poorer fits as the distance travelled by the corresponding seismic wave increases and/or as it samples deeper in the mantle. In particular, models 1 and 2 are biased towards fast shear velocities, on average. Our results suggest that, given a comparable frequency band, models constructed using fundamental mode surface wave data alone and the path average approximation, fail to provide acceptable fits to the corresponding waveforms. However, the shallow part of such a 3D radially anisotropic model can be a good starting model for further full waveform inversion using numerical wavefield computations. Moreover, the shallow part of such a model, including its smooth crustal model, and down to a maximum depth that depends on the frequency band considered, can be fixed in FWI iterations for deeper structure. This can save significant computational time when higher resolution is sought in the deeper mantle. In the future, additional constraints for the construction of the homogenized model of the crust can be implemented from independent short period studies, either globally or regionally.

Full‐waveform inversion as a tool to predict fault zone acoustic properties

AbstractUnderstanding the physical properties of fault zones is essential for various subsurface applications, including carbon capture and geologic storage, geothermal energy and seismic hazard assessment. Although three‐dimensional seismic reflection data can image the geometries of faults in the sub‐surface, it does not provide any direct information on the physical properties of fault zones. We currently cannot use seismic reflection data to infer directly which faults may be leaking or sealing and are reliant instead on shale‐gauge ratio type calculations, which are fraught with uncertainties. In this paper, we propose that full‐waveform inversion P‐wave velocity models can be used to extract information on fault zone acoustic properties directly, which may be a proxy for subsurface fault transmissibility. In this study, we use high‐quality post‐stack depth–migrated seismic reflection and full‐waveform inversion velocity data to investigate the characteristics of fault zones in the Samson Dome in the SW Barents Sea. We analyse the variance attribute of the post‐stack depth migrated and full‐waveform inversion volumes, revealing linear features that consistently appear in both datasets. These features correspond to locations of rapid velocity changes and seismic trace distortions, which we interpret as faults. These observations demonstrate the capability of full‐waveform inversion to recover fault zone velocity structures. Our findings also reveal the natural heterogeneity and complexity of fault zones, with varying P‐wave velocity anomalies within the studied fault network and along individual faults. Our results indicate a correlation between P‐wave velocity anomalies within fault zones and the modern‐day stress orientation. Faults with high P‐wave velocity are the ones that are perpendicular to the present‐day maximum horizontal stress orientation and are likely under compression. Faults with lower P‐wave velocity are the ones more parallel to the present‐day maximum horizontal stress orientation and are likely in extension. We propose that these P‐wave velocity anomalies may indicate differences in how ‘open’ and fluid filled the fault zones are (i.e. faults in extension are more open, more fluid filled and have lower VP) and therefore may provide a promising proxy for fault transmissibility.

A fast and robust two‐point ray tracing method in layered vertical transversely isotropic media with strong anisotropy

AbstractA fast and robust two‐point ray tracing method was developed for layered vertical transversely isotropic media with strong anisotropy. Utilizing the Christoffel slowness equation, a novel generalized dimensionless ray parameter, , modified from the ray parameter (horizontal slowness), was proposed to efficiently and simultaneously determine the ray paths and travel times for direct and reflected quasi‐P, quasi‐SV and quasi‐SH waves. The Newton optimization algorithm was employed to solve the nonlinear offset equation accurately, resulting in rapid convergence to the true value. The inferred analytical equations show that the generalized ray parameter stabilizes the inversion process at large offsets. Additionally, a piecewise function was introduced to enhance the initial value estimation and calculation efficiency. The numerical results demonstrate that this novel approach can reduce the iteration error to 10−10 m in less than three iterations. Monte Carlo simulations further validated the effectiveness of the method for inferring the true ray paths at various offsets within complex velocity models. Furthermore, the method can address the triplication issue in quasi‐SV waves and exhibit robustness in strong‐layered vertical transversely isotropic media.

Deep low velocity layer in the sublithospheric mantle beneath India

SUMMARY Globally, there is now a growing evidence for a low velocity layer in the deeper parts of the upper mantle, above the 410 km discontinuity (hereafter called LVL-410). The origin of this layer is primarily attributed to interaction of slabs or plumes with a hydrous mantle transition zone (MTZ) that results in dehydration melting induced by water transport upward out of the MTZ. However, the ubiquitous nature of this layer and its causative remain contentious. In this study, we use high quality receiver functions (RFs) sampling diverse tectonic units of the Indian subcontinent to identify Ps conversions from the LVL-410. Bootstrap and differential slowness stacking of RFs migrated to depth using a 3-D velocity model reveal unequivocal presence of a deep low velocity layer at depths varying from 290 to 400 km. This layer appears more pervasive and deeper beneath the Himalaya, where detached subducted slabs in the MTZ have been previously reported. Interestingly, the layer is shallower in plume affected regions like the Deccan Volcanic Province and Southern Granulite Terrane. Even though a common explanation does not appear currently feasible, our observations reaffirm deep low velocity layers in the bottom part of the upper mantle and add to the list of regions that show strong presence of such layers above the 410 km discontinuity.

Structure of the crust and upper mantle in Greenland and northeastern Canada: Insights from anisotropic rayleigh-wave tomography

Summary Seismic velocity models provide important constraints on Greenland’s deep structure, which, in turn, has profound implications for our understanding of the tectonic history of this region. However, the resolution of seismic models has been limited by a sparse network, particularly in northern and central Greenland. We address these limitations by generating new high-resolution Rayleigh-wave phase velocity maps encompassing Greenland and northeastern Canada by processing over three decades of teleseismic earthquake records and incorporating recently added stations in Greenland and Arctic Canada. These phase velocity maps are sensitive to structure from the lower crust down to the sub-lithospheric mantle (25-185 s period). We find significant heterogeneity and a strong correlation between isotropic and anisotropic seismic velocities with inferred geological structure. High seismic velocities associated with cratonic lithosphere are broadly divided into two regions, with a belt of reduced velocity spanning central Greenland, which we interpret as lithospheric erosion resulting from interaction between the Greenland continental keel and the Iceland plume. Within each region, we identify tectonic subdivisions that suggest fundamental differences between the blocks that make up Precambrian Greenland. In the south, the North Atlantic craton (NAC) has a high-velocity keel exhibiting anisotropic stratification. Between the NAC and the cratonic lithosphere further north, the Proterozoic Nagssugtoqidian orogenic belt shows a distinct signature of reduced seismic velocity to ∼75 s period, but then appears to pinch out at depth. The northern Greenland lithosphere exhibits significant isotropic heterogeneity, with a distinct core of high velocities in the northwest (∼55-75 s period) giving way to a set of distinct east-west trending high-velocity belts at longer periods. At all periods sensitive to the lithospheric mantle in this region, anisotropic fast orientations are E-W, consistent with a north-south Precambrian assembly of the Greenland shield. In contrast to the NAC, there is no evidence of anisotropic stratification in the northern part of the cratonic keel. Based on both isotropic and anisotropic phase-velocity anomalies, we suggest that the Phanerozoic Caledonian and Ellesmerian-Franklinian fold belts are relatively thin-skinned features onshore Greenland, though the Caledonian belt may have a stronger signature off the east coast. At the longest periods, a prominent low-velocity anomaly initially centred on Iceland migrates northwards and spreads beneath central-eastern Greenland. Coupled with NW-SE trending anisotropy, this feature is interpreted as the effect of mantle flow radiating outward from the Iceland plume and interacting with the eroded Greenland lithosphere.

Full-waveform automatic location of small seismic events in an underground mine using synthetic strain Greens tensors

A reliable method to locate microseismic evets is an important objective for routine mine seismic monitoring. We developed an algorithm that can automatically and reliably locate an event from a single uniaxial trace. Signals created by microseismic events are recorded as seismograms that consist of direct P- and S-waves as well as the waves reflected off the underground excavations. We take advantage of the extra complexity these reflections introduce as they make the signal more unique. In our investigation we cross-correlated unaltered synthetic seismograms against trial waveforms created from a linear combination of appropriate waveforms from a library of numerically calculated strain Green's tensors. We use a single uniaxial synthetic seismogram to invert for the source location, six moment tensor components, and frequency. The resulting cost functions are constructed as maximal correlations at each node in the mine grid through the true source x-planes. This function is maximized by differential evolution to simultaneously yield the source parameters. Sensitivity tests such as shifting the mine plans and altering the background velocities of the media are investigated to test the robustness of the method. We contaminate the data with 40% white noise to test the resolution against simulated real recorded seismograms. The method is demonstrated using synthetic data calculated with a realistic underground 3D velocity model.

A fine‐tuning workflow for automatic first‐break picking with deep learning

AbstractFirst‐break picking is an essential step in seismic data processing. For reliable results, first arrivals should be picked by an expert. This is a time‐consuming procedure and subjective to a certain degree, leading to different results for different operators. In this study, we have used a U‐Net architecture with residual blocks to perform automatic first‐break picking based on deep learning. Focusing on the effects of weight initialization on first‐break picking, we conduct this research by using the weights of a pre‐trained network that is used for object detection on the ImageNet dataset. The efficiency of the proposed method is tested on two real datasets. For both datasets, we pick manually the first breaks for less than 10 of the seismic shots. The pre‐trained network is fine‐tuned on the picked shots, and the rest of the shots are automatically picked by the neural network. It is shown that this strategy allows to reduce the size of the training set, requiring fine‐tuning with only a few picked shots per survey. Using random weights and more training epochs can lead to a lower training loss, but such a strategy leads to overfitting as the test error is higher than the one of the pre‐trained network. We also assess the possibility of using a general dataset by training a network with data from three different projects that are acquired with different equipment and at different locations. This study shows that if the general dataset is created carefully it can lead to more accurate first‐break picking; otherwise, the general dataset can decrease the accuracy. Focusing on near‐surface geophysics, we perform traveltime tomography and compare the inverted velocity models based on different first‐break picking methodologies. The results of the inversion show that the first breaks obtained by the pre‐trained network lead to a velocity model that is closer to the one obtained from the inversion of expert‐picked first breaks.

Ultrasonic measurement method for three-dimensional assembly stress of aero-engine rotors

The precise measurement of three-dimensional stress in the rotor is a key means to ensure the assembly accuracy and operational safety of aero-engines. The unclear coupling influence mechanism of different directions of stress on ultrasonic propagation makes it challenging to establish a three-dimensional stress measurement model. This paper proposes a three-dimensional assembly stress measurement method for aero-engine rotors based on ultrasonic propagation time difference decoupling. The three-dimensional assembly stress is converted into the stress component parallel and perpendicular to the ultrasonic wave propagation direction in this method. The acoustoelastic equation including three-dimensional stress parameters is established based on the function model of different direction stress and ultrasonic wave propagation velocity. An innovative three-transducer stress measurement scheme is designed, which integrates ultrasonic wave transmitting and receiving transducers and time difference feedback transducers. The acoustoelastic equations in three propagation directions are constructed by rotating the transducer group to change the propagation direction of ultrasonic waves. Finally, the acoustoelastic equations are solved simultaneously to decouple the ultrasonic propagation time difference, thus achieving a three-dimensional assembly stress measurement of the aero-engine rotor. The experimental results show that the maximum deviation between the measured values and the simulation results in the x ,y, and z directions is 15.885 MPa, 14.113 MPa, and 17.093 MPa, respectively. This study provides effective theoretical and technical support for precisely measuring three-dimensional assembly stress in high-end equipment represented by aero-engines.

Estimation of antigorite wave velocities in subduction conditions based on first-principles thermoelasticity

The most abundant serpentine mineral in subduction settings, antigorite has one of the highest water storage capacities and is involved in seismicity. Seismic wave velocities of antigorite are important for detecting and quantifying serpentinization within the mantle wedge and the subducting oceanic plate. At present, the elastic properties of antigorite at high pressures and temperatures are unclear. In this study, we have investigated pressure-volume-temperature (P-V-T) data and thermodynamic properties of antigorite using first-principles molecular dynamics (FPMD) simulations. Using these simulations results, we computed the relevant thermoelastic parameters and estimated compressional and shear wave velocities (vP and vS) of antigorite in subduction conditions. A simplified velocity model of antigorite with its coexisting mantle anhydrous phases was introduced to help us understand the potential effect of serpentinization on the seismic velocity of mantle rocks. Combined with seismic observations, we re-evaluated some velocity anomalies within forearc mantle wedges and established reliable serpentinization budgets. These results can provide preliminary evaluations and reliable constraints on serpentinization and water content in mantle rocks, which has important implications for understanding global plate dynamics and the deep water cycle.

Kinematic Modeling of Pitch Velocity in High School and Professional Baseball Pitchers: Comparisons With the Literature.

Kinematic parameters predictive of pitch velocity have been evaluated in adolescent and collegiate baseball pitchers; however, they have not been established for high school or professional pitchers. To create multiregression models using anthropometric and kinematics features most predictive for pitch velocity in high school and professional pitchers and compare them with prior multiregression models evaluating other playing levels. Descriptive laboratory study. High school (n = 59) and professional (n = 337) baseball pitchers threw 8 to 12 fastballs while being evaluated with 3-dimensional motion capture (480 Hz). Using anthropometric and kinematic variables, multiregression models for pitch velocity were created for each group. A systematic review was conducted to determine previous studies that established kinematic models for ball velocity in youth, high school, and collegiate pitchers. Leg length was predictive of pitch velocity for high school and professional pitchers (P < .001 for both). When compared with previously established models for pitch velocity, almost all groups were distinct from one another when assessing age (P maximum < .001), weight (P max = .0095), and pitch velocity (P max < .001). Stride length was a significant predictor for the youth/high school pitchers, as well as the current study's high school and professional pitchers (P < .001 for all). Maximal shoulder external rotation (collegiate: P = .001; professional: P < .001) and maximal elbow extension velocity (high school/collegiate: P = .024; collegiate: P < .001; professional: P = .006) were shared predictors for the collegiate and current study's professional group multiregression models. Trunk flexion at ball release was a commonly shared predictor in the youth/high school (P = .04), high school/collegiate (P = .003), collegiate (P < .001), and the current study's professional group (P < .001). Youth, high school, collegiate, and professional pitchers had unique, predictive kinematic and anthropometric features predictive of pitch velocity. Leg length, stride length, trunk flexion at ball release, and maximal shoulder external rotation were predictive features that were shared between playing levels. Coaches, clinicians, scouts, and pitchers can consider both the unique and the shared predictive features at each playing level when attempting to maximize pitch velocity.

TOI-1408: Discovery and Photodynamical Modeling of a Small Inner Companion to a Hot Jupiter Revealed by Transit Timing Variations

We report the discovery and characterization of a small planet, TOI-1408 c, on a 2.2 day orbit located interior to a previously known hot Jupiter, TOI-1408 b (P = 4.42 days, M = 1.86 ± 0.02 M Jup, R = 2.4 ± 0.5 R Jup) that exhibits grazing transits. The two planets are near 2:1 period commensurability, resulting in significant transit timing variations (TTVs) for both planets and transit duration variations for the inner planet. The TTV amplitude for TOI-1408 c is 15% of the planet’s orbital period, marking the largest TTV amplitude relative to the orbital period measured to date. Photodynamical modeling of ground-based radial velocity (RV) observations and transit light curves obtained with the Transiting Exoplanet Survey Satellite and ground-based facilities leads to an inner planet radius of 2.22 ± 0.06 R ⊕ and mass of 7.6 ± 0.2 M ⊕ that locates the planet into the sub-Neptune regime. The proximity to the 2:1 period commensurability leads to the libration of the resonant argument of the inner planet. The RV measurements support the existence of a third body with an orbital period of several thousand days. This discovery places the system among the rare systems featuring a hot Jupiter accompanied by an inner low-mass planet.

Spatial Distribution of Tremor Episodes From Long‐Term Monitoring in the Northern Cascadia Subduction Zone

AbstractLarge bursts of non‐volcanic tremor (“major” tremor episodes) correlated with geodetic deformation recur regularly in the Cascadia subduction zone and are often called episodic tremor and slip (ETS). Minor episodes of tremor between ETS are ubiquitous but have been understudied. This paper assesses time‐invariant characteristics of tremor episodes of all sizes within northern Cascadia. We derive a catalog of tremor episodes ranging in size from 10 to &gt;13,000 tremor events using the results of 17 years of tremor monitoring. Minor episodes represent ∼96% of all 896 tremor episodes and their occurrence varies on 10‐km scales. Using estimates for the depth of the forearc Moho and subducting slab, we observe an association between the location of the forearc mantle corner (FMC) and tremor occurrence that leads to along‐dip modality. Bimodality, present in southern Washington and Vancouver Island, represents the segmentation of major and minor episodes up‐dip and down‐dip of the FMC, respectively. Unimodality, present in Puget Sound, results when the FMC is located near the down‐dip edge of the ETS zone and no segmentation occurs. We also use our extensive tremor episode catalog alongside three‐dimensional regional tomographic velocity models to reassess the relationship between tremor activity and low Vp/Vs signatures in the forearc. We do not find a correlation between tremor episode recurrence intervals and Vp/Vs, contrary to some previous work, suggesting that controls on silica precipitation in the forearc crust are not dominant controls of tremor episode recurrence, or that the association is not widely observable.

Compressional and shear wave velocities of Fe-bearing silicate post-perovskite in Earth’s lowermost mantle

The bridgmanite (Bgm) to silicate post-perovskite (PPv) phase transition is believed to be the main cause for the distinct seismic features observed in the D'' layer, the lowermost region of the Earth’s mantle. However, the transition depth and elasticity of the PPv phase have been highly debated, as the chemical complexity within the D'' layer can largely affect the Bgm-PPv transition pressure and the associated velocity contrast. Experimental measurements of sound velocities for PPv with different chemical compositions under relevant lowermost-mantle conditions are essential but remain limited. In this study, we have reliably measured both compressional wave velocity (VP), shear wave velocity (VS), and density, for two Fe-bearing PPv compositions [(Mg0.85Fe0.15)SiO3 and (Mg0.75Fe0.25)SiO3] at lowermost mantle pressures using Impulsive Stimulated Light Scattering (ISS), Brillouin Light Scattering (BLS), and X-ray Diffraction (XRD) in diamond anvil cells. Our results indicate that the velocities of Fe-bearing PPv at 120 GPa can be described by the following relationships: VS (km/s) = 7.65–2.8XFe and VP (km/s) = 14.11–3.8XFe, where XFe represents mole fraction of the Fe content. The variations in the Fe content of PPv may provide one of the explanations for the seismic lateral variations observed at the Earth’s core mantle boundary. By comparing our results with the high-pressure velocities of Bgm, our velocity model suggests significant discontinuities across the Bgm-PPv transition, characterized by a reduction in both VP and VΦ, and an increase in VS. These findings highlight the importance of considering the influence of chemical composition, particularly Fe content which could vary significantly at the D'' region, on the seismic properties of the PPv phase. The observed velocity contrasts across the Bgm-PPv transition may contribute to the complex seismic signatures observed in the D'' layer, underscoring the potential role of this phase transition in interpreting the seismic features of the lowermost mantle region.