Source model of the 2018 M_(w)7.5 Papua New Guinea earthquake constrained by InSAR measurements: Insight into a blind rupture beneath the New Guinea Highlands

The Southern Highland Fold and Thrust Belt(SHFTB),the boundary of the Australian plate and the New Guinea Highland block,significantly contributes to the convergent deformation along the plate bound-ary.However,due to the lack of observation data,the detailed slip pattern of the SHFTB and the orogenic mechanism beneath the New Guinea Highlands remains controversial.On 25 February 2018,the M_(w)7.5 Papua New Guinea(PNG)earthquake struck the southeastern segment of the SHFTB.The detailed rupture characteristics of this event is significant for further clarifying the inter-seismic slip pattern along the SHFTB.Here,the coseismic deformation field of this earthquake was obtained using high-resolution ALOS-2 satellite images.We find that the 2018 M_(w)7.5 PNG earthquake ruptured a large-scaled fault(SHFTB)extending to the lower crust(deeper than 20 km)beneath the New Guinea Highlands,with a dip angle of 24°.The slips on the fault plane are equivalent to moment magnitudes of M_(w)7.51.Three major asperities with thrust-dominated slip of up to 3.94 m are detected on the fault plane.This finding implies that the slip pattern on the eastern segment of the SHFTB is dominated by thrust,rather than with significant sinistral movement,as previously reported.The tectonic deformation across the New Guinea Highlands is possibly concentrated on the large-scale fault SHFTB and primarily controls the intra-continental orogeny in the central Papua New Guinea.

Present-day movement characteristics of the Qinghai Nanshan fault and its surrounding area from GPS observation

The Qinghai Nanshan fault is a larger fault in the Northeastern Xizang Plateau.In previous studies,its movement characteristics are mainly investigated with geological and seismic observations,and the tectonic transformation role of the fault on its east is not yet clear.This study uses data fusion to obtain denser GPS observations near the Qinghai Nanshan fault.Based on tectonic characteristics,we establish a block model to investigate the fault slip rate,locking degree,and slip deficit.The results show that the Qinghai Nanshan fault slip rate is characterized by sinistral and convergent movement.Both the sinistral and convergent rates display a decreasing trend from west to east.The locking degree and slip deficit are higher in the western segment(with an average of about 0.74 and 1.1 mm/a)and lower in the eastern segment.Then,we construct a strain rate field using GPS observations to analyze the regional strain characteristics.The results indicate that along the fault,the western segment shows a larger shear strain rate and negative dilation rate.Regional earthquake records show that the frequency of earthquakes is lower near the fault.The joint results suggest that the western segment may have a higher earthquake risk.In addition,the insignificant fault slip rate in the eastern segment may indicate that it does not participate in the tectonic transformation among the Riyueshan,Lajishan,and West Qinling faults.

Crustal and uppermost mantle structure of the northeastern Qinghai-Xizang Plateau from joint inversion of surface wave dispersions and receiver functions with P velocity constraints

Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.

Review of the Research Progress on Static Earth Gravity Field and Vertical Datum in China during 2019—2023

The contribution presents the representative research progress on global static gravity field modeling,regional geoid/quasigeoid determination,vertical datum study,as well as the theory,algorithm and software for gravity field study in China from 2019 to 2023,which are the highlights of the chapter 6“Progress in Earth Gravity Model and Vertical Datum”in the“2019—2023 China National Report on Geodesy”that submitted to the International Association of Geodesy(IAG).In addition,suggestions are proposed to promote the research in the fields of earth gravity field,geoid/quasigeoid and vertical datumin China according to trends of international geodesy and related disciplines.

Long-term and short-term stability characteristics of receiver inter system bias for BDS3/BDS2

The study of inter-system bias(ISB)is important for multi-system fusion and the performance of different signal compatibility.In this paper,the stability of ISB at the BDS3/BDS2 receiver end is calculated and analyzed for different time spans(DOY 060~090 in 2021)from a total of 31 MGEX and iGMAS stations.We adopted two estimation strategies,random walk and constant approach,using the precision products of orbit and clock bias provided by WUM,the influence of which on ISB was also analyzed.Our results showed that the ISB value varied little within a day,and the mean of daily ISB standard deviation was only 0.037 m when the observation condition was good.The signal reception was continuous,indicating a high ISB stability for one day.If extending the time series to one month,however,the ISB standard deviation calculated by constant approach,in which a constant ISB is estimated on a daily basis was about 0.1 m,and the results of adjacent days were not continuous,with no apparent pattern.Concerning the random walk approach,the obtained ISB time series also had a jump,and the conclusion was the same as that of the constant strategy.Besides,receiver types showed a strong regularity in ISB numerical situation,and the distribution of ISB values corresponding to the same receiver type was relatively close.Therefore,we conclude that the ISB parameters remain stable in the short term(one day)and less stable in the long-term period.It is recommended that the ISB term should be set as a constant estimate every day in BDS3/BDS2 solutions,regardless of receiver type consistency.

Extraction and analysis of tree canopy height information in high-voltage transmission-line corridors by using integrated optical remote sensing and LiDAR

Traditional inspection methods cannot quickly and accurately monitor tree barriers and safeguard the transmission lines.To solve these problems,in this study,we proposed a rapid canopy height information extraction method using optical remote sensing and LiDAR,and used UAV optical imagery with LiDAR to monitor the height of trees in a university and a high-voltage transmission line corridor in the Ningxia region.The results showed that the relative error of tree height extraction using UAV optical images was less than 5%,and the lowest relative error was 0.11%.The determination coefficient R^(2) between the optical image tree height extraction results and the measured tree height was 0.97,thus indicating a high correlation for both.In the field of tree barrier monitoring,the determination coefficient R^(2) of tree height extracted using airborne LiDAR point cloud,and canopy height model(CHM)and of the measured tree height were 0.947 and 0.931,respectively.The maximum and minimum relative error in tree height extraction performed using point cloud was 2.91%and 0.2%,respectively,with an extraction accuracy of over 95%.The experimental results demonstrated that it is feasible to use UAV optical remote sensing and LiDAR in monitoring tree barriers and tree height information extraction quickly and accurately,which is of great significance for the risk assessment and early warning of tree barriers in transmission-line corridors.

Progress of Geodesy Related Ionosphere from Chinese Scientists in the Period of 2019—2023

The ionosphere is the ionized part of the upper atmosphere of the Earth,which plays an important role in atmospheric electricity and forms the inner edge of the magnetosphere.It influences radio propagation significantly,such as the Global Navigation Satellite System(GNSS).Meanwhile,the GNSS is also an essential technique for sensing the variation of ionosphere.During the years of 2019—2023,a large number of Chinese geodesy scientists devoted much efforts to the geodesy related ionosphere.Due to the very limited length,the achievements are carried out from the following six aspects,including:①The ionospheric correction models for BDS and BDSBAS;②Real-time global ionospheric monitoring and modeling;③The ionospheric 2D and 3D modeling based on GNSS and LEO satellites;④The ionospheric prediction based on artificial intelligence;⑤The monitoring and mitigation of ionospheric disturbances for GNSS users;⑥The ionospheric related data products and classical applications.

Coseismic Coulomb stress changes induced by a 2020-2021 M_(W)>7.0 Alaska earthquake sequence in and around the Shumagin gap and its influence on the Alaska-Aleutian subduction interface

Three M_(W)>7.0 earthquakes in 2020-2021 occurred in the Shumagin seismic gap and its adjacent area of the Alaska-Aleutian subduction zone,including the Mw7.8 Simeonof thrust earthquake on July 22,2020,the M_(W)7.6 Sand Point strike-slip earthquake on October 19,2020,and the M_(W)8.2 Chignik thrust earthquake on July 29,2021.The spatial and temporal proximity of these three earthquakes prompts us to probe stress-triggering effects among them.Here we examine the coseismic Coulomb stress change imparted by the three earthquakes and their influence on the subduction interface.Our results show that:(1)The Simeonof earthquake has strong loading effects on the subsequent Sand Point and Chignik earthquakes,with the Coulomb stress changes of 3.95 bars and 2.89 bars,respectively.The Coulomb stress change caused by the Sand Point earthquake at the hypocenter of the Chignik earthquake is merely around 0.01 bars,suggesting the negligible triggering effect on the latter earthquake;(2)The triggering effects of the Simeonof,Sand Point,and Chignik earthquakes on aftershocks within three months are not well pronounced because of the triggering rates of 38%,14%,and 43%respectively.Other factors may have played an important role in promoting the occurrence of these aftershocks,such as the roughness of the subduction interface,the complicated velocity structure of the lithosphere,and the heterogeneous prestress therein;(3)The three earthquakes caused remarkable coseismic Coulomb stress changes at the subduction interface nearby these mainshocks,with an average Coulomb stress change of 3.2 bars in the shallow region directly inwards the trench.

Suitable region of dynamic optimal interpolation for efficiently altimetry sea surface height mapping

The dynamic optimal interpolation(DOI)method is a technique based on quasi-geostrophic dynamics for merging multi-satellite altimeter along-track observations to generate gridded absolute dynamic topography(ADT).Compared with the linear optimal interpolation(LOI)method,the DOI method can improve the accuracy of gridded ADT locally but with low computational efficiency.Consequently,considering both computational efficiency and accuracy,the DOI method is more suitable to be used only for regional applications.In this study,we propose to evaluate the suitable region for applying the DOI method based on the correlation between the absolute value of the Jacobian operator of the geostrophic stream function and the improvement achieved by the DOI method.After verifying the LOI and DOI methods,the suitable region was investigated in three typical areas:the Gulf Stream(25°N-50°N,55°W-80°W),the Japanese Kuroshio(25°N-45°N,135°E-155°E),and the South China Sea(5°N-25°N,100°E-125°E).We propose to use the DOI method only in regions outside the equatorial region and where the absolute value of the Jacobian operator of the geostrophic stream function is higher than1×10^(-11).

A new method for deriving broad-band polar motion geodetic excitations

While the geodetic excitationχ(t)of polar motion p(t)is essential to improve our understanding of global mass redistributions and relative motions with respect to the terrestrial frame,the widely adopted method to deriveχ(t)from p(t)has biases in both amplitude and phase responses.This study has developed a new simple but more accurate method based on the combination of the frequency-and time-domain Liouville's equation(FTLE).The FTLE method has been validated not only with 6-h sampled synthetic excitation series but also with daily and 6-h sampled polar motion measurements as well asχ(t)produced by the interactive webpage tool of the International Earth Rotation and Reference Systems Service(IERS).Numerical comparisons demonstrate thatχ(t)derived from the FTLE method has superior performances in both the time and frequency domains with respect to that obtained from the widely adopted method or the IERS webpage tool,provided that the input p(t)series has a length around or more than 25 years,which presents no practical limitations since the necessary polar motion data are readily available.The FTLE code is provided in the form of Mat Lab function.

Deformation monitoring of long-span railway bridges based on SBAS-InSAR technology

The deformation monitoring of long-span railway bridges is significant to ensure the safety of human life and property.The interferometric synthetic aperture radar(In SAR)technology has the advantage of high accuracy in bridge deformation monitoring.This study monitored the deformation of the Ganjiang Super Bridge based on the small baseline subsets(SBAS)In SAR technology and Sentinel-1A data.We analyzed the deformation results combined with bridge structure,temperature,and riverbed sediment scouring.The results are as follows:(1)The Ganjiang Super Bridge area is stable overall,with deformation rates ranging from-15.6 mm/yr to 10.7 mm/yr(2)The settlement of the Ganjiang Super Bridge deck gradually increases from the bridge tower toward the main span,which conforms to the typical deformation pattern of a cable-stayed bridge.(3)The sediment scouring from the riverbed cause the serious settlement on the bridge’s east side compared with that on the west side.(4)The bridge deformation negatively correlates with temperature,with a faster settlement at a higher temperature and a slow rebound trend at a lower temperature.The study findings can provide scientific data support for the health monitoring of long-span railway bridges.

An improved typhoon monitoring model based on precipitable water vapor and pressure

The potential of monitoring the movement of typhoons using the precipitable water vapor(PWV) has been confirmed. However, monitoring the movement of typhoon is focused on PWV, making it difficult to describe the movement of a typhoon in detail minutely and resulting in insufficient accuracy. Hence,based on PWV and meteorological data, we propose an improved typhoon monitoring mode. First, the European Centre for Medium-Range Weather Forecasts Reanalysis 5-derived PWV(ERA5-PWV) and the Global Navigation Satellite System-derived PWV(GNSS-PWV) were compared with the reference radiosonde PWV(RS-PWV). Then, using the PWV and atmospheric parameters derived from ERA5, we discussed the anomalous variations of PWV, pressure(P), precipitation, and wind speed during different typhoons. Finally, we compiled a list of critical factors related to typhoon movement, PWV and P. We developed an improved multi-factor typhoon monitoring mode(IMTM) with different models(i.e.,IMTM-I and IMTM-II) in different cases with a higher density of GNSS observation or only Numerical Weather Prediction(NWP) data. The IMTM was evaluated through the reference movement speeds of HATO and Mangkhut from the China Meteorological Observatory Typhoon Network(CMOTN). The results show that the root mean square(RMS) of the IMTM-I is 1.26 km/h based on ERA5-P and ERA5-PWV,and the absolute bias values are mostly within 2 km/h. Compared with the models considering the single factor ERA5-P/ERA5-PWV, the RMS of the IMTM-I is improved by 26.3% and 38.5%, respectively. The IMTM-II model manifests a residual of only 0.35 km/h. Compared with the single-factor model based on GNSS-PWV/P, the residual of the IMTM-II model is reduced by 90.8% and 84.1%, respectively. These results propose that the typhoon movement monitoring approach combining PWV and P has evident advantages over the single-factor model and is expected to supplement traditional typhoon monitoring.

Constraint on the focal mechanism of the 2011 Tohoku earthquake from the radial modes

Different from other normal modes of the Earth’s free oscillation that depend on all the six components(M_(rr),M_(tt),M_(pp),M_(rt),M_(rp),and M_(tp))of the centroid moment tensor,the amplitudes of the radial modes depend on the M_(rr)component(e.g.,scalar moment(M_(0)),dip(δ),and slip(λ))and hypocenter depth of the focal mechanism,and hence can be easily used to constrain these parameters of the focal mechanism.In this study,we use the superconducting gravimeter(SG)records after the 2011 Tohoku earthquake to analyze the radial modes_(0)S_(0)and_(1)S_(0).Based on the solutions of the focal mechanism provided by the GCMT and USGS,we can obtain the theoretical amplitudes of these two radial modes.Comparing the theoretical amplitudes with the observation amplitudes,it is found that there are obvious differences between the former and the latter,which means that the GCMT and USGS focal mechanisms cannot well represent the real focal mechanism of the 2011 event.Taking the GCMT solution as a reference and changing the depth and the three parameters of the M_(rr)moment,the scalar moment(M_(0))and the dip(δ)have significant influences,but the effects of the slip(λ)and the depth are minor.After comparisons,we provide a new constraint(M_(0)=5.8±0.09×10^(22)N·m,δ=10.1±0.08°,λ=88°,and depth=20 km)for the focal mechanism of the 2011 event.In addition,we further determine the center frequency(1.631567±2.6e^(-6)mHz)and quality factor(2046.4±50.1)of the_(1)S_(0)mode.

The GNSS PWV retrieval using non-observation meteorological parameters based on ERA5 and its relation with precipitation

The pressure and temperature significantly influence precipitable water vapor(PWV) retrieval. Global Navigation Satellite System(GNSS) PWV retrieval is limited because the GNSS stations lack meteorological sensors. First, this article evaluated the accuracy of pressure and temperature in 68 radiosonde stations in China based on ERA5 Reanalysis data from 2015 to 2019 and compared them with GPT3model. Then, the accuracy of pressure and temperature calculated by ERA5 were estimated in 5 representative IGS stations in China. And the PWV calculated by these meteorological parameters from ERA5(ERA5-PWV) were analyzed. Finally, the relation between ERA5-PWV and precipitation was deeply explored using wavelet coherence analysis in IGS stations. These results indicate that the accuracy of pressure and temperature of ERA5 is better than the GPT3 model. In radiosonde stations, the mean BIAS and MAE of pressure and temperature in ERA5 are-0.41/1.15 hpa and-0.97/2.12 K. And the mean RMSEs are 1.35 hpa and 2.87 K, which improve 74.77% and 40.58% compared with GPT3 model. The errors of pressure and temperature of ERA5 are smaller than the GPT3 model in bjfs, hksl and wuh2, and the accuracy of ERA5-PWV is improved by 18.77% compared with the GPT3 model. In addition, there is a significant positive correlation between ERA5-PWV and precipitation. And precipitation is always associated with the sharp rise of ERA5-PWV, which provides important references for rainfall prediction.

A Web-Based Approach for the Efficient Management of Massive Multi-source 3D Models

Effectively managing extensive,multi-source,and multi-level real-scene 3D models for responsive retrieval scheduling and rapid visualization in the Web environment is a significant challenge in the current development of real-scene 3D applications in China.In this paper,we address this challenge by reorganizing spatial and temporal information into a 3D geospatial grid.It introduces the Global 3D Geocoding System(G_(3)DGS),leveraging neighborhood similarity and uniqueness for efficient storage,retrieval,updating,and scheduling of these models.A combination of G_(3)DGS and non-relational databases is implemented,enhancing data storage scalability and flexibility.Additionally,a model detail management scheduling strategy(TLOD)based on G_(3)DGS and an importance factor T is designed.Compared with mainstream commercial and open-source platforms,this method significantly enhances the loadable capacity of massive multi-source real-scene 3D models in the Web environment by 33%,improves browsing efficiency by 48%,and accelerates invocation speed by 40%.

Numerical simulations of earthquake rupture induced by pressure perturbation

The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.

Error Analysis and Accuracy Assessment of GPS Absolute Velocity Determination without SA

Error sources which decrease the accuracy of GPS in absolute velocity determination have been changed since SA was turned off. Firstly, quantities of all kinds of error sources that influence velocity deter-mination are analyzed. The potential accuracy of GPS absolute velocity determination is derived from both theory and field GPS data simulation. After that, two tests were carried out to evaluate the performance of GPS absolute velocity determination in the case of a static and an airborne GPS receiver and INS (Inertial Navigation System) instrument in kinematic mode. In static mode, the receiver velocity has been estimated to be several mm/s with the carrier-phase derived Doppler measurements, and several cm/s with the receiver generated Doppler measurements. In kinematic mode, GPS absolute velocity estimates are compared with the synchronized measurements from the high accuracy INS. The root mean square statistics of the velocity discrepancies between GPS and INS come up to dm/s. Moreover, it has a strong correlation with the accel-eration or jerk of the aircraft.

Earth＇s Temporal Principal Moments of Inertia and Variable Rotation

Based on the gravity field models EGM96 and EIGEN-GL04C, the Earth＇s time-dependent principal moments of inertia A, B, C are obtained, and the variable rotation of the Earth is determined. Numerical results show that A, B, and C have increasing tendencies; the tilt of the rotation axis increases 2.1×10^ 8 mas/yr; the third component of the rotational angular velocity, ω3 , has a decrease of 1.0×10^ 22 rad/s^2, which is around 23% of the present observed value. Studies show in detail that both 0 and ω3 experience complex fluctuations at various time scales due to the variations of A, B and C.

Measurements of ocean wave and current field using dual polarized X-band radar

A new ocean wave and sea surface current monitoring system with horizontally-(HH) and vertically-(VV) polarized X-band radar was developed.Two experiments into the use of the radar system were carried out at two sites,respectively,for calibration process in Zhangzi Island of the Yellow Sea,and for validation in the Yellow Sea and South China Sea.Ocean wave parameters and sea surface current velocities were retrieved from the dual polarized radar image sequences based on an inverse method.The results obtained from dual-polarized radar data sets acquired in Zhangzi Island are compared with those from an ocean directional buoy.The results show that ocean wave parameters and sea surface current velocities retrieved from radar image sets are in a good agreement with those observed by the buoy.In particular,it has been found that the vertically-polarized radar is better than the horizontally-polarized radar in retrieving ocean wave parameters,especially in detecting the significant wave height below 1.0 m.

Crust and upper mantle structure of East Asia from ambient noise and earthquake surface wave tomography

The complex tectonic background of East Asia makes it an ideal region for investigating the evolution of the continental lithosphere,for which high-resolution lithospheric structural models are essential.In this study,we measured Rayleigh-wave phase-velocity dispersion curves at periods of 10-120 s and group velocity dispersion curves at periods of 10-140 s using event records from more than 1,000 seismic stations in and around China.By jointly inverting new and previously published dispersion data from ambient noise and earthquakes,we developed a high-resolution shear-wave velocity model down to a depth of~300 km beneath East Asia.Our model revealed heterogeneous lithospheric structures beneath East Asia,and allowed us to investigate the velocity structure of the entire lithosphere.We also derived crustal and lithospheric thickness models from the three-dimensional(3D)shear-wave model,revealing strong spatial heterogeneity and a general thinning trend of lithospheric thickness from west to east across the study region.Overall,our models reveal important lithospheric features beneath East Asia and provide a valuable baseline dataset for understanding continental-scale dynamics and evolution.