The 2023 M_(w)6.8 Adassil Earthquake(Chichaoua,Morocco)on a steep reverse fault in the deep crust and its geodynamic implications

The Mw 6.8 Adassil earthquake that occurred in the High Atlas on September 8,2023,was a catastrophic event that provided a rare opportunity to study the mechanics of deep crustal seismicity.This research aimed to decipher the rupture characteristics of the Adassil earthquake by analyzing teleseismic waveform data in conjunction with interferometric synthetic aperture radar(InSAR)observations from both ascending and descending orbits.Our analysis revealed a reverse fault mechanism with a centroid depth of approximately 28 km,exceeding the typical range for crustal earthquakes.This result suggests the presence of cooler temperatures in the lower crust,which facilitates the accumulation of tectonic stress.The earthquake exhibited a steep reverse mechanism,dipping at 70°,accompanied by minor strike-slip motion.Within the geotectonic framework of the High Atlas,known for its volcanic legacy and resulting thermal irregularities,we investigated the potential contributions of these factors to the initiation of the Adassil earthquake.Deep seismicity within the lower crust,away from plate boundaries,calls for extensive research to elucidate its implications for regional seismic hazard assessment.Our findings highlight the critical importance of studying and preparing for significant seismic events in similar geological settings,which would provide valuable insights into regional seismic hazard assessments and geodynamic paradigms.

Assessment of natural and anthropogenic impacts on terrestrial water storage in the Loess Plateau based on different types of GRACE/GRACE-FO solutions

Changes in water resource storage are inevitable due to climate change and human activities,thus understanding alterations in water storage within a specific region is imperative for the planning and management of water resources.Data from the Gravity Recovery and Climate Experiment(GRACE)satellite mission are extensively employed to analyze large-scale total terrestrial water storage anomalies(TWSA).In this study,we derived a more reliable TWSA using different types of GRACE gravity models,which served as the basis for evaluating spatial and temporal variations in total terrestrial water storage and its hydrological components(soil moisture and groundwater)across the Loess Plateau.Additionally,we analyzed the impact of natural and anthropogenic influences on water storage in the Loess Plateau,categorizing them into primary and secondary influences,utilizing data on climate and human activities.The findings revealed a declining trend in the overall TWSA of the Loess Plateau,with a rate of decrease at-0.65±0.05 cm/yr from 2003 to 2020(P<0.01).As the direct factors affecting TWSA,soil moisture dominated the change of TWSA before 2009,and groundwater dominated the change of TWSA after 2009.Spatially,there was variability in the changes of TWSA in the Loess Plateau.More in-depth studies showed that soil moisture changes in the study area were primarily driven by evapotranspiration and temperature,with precipitation and vegetation cover status playing a secondary role.Human activities had a secondary effect on soil moisture in some sub-regions.Population change and agricultural development were major factors in altering groundwater storage in the study area.Other than that,groundwater was influenced by natural factors to a limited extent.These findings provided valuable insights for local governments to implement proactive water management policies.

Physical analysis on improving the recovery accuracy of the Earth's gravity field by a combination of satellite observations in along-track and cross-track directions

The physical investigations on the accuracy improvement to the measurement of the Earth＇s gravity field recovery are carried out based on the next-generation Pendulum-A/B out-of-plane twin-satellite formation in this paper. Firstly, the Earth＇s gravity field complete up to degree and order 100 is, respectively, recovered by the collinear and pendulum satellite formations using the orbital parameters of the satellite and the matching accuracies of key payloads from the twin GRACE satellites. The research results show that the accuracy of the Earth＇s gravity field model from the Pendulum-A/B satellite formation is about two times higher than from the collinear satellite formation, and the further improvement of the determination accuracy of the Earth＇s gravity field model is feasible by the next-generation Pendulum-A/B out-of-plane twin-satellite formation. Secondly, the Earth＇s gravity field from Pendulum-A/B complete up to degree and order 100 is accurately recovered based on the orbital parameters of the satellite （e.g., an orbital altitude of 400 km, an intersatellite range of 100 km, an orbital inclination of 89° and an orbital eccentricity of 0.001）, the matching accuracies of space- borne instruments （e.g. 10-6 m in the intersatellite range, 10-3 m in the orbital position, 10-6 m/s in orbital velocity, and 10-11 m/s2 in non-conservative force）, an observation time of 30 days and a sampling interval of 10 s. The measurement accuracy of the Earth＇s gravity field from the next-generation Pendulum-A/B out-of-plane twin-satellite formation is full of promise for being improved by about l0 times compared with that from the current GRACE satellite formation. Finally, the physical requirements for the next-generation Pendulum-A/B out-of-plane twin-satellite formation are analyzed, and it is proposed that the satellite orbital altitude be preferably designed to be close to 400±50 km and the matching precision of key sensors from the Pendulum-A/B mission be about one order of magnitude higher than from the GRACE program.

High-resolution seismic tomography of the upper crust in the source region of the April 20,2013 Lushan earthquake and its seismotectonic implications

Both P- and S-wave arrivals were collected for imaging upper crustal structures in the source region of the April 20, 2013 Lushan earthquake. High-resolution, three- dimensional P and S velocity models were constructed by travel-time tomography. Moreover, more than 3700 after- shocks of the Lushan earthquake were relocated via a grid search method. The P- and S-wave velocity images of the upper crust show largely similar characters, with high and low velocity anomalies, which mark the presence of sig- nificant lateral and vertical heterogeneity at the source region of the Lushan earthquake. The characteristics of the velocity anomalies also reflect the associated surface geo- logical tectonics in this region. The distributions of high velocity anomalies of both P- and S-waves to 18 km depth are consistent with the distributions of relocated after- shocks, suggesting that most of the ruptures were localized inside the high velocity region. In contrast, low P and S velocities were found in the surrounding regions without aftershocks, especially in the region to the northeast of the Lushan earthquake. For the relocated aftershocks of the Lushan earthquake from this study, we found that mostaftershocks were concentrated in a zone of about 40 km long and 20 km wide, and were located in the hanging wall of Dayi-Mingshan fault. The focal depths of aftershocks increase from the southeast to the northwest region in the direction perpendicular to the fault strike, suggesting that the fault ruptured at an approximate dip angle of 45°. The main depths of the aftershocks in the northwest of the main shock are significantly shallower than expected, revealing the different seismogenic conditions in the source region.

Theories and applications of earthquake-induced gravity variation:Advances and perspectives

Earthquake-induced gravity variation refers to changes in the earth’s gravity field associated with seismic activities.In recent years,development in the theories has greatly promoted seismic deformation research,laying a solid theoretical foundation for the interpretation and application of seismological gravity monitoring.Traditional terrestrial gravity measurements continue to play a significant role in studies of interseismic,co-seismic,and post-seismic gravity field variations.For instance,superconducting gravimeter networks can detect co-seismic gravity change at the sub-micro Gal level.At the same time,the successful launch of satellite gravity missions(e.g.,the Gravity Recovery and Climate Experiment or GRACE)has also facilitated applied studies of the gravity variation associated with large earthquakes,and several remarkable breakthroughs have been achieved.The progress in gravity observation technologies(e.g.,GRACE and superconducting gravimetry)and advances in the theories have jointly promoted seismic deformation studies and raised many new research topics.For example,superconducting gravimetry has played an important role in analyses of episodic tremor,slow-slip events,and interseismic strain patterns;the monitoring of transient gravity signals and related theories have provided a new perspective on earthquake early warning systems;the mass transport detected by the GRACE satellites several months before an earthquake has brought new insights into earthquake prediction methods;the use of artificial intelligence to automatically identify tiny gravity change signals is a new approach to accurate and rapid determination of earthquake magnitude and location.Overall,many significant breakthroughs have been made in recent years,in terms of the theory,application,and observation measures.This article summarizes the progress,with the aim of providing a reference for seismologists and geodetic researchers studying the phenomenon of gravity variation,advances in related theories and applications,and future research directions in this discipline.

Co-seismic deformation and gravity changes of the 2011 India-Nepal and Myanmar earthquakes

Co-seismic deformation and gravity field changes caused by the 2011 Mw6. 8 Myanmar and Mw6. 9 India-Nepal earthquakes are calculated with a finite-element model and an average-slip model, respectively, based on the multi-layered elastic half-space dislocation theory. The calculated maximum horizontal displace- ment of the Myanmar earthquake is 36 era, which is larger than the value of 9. 5 cm for the India-Nepal earth- quake. This difference is attributed to their different focal depths and our use of different models. Except cer- tain differences in the near field, both models give similar deformation and gravity results for the Myanmar event.

Machine-learning-facilitated earthquake and anthropogenic source detections near the Weiyuan Shale Gas Blocks,Sichuan,China

Seismic hazard assessment and risk mitigation depend critically on rapid analysis and characterization of earthquake sequences.Increasing seismicity in shale gas blocks of the Sichuan Basin,China,has presented a serious challenge to monitoring and managing the seismicity itself.In this study,to detect events we apply a machine-learning-based phase picker(PhaseNet)to continuous seismic data collected between November 2015 and November 2016 from a temporary network covering the Weiyuan Shale Gas Blocks(SGB).Both P-and S-phases are picked and associated for location.We refine the velocity model by using detected explosions and earthquakes and then relocate the detected events using our new velocity model.Our detections and absolute relocations provide the basis for building a high-precision earthquake catalog.Our primary catalog contains about 60 times as many earthquakes as those in the catalog of the Chinese Earthquake Network Center(CENC),which used only the sparsely distributed permanent stations.We also measure the local magnitude and achieve magnitude completeness of ML0.We relocate clusters of events,showing sequential migration patterns overlapping with horizontal well branches around several well pads in the Wei202 and Wei204 blocks.Our results demonstrate the applicability of a machine-learning phase picker to a dense seismic network.The algorithms can facilitate rapid characterization of earthquake sequences.

Effect of different inter-satellite range on measurement precision of Earth’s gravitational field from GRACE

The precision of Earth＇s gravitational field from GRACE up to degree and order 120 was studied for different inter-satellite ranges using the improved energy conservation principle. Our simulated result shows that： For long wavelength （L≤20） at degree 20, the cumulative geoid-height error gradually decreased with increasing range, from 0. 052 cm for 110 km to 1. 156 times and 1. 209 times as large for 220 km and 330 kin, respectively. For medium-wavelength （ 100 ≤ L ≤ 120） at degree 120, the cumulative geoid-height error de- creased from 13. 052 cm for 110 km, to 1. 327 times and 1. 970 times as large for the ranges of 220 km and 330 km, respectively; By adopting an optimal range of 220 ± 50 km, we can suppress considerably the loss of precision in the measurement of the Earth＇ s long-wavelength and medium-wavelength gravitational field.

Far-field coseismic gravity changes related to the 2015 MW7.8 Nepal(Gorkha)earthquake observed by superconducting gravimeters in Chinese mainland

Using data from five SGs at four stations in Chinese mainland,obvious permanent gravity changes caused by the 2015 MW7.8 Nepal(Gorkha)earthquake were detected.We analyzed the gravity effects from ground vertical deformation(VD)using co-site continuous GPS(cGPS)data collocated at the Lijiang and the Wuhan station,and hydrological effects using GLDAS models and groundwater level records.After removing these effects,SG observations before and after the earthquake revealed obvious permanent gravity changes:−3.0μGal,7.3μGal and 8.0μGal at Lhasa,Lijiang and Wuhan station,respectively.We found that the gravity changes cannot be explained by the results of dislocation theory.

Preliminary study on aftershock decay rate of the 2013 Ms 7.0 Lushan earthquake

We obtained a catalog of early aftershocks of the 2013 Lushan earthquake by examining waveform from a nearby station MDS which is 30.2 km far away from the epicenter, and then we analyzed the relation between aftershock rate and time. We used time-window ratio method to identify aftershocks from continuous waveform data and compare the result with the catalog provided by China Earthquake Networks Center （CENC）. As expected, a significant amount of earthquakes is missing in CENC catalog in the 24 h after the main shock. Moreover, we observed a steady seismicity rate of aftershocks nearly in the first 10,000 s before an obvious power-law decay of aftershock activity. We consider this distinct early stage which does not fit the Omori law with a constant p （p - 1） value as early aftershock deficiency （EAD）, as proposed by previous studies. Our study suggests that the main shock rupture process is different from aftershocks＇ processes, and EAD can vary in different cases as compared to earthquakes of strike-slip mechanism in California.

Precise geoid computation using Stokes-Helmert’s scheme and strict integrals of topographic effects

The high-precision local geoid model was computed based on the improved Stokes-Helmert0 s boundary value problem and strict integrals of topographic effects. This proposed method involves three steps.First, the mathematical form of Stokes-Helmert0 s boundary value problem was derived, and strict computational formulas regarding topographic effects were provided to overcome the disadvantage of planar approximations. Second, a gravimetric geoid model was constructed using the proposed StokesHelmert0 s scheme with a heterogeneous data set. Third, a least squares adjustment method combined with a multi-surface function model was employed to remove the bias between the gravimetric geoid model and the GNSS/leveling data and to refine the final local geoid model. The accuracy of the final geoid model was evaluated using independent GNSS/leveling data. Numerical results show that an external precision of 1.45 cm is achievable.

Multi-scale Analysis of the Relationships between Solar Activity, CO_(2) and Global Surface Temperature

To reveal whether the dynamics of solar activity precede those of global temperature, especially in terms of global warming, the relationship between total solar irradiance(TSI), which is treated as a proxy of solar activity, and global surface temperature(GST) is investigated in the frequency domain using wavelet coherence. The results suggest that the effect of TSI on GST is mainly reflected on the characteristic scale around 22 yr, and variations in TSI lead to changes in GST with some delay effect as shown by the phase difference arrows. However, this implicated relationship has been perturbed by excessive CO_(2) emissions since 1960. Through the combination of co-integration analysis and wavelet coherence, the hidden relationship between TSI and GST has been uncovered without the CO_(2) effect and the results further indicate that TSI has a positive effect on GST at the characteristic scale around 22 yr with a 3 yr lead.

Marginal weakening promoted rift propagation of the Petermann Ice Shelf in northwestern Greenland from 2016 to 2022

Ice-shelf rifts are precursors of glacier calving, and thus they serve as indicators of ice shelf instability, especially under oceanic and atmospheric warming conditions. Therefore, understanding the dynamic processes underlying rift propagation and the associated damage mechanisms is essential to evaluate ice-shelf instability and to predict glacier calving. In this study, we investigated the effect of marginal weakening on rift propagation on the ice shelf of the Petermann Glacier, among the largest in Greenland, during 2016-2022. First, we analyzed satellite optical images to monitor rift growth(length and width) by tracking the tip trajectory of three large rifts identified on the Petermann Ice Shelf. Then, we estimated rift depth using ArcticDEM and ICESat-2 data. Our results indicated consistent increases of the rift widths and depths over the study period, with mean values of 133 m·a^(-1) and 0.3 m·a^(-1), respectively. We also combined remote-sensing observations with an ice-sheet numerical model to calculate the stress and damage fields on the Petermann Ice Shelf and to assess the ice shelf margin stability and strength. We determined that damage and lateral shear in the fracture zone degraded ice shelf integrity by decreasing the contact length with the fjord wall. In conclusion, marginal weakening effectively promoted rift propagation on the Petermann Ice Shelf, increasing the risk of future glacier calving.

Status of UnDifferenced and Uncombined GNSS Data Processing Activities in China

With the continued development of multiple Global Navigation Satellite Systems(GNSS)and the emergence of various frequencies,UnDifferenced and UnCombined(UDUC)data processing has become an increasingly attractive option.In this contribution,we provide an overview of the current status of UDUC GNSS data processing activities in China.These activities encompass the formulation of Precise Point Positioning(PPP)models and PPP-Real-Time Kinematic(PPP-RTK)models for processing single-station and multi-station GNSS data,respectively.Regarding single-station data processing,we discuss the advancements in PPP models,particularly the extension from a single system to multiple systems,and from dual frequencies to single and multiple frequencies.Additionally,we introduce the modified PPP model,which accounts for the time variation of receiver code biases,a departure from the conventional PPP model that typically assumes these biases to be time-constant.In the realm of multi-station PPP-RTK data processing,we introduce the ionosphere-weighted PPP-RTK model,which enhances the model strength by considering the spatial correlation of ionospheric delays.We also review the phase-only PPP-RTK model,designed to mitigate the impact of unmodelled code-related errors.Furthermore,we explore GLONASS PPP-RTK,achieved through the application of the integer-estimable model.For large-scale network data processing,we introduce the all-in-view PPP-RTK model,which alleviates the strict common-view requirement at all receivers.Moreover,we present the decentralized PPP-RTK data processing strategy,designed to improve computational efficiency.Overall,this work highlights the various advancements in UDUC GNSS data processing,providing insights into the state-of-the-art techniques employed in China to achieve precise GNSS applications.

A new approach for separating mixed model parameters:application to simultaneous inversion of earthquake source parameters

A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single model space shows instabilities and trade-offs of the solutions. We separate the model space into N-subspaces based on their physical properties or computational convenience and solve the N-subspaces systems by damped least-squares and singular-value decomposition. Since the condition number of each subsystem is smaller than that of the single global system,the approach can greatly increase the stability of the inversion. We also introduce different damping factors into the subsystems to reduce the tradeoffs between the different parameters. The damping factors depend on the conditioning of the subsystems and may be adequately chosen in a range from 0.1 % to 10 % of the largest singular value. We illustrate the method with an example of simultaneous determination of source history,source geometry,and hypocentral location from regional seismograms,although it is applicable to any geophysical inversion.

Upper-mantle velocity structures beneath the Tibetan Plateau and surrounding areas inferred from triplicated P waveforms

P-wave waveforms in the distance range between 12°and 30°were analyzed to investigate upper-mantle P velocity structures beneath the Tibetan Plateau and surrounding areas.The waveform data from 504 earthquakes with magnitudes larger than 5.0 between 1990 and 2005 that occurred within 30°from the center of the Plateau were modelled.We divided the study area into 6 regions and modeled upper-mantle-distance P waveforms with turning points beneath each region separately.The results show that the uppermantle P-wave velocity structures beneath India,the Himalayas,and the Lhasa Terrane are similar and contain a high-velocity lid about 250 km thick.The upper-mantle velocities down to 200 km beneath the Qiangtang Terrane,the Tarim Basin,and especially the Songpan-GarzêTerrane are lower than those in the south.The 410-km discontinuity beneath these two terranes is elevated by about 20 km.Highvelocity anomalies are found in the transition zone below 500 km under the Lhasa and Qiangtang Terranes.The results suggest that the Tibetan Plateau was generated by thrusting of the Indian mantle lithosphere under the southern part of Tibet.Portions of the thickened Eurasian mantle lithosphere were delaminated;they are now sitting in the transition zone beneath southern Tibet and atop of the 410-km discontinuity underneath northern Tibet.

Vertical crustal deformation velocity and its influencing factors over the Qinghai-Tibet Plateau based on satellite gravity data

The uplift of the Qinghai-Tibet Plateau(TP)strongly influences climate change,both regionally and globally.Surface observation data from this region have limited coverage and are difficult to obtain.Consequently,the vertical crustal deformation velocity(VCDV)distribution of the TP is poorly constrained.In this study,the VCDV from the TP was inverted by using data from the gravity recovery and climate experiment(GRACE).We were able to obtain the vertical crustal movement by deducting the hydrological factors,based on the assumption that the gravity signal detected by GRACE is mainly composed of hydrological factors and vertical crustal movement.From the vertical crustal movement,we inverted the distribution of the VCDV across the TP.The results showed that the VCDV of the southern,eastern,and northern TP is~1.1 mm/a,~0.5 mm/a,and−0.1 mm/a,respectively,whereas that of the region between the Qilian Haiyuan Fault and the Kunlun Fault is~0.0 mm/a.These results are consistent with the distribution of crustal deformation,thrust earthquakes and faults,and regional lithospheric activity.The hydrology,crustal thickness,and topographic factors did not change the overall distribution of the VCDV across the TP.The influence of hydrological factors is marked,with the maximum differences being approximately−0.4 mm/a in the northwest and 1.0 mm/a in the central area.The results of this study are significant for understanding the kinematics of the TP.

Template matching for simple waveforms with low signal-to-noise ratio and its application to icequake detection

Template matching is a useful method to detect seismic events through waveform similarity between two signals.The traditional template matching method works well in detecting small tectonic earthquakes.However,the method has some difficulty when the signals have relatively low signal-to-noise ratios(SNRs)and simple shapes,e.g.a sinusoidal function.In this study,we modify the traditional template matching approach for this situation.We first construct a virtual three-component seismic station using vertical-component waveforms recorded by three stations.Next,we select a template event from the virtual station,and apply the traditional template matching.We then verify this method by detecting icequakes with simple waveforms on the Urumqi Glacier No.1 and compare the results with those from the short-term-averages over long-term-average(STA/LTA),the REST method,and traditional template matching method.It can be concluded that the modified template matching method using virtual stations has some advantages for seismic data with low SNRs.

Crustal S-wave velocity structure of the Yellowstone region using a seismic ambient noise method

The Yellowstone volcano is one of the largest active volcanoes in the world, and its potential hazards demand detailed seismological and geodetic studies. Previous studies with travel time tomography and receiver functions have revealed a low-velocity layer in the crust beneath the Yellowstone volcano, suggesting the presence of a magma chamber at depth. We use ambient seismic noise from regional seismic stations to retrieve short-period surface waves and then study the shallow shear velocity structure of the Yellowstone region by surface wave dispersion analysis. We first obtained a crustal model of the area outside of the Yellowstone volcano and then constructed an absolute shear wave velocity structure in combination with receiver function results for the crust beneath the Yellowstone volcano. The velocity model shows a low-velocity layer with shear velocity at around 1.3 km/s, suggesting that a large-scale magma chamber exists at shallow levels within the crust of the Yellowstone volcanic region.

On Inversion of Continuous Wavelet Transform

This study deduces a general inversion of continuous wavelet transform (CWT) with timescale being real rather than positive. In conventional CWT inversion, wavelet’s dual is assumed to be a reconstruction wavelet or a localized function. This study finds that wavelet’s dual can be a harmonic which is not local. This finding leads to new CWT inversion formulas. It also justifies the concept of normal wavelet transform which is useful in time-frequency analysis and time-frequency filtering. This study also proves a law for CWT inversion: either wavelet or its dual must integrate to zero.