Permanent Stations Research Articles

Directional amplification across the San Jacinto fault zone, CA

SUMMARY The amplitude, frequency and polarization of ground motion at the surface can be affected by the local geology. While low-velocity sediments and fill can amplify ground motions in certain frequency ranges, the low velocities found in fault zones can also produce prominent wavelets. In this paper, we provide further evidence that polarization of ground motion can be affected by the geological fabric in fault zones that have sustained significant brittle deformation. Aside from the well-known effect of fault-trapped waves in the low-velocity zone with polarization azimuths parallel to the fault strike, the effect of stiffness anisotropy was recently recognized with polarization azimuths at high-angle to the fault strike and orthogonal to the locally predominant fracture field in the fault damage zone. To clarify further such features, we investigate directional amplification effects across the San Jacinto fault zone in Southern California using seismic data recorded by permanent seismic stations and dense across-fault arrays. We observe three main polarization trends. The first trend parallel to the fault strike is ascribed to fault-trapped waves along the low-velocity zone, in agreement with several studies in the last decade in the same region. The second and third trends are orthogonal to the orientation of R and T Riedel planes, respectively. They are related to the stiffness anisotropy in densely fractured rocks in the damage zone, which are more compliant orthogonal to their fractures. At some locations the two effects are superimposed, occurring in different and distinct frequency ranges. Directional amplification at rock sites can be important for expected ground motion and seismic hazard. However, in seismic engineering the current prescriptions of seismic codes do not account for amplification effects at rock sites at frequencies of engineering interest.

Влияние смены оборудования пулковской ГНСС-станции PULK на результаты наблюдений

Pulkovo Observatory has been operating a permanent GNSS station PULK since 2002. This paper presents the results of preliminary analysis of the changes in the errors in the coordinates of the station over time, primarily after replacement of the receiver or antenna. Two accuracy estimates were applied: the formal error of the daily coordinate estimates and the noise component of the PULK coordinate series. In the first case, the median estimate was used, in the second, the Allan variance. As a result, it turned out that two receiver replacements in the first years of operation led to a noticeable decrease in the formal coordinate error, whereas subsequent receiver replacements, as well as both antenna replacements, showed no noticeable changes in this error. The magnitude of the noise component of the coordinate series did not change much during the entire observation period.

The KISS Network in 2015–2016: Catalogs and Comparison of the Processing Results with Operational Estimates from the Permanent Network

As part of the international collaboration of several research groups from Russia, France, and Germany, 77 temporary seismic stations were installed in the summer of 2015 for one-year period to conduct a detailed study of the deep structure of the Earth’s crust and upper mantle in the region of the Klyuchevskoi Volcano Group (KGV) in the Kamchatka Peninsula. One of the results of the KISS experiment (Klyuchevskoi Investigation – Seismic Structure of an extraordinary volcanic system) was the final catalog of the joint data from the temporary stations and the permanent network of the Kamchatka Branch of the Geophysical Survey of the Russian Academy of Sciences (KB GS RAS). The catalog comprises 2136 events, including 560 for which the permanent network catalog lacked sufficient data for correct processing. The catalog in .xlsx format and the station bulletin in .isf format are presented in the supplementary material to the paper. A comparative analysis was conducted on the joint solutions of two catalogs: one obtained solely from the data of the KB GS RAS permanent network stations and another from a denser seismic network integrated with KISS stations.

An Evaluation of Optimization Algorithms for the Optimal Selection of GNSS Satellite Subsets

Continuous advancements in GNSS systems have led, apart from the broadly used GPS, to the development of other satellite systems (Galileo, BeiDou, GLONASS), which have significantly increased the number of available satellites for GNSS positioning applications. However, despite GNSS satellites’ redundancy, a potential poor GNSS satellite signal (i.e., low signal-to-noise ratio) can negatively affect the GNSS’s performance and positioning accuracy. On the other hand, selecting high-quality GNSS satellite signals by retaining a sufficient number of GNSS satellites can enhance the GNSS’s positioning performance. Various methods, including optimization algorithms, which are also commonly adopted in artificial intelligence (AI) methods, have been applied for satellite selection. In this study, five optimization algorithms were investigated and assessed in terms of their ability to determine the optimal GNSS satellite constellation, such as Artificial Bee Colony optimization (ABC), Ant Colony Optimization (ACO), Genetic Algorithm (GA), Particle Swarm Optimization (PSO), and Simulated Annealing (SA). The assessment of the optimization algorithms was based on two criteria, such as the robustness of the solution for the optimal satellite constellation and the time required to find the solution. The selection of the GNSS satellites was based on the weighted geometric dilution of precision (WGDOP) parameter, where the geometric dilution of precision (GDOP) is modified by applying weights based on the quality of the satellites’ signal. The optimization algorithms were tested on the basis of 24 h of tracking data gathered from a permanent GNSS station, for GPS-only and multi-GNSS data (GPS, GLONASS, and Galileo). According to the comparison results, the ABC, ACO, and PSO algorithms were equivalent in terms of selection accuracy and speed. However, ABC was determined to be the most suitable algorithm due it requiring the fewest number of parameters to be set. To further investigate ABC’s performance, the method was applied for the selection of an optimal GNSS satellite subset according to the number of total available tracked GNSS satellites (up to 31 satellites), leading to more than 300 million possible combinations of 15 GNSS satellites. ABC was able to select the optimal satellite subsets with 100% accuracy.

Seismicity Migration and the Upper Crustal Structure in the Xinfengjiang Reservoir

Abstract After the impoundment of the Xinfengjiang Reservoir (XFJR) in Guangdong, China, numerous earthquakes occurred in the area, including a magnitude 6.1 event in 1962. Analysis of historical earthquakes indicates that M ≥ 4 earthquakes began occurring in the northwestern XFJR in 2012, and seismicity has gradually migrated from the southeastern to the northwestern reservoir (NWR). However, the mechanisms governing the migration of seismicity and the current upper-crustal structure beneath the reservoir area remain unclear. In our study, we conducted tomographic imaging by combining waveform data from short-period and permanent stations to construct a 3D velocity model. Our high-resolution velocity models revealed a horizontal fractured zone at ∼5 km depth that extends from the southeastern to northwestern XFJR, and a steep fault that extends to about 9 km depth. These two fractured zones may interact with each other, allowing for fluid infiltration and contributing to earthquake triggering via pore pressure diffusion in the XFJR areas. Furthermore, the calculation of Coulomb stress changes indicated that microearthquakes in the southeastern XFJR may contribute to the seismicity in the NWR. However, the influences of M ≥ 4 earthquakes in the northwestern XFJR on subsequent M ≥ 4 earthquakes in the southeastern XFJR vary differently. Our results provide crucial insights for understanding the migration of microearthquakes in the XFJR area.

Distribution, source identification, and contamination level of trace metals in the sediment of the Shari-Goyain River in Bangladesh: Implications for ecological health risks

The major issue associated with coal mine drainage is its contribution to river pollution, which occurs at the local, regional, and global levels. The pollution of surface sediment by heavy metals is a major environmental and health concern in coal mining and downstream areas. This study explored the concentrations of eight metal components in the sediments from the Shari-Goyain River of Bangladesh. The trend of decreasing metal concentration was identified as Fe > Mn > Ni > Zn > Cu > Cr > Pb > Cd in sediment. The level of metal pollution in the study area was assessed by using various indicators like geo-accumulation index (Igeo), pollution load index (PLI), contamination factor (CF), potential ecological risk factor, and risk index (RI). The sediment in the river exhibited PLI values ranging from 0.133 to 0.543, suggesting a low level of pollution from the evaluated heavy metals (PLI < 1). The mean Igeo showed that the study area was unpolluted whereas Ni showed unpolluted to moderately polluted status. For most of the metals, the sediment samples recorded a low degree of contamination (CF < 1) except Pb, Cd, and Ni which exhibited moderate degree of contamination (1 ≤ CF < 3). Through the application of various statistical analyses, coal mine drainage has been identified as the possible source of pollution of the analyzed metals in the Shari-Goyain River. However, the risk index and RI suggested a low risk of metal pollution in the studied areas. To improve the environmental conditions of the Shari-Goyain River, it is crucial to construct permanent sediment quality monitoring stations and conduct extensive ecological investigations.

Analysis of Multipath Changes in the Polish Permanent GNSS Stations Network

This study examines the influence of multipath errors on Global Navigation Satellite System (GNSS) measurements collected at ASG-EUPOS reference stations between 2010 and 2021. Multipath occurs when GNSS signals reflect off surrounding objects before reaching the receiver antenna, leading to positioning errors. In the case of reference stations, all available mitigation techniques were used to minimize the impact of multipath. However, it is still detectable and affects the measurement results. For carrier phase differential positioning, it increases the ambiguous search space, which results in a decrease in determining rover—reference station vector accuracy. The study employs two linear combinations (Code-Minus-Carrier and Multipath Pseudorange Observable) to quantify the multipath effect on both pseudorange and carrier phase measurements. Based on the research, it was found that the multipath values changed depending on the change of the receiver and the terrain around the reference stations. The study observed a gradual decrease in multipath errors from 2010 to 2021, likely due to technological advancements in receiver design. No significant increase in multipath errors was observed due to environmental changes around the stations, suggesting a minimal influence from new reflecting objects nearby. Based on the analyses conducted, it is also recommended to perform periodic tests to detect incorrect receiver configuration or operation.

GNSS and Sentinel-1 InSAR Integrated Long-Term Subsidence Monitoring in Quetta and Mastung Districts, Balochistan, Pakistan

Land subsidence (LS) is a global phenomenon that has affected several urban centres around the world such as Jakarta (Indonesia), Mexico City (Mexico), Xi’an (China), and Iron County (US). It has mainly been attributed to anthropogenic activities such as groundwater exploitation, especially in unconsolidated aquifer systems rich in highly compressible clay and silt. The platy clay minerals rearrange into horizontal stacks after dewatering, leading to a volume change due to overburden. In this study, land subsidence is investigated in the Quetta and Mastung districts, Balochistan, Pakistan, by employing Small Baseline Subset (SBAS) Interferometric Synthetic Aperture Radar (InSAR), Global Navigation Satellite System (GNSS), and groundwater level (GWL) variations. This study represents the first attempt in Pakistan to measure the long-term land subsidence by fusing GNSS and InSAR data for improved validity. InSAR data from the Sentinel-1 satellite in the Ascending (195 scenes) and Descending (183 scenes) tracks were used to analyse LS from December 2015 to December 2022. High-accuracy Trimble NetRS GNSS receivers were used in five locations from October 2006 to December 2022. An average subsidence ranging from 3.2 cm/y to 16 cm/y was recorded in the valley mainly due to the GWL decline and clay-rich sediments, which are prone to compaction due to dewatering. An accumulative LS of 2 m was recorded by the permanent GNSS station in central Quetta from October 2008 to January 2023 (14.2 years). An acceleration in the subsidence from 12 cm/y to 16.6 cm/y after 2016 was recorded by the continuous GNSS. Additionally, the InSAR and GNSS values were compared for validation, resulting in a good correlation between both techniques. A GWL decline ranging from 1.7 m to 6 m was recorded by the piezometers in Quetta during the period 1987–2022. Large- and small-scale fissures were observed in the study area during the surveys. These fissures are responsible for damage to the city’s infrastructure and aquifer contamination. The subsidence profile also agrees with the subsurface lithology. Our assessment concludes that Quetta may be the fastest-sinking metropolitan city in Pakistan. The overexploitation of groundwater and the population explosion may be the main contributing factors for the land subsidence.

Identifying and localizing of seismogenic electromagnetic anomalies from data observed by permanent MT stations

Identifying and localizing of seismogenic electromagnetic anomalies from data observed by permanent MT stations

Crustal structure beneath Central Kamchatka inferred from ambient noise tomography

In this study, we investigated the crustal structure beneath the area of Central Kamchatka (Kamchatka Peninsula, Russian Far East) limited by latitudes 52.8 and 54.1 degrees. This area comprises parts of the active Eastern Volcanic Front and two mountain ridges Sredinny and Ganalsky, as well as several backarc volcanoes. We used the continuous data from the permanent stations and the temporary network installed in 2019–2020 to perform ambient noise tomography. We implemented a two-step tomographic inversion algorithm that initially generated the two-dimensional maps of group velocities for periods ranging from 1 s to 30 s, then found the optimal 1D model of shear wave velocities and derived the three-dimensional S-wave velocity distribution. In the resulting model, we observe a prominent low-velocity anomaly below the active volcanoes of Avacha and Koryaksky, which is traced to depths of ∼40 km and may be associated with a zone of heated rocks suitable for petrothermal energy extraction. We propose that this anomaly represents a common magma pathway through the crust feeding both volcanoes. Beneath another active Zhupanovsky volcano, we observe three low-velocity layers at depths of 2 km, 8 km and 20 km indicative of a multilevel magma feeding system. We believe our findings will further the understanding of the complex processes occurring within these active volcanoes, which pose a tangible hazard to the relatively populated Petropavlovsk-Kamchatsky urban district. Beneath the Sredinny Ridge we reveal a high-velocity pattern at shallow depths, likely associated with consolidated magmatic structures. The low-velocity in the lower crust is interpreted as a high-temperature zone conserved here since the time when this area hosted an active volcanic arc. Upon further exploration, this high-temperature zone could prove to be a prospective source of geothermal energy. Below the Ganalsky Ridge, a prominent high-velocity anomaly extending to a depth of ∼15 km can be attributed to the Precambrian to Carboniferous metamorphic rocks of predominantly mafic composition.

The 2020 Mw 6.5 Stanley, Idaho, Earthquake and Aftershock Sequence: Complex Faulting at the Northern End of the Basin and Range Province

ABSTRACT The magnitude 6.5 Stanley, Idaho, earthquake occurred on 31 March 2020 in a sparsely populated region north of the Sawtooth normal fault. We used seismic data from temporary and permanent stations to derive a 1D velocity model and relocate 1401 M ≥ 2.4 earthquakes with hypoDD, including a foreshock, the mainshock, and 3 yr of aftershocks. We used broadband data to determine seismic moment tensors for 173 Mw≥3.1 earthquakes. Combining locations and mechanisms shows the mainshock ruptured an unmapped north-trending, steeply west-dipping, left-lateral strike-slip fault. Rupture initiated near the bottom of the seismogenic zone and propagated upward and bilaterally for ∼20 km north and ∼3–5 km south, where the fault likely ends and deformation changes to extension. There, the rupture may have jumped west to another unmapped blind fault accommodating oblique extension. Support for a late mainshock rupture on a northwest-trending, likely east-dipping fault comes from normal faulting aftershocks. The total rupture length is ∼25–30 km, because oblique fault activity ends at the latitude of the northern terminus of the Sawtooth fault, but its trace, if it reached the surface, would be ∼6 km to the west. The Sawtooth fault was not active, even though aftershock clusters indicate that short strike-slip and normal faults are abundant in its footwall and hanging wall. Extension in the northern Basin and Range seems to terminate where major normal faults reach the inactive Eocene Trans-Challis fault system (TCFS), suggesting the TCFS exerts structural control. The deformation north of the TCFS is low, and the strike-slip character was unknown before the Stanley earthquake. Faults rupturing in the Stanley earthquake lack surface expression and are immature with low cumulative displacement. Complex transitions between tectonic regimes are common and may result in blind ruptures on unknown, immature faults, posing an underrated hazard.

Comparison of anomalies in the VLF spectrum of the natural electromagnetic field with data from the seismometer in a landslide-affected area

The data from a permanent monitoring station, based on a sound card as an AD/DA converter and a magnetic loop antenna for continuous recording of electromagnetic field intensities in the Very Low Frequency (VLF) range, were recorded and compared with data from a seismic station based on Raspberry Pi located in a landslide-affected area. The stations operated 24/7 and were placed in a room below ground level. Correlations were observed between seismic signals in the frequency spectrum from 0.1 Hz to 10 Hz (some extending beyond 20 Hz) and anomalies in the electromagnetic field in the VLF range. These anomalies are likely associated with micro-fracturing, piezo-electromagnetic, and triboelectric phenomena within the landslide body, producing relatively weak VLF emissions and a faint seismic signal. A single-component 4.5 Hz 395 Ohm vertical Racotech RGI-20DX geophone with electronic extension to lower frequencies (&lt; 1 Hz) was used with a sampling rate of 100 sps (samples per second). VLF emissions cover almost the entire spectrum from 6 kHz with a peak at 14 kHz to 18 kHz. The received spectrum was divided into following sections: VLF band 4 kHz – 6 kHz; 6 kHz – 8 kHz; 8 kHz – 10 kHz; 10 kHz – 12 kHz; 12 kHz – 14 kHz; 14 kHz – 16 kHz; 16 kHz – 18 kHz. Simultaneously, there were changes in these sections analysed and compared with the seismic record within the same time interval.

Altimeter Calibrations in the Preliminary Four Years’ Operation of Wanshan Calibration Site

In order to accomplish the calibration and validation (Cal/Val) of altimeters, the Wanshan calibration site (WSCS) has been used as a calibration site for satellite altimeters since its completion in August 2019. In this paper, we introduced the WSCS and the dedicated equipment including permanent GNSS reference stations (PGSs), acoustic tide gauges (ATGs), and dedicated GNSS buoys (DGB), etc. placed on Zhi’wan, Wai’ling’ding, Dan’gan, and Miao’Wan islands of the WSCS. The PGSs data of Zhi’wan and Wai’ling’ding islands were processed and analyzed using the GAMIT/GLOBK (Version 10.7) and Hector (Version 1.9) software to define the datum for Cal/Val of altimeters in WSCS. The DGB was used to transfer the datum from the PGSs to the ATGs of Zhi’wan, Wai’ling’ding, and Dan’gan islands. Separately, the tidal and mean sea surface (MSS) corrections are needed in the Cal/Val of altimeters. We evaluated the global/regional tide models of FES2014, HAMTIDE12, DTU16, NAO99jb, GOT4.10, and EOT20 using the three in situ tide gauge data of WSCS and Hong Kong tide gauge data (No. B329) derived from the Global Sea Level Observing System. The HAMTIDE12 tide model was chosen to be the most accurate one to maintain the tidal difference between the locations of the ATGs and the altimeter footprints. To establish the sea surface connections between the ATGs and the altimeter footprints, a GPS towing body and a highly accurate ship-based SSH measurement system (HASMS) were used to measure the sea surface of this area in 2018 and 2022, respectively. The global/regional mean sea surface (MSS) models of DTU 2021, EGM 2008 (mean dynamic topography minus by CLS_MDT_2018), and CLS2015 were accurately evaluated using the in situ measured data and HY-2A altimeter, and the CLS2015 MSS model was used for Cal/Val of altimeters in WSCS. The data collected by the equipment of WSCS, related auxiliary models mentioned above, and the sea level data of the hydrological station placed on Dan’gan island were used to accomplish the Cal/Val of HY-2B, HY-2C, Jason-3, and Sentinel-3A (S3A) altimeters. The bias of HY-2B (Pass No. 375) was −16.7 ± 45.2 mm, with a drift of 0.5 mm/year. The HY-2C biases were −18.9 ± 48.0 mm with drifts of 0.0 mm/year and −5.6 ± 49.3 mm with −0.3 mm/year drifts for Pass No. 170 and 185, respectively. The Jason-3 bias was −4.1 ± 78.7 mm for Pass No. 153 and −25.8 ± 85.5 mm for Pass No. 012 after it has changed its orbits since April 2022, respectively. The biases of S3A were determined to be −16.5 ± 46.3 mm with a drift of −0.6 mm/year and −9.8 ± 30.1 mm with a drift of 0.5 mm/year for Pass No. 260 and 309, respectively. The calibration results show that the WSCS can commercialize the satellite altimeter calibration. We also discussed the calibration potential for a wide swath satellite altimeter of WSCS.

Post-collisional lithospheric delamination in eastern Iran, revealed by non-linear teleseismic tomography and residual topography

The Eastern Iranian Mountain Ranges (EIR) emerged as a consequence of the Late Cretaceous collision between the Afghan and Lut blocks. However, the response of the uppermost mantle to this collision remains enigmatic. Additionally, although petrological evidence suggests that post-collisional delamination is possible, it has not been conclusively identified in prior regional seismic imagery. This observation leads us to further explore this possibility using a dense seismic network. To gain insight into the geodynamic implications for eastern Iran and address knowledge gaps, we extensively investigated the seismic structure of the uppermost mantle beneath the EIR using a dense seismic network of 34 temporary stations, complemented by data from nine additional local permanent stations. By meticulously analyzing 6589 relative arrival time residuals from teleseismic records with favorable signal-to-noise ratios, we applied a non-linear tomography method to map P-wave velocity perturbations in a relative sense. Our tomographic images unveiled distinct instances of rapid high-velocity anomalies beneath low-velocity regions in the shallow mantle, suggesting the potential occurrence of lithospheric dripping, followed by subsequent asthenospheric upwelling. This observation offers a plausible explanation for the observed post-collisional magmatism over the Lut Block. Furthermore, to maintain the approximately 1.5-km positive residual topography across the EIR, beyond the influence of crustal properties, additional support from the hot and buoyant asthenosphere becomes crucial, particularly in the absence of a substantial lithospheric mantle.

Multiscale structures of crust-mantle beneath the South China block and their geodynamic implication

The South China block is an ideal location in which to study the impact of paleo-Pacific plate subduction on the upper mantle and crustal structures of the overriding plate, as well as the mechanism of lithospheric thinning. By integrating data from permanent seismic stations of the China Seismic Network and a dense seismic array at Xuefeng Mountain, this study employed P- and S-wave receiver functions (RFs) to image multiscale structures of crust-mantle across the South China block, thereby providing constraints for a comprehensive understanding of its tectonic evolution. The main findings are as follows: (1) The thickness of lithosphere and crust gradually decreases from northwest (NW) to southeast (SE). (2) Beneath the Xuefeng Mountain uplift zone, the P-to-S conversion from the 660 km discontinuity arrives earlier than that predicted from a one-dimensional model, other weak discontinuities can be traced in the upper mantle, including the lithosphere-asthenosphere boundary, the crust becomes thinner, and the Conrad discontinuity becomes deeper. Referring to previous results from tomography and geological reconstructions, we deduce that lithospheric delamination might have occurred beneath the Xuefeng Mountain uplift zone, and the delaminated lithosphere penetrated the upper mantle and reached the mantle transition zone.

Improving image accuracy of ambient noise data by temporary seismic arrays at different observation periods

When using ambient noise data to invert velocity and anisotropic structures, the two-station inter-correlation method requires synchronous stations. If there are multiple temporary seismic arrays with different observation periods in the study area, the seismic arrays are usually used selectively. This paper takes the Sanjiang lateral collision zone as an example, and utilizes the ambient noise data of multiple temporary seismic arrays at different observation periods to improve the accuracy of regional velocity structure and anisotropy by anchoring permanent seismic stations. In this paper, notable enhancements in S-wave velocity and azimuthal anisotropy imaging accuracy are achieved by integrating data from three temporary seismic arrays (SJ-Array, SL-Array, and ChinArray-I) with the permanent seismic network. The imaging resolutions for the S-wave velocity and azimuthal anisotropy above 40 ​km are 0.4° ​× ​0.4° and 0.5° ​× ​0.5°, respectively. In the region of the most concentrated array coverage, the imaging resolution of S-wave velocity can reach 0.33° ​× ​0.33° at depths of less than 30 ​km. These findings underscore the significant improvement in deep structure imaging accuracy by the synergistic integration of ambient noise data from multiple temporary seismic arrays.

AI based 1-D P- and S-wave velocity models for the greater alpine region from local earthquake data

SUMMARY The recent rapid improvement of machine learning techniques had a large impact on the way seismological data can be processed. During the last years several machine learning algorithms determining seismic onset times have been published facilitating the automatic picking of large data sets. Here we apply the deep neural network PhaseNet to a network of over 900 permanent and temporal broad-band stations that were deployed as part of the AlpArray research initiative in the Greater Alpine Region (GAR) during 2016–2020. We selected 384 well distributed earthquakes with ML ≥ 2.5 for our study and developed a purely data-driven pre-inversion pick selection method to consistently remove outliers from the automatic pick catalogue. This allows us to include observations throughout the crustal triplication zone resulting in 39 599 P and 13 188 S observations. Using the established VELEST and the recently developed McMC codes we invert for the 1-D P- and S-wave velocity structure including station correction terms while simultaneously relocating the events. As a result we present two separate models differing in the maximum included observation distance and therefore their suggested usage. The model AlpsLocPS is based on arrivals from ≤130 km and therefore should be used to consistently (re)locate seismicity based on P and S observations. The model GAR1D_PS includes the entire observable distance range of up to 1000 km and for the first time provides consistent P- and S-phase synthetic traveltimes for the entire Alpine orogen. Comparing our relocated seismicity with hypocentral parameters from other studies in the area we quantify the absolute horizontal and vertical accuracy of event locations as ≈2.0 and ≈6.0 km, respectively.

Layered Mantle Flow Field Associated With Plate Kinematics and Slab Modulation Underneath the Horseshoe‐Shaped Banda Arc‐Islands

AbstractThe Banda arc‐continent collision zone signifies one of the most seismically active and tectonically intricate zones. The high convergence rate across the region, coupled with the exceptionally arcuate arc and subducted slab, makes it an ideal locale for investigating interactions between plate (slab) kinematics and plastic flow in the asthenosphere, which can be diagnosed by seismic anisotropy from shear wave splitting analyses. In total, 206 pairs of splitting measurements using teleseismic SKS, SKKS, and PKS, along with 43 pairs using local S phases, are obtained by utilizing broadband seismic data from five permanent seismic stations. To reduce the ambiguity in determining the origin of anisotropy leading to the teleseismic splittings, which lack vertical resolution, crustal anisotropy is constrained according to the sinusoidal moveout of converted S phases at the Moho using receiver functions. A layered anisotropic structure based on joint analyses of the anisotropy measurements characterizing different depth layers suggests the presence of trench‐parallel flow both in the mantle wedge and the sub‐slab region. The northeastward motion of the slab, entrained by the fast‐moving Australian Plate, deflects asthenospheric materials. The modulation results in trench‐parallel plastic mantle flows and leads to the steepening of the southern portion of the asymmetric spoon‐shaped Banda slab. In the shallower part of the sub‐slab region, the northeastward Australian Plate motion produces simple shear in the transitional layer between the rigid lithosphere and the viscous asthenosphere. The shear deformation induces seismic anisotropy with resulting fast orientations in accordance with the plate motion direction.

An Accurate and Automated Convective Vortex Detection Method for Long-Duration Infrasound Microbarometer Data

Abstract Heating of the surficial layer of the atmosphere often generates convective vortices, known as “dust devils” when they entrain visible debris. Convective vortices are common on both Earth and Mars, where they affect the climate via dust loading, contribute to wind erosion, impact the efficiency of photovoltaic systems, and potentially result in injury and property damage. However, long-duration terrestrial convective vortex activity records are rare. We have developed a high-precision and high-recall method to extract convective vortex signatures from infrasound microbarometer data streams. The techniques utilizes a wavelet-based detector to capture potential events and then a template matching system to extract the duration of the vortex. Since permanent and temporary infrasound sensors networks are present throughout the globe (many with open data), our method unlocks a vast new convective vortex dataset without requiring the deployment of specialized instrumentation. Significance Statement Convective vortices, or “dust devils,” contribute to regional dust loading in Earth’s atmosphere. However, long-duration convective vortex activity records are rare. We came up with a way to autonomously detect the pressure signatures left by convective vortices striking low-frequency sound, or “infrasound,” sensors. Since permanent infrasound stations have been active for decades, our method has the potential to add orders-of-magnitude more events than previously catalogued.

Upper‐Mantle Anisotropy in the Southeastern Margin of Tibetan Plateau Revealed by Fullwave SKS Splitting Intensity Tomography

AbstractThe southeastern margin of the Tibetan Plateau has undergone complex deformation since the Cenozoic, resulting in a high level of seismicity and seismic hazard. Knowledge about the seismic anisotropy provides important insight about the deformation mechanism and the regional seismotectonics beneath this tectonically active region. In this study, we conduct fullwave multi‐scale tomography to investigate the seismic anisotropy in the southeastern margin of the Tibetan Plateau. Broadband records at 111 permanent stations in the region from 470 teleseismic events are used to obtain 5,216 high‐quality SKS splitting intensity (SI) measurements, which are then inverted in conjunction with 3D sensitivity kernels to obtain an anisotropic model with multi‐scale resolution. Resolution tests show that our data set recovers anisotropy anomalies reasonably well on the scale of 1° × 1° horizontally and ∼100 km vertically. Our result suggests that in the southeastern margin of the Tibetan Plateau the deformation in the lithosphere and asthenosphere are decoupled. The anisotropy in the lithosphere varies both laterally and vertically as a result of dynamic interactions of neighboring blocks as well as lithospheric reactivation. The anisotropy in the asthenosphere largely follows the direction of regional absolute plate motion. The SKS splittings observed at the surface are shown to be consistent with the vertical integral of our depth‐dependent anisotropy model over lithospheric and asthenospheric depths.