Velocity Anomalies Research Articles

Using time-lapse seismic velocity changes to monitor the Domo de San Pedro Geothermal field, Mexico

The Domo de San Pedro Geothermal field, located in the westernmost part of the trans-Mexican volcanic belt, is part of a silicic volcanic complex. Geothermal energy is produced since 2015 from 9 wells as deep as 3 km. We monitor geothermal operations using a continuous seismic dataset from a temporary network of 20 broadband stations from March 2021 to January 2022. The well logs describing the geothermal operations allow us to investigate the seismic response of the geothermal field to fluid extraction and re-injection across the geothermal wells. We produced a catalog of 217 local earthquakes and derive an optimized 1D velocity model for the geothermal field between 1.2 m above the sea level (masl) and 12 km depth. Apparent velocity changes (dv/v) are measured using seismic ambient noise interferometry from the coda of noise cross-correlations. Next, we identify and locate subsurface velocity anomalies, combining the methods of station collapse sensitivity and ambient noise imaging applied to the continuous seismic noise records. We distinguish three periods of seismic activity where prominent changes in velocity coincide with the re-injection of fluids and significant variation in the production rate of the geothermal wells. We determine the distribution of these velocity transients inside the geothermal field and propose that they are associated with variations in reservoir fluid content and effective stress changes beneath the volcanic system and are responsible for induced microseismicity. Our study shows that implementing dv/v quantification as an additional tool for geothermal seismic monitoring can lead to an improved understanding of the behavior of the geothermal system over time and ultimately de-risk geothermal operations.

Fluids control along-strike variations in the Alaska megathrust slip

Interseismic and coseismic slip on the subduction zone megathrust show complex along-strike variations and the controlling factors remain debated. Here we image seismic velocity anomalies to infer fluid distribution in the Alaska subduction zone (Alaska Peninsula section) with seismic tomography. The weakly locked Shumagin segment is characterized by abundant fluids on the plate interface and in the overriding plate, whereas the moderately-to-highly locked Chignik and Chirikof segments that hosted M > 8 earthquakes are relatively dry. The Kodiak segment, which was ruptured by the 1964 M9.2 Great Alaska earthquake and is presently fully locked near the trench, is characterized by a fluid-rich plate interface but a fluid-poor and inferred low-porosity overriding plate. Multiple large earthquakes have nucleated in fluid-rich regions at the plate interface. Our study highlights the important roles of fluids, particularly in the overriding plate, in controlling interseismic slip deficit and earthquake rupture. We propose that a fluid-rich overriding plate means that there has been a high sustained flux of fluids across the plate interface, which enhances the formation of clay minerals that promote fault creep.

Crustal shear wave velocity and radial anisotropy beneath the Mississippi embayment from ambient noise tomography

Imaging of shallow multiple basins and middle to lower crustal structure beneath the Mississippi embayment is still limited. No tomography studies have delineated multiple basins as effectively as those from reflection and drilling experiments. Several studies have constructed different velocity models for the middle to lower crust. We take advantage of the Northern Embayment Lithosphere Experiment dataset and long-term deployment of 277 3-component broadband stations installed from 1990 to 2018 to image crustal shear wave velocity and radial anisotropy in the Mississippi embayment. We first construct 5–40s group velocity maps for Rayleigh wave and 5–30s ones for Love wave. In the 5 s group velocity map, we observe that low velocity features correlate spatially with the sediment thickness from the Crust 1.0 model. High velocity anomalies are associated with the Ozark uplift and the Nashville Dome. In the 32 s group velocity map, high-low velocity boundaries align with 40 km crustal thickness contours from the Crust 1.0 model. We then invert Rayleigh and Love group velocity dispersion for shear wave velocity and radial anisotropy at each grid point and construct a 3D model accordingly. At the depth of 0–5 km of the model, we find that low velocity features characterize the Mississippi Valley Graben, Rough Creek Graben, Black Warrior Basin, and coastal flood plain. High velocity features mark thinly sedimented regions, like the Ozark uplift and the Nashville Dome. Regions with high gravity anomalies are characterized with negative anisotropies. At the depth of 5–12 km, we observe low velocity features in the Missouri Batholith, centering around 10 km. This region is characterized with negative radial anisotropy (Vsh < Vsv), inferring possible intrusions. High-low velocity boundaries correlate with the Appalachian-Ouachita thrust front. At the depth of 12–24 km, the Mississippi Valley Graben is filled with high velocity material, which may be related to the formation of the graben due to isostatic equilibrium. At the depth of 24–43 km, high velocity anomalies characterize the lower crust of the Mississippi embayment, consistent with high density features from the inversion of free air anomaly. The high velocity anomaly might be responsible for the formation of the embayment due to isostatic adjustment. The average velocity in the lower crust is approximately 4.3 km/s. Mafic Pikwitonei granulite is a possible material, indicating relative strong lower crust. This work provides useful information to explain earthquake occurrence in the New Madrid Seismic Zone and regional tectonic settings.

Computation and Modelling of Seismic Refraction Tomography for Anisotropic Medium Based on Pseudo Bending Method

Seismic refraction tomography is one of the imaging techniques in geophysical methods used to remodel the near-surface velocity layer structure of the Earth. In this study, we carried out a new computational approach and modelling of seismic refraction tomography using the pseudo-bending method. The true model of the near-surface is designed to be anisotropic medium which is having a low velocity anomaly distribution. This anomaly is constructed in such a way as to be similar to a model of liquid waste away spreading, which exhibits seismic velocities ranging from 1600 m/s to 1800 m/s. Based on our computations and numerical modelling results, it was found that the ray tracing path using pseudo-bending method displays an asymmetrical trajectory when the positions of the source and geophone are exchanged. Altering the shooting configuration from direct shoot (DS) to reversed shoot (RS) also reveals a significant difference in travel time values. The results of delay time tomography inversion, which represents the difference between travel times in the true model and the initial model using the SIRT method, indicate the presence of a low velocity anomaly that can be interpreted as the distribution of liquid waste.

Uppermost Mantle Pn Velocity and Anisotropy Structures Beneath the Sakhalin–Kuril–Kamchatka Region

AbstractIn this study, we used the Pn tomography method to obtain detailed velocity and anisotropy structures of the uppermost mantle beneath Sakhalin–Kuril–Kamchatka region for improving the understanding of plate subduction, arc–arc collision, and volcanism. We found low Pn velocities beneath volcanoes and areas characterized by pronounced tectonic activity and high Pn velocities with strong anisotropy in the subducting plate. Low Pn velocity anomalies beneath southern Sakhalin connected the low‐velocity anomalies in the mantle and crust, indicating the ascent of fluid or melt, and may provide a magmatic source for the rear‐arc Rishiri volcano. In the absence of plate subduction, low‐velocity anomalies and north–south Pn anisotropy manifested beneath northern Kamchatka, revealing the lateral propagation of mantle flow beneath this northern region. We suggest that the eastern boundary of the slab window at the Kamchatka–Aleutian junction is likely to be located near the Komandorsky Islands.

Processes behind the decrease in Congo Basin precipitation during the rainy seasons inferred from ERA‐5 reanalysis

AbstractThe ongoing Congo Basin drying is linked to decreasing atmospheric water vapour. To comprehend the processes reducing precipitation, we analysed water and Moist Static Energy components, associated with heating sources using ERA5 reanalysis data over the period 1981–2022. The findings reveal that the precipitation deficit in the Congo Basin during March–April–May and September–October–November is primarily attributed to a decreasing trend in the vertical moisture advection anomaly induced by vertical velocity anomalies. Following this, a horizontal moisture advection anomaly induced by specific humidity anomalies plays a significant role. However, evaporation, horizontal moisture advection anomalies induced by wind anomalies, and vertical moisture advection anomalies induced by specific humidity anomalies do not appear to contribute significantly to the drying trend in the region. The horizontal propagation of MSE advection anomalies exhibits a similar pattern to that of moisture flux advection anomalies, both at different boundaries and within the basin itself. Conversely, the contribution of vertical advection of MSE anomalies remains weak. By integrating moisture and MSE budgets, we interpret the declining precipitation trend as a consequence of uplift driven by horizontal MSE advection. Furthermore, our analysis indicates a decreasing trend in the atmospheric heating source within the basin in association with reduced precipitation. This reduction could signify a decrease in moisture influx into the basin, thereby providing more explanations for the observed drying.

High P‐T Sound Velocities of Amphiboles: Implications for Low‐Velocity Anomalies in Metasomatized Upper Mantle

AbstractMetasomatized mantle xenoliths containing hydrous minerals, such as amphiboles, serpentine, and phlogopite, likely represent the potential mineralogical compositions of the metasomatized upper mantle, where low seismic velocities are commonly observed. This study presents the first experimentally determined single‐crystal elasticity model of an Fe‐free near Ca, Mg‐endmember amphibole tremolite at high pressure and/or temperature conditions (maximum pressure 7.3(1) GPa, maximum temperature 700 K) using Brillouin spectroscopy. We found that sound velocities of amphiboles strongly depend on the Fe content. We then calculated the sound velocities of 441 hydrous‐mineral‐bearing mantle xenoliths collected around the globe, and quantitatively evaluated the roles that amphiboles, phlogopite and serpentine played in producing the low velocity anomalies in the metasomatized upper mantle.

The Magmatic Patterns Formed by the Interaction of the Hainan Mantle Plume and Lei–Qiong Crust Revealed through Seismic Ambient Noise Imaging

Magmatism on continental lithospheres induced by mantle plumes is more complex compared to oceanic intraplate volcanism owing to the heterogeneous nature of continental crustal and lithospheric structures. Substantial evidence points to the deep-oriented Hainan mantle plume beneath the Lei–Qiong region, the southernmost of the South China block. In this study, we present a detailed shear wave velocity model of the crust and uppermost mantle in the Lei–Qiong volcanic region, derived from 3-year seismic data (2016–2018) from 34 stations and the use of the ambient noise tomography method. An evident columnar low-velocity anomaly was imaged in the crust and uppermost mantle beneath the Wushi Sag (WSS), Beibu Gulf, potentially suggesting that the center of either one branch or the entirety of the Hainan mantle plume impacts the crust here. This low-velocity anomaly is overlaid by a local Moho deepening, indicative of underplating beneath the existing crust. The Maanling–Leihuling Volcanic Field (MLVF) in northern Hainan Island, previously considered the center of the hotspot, does not exhibit such distinct velocity anomalies. Instead, subtle lower crustal anomalies beneath the MLVF are linked with the upper mantle low-velocity zone beneath the WSS. Additionally, the high-conductivity bodies beneath the MLVF indicate lateral magma transport. Earthquake swarms and deep-seated seismic events beneath the WSS further support the presence of magmatic processes. This study indicates that in the Lei–Qiong region, the interaction of the continental crust with the mantle plume centered in the WSS results in magma exhibiting both vertical ascent and lateral migration, leading to a dual low-velocity shear wave pattern in the upper crust, which significantly influences the surface volcanic activity.

Subsurface anatomy of the Irazú–Turrialba volcanic complex, inferred from the integration of local and ambient seismic tomographic methods

SUMMARY Irazú and Turrialba are a twin volcanic complex that marks a distinct stop in volcanism along the Central America volcanic arc. We present a new traveltime velocity model of the crust beneath Irazú and Turrialba volcanoes, Costa Rica, and interpret it considering the results of previous ambient noise tomographic inversions. Data were acquired by a temporary seismic network during a period of low activity of the Irazú–Turrialba volcanic complex in 2018–2019. Beneath the Irazú volcano, we observe low P-wave velocities (VP = 5 km s−1) and low velocity ratios (VP/VS = 1.6). In contrast, below the Turrialba volcano, we observe low S-wave velocities (VS = 3 km s−1) and a high VP/VS (= 1.85) anomaly. We found that locations of low VP and VS anomalies (−15 %) correspond well with shear wave velocity anomalies retrieved from ambient noise tomography. At shallower depths, we observe high VP and VS anomalies (+15 %) located between the summits of the volcanoes. Subvertical velocity anomalies are also observed at greater depths, with high VP and VS anomalies appearing at the lower limits of our models. We propose a complex structure of an intermediate magmatic reservoir, presenting multiphase fluid states of a liquid-to-gas transition beneath Irazú and a juvenile store of magmatic fluid beneath Turrialba, while shallow fluid transport provides evidence of magmatic–hydrothermal interactions.

Topography effect on seismic waveform tomography: a quantitative study

SUMMARY In seismic tomography practices the Earth's surface is sometimes assumed to be either spherical or flat for convenience in forward modelling calculations. The effect of irregular surface topography on seismic wave propagation is thus ignored, resulting in biases in the phases and amplitudes of synthetic seismograms, which contribute to the residuals that are mapped into velocity structures in tomography inversions. In this study, we conduct a series of inversion experiments based on the adjoint waveform tomography method to quantitatively assess the topography effect on waveform-based inversion results. We first employ models with simplified topography to better highlight and quantify the topography effect. Results show that when topography effect is ignored in the forward modelling, it is mapped into velocity perturbations, leading to spurious velocity anomalies in tomography models. The strength of the spurious velocity anomaly is quasi-linearly related to locally averaged topography gradient. Our inversion experiments demonstrate that in places of strong topography variation, such as the Longmenshan Fault Zone region where the 2008 Mw 7.9 Wenchuan earthquake occurred, topography effect can lead to spurious relative velocity anomalies of up to 10 per cent, which cannot be ignored in waveform-based tomography inversions.

The structure of the mantle transition zone and its implications for mantle upwelling beneath the central Tibetan plateau from triplicated P-waveforms modeling

The Tibetan plateau is considered as an ideal place for studying the continental collisions, tectonic movements and regional geodynamics. Although multiple episodes of igneous activities in Tibetan plateau have been suggested derived from mantle upwelling, their origins and corresponding geodynamic models remain vigorously debate. This study presents a high-resolution upper mantle velocity model of P-wave and its lateral variation beneath the central Tibetan plateau based on triplicated waveform modeling. The results show that, compared with IASP91 model, the 660-km discontinuity is deepened approximately 17–28 km from the Qiangtang terrane to the Beishan block. The lower mantle transition zone (MTZ) is characterized by high P-wave velocity anomalies (HVA) in range of 0.3%–0.9%. We infer that they are mainly related to the hydrous materials with low temperature in the lower MTZ. Besides, the depression of the 410-km discontinuity is about 6–15 km. A remarkable low velocity layer (LVL) is detected atop the 410-km discontinuity, exhibiting the thickness of 31–40 km and low P-wave velocity anomaly of −2.5% to −1.6%. This is against completely temperature-dominated origin of the LVL, and suggests that the material composition has also contributes. Additionally, we propose a new geodynamic model to demonstrate the asthenospheric upwelling beneath the Qiangtang and Songpan-Ganzi terranes, which suggests that the passively ascending water-enriched materials from the lower MTZ and their dehydration are responsible for the observed LVL atop the 410-km discontinuity, as well as the origin of the Cenozoic magmatic intrusion and volcanic eruption of the northern Tibet volcano in Pleistocene.

Ambient noise tomography of El Hierro island (Canary Islands)

El Hierro island is one of the most active islands in the Canary Islands from a volcanological point of view. This is the reason why the imaging of the internal crustal structure is of huge importance. The geophysical exploration methods employed on El Hierro Island, such as gravimetry and seismic tomography, allowed obtaining the high-resolution characterization of the crust’s deep part. However, these methods did not yield significant information about the surface and the shallower part of the crust. To gain a deeper insight into the shallow geological structure of El Hierro island, we employed Ambient Noise Tomography to construct a 3D S-wave velocity model. Our investigation revealed the presence of seven significant seismic velocity anomalies, partly identified by previous studies. We identified two high-velocity anomalies located in the eastern and western parts of the island at a depth between 0 and 3 km below sea level (b.s.l.). We interpreted these anomalies as dense intrusive complexes of dikes, possibly linked to the Tanganasoga volcano and the formation of the Tiñor edifice. Additionally, we observed two high-velocity anomalies in the northern and southern parts of the island at a depth between 3 and 4 km b.s.l., which we related to the accumulation of solidified igneous rocks. On the other hand, a low-velocity anomaly was observed in the Golfo valley, between 0 and 0.5 km b.s.l., and we interpreted it as megalandslide deposits. This anomaly was evidenced for the first time in the present study. Finally, two low-velocity anomalies were observed in the southern part of the island at different depths, between 0–0.5 km b.s.l. and 0–2 km b.s.l. These were interpreted as fractures generated by Quaternary volcanism along the SSE Rift. Also, one of them was evidenced for the first time in this study, corresponding to the zone of the fractures produced during the Quaternary volcanism. This study has allowed us to gain a more detailed understanding of the shallow geological structure of the island. Even if most of the anomalies had been evidenced previously, we could observe the existence of two low-velocity zones in the shallow crust that have not been observed before.

Three-dimensional fine crust-mantle structure imaging and structural characteristics of Bohai Sea

This study used the consistency-constrained double-difference tomography method to invert 3D fine structure models of Vp, Vs, Vp/Vs at depths above 60 km, and precise relocation parameters of earthquakes in the Bohai Sea. According to the results, the velocity structures of P-wave and S-wave in the Bohai Sea area were highly similar and demonstrated noticeable lateral non-uniformity. The crust beneath the Tanlu Fault displayed a clear stratification structure, with a continuous velocity transition in the middle of the crust. The intricate crustal structure beneath the Zhangpeng Fault displayed high-velocity bodies in the crust and low-velocity anomalous zones connected to the top of the mantle in the lower section of the crust. The structural pattern in the deep crust of the Bohai Sea controls the occurrence characteristics of the Zhangpeng Fault and the Tanlu Fault. The earthquakes in the Bohai Sea area were concentrated mainly in the southern part near the Zhangpeng Fault. There is a good correspondence between the relocated earthquakes and velocity structure. There are many significant differences in crustal structure between the north and south of the Bohai Strait, and there are obvious velocity anomalies in the middle and upper crust. The distribution shape of high Vp/Vs value indicates that mantle material migration has occurred at the bottom of the crust. This paper provided important reference for further research on the relationship between deep tectonic features and tectonic activity in the North China Craton.

An updated 3D crustal velocity model in Cyprus by non-linear traveltime tomography

We have used >6000 regional earthquakes in the broader Cyprus island area, with 100,000+ P and S-phases recorded by 26 local seismological stations to obtain a 3D velocity model for the Cyprus area using a non-linear seismic tomography approach. Due to the limited aperture of the recording network, the tomographic model mainly recovers the dominant geophysical features of the upper crust of the Cyprus island. The obtained results reveal new information on the 3D structure of the Troodos ophiolithic complex, identifying a high velocity anomaly that coincides with the NW-SE trending ophiolitic units of the Troodos mountain. Furthermore, they clearly show that the anomaly associated with the ophiolite zone of Troodos is well preserved up to a depth of at least 10 Km, dipping towards the Northeast. The Mesaoria basin in central Cyprus is recognized as a large low-velocity feature, with the results suggesting a significant thickness of the basin, exceeding 5 km. A similar low-velocity pattern can be observed for the Circum Troodos formations (Pakhna and Lefkara) to the southwest of the Troodos mountain, while the Kyrenia range is not associated with a crustal velocity anomaly. To assess the reliability of the results we have conducted checkerboard resolution tests, which suggest that the uppermost crustal features are well reconstructed, while for larger depths (>10 km) we cannot recover mid-crustal anomalies due to the insufficient traveltime data coverage and quality (e.g., traveltime picking errors).

Impacts of compound extreme weather events on summer ozone in the Beijing-Tianjin-Hebei region

The Beijing-Tianjin-Hebei (BTH) region faces severe summer ozone (O3) pollution. In this study, we conducted a comprehensive analysis of O3 concentration variations during compound extreme (CE) weather events that occur simultaneously with atmospheric stagnation and heat wave during the summer period from 2015 to 2021. Our results showed that during CE events, the rate of O3 exceedance was nearly three times higher compared to non-CE days. O3 concentrations remained relatively stable during CE events lasting five days or less in regions other than Beijing, but when CE events extended beyond five days, both O3 concentration and the exceedance rate increased rapidly. Under CE conditions, daily O3 pollution persisted for a longer duration with 70.0% of daily O3 peaks surpassing 200 μg m−3, especially between 7 a.m. and 12 p.m., the time when CE events had a significant impact on O3 concentration characterized by a rapid rise in O3 concentration. The movement of the western Pacific subtropical high towards the east and its northward shift may contribute to the occurrence of favorable conditions for CE events in the BTH region, but meteorologically these results can be attributed to anticyclonic and high-pressure system anomalies during CE events in the BTH region, which are characterized by low total cloud cover, high temperature, low relative humidity, and high radiation. Positive anomalies in vertical velocity further contribute to favorable conditions for photochemical O3 formation, accumulation, and persistent O3 exceedances.

Robust estimates of the ratio between S- and P-wave velocity anomalies in the Earth's mantle using normal modes

Seismic tomography allows us to image the interior of the Earth. In general, to determine the nature of seismic anomalies, constraints on more than one seismic parameter are required. For example, the ratio R between perturbations in vs and vp (dlnvs and dlnvp, respectively) is studied extensively in the lowermost mantle and interpreted in terms of thermal and/or chemical anomalies. However, to jointly interpret tomographic models of variations in vs and vp or their ratio R, it is essential for them to share the same local resolution. Most existing models do not provide resolution information, and thus cannot guarantee to honour this condition. In addition, uncertainties are typically not provided, making it difficult to robustly interpret the ratio R=dlnvs/dlnvp. To overcome these issues, we utilise the recently developed SOLA tomographic method, a variant of the linear Backus–Gilbert inversion scheme. SOLA retrieves local-average model estimates, together with information on their uncertainties, whilst it also provides direct control on model resolution through target kernels. In this contribution, we apply SOLA to normal-mode data with sensitivity to both vs and vp, as well as density throughout the mantle. Specifically, we aim to develop models of both vs and vp with the same local resolution. We test our methodology and approach using synthetic tests for various noise cases (random noise, data noise or also additional 3D Earth noise due to variations in other physical parameters than the one of interest). We find that the addition of the 3D noise increases the uncertainties in our model estimates significantly, only allowing us to find model estimates in six or four layers for vs and vp, respectively. While the synthetic tests indicate that no satisfactory density models can be obtained, we easily manage to construct models of dlnvs and dlnvp with almost identical resolution, from which the ratio R can be robustly inferred. The obtained values of R in our synthetic experiments significantly depend on the noise case considered and the method used to calculate it, with the addition of 3D noise always leading to an overestimate of R. When applying our approach to real data, we obtain values of R in the range of 2.5–4.0 in the lowest 600 km of the mantle, which are consistent with previous studies. Our model estimates with related resolving kernels and uncertainties can be used to test geodynamic model predictions to provide further insights into the temperature and composition of the mantle.

Spatial Correlations between the Terrain Features, Gravitational Field and Seismic Velocity Anomalies in the Central Kamchatka Region

Abstract ––Different approaches to tectonic zoning and different types (complexes) of used data, on the basis of which the zoning is performed, are the plurality of tectonic schemes of Kamchatka and the Kamchatka region previously obtained by different authors. In this work, a joint analysis of gravity field anomalies, gravity field and relief transforms, magnetic field and modern seismotomographic models of velocity anomalies was used to study spatial relationships of tectonic manifestations of geological processes. As a result of this study, a linear-block type tectonic scheme was constructed. The obtained linear structural elements are in good agreement with the data of potential fields, seismic constructions and previously published materials. The productivity of the presented method of complex analysis of geophysical fields for revealing their spatial relations and block structures of the environment is shown. The analysis outcomes can be used to establish spatial constraints for gravity field sources when solving inverse problems, potential field separation problems and their detailed geological interpretation.

Imaging high-resolution 3D shallow crustal structure by ambient noise method in the Shiziyang Waterway area, Guangzhou, South China

The Shiziyang Waterway, located in the northern continental margin of the South China Sea, has developed a complex fault system owing to the coupling action of the oceanic crust and continental crust. However, the geological structure in the shallow crust beneath the Shiziyang Waterway and the relationship between faults and waterway development are still unclear, resulting in difficulty in assessing the earthquake risk to the surrounding cities. In this study, based on a high-density short-period seismic array that was recently deployed around the Shiziyang Waterway, a high-resolution shear wave velocity model is inverted to as shallow as 4 km by the ambient noise tomography method. The results showed that at depths ranging from 0.2 to 4 km, three sets of low-velocity zones (LVZs), trending east–west, northeast, and northwest, are revealed, which are speculated to be the Shougouling fault (SGLF), Xinhui-Shiqiao fault (XSF), and Shiziyang fault bundle (SZYFb), respectively. The shear-wave velocity anomalies illustrate that the SGLF stopped extending east at the north endpoint of the Wenchong-Zhujiangkou fault (WZF). The two-stage differential vertical movements of the WZF, Nangang-Humen fault (NHF), and Hualong-Haioudao fault (HHF) dominated the opening of the Shiziyang Waterway. This study provides new seismological evidence of the opening process for the Shiziyang Waterway and reveals that internal fault activity dominated the development and evolution of the Pearl River Delta.

Near-Surface Radial Anisotropy Tomography of Geothermal Reservoir Using Dense Seismic Nodal Array

The development and utilization of geothermal energy has become an indispensable part in the process of urbanization in China. Near-surface tomographic imaging based on ambient noise is increasingly widely used in geothermal exploration. To construct the subsurface structure of Shear-wave velocity in Andi town, Zhejiang Province, recordings of ambient noise data of 4.7 days were recorded by a dense seismic array composed of 192 stations with an aperture of about 5 km. Analysis of array beamforming in high-frequency bands indicated the approximately uniformly distributed of noise sources. We applied seismic interferometry to retrieve both Rayleigh and Love waves. We used two-station analysis to obtain dispersion images and the corresponding high-quality dispersion curves were extracted by automated codes. 3-D tomographic results down to a depth of 1.6 km of S-wave velocity structure and radial anisotropy parameters were obtained using direct surface wave radial anisotropy tomography method. The anomalies of S-wave velocity and radial anisotropy parameters were combined to infer the storage of geothermal reservoir in Andi town. Furthermore, the accuracy of our results was verified by the borehole data.

Low‐Velocity Structure of Subducted Oceanic Crust in the Upper Mantle: Insights From High Pressure and Temperature Elasticity Measurements of Aragonite

AbstractThis study examines the effects of water and carbon on the velocity profiles of subducted oceanic crust in the upper mantle. High pressure and temperature Brillouin measurements were conducted to determine the single‐crystal elasticity of aragonite (CaCO3) up to 20 GPa and 600 K, respectively. Using the finite‐strain method, we determined the elastic parameters: KS0 = 70.7 GPa, G0 = 36.1(3) GPa with KS′ = 5.0(1), G′ = 1.3(1), ∂KS/∂T = −0.020(2) GPa/K, and ∂G/∂T = −0.015(1) GPa/K. When combined with literature results, our findings reveal that neither hydrous minerals nor carbonate alone can explain the observed 3%–4% and 3%–7% low velocity anomalies in the compressional (VP) and shear‐wave (VS) velocities of subducted oceanic crust in the circum‐Pacific region at 150–250 km depth. Considering the combined effect of water and carbonate, the addition of 5.9–6.9 vol.% aragonite together with 10 vol.% lawsonite in the oceanic crust can produce a 3%–4% and 6.1%–8.1% low velocity anomalies in the VP and VS at 150–250 km depth, respectively, consistent with the seismic observations in the region. Complete dehydration of lawsonite accompanied by the gradual decarbonization of the subducted oceanic crust explains the absence of low‐velocity anomalies below 300‐km depth. Our findings help for a better understanding of the possible influence of carbon and water on seismic velocities of the mantle. As a result, the circulation of carbon and water may be better understood while taking the complex velocity structure of subduction zones and mineral physics findings into account.