Nonlinear variations in the coordinate time series of global navigation satellite system(GNSS) reference stations are strongly correlated with surface displacements caused by environmental loading effects,including at...Nonlinear variations in the coordinate time series of global navigation satellite system(GNSS) reference stations are strongly correlated with surface displacements caused by environmental loading effects,including atmospheric, hydrological, and nontidal ocean loading. Continuous improvements in the accuracy of surface mass loading products, performance of Earth models, and precise data-processing technologies have significantly advanced research on the effects of environmental loading on nonlinear variations in GNSS coordinate time series. However, owing to theoretical limitations, the lack of high spatiotemporal resolution surface mass observations, and the coupling of GNSS technology-related systematic errors, environmental loading and nonlinear GNSS reference station displacements remain inconsistent. The applicability and capability of these loading products across different regions also require further evaluation. This paper outlines methods for modeling environmental loading, surface mass loading products, and service organizations. In addition, it summarizes recent advances in applying environmental loading to address nonlinear variations in global and regional GNSS coordinate time series. Moreover, the scientific questions of existing studies are summarized, and insights into future research directions are provided. The complex nonlinear motion of reference stations is a major factor limiting the accuracy of the current terrestrial reference frame. Further refining the environmental load modeling method, establishing a surface mass distribution model with high spatiotemporal resolution and reliability, exploring other environmental load factors such as ice sheet and artificial mass-change effects, and developing an optimal data-processing model and strategy for reprocessing global reference station data consistently could contribute to the development of a millimeter-level nonlinear motion model for GNSS reference stations with actual physical significance and provide theoretical support for establishing a terrestrial reference frame with 1 mm accuracy by 2050.展开更多
Measurements from geomagnetic satellites continue to underpin advances in geomagnetic field models that describe Earth's internally generated magnetic field.Here,we present a new field model,MSCM,that integrates v...Measurements from geomagnetic satellites continue to underpin advances in geomagnetic field models that describe Earth's internally generated magnetic field.Here,we present a new field model,MSCM,that integrates vector and scalar data from the Swarm,China Seismo-Electromagnetic Satellite(CSES),and Macao Science Satellite-1(MSS-1)missions.The model spans from 2014.0 to 2024.5,incorporating the core,lithospheric,and magnetospheric fields,and it shows characteristics similar to other published models based on different data.For the first time,we demonstrate that it is possible to successfully construct a geomagnetic field model that incorporates CSES vector data,albeit one in which the radial and azimuthal CSES vector components are Huber downweighted.We further show that data from the MSS-1 can be integrated within an explicitly smoothed,fully time-dependent model description.Using the MSCM,we identify new behavior of the South Atlantic Anomaly,the broad region of low magnetic field intensity over the southern Atlantic.This prominent feature appears split into a western part and an eastern part,each with its own intensity minimum.Since 2015,the principal western minimum has undergone only modest intensity decreases of 290 nT and westward motion of 20 km per year,whereas the recently formed eastern minimum has shown a 2–3 times greater intensity drop of 730 nT with no apparent east-west motion.展开更多
The mathematical method of ZTD(zenith tropospheric delay)spatial prediction is important for precise ZTD derivation and real-time precise point positioning(PPP)augmentation.This paper analyses the performance of the p...The mathematical method of ZTD(zenith tropospheric delay)spatial prediction is important for precise ZTD derivation and real-time precise point positioning(PPP)augmentation.This paper analyses the performance of the popular optimal function coefficient(OFC),sphere cap harmonic analysis(SCHA),kriging and inverse distance weighting(IDW)interpolation in ZTD spatial prediction and Beidou satellite navigation system(BDS)-PPP augmentation over China.For ZTD spatial prediction,the average time consumption of the OFC,kriging,and IDW methods is less than 0.1 s,which is significantly better than that of the SCHA method(63.157 s).The overall ZTD precision of the OFC is 3.44 cm,which outperforms those of the SCHA(9.65 cm),Kriging(10.6 cm),and IDW(11.8 cm)methods.We confirmed that the low performance of kriging and IDW is caused by their weakness in modelling ZTD variation in the vertical direction.To mitigate such deficiencies,an elevation normalization factor(ENF)is introduced into the kriging and IDW models(kriging-ENF and IDW-ENF).The overall ZTD spatial prediction accuracies of IDW-ENF and kriging-ENF are 2.80 cm and 2.01 cm,respectively,which are both superior to those of the OFC and the widely used empirical model GPT3(4.92 cm).For BDS-PPP enhancement,the ZTD provided by the kriging-ENF,IDW-ENF and OFC as prior constraints can effectively reduce the convergence time.Compared with unconstrained BDS-PPP,our proposed kriging-ENF outperforms IDW-ENF and OFC by reducing the horizontal and vertical convergence times by approximately 13.2%and 5.8%in Ningxia and 30.4%and 7.84%in Guangdong,respectively.These results indicate that kriging-ENF is a promising method for ZTD spatial prediction and BDS-PPP enhancement over China.展开更多
The Main Himalayan Thrust(MHT),where the 2015 MW7.8 Gorkha earthquake occurred,features the most seismicity of any structure in Nepal.The structural complexity of the MHT makes it difficult to obtain a definitive inte...The Main Himalayan Thrust(MHT),where the 2015 MW7.8 Gorkha earthquake occurred,features the most seismicity of any structure in Nepal.The structural complexity of the MHT makes it difficult to obtain a definitive interpretation of deep seismogenic structures.The application of new methods and data in this region is necessary to enhance local seismic hazard analyses.In this study,we used a well-designed machine learning-based earthquake location workflow(LOC-FLOW),which incorporates machine learning phase picking,phase association,absolute location,and double-difference relative location,to process seismic data collected by the Hi-CLIMB and NAMASTE seismic networks.We built a high-precision earthquake catalog of both the quiet-period and aftershock seismicity in this region.The seismicity distribution suggests that the quietperiod seismicity(388 events)was controlled by a mid-crustal ramp and the aftershock seismicity(12,669 events)was controlled by several geological structures of the MHT.The higher-level detail of the catalogs derived from this machine learning method reveal clearer structural characteristics,showing how the flat-ramp geometry and a possible duplex structure affect the depth distribution of the seismic events,and how a tear fault changes this distribution along strike.展开更多
Since April 2002,the Gravity Recovery and Climate Experiment Satellite(GRACE)has provided monthly total water storage anomalies(TWSAs)on a global scale.However,these TWSAs are discontinuous because some GRACE observat...Since April 2002,the Gravity Recovery and Climate Experiment Satellite(GRACE)has provided monthly total water storage anomalies(TWSAs)on a global scale.However,these TWSAs are discontinuous because some GRACE observation data are missing.This study presents a combined machine learning-based modeling algorithm without hydrological model data.The TWSA time-series data for 11 large regions worldwide were divided into training and test sets.Autoregressive integrated moving average(ARIMA),long short-term memory(LSTM),and an ARIMA-LSTM combined model were used.The model predictions were compared with GRACE observations,and the model accuracy was evaluated using fi ve metrics:the Nash-Sutcliff e effi ciency coeffi cient(NSE),Pearson correlation coeffi cient(CC),root mean square error(RMSE),normalized RMSE(NRMSE),and mean absolute percentage error.The results show that at the basin scale,the mean CC,NSE,and NRMSE for the ARIMA-LSTM model were 0.93,0.83,and 0.12,respectively.At the grid scale,this study compared the spatial distribution and cumulative distribution function curves of the metrics in the Amazon and Volga River basins.The ARIMA-LSTM model had mean CC and NSE values of 0.89 and 0.61 and 0.92 and 0.61 in the Amazon and Volga River basins,respectively,which are superior to those of the ARIMA model(0.86 and 0.48 and 0.88 and 0.46,respectively)and the LSTM model(0.80 and 0.41 and 0.89 and 0.31,respectively).In the ARIMA-LSTM model,the proportions of grid cells with NSE>0.50 for the two basins were 63.3%and 80.8%,while they were 54.3%and 51.3%in the ARIMA model and 53.7%and 43.2%in the LSTM model.The ARIMA-LSTM model significantly improved the NSE values of the predictions while guaranteeing high CC values in the GRACE data reconstruction at both scales,which can aid in fi lling in discontinuous data in temporal gravity fi eld models..展开更多
Long-period free oscillations provide robust information for the spatio-temporal characteristics of large earthquakes.In this study,we employ a normal-mode summation algorithm to generate threecomponent seismograms wi...Long-period free oscillations provide robust information for the spatio-temporal characteristics of large earthquakes.In this study,we employ a normal-mode summation algorithm to generate threecomponent seismograms within an aspherical,anelastic,and rotating Earth model,focusing on the excitation of seismic normal modes by the 2011 Tohoku megathrust earthquake.Specifically,we analyze the effects of seismic source parameters,including fault geometry,focal depth,and rupture duration.By comparing synthetic free oscillation spectra with observed data,we validate several source mechanisms and emphasize the need for horizontal observations to improve seismic mechanism inversions.Our quantitative analyses reveal that among fault geometry parameters,dip and slip angles have a more pronounced impact on excitation amplitudes than fault strikes.Certain fault configurations enhance the detectability of specific modes,underscoring the relationship between fault geometry and mode excitation.Normal modes also exhibit varying sensitivity to focal depth,with significant excitation amplitude changes at discontinuity boundaries.Additionally,we demonstrate that while rupture duration can be inferred by minimizing differences between observed and synthetic spectra,more extensive records and modes should be included.展开更多
As a major fault in the northeastern Qinghai-Xizang Plateau,the Haiyuan fault zone is important for understanding the regional deformation.Aiming at the differences in the slip rate and locking degree obtained from di...As a major fault in the northeastern Qinghai-Xizang Plateau,the Haiyuan fault zone is important for understanding the regional deformation.Aiming at the differences in the slip rate and locking degree obtained from different studies,this study constructs a refined block model(including Qilian,Alxa,Ordos,Xining,Haiyuan,and Lanzhou blocks)and uses the grid search and simulated annealing methods to invert GPS data for slip rate and locking degree of the Haiyuan fault zone.The results are as follows:(1)The sinistral slip rates in the western,middle,and eastern segments are 4.93-5.22 mm/a,1.52-4.94 mm/a,and 0.43-1.18 mm/a,decreasing eastward on the whole,while the compression rates are 0.45-1.26 mm/a,0.58-2.62 mm/a,and3.52-4.48 mm/a,increasing eastward on the whole.(2)The locking depth of the western segment increases from about 5 km to about 20 km eastward;the middle segment decreases and then increases eastward;the eastern segment concentrates at about 20 km(PHI is about 0.86).(3)The slip deficit is relatively higher in the Lenglongling,Jinqianghe,Maomaoshan,and Liupanshan faults(averaging about 3.42 mm/a,4.16 mm/a,4.23 mm/a,and 3.43 mm/a within 20 km).(4)The Qilian,Alxa,Xining,Lanzhou,and Haiyuan blocks rotate clockwise,while the Ordos block rotates counterclockwise.Additionally,by comparing different block models,the Haiyuan block should be considered independently.The Haiyuan fault zone adjusts surrounding block movements and uplifts Liupanshan mountain tectonically.The results can provide important references for understanding the regional earthquake risk and deformation mechanism.展开更多
The early excessive groundwater exploitation caused extensive land subsidence in the Hangjiahu Plain(HJHP),threatening geological stability and infrastructure.Government restrictions on groundwater exploitation led to...The early excessive groundwater exploitation caused extensive land subsidence in the Hangjiahu Plain(HJHP),threatening geological stability and infrastructure.Government restrictions on groundwater exploitation led to considerable changes in surface subsidence.To explore the distribution and causes of surface deformation after restrictions on groundwater exploitation,we used 46 Sentinel-1A TOPS images to detect the spatiotemporal distribution of surface deformation in the HJHP area from January 2019 to October 2022 by SBAS-InSAR technique and verified it with leveling results.Results indicated the following:(1)The uneven surface subsidence in the HJHP was relatively apparent,with the localized subsidence concentrated in a few areas.During InSAR monitoring,the maximum subsidence velocity reached-21.0 mm/yr.The most severe subsidence area was in the Shengze town industrial park,with remarkable subsidence bowls and a maximum subsidence of-85.8 mm.(2)Results were verified by comparing the InSAR results with leveling,and the results were in good agreement,indicating the reliability of the InSAR results.(3)Comprehensive analysis of the spatial distribution characteristics of surface subsidence,geological data and rainfall changes revealed that fault activity and rainfall changes had varying degrees of effect on surface subsidence in the HJHP area.(4)The comparison with land use data indicated that the built area was the dominant mechanism of surface subsidence in the study area.In addition,the effect of industrial activities on the surface cannot be ignored.The study results provide scientific data support for emergency disaster prevention work in the HJHP.展开更多
The time-varying periodic variations in Global Navigation Satellite System(GNSS)stations affect the reliable time series analysis and appropriate geophysical interpretation.In this study,we apply the singular spectrum...The time-varying periodic variations in Global Navigation Satellite System(GNSS)stations affect the reliable time series analysis and appropriate geophysical interpretation.In this study,we apply the singular spectrum analysis(SSA)method to characterize and interpret the periodic patterns of GNSS deformations in China using multiple geodetic datasets.These include 23-year observations from the Crustal Movement Observation Network of China(CMONOC),displacements inferred from the Gravity Recovery and Climate Experiment(GRACE),and loadings derived from Geophysical models(GM).The results reveal that all CMONOC time series exhibit seasonal signals characterized by amplitude and phase modulations,and the SSA method outperforms the traditional least squares fitting(LSF)method in extracting and interpreting the time-varying seasonal signals from the original time series.The decrease in the root mean square(RMS)correlates well with the annual cycle variance estimated by the SSA method,and the average reduction in noise amplitudes is nearly twice as much for SSA filtered results compared with those from the LSF method.With SSA analysis,the time-varying seasonal signals for all the selected stations can be identified in the reconstructed components corresponding to the first ten eigenvalues.Moreover,both RMS reduction and correlation analysis imply the advantages of GRACE solutions in explaining the GNSS periodic variations,and the geophysical effects can account for 71%of the GNSS annual amplitudes,and the average RMS reduction is 15%.The SSA method has proved to be useful for investigating the GNSS timevarying seasonal signals.It could be applicable as an auxiliary tool in the improvement of nonlinear variations investigations.展开更多
This paper introduces the comparison of absolute gravity measurements in China during 2006-2008,and analyzing the survey and comparing the results,it shows that there is no obvious system error between the gravimeter ...This paper introduces the comparison of absolute gravity measurements in China during 2006-2008,and analyzing the survey and comparing the results,it shows that there is no obvious system error between the gravimeter of FG5/214 and FG5/232,the surveying accuracy is very high and repetition is good,and their inner surveying accuracy is about 2-3 microGal.展开更多
Gravity Recovery and Climate Experiment（GRACE） observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 2...Gravity Recovery and Climate Experiment（GRACE） observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 26,2004.This article focuses on investigating some gravita-tional effects caused by this huge earthquake.We computed the geoid height changes,the equivalent water height（EWH） changes,and the gravity changes using the GRACE Level-2 monthly spherical harmonic（SH） solutions released by University of Texas Center for Space Research（UTCSR）.The GRACE results agree well with the prediction by a dislocation model and are consistent with the results obtained by some previous scholars.In particular,we calculated the three components of the gravity gradient variations and found that they can recover the seismic-related signature more sensitively due to a certain degree of amplification of the signals.A positive-negative-positive mode predominates in the spatial distribution of the horizontal components of the gravity gradient variations,which is possibly attributed to the anomalies in the crustal density distribution caused by the uplift-subduction effect of the dip-slip earthquake.Moreover,the latitude components of the gravity gradient changes show strong suppression of the north-south stripes,which is due to the along-orbit measurements of the two GRACE satellites.We conclude that the posi-tive-negative-positive mode in latitude gravity gradient changes would be a more sensitive fea-ture to detect the deformations of some major dip-slip earthquakes by GRACE data.展开更多
With the undergoing and planned implementations of mega constellations of thousands of Low Earth Orbiting(LEO)satellites,space will become even more congested for satellite operations.The enduring effects on the long-...With the undergoing and planned implementations of mega constellations of thousands of Low Earth Orbiting(LEO)satellites,space will become even more congested for satellite operations.The enduring effects on the long-term space environment have been investigated by various researchers using debris environment models.This paper is focused on the imminent short-term effects of LEO mega constellations on the space operation environment concerned by satellite owners and operators.The effects are measured in terms of the Close Approaches(CAs)and overall collision probability.Instead of using debris environment models,the CAs are determined from integrated orbit positions,and the collision probability is computed for each CA considering the sizes and position covariance of the involving objects.The obtained results thus present a clearer picture of the space operation safety environment when LEO mega constellations are deployed.Many mega constellations are simulated,including a Starlink-like constellation of 1584 satellites,four possible generic constellations at altitudes between 1110 km and 1325 km,and three constellations of 1584 satellites each at the altitudes of 650 km,800 km,and 950 km,respectively,where the Resident Space Object(RSO)spatial density is the highest.The increases in the number of CAs and overall collision probability caused by them are really alarming.The results suggest that highly frequent orbital maneuvers are required to avoid collisions between existing RSOs and constellation satellites,and between satellites from two constellations at a close altitude,as such the constellation operation burden would be very heavy.The study is not only useful for satellite operators but a powerful signal for various stakeholders to pay serious attention to the development of LEO mega constellations.展开更多
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 a...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.展开更多
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...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.展开更多
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...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.展开更多
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...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.展开更多
Urban faults in Shenzhen are potential threats to city security and sustainable development. In consideration of the importance of the Shenzhen fault zone, the author provide a detailed interpretation on gravity data ...Urban faults in Shenzhen are potential threats to city security and sustainable development. In consideration of the importance of the Shenzhen fault zone, the author provide a detailed interpretation on gravity data model. Bouguer gravity covering the whole Shenzhen City was calculated with a 1-km resolution. Wavelet multi-scale analysis（MSA） was applied to the Bouguer gravity data to obtain the multilayer residual anomalies corresponding to different depths. In addition, 2D gravity models were constructed along three profiles. The Bouguer gravity anomaly shows an NE-striking high-low-high pattern from northwest to southeast, strongly related to the main faults. According to the results of MSA, the correlation between gravity anomaly and faults is particularly significant from 4 to 12 km depth. The residual gravity with small amplitude in each layer indicates weak tectonic activity in the crust. In the upper layers, positive anomalies along most of faults reveal the upwelling of high-density materials during the past tectonic movements. The multilayer residual anomalies also yield important information about the faults, such as the vertical extension and the dip direction. The maximum depth of the faults is about 20 km. In general, NE-striking faults extend deeper than NW-striking faults and have a larger dip angle.展开更多
Both the Global Positioning System(GPS)and Gravity Recovery and Climate Experiment(GRACE)/GRACE Follow-On(GFO)provide effective tools to infer surface mass changes.In this paper,we combined GPS,GRACE/GFO spherical har...Both the Global Positioning System(GPS)and Gravity Recovery and Climate Experiment(GRACE)/GRACE Follow-On(GFO)provide effective tools to infer surface mass changes.In this paper,we combined GPS,GRACE/GFO spherical harmonic(SH)solutions and GRACE/GFO mascon solutions to analyze the total surface mass changes and terrestrial water storage(TWS)changes in the Shaan-Gan-Ning Region(SGNR)over the period from December 2010 to February 2021.To improve the reliability of GPS inversion results,an improved regularization Laplace matrix and monthly optimal regularization parameter estimation strategy were employed to solve the ill-posed problem.The results show that the improved Laplace matrix can suppress the edge effects better than that of the traditional Laplace matrix,and the corre-lation coefficient and standard deviation(STD)between the original signal and inversion results from the traditional and improved Laplace matrix are 0.84 and 0.88,and 17.49 mm and 15.16 mm,respectively.The spatial distributions of annual amplitudes and time series changes for total surface mass changes derived from GPS agree well with GRACE/GFO SH solutions and mascon solutions,and the correlation coefficients of total surface mass change time series between GPS and GRACE/GFO SH solutions,GPS and GRACE/GFO mascon solutions are 0.80 and 0.77.However,the obvious differences still exist in local regions.In addition,the seasonal characteristics,increasing and decreasing rate of TWS change time series derived from GPS,GRACE/GFO SH and mascon solutions agree well with the Global Land Data Assimilation System(GLDAS)hydrological model in the studied area,and generally consistent with the precipitation data.Meanwhile,TWS changes derived from GPS and GRACE mascon solutions in the SGNR are more reliable than those of GRACE SH solutions over the period from January 2016 to June 2017(the final operation phase of the GRACE mission).展开更多
On November 13, 2016, an MW7.8 earthquake struck Kaikoura in South Island of New Zealand. By means of back-projection of array recordings, ASTFs-analysis of global seismic recordings, and joint inversion of global sei...On November 13, 2016, an MW7.8 earthquake struck Kaikoura in South Island of New Zealand. By means of back-projection of array recordings, ASTFs-analysis of global seismic recordings, and joint inversion of global seismic data and co-seismic In SAR data, we investigated complexity of the earthquake source. The result shows that the 2016 MW7.8 Kaikoura earthquake ruptured about 100 s unilaterally from south to northeast(~N28°–33°E), producing a rupture area about 160 km long and about 50 km wide and releasing scalar moment 1.01×1021 Nm. In particular, the rupture area consisted of two slip asperities, with one close to the initial rupture point having a maximal slip value ~6.9 m while the other far away in the northeast having a maximal slip value ~9.3 m. The first asperity slipped for about 65 s and the second one started 40 s after the first one had initiated. The two slipped simultaneously for about 25 s.Furthermore, the first had a nearly thrust slip while the second had both thrust and strike slip. It is interesting that the rupture velocity was not constant, and the whole process may be divided into 5 stages in which the velocities were estimated to be 1.4 km/s, 0 km/s, 2.1 km/s, 0 km/s and 1.1 km/s, respectively. The high-frequency sources distributed nearly along the lower edge of the rupture area, the highfrequency radiating mainly occurred at launching of the asperities, and it seemed that no high-frequency energy was radiated when the rupturing was going to stop.展开更多
A GIS for ocean applications called ＂the Xiamen Atmospheric and Oceanographic Data Management and Display System （AODMDS）＂ has been designed and developed. The system is based on ArcObjects （AO）, a component-bas...A GIS for ocean applications called ＂the Xiamen Atmospheric and Oceanographic Data Management and Display System （AODMDS）＂ has been designed and developed. The system is based on ArcObjects （AO）, a component-based GIS de- velopment tool. The paper discusses in detail the storage and organization of the atmospheric and oceanographic data, the strategy and methods for the visualization and mapping of oceanographic and atmospheric data, and the implementation of the methods in AODMDS. It also discusses some advanced display control techniques that expand the functions of ArcObjects One of the techniques is ＂gradient-fill-style color-map control,＂ which provides a feasible color-rich display control for all types of raster maps. As a stand-alone desktop GIS system built on AO, AODMDS provides effective data management and powerful mapping and visualization functions for atmospheric and oceanographic data.展开更多
基金supported by the Basic Science Center Project of the National Natural Science Foundation of China(42388102)the National Natural Science Foundation of China(42174030)+2 种基金the Special Fund of Hubei Luojia Laboratory(220100020)the Major Science and Technology Program for Hubei Province(2022AAA002)the Fundamental Research Funds for the Central Universities of China(2042022dx0001 and 2042023kfyq01)。
文摘Nonlinear variations in the coordinate time series of global navigation satellite system(GNSS) reference stations are strongly correlated with surface displacements caused by environmental loading effects,including atmospheric, hydrological, and nontidal ocean loading. Continuous improvements in the accuracy of surface mass loading products, performance of Earth models, and precise data-processing technologies have significantly advanced research on the effects of environmental loading on nonlinear variations in GNSS coordinate time series. However, owing to theoretical limitations, the lack of high spatiotemporal resolution surface mass observations, and the coupling of GNSS technology-related systematic errors, environmental loading and nonlinear GNSS reference station displacements remain inconsistent. The applicability and capability of these loading products across different regions also require further evaluation. This paper outlines methods for modeling environmental loading, surface mass loading products, and service organizations. In addition, it summarizes recent advances in applying environmental loading to address nonlinear variations in global and regional GNSS coordinate time series. Moreover, the scientific questions of existing studies are summarized, and insights into future research directions are provided. The complex nonlinear motion of reference stations is a major factor limiting the accuracy of the current terrestrial reference frame. Further refining the environmental load modeling method, establishing a surface mass distribution model with high spatiotemporal resolution and reliability, exploring other environmental load factors such as ice sheet and artificial mass-change effects, and developing an optimal data-processing model and strategy for reprocessing global reference station data consistently could contribute to the development of a millimeter-level nonlinear motion model for GNSS reference stations with actual physical significance and provide theoretical support for establishing a terrestrial reference frame with 1 mm accuracy by 2050.
基金supported by the National Natural Science Foundation of China(Grant No.42274003)PWL was supported by Swarm DISC(Swarm Data,Innovation,and Science Cluster)+2 种基金funded by the European Space Agency(ESAContract No.4000109587)HFR acknowledges funding from the UK Natural Environment Research Council(Grant No.NE/V010867/1)。
文摘Measurements from geomagnetic satellites continue to underpin advances in geomagnetic field models that describe Earth's internally generated magnetic field.Here,we present a new field model,MSCM,that integrates vector and scalar data from the Swarm,China Seismo-Electromagnetic Satellite(CSES),and Macao Science Satellite-1(MSS-1)missions.The model spans from 2014.0 to 2024.5,incorporating the core,lithospheric,and magnetospheric fields,and it shows characteristics similar to other published models based on different data.For the first time,we demonstrate that it is possible to successfully construct a geomagnetic field model that incorporates CSES vector data,albeit one in which the radial and azimuthal CSES vector components are Huber downweighted.We further show that data from the MSS-1 can be integrated within an explicitly smoothed,fully time-dependent model description.Using the MSCM,we identify new behavior of the South Atlantic Anomaly,the broad region of low magnetic field intensity over the southern Atlantic.This prominent feature appears split into a western part and an eastern part,each with its own intensity minimum.Since 2015,the principal western minimum has undergone only modest intensity decreases of 290 nT and westward motion of 20 km per year,whereas the recently formed eastern minimum has shown a 2–3 times greater intensity drop of 730 nT with no apparent east-west motion.
基金co-supported by the National Nature Science Foundation of China(No.12303071)the Shanghai Science and Technology Plan Project,China(No.23YF1455500)+1 种基金the China Postdoctoral Science Foundation(No.2023M743653)Ministry of Industry and Information Technology of China through the High Precision Timing Service Project(No.TC220A04A-80)。
文摘The mathematical method of ZTD(zenith tropospheric delay)spatial prediction is important for precise ZTD derivation and real-time precise point positioning(PPP)augmentation.This paper analyses the performance of the popular optimal function coefficient(OFC),sphere cap harmonic analysis(SCHA),kriging and inverse distance weighting(IDW)interpolation in ZTD spatial prediction and Beidou satellite navigation system(BDS)-PPP augmentation over China.For ZTD spatial prediction,the average time consumption of the OFC,kriging,and IDW methods is less than 0.1 s,which is significantly better than that of the SCHA method(63.157 s).The overall ZTD precision of the OFC is 3.44 cm,which outperforms those of the SCHA(9.65 cm),Kriging(10.6 cm),and IDW(11.8 cm)methods.We confirmed that the low performance of kriging and IDW is caused by their weakness in modelling ZTD variation in the vertical direction.To mitigate such deficiencies,an elevation normalization factor(ENF)is introduced into the kriging and IDW models(kriging-ENF and IDW-ENF).The overall ZTD spatial prediction accuracies of IDW-ENF and kriging-ENF are 2.80 cm and 2.01 cm,respectively,which are both superior to those of the OFC and the widely used empirical model GPT3(4.92 cm).For BDS-PPP enhancement,the ZTD provided by the kriging-ENF,IDW-ENF and OFC as prior constraints can effectively reduce the convergence time.Compared with unconstrained BDS-PPP,our proposed kriging-ENF outperforms IDW-ENF and OFC by reducing the horizontal and vertical convergence times by approximately 13.2%and 5.8%in Ningxia and 30.4%and 7.84%in Guangdong,respectively.These results indicate that kriging-ENF is a promising method for ZTD spatial prediction and BDS-PPP enhancement over China.
基金funded by the National Key R&D Program of China(2022YFF0800601)National Natural Science Foundation of China(42174069,U1939204).
文摘The Main Himalayan Thrust(MHT),where the 2015 MW7.8 Gorkha earthquake occurred,features the most seismicity of any structure in Nepal.The structural complexity of the MHT makes it difficult to obtain a definitive interpretation of deep seismogenic structures.The application of new methods and data in this region is necessary to enhance local seismic hazard analyses.In this study,we used a well-designed machine learning-based earthquake location workflow(LOC-FLOW),which incorporates machine learning phase picking,phase association,absolute location,and double-difference relative location,to process seismic data collected by the Hi-CLIMB and NAMASTE seismic networks.We built a high-precision earthquake catalog of both the quiet-period and aftershock seismicity in this region.The seismicity distribution suggests that the quietperiod seismicity(388 events)was controlled by a mid-crustal ramp and the aftershock seismicity(12,669 events)was controlled by several geological structures of the MHT.The higher-level detail of the catalogs derived from this machine learning method reveal clearer structural characteristics,showing how the flat-ramp geometry and a possible duplex structure affect the depth distribution of the seismic events,and how a tear fault changes this distribution along strike.
基金financially supported by The National Natural Science Foundation of China (42374004)the Open Fund of Hubei Luojia Laboratory (220100045)the Natural Science Foundation of Sichuan Province (2022NSFSC1047)。
文摘Since April 2002,the Gravity Recovery and Climate Experiment Satellite(GRACE)has provided monthly total water storage anomalies(TWSAs)on a global scale.However,these TWSAs are discontinuous because some GRACE observation data are missing.This study presents a combined machine learning-based modeling algorithm without hydrological model data.The TWSA time-series data for 11 large regions worldwide were divided into training and test sets.Autoregressive integrated moving average(ARIMA),long short-term memory(LSTM),and an ARIMA-LSTM combined model were used.The model predictions were compared with GRACE observations,and the model accuracy was evaluated using fi ve metrics:the Nash-Sutcliff e effi ciency coeffi cient(NSE),Pearson correlation coeffi cient(CC),root mean square error(RMSE),normalized RMSE(NRMSE),and mean absolute percentage error.The results show that at the basin scale,the mean CC,NSE,and NRMSE for the ARIMA-LSTM model were 0.93,0.83,and 0.12,respectively.At the grid scale,this study compared the spatial distribution and cumulative distribution function curves of the metrics in the Amazon and Volga River basins.The ARIMA-LSTM model had mean CC and NSE values of 0.89 and 0.61 and 0.92 and 0.61 in the Amazon and Volga River basins,respectively,which are superior to those of the ARIMA model(0.86 and 0.48 and 0.88 and 0.46,respectively)and the LSTM model(0.80 and 0.41 and 0.89 and 0.31,respectively).In the ARIMA-LSTM model,the proportions of grid cells with NSE>0.50 for the two basins were 63.3%and 80.8%,while they were 54.3%and 51.3%in the ARIMA model and 53.7%and 43.2%in the LSTM model.The ARIMA-LSTM model significantly improved the NSE values of the predictions while guaranteeing high CC values in the GRACE data reconstruction at both scales,which can aid in fi lling in discontinuous data in temporal gravity fi eld models..
基金supported by the National Natural ScienceFoundation of China(42204003,42274011,42388102,42192533,42192531)the Natural Science Foundation of Wuhan(2024040701010027)+1 种基金the Open Fund Supported by the State KeyLaboratory of Precision Geodesy,Innovation Academy for PrecisionMeasurement Science and Technology,Chinese Academy of Sci-ences(SKLGED2024-1-1)the Open Fund Supported by KeyLaboratory of Polar Environment Monitoring and Public Gover-nance(Wuhan University),Ministry of Education(202401)。
文摘Long-period free oscillations provide robust information for the spatio-temporal characteristics of large earthquakes.In this study,we employ a normal-mode summation algorithm to generate threecomponent seismograms within an aspherical,anelastic,and rotating Earth model,focusing on the excitation of seismic normal modes by the 2011 Tohoku megathrust earthquake.Specifically,we analyze the effects of seismic source parameters,including fault geometry,focal depth,and rupture duration.By comparing synthetic free oscillation spectra with observed data,we validate several source mechanisms and emphasize the need for horizontal observations to improve seismic mechanism inversions.Our quantitative analyses reveal that among fault geometry parameters,dip and slip angles have a more pronounced impact on excitation amplitudes than fault strikes.Certain fault configurations enhance the detectability of specific modes,underscoring the relationship between fault geometry and mode excitation.Normal modes also exhibit varying sensitivity to focal depth,with significant excitation amplitude changes at discontinuity boundaries.Additionally,we demonstrate that while rupture duration can be inferred by minimizing differences between observed and synthetic spectra,more extensive records and modes should be included.
基金supported by the National Natural Science Foundation of China(42474003,42074007)the Fundamental Research Funds for the Central Universities(2042023kfyq01)。
文摘As a major fault in the northeastern Qinghai-Xizang Plateau,the Haiyuan fault zone is important for understanding the regional deformation.Aiming at the differences in the slip rate and locking degree obtained from different studies,this study constructs a refined block model(including Qilian,Alxa,Ordos,Xining,Haiyuan,and Lanzhou blocks)and uses the grid search and simulated annealing methods to invert GPS data for slip rate and locking degree of the Haiyuan fault zone.The results are as follows:(1)The sinistral slip rates in the western,middle,and eastern segments are 4.93-5.22 mm/a,1.52-4.94 mm/a,and 0.43-1.18 mm/a,decreasing eastward on the whole,while the compression rates are 0.45-1.26 mm/a,0.58-2.62 mm/a,and3.52-4.48 mm/a,increasing eastward on the whole.(2)The locking depth of the western segment increases from about 5 km to about 20 km eastward;the middle segment decreases and then increases eastward;the eastern segment concentrates at about 20 km(PHI is about 0.86).(3)The slip deficit is relatively higher in the Lenglongling,Jinqianghe,Maomaoshan,and Liupanshan faults(averaging about 3.42 mm/a,4.16 mm/a,4.23 mm/a,and 3.43 mm/a within 20 km).(4)The Qilian,Alxa,Xining,Lanzhou,and Haiyuan blocks rotate clockwise,while the Ordos block rotates counterclockwise.Additionally,by comparing different block models,the Haiyuan block should be considered independently.The Haiyuan fault zone adjusts surrounding block movements and uplifts Liupanshan mountain tectonically.The results can provide important references for understanding the regional earthquake risk and deformation mechanism.
基金funded by the National Natural Science Foundation of China under Grant 42264004the Open Fund of Hubei Luojia Laboratory under Grant 230100018the Guangxi Universities'1000 Young and Middleaged Backbone Teachers Training Program。
文摘The early excessive groundwater exploitation caused extensive land subsidence in the Hangjiahu Plain(HJHP),threatening geological stability and infrastructure.Government restrictions on groundwater exploitation led to considerable changes in surface subsidence.To explore the distribution and causes of surface deformation after restrictions on groundwater exploitation,we used 46 Sentinel-1A TOPS images to detect the spatiotemporal distribution of surface deformation in the HJHP area from January 2019 to October 2022 by SBAS-InSAR technique and verified it with leveling results.Results indicated the following:(1)The uneven surface subsidence in the HJHP was relatively apparent,with the localized subsidence concentrated in a few areas.During InSAR monitoring,the maximum subsidence velocity reached-21.0 mm/yr.The most severe subsidence area was in the Shengze town industrial park,with remarkable subsidence bowls and a maximum subsidence of-85.8 mm.(2)Results were verified by comparing the InSAR results with leveling,and the results were in good agreement,indicating the reliability of the InSAR results.(3)Comprehensive analysis of the spatial distribution characteristics of surface subsidence,geological data and rainfall changes revealed that fault activity and rainfall changes had varying degrees of effect on surface subsidence in the HJHP area.(4)The comparison with land use data indicated that the built area was the dominant mechanism of surface subsidence in the study area.In addition,the effect of industrial activities on the surface cannot be ignored.The study results provide scientific data support for emergency disaster prevention work in the HJHP.
基金supported by the National Natural Science Foundation of China(NO.42104028,42174030 and 42004017)the Open Fund of Hubei Luojia Laboratory(No.220100048 and 230100021)the Scientific Research Project of Hubei Provincial Department of Education,and Research Foundation of the Department of Natural Resources of Hunan Province(No.20230104CH)。
文摘The time-varying periodic variations in Global Navigation Satellite System(GNSS)stations affect the reliable time series analysis and appropriate geophysical interpretation.In this study,we apply the singular spectrum analysis(SSA)method to characterize and interpret the periodic patterns of GNSS deformations in China using multiple geodetic datasets.These include 23-year observations from the Crustal Movement Observation Network of China(CMONOC),displacements inferred from the Gravity Recovery and Climate Experiment(GRACE),and loadings derived from Geophysical models(GM).The results reveal that all CMONOC time series exhibit seasonal signals characterized by amplitude and phase modulations,and the SSA method outperforms the traditional least squares fitting(LSF)method in extracting and interpreting the time-varying seasonal signals from the original time series.The decrease in the root mean square(RMS)correlates well with the annual cycle variance estimated by the SSA method,and the average reduction in noise amplitudes is nearly twice as much for SSA filtered results compared with those from the LSF method.With SSA analysis,the time-varying seasonal signals for all the selected stations can be identified in the reconstructed components corresponding to the first ten eigenvalues.Moreover,both RMS reduction and correlation analysis imply the advantages of GRACE solutions in explaining the GNSS periodic variations,and the geophysical effects can account for 71%of the GNSS annual amplitudes,and the average RMS reduction is 15%.The SSA method has proved to be useful for investigating the GNSS timevarying seasonal signals.It could be applicable as an auxiliary tool in the improvement of nonlinear variations investigations.
基金Supported by the Fund of Institute of Seismology,CEA (No.IS200736026)
文摘This paper introduces the comparison of absolute gravity measurements in China during 2006-2008,and analyzing the survey and comparing the results,it shows that there is no obvious system error between the gravimeter of FG5/214 and FG5/232,the surveying accuracy is very high and repetition is good,and their inner surveying accuracy is about 2-3 microGal.
基金supported by the National Natural Science Foundation of China (Nos. 40974015,40637034)the Fund of Key Laboratory of Geodynamic Geodesy, Chinese Academy of Sciences (No. 09-18)the Fund of Key Laboratory of Geospace Environment and Geodesy,Ministry of Education,China (No. 07-12)
文摘Gravity Recovery and Climate Experiment（GRACE） observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 26,2004.This article focuses on investigating some gravita-tional effects caused by this huge earthquake.We computed the geoid height changes,the equivalent water height（EWH） changes,and the gravity changes using the GRACE Level-2 monthly spherical harmonic（SH） solutions released by University of Texas Center for Space Research（UTCSR）.The GRACE results agree well with the prediction by a dislocation model and are consistent with the results obtained by some previous scholars.In particular,we calculated the three components of the gravity gradient variations and found that they can recover the seismic-related signature more sensitively due to a certain degree of amplification of the signals.A positive-negative-positive mode predominates in the spatial distribution of the horizontal components of the gravity gradient variations,which is possibly attributed to the anomalies in the crustal density distribution caused by the uplift-subduction effect of the dip-slip earthquake.Moreover,the latitude components of the gravity gradient changes show strong suppression of the north-south stripes,which is due to the along-orbit measurements of the two GRACE satellites.We conclude that the posi-tive-negative-positive mode in latitude gravity gradient changes would be a more sensitive fea-ture to detect the deformations of some major dip-slip earthquakes by GRACE data.
基金the National Natural Science Foundation of China[Grant No.41874035]the National Natural Science Foundation of China[Grant No.12103035]+1 种基金the Natural Science Foundation of Hubei province,China[Grant No.2020CFB396]the Fundamental Research Funds for the Central Universities,China[Grant No.2042021kf0001].
文摘With the undergoing and planned implementations of mega constellations of thousands of Low Earth Orbiting(LEO)satellites,space will become even more congested for satellite operations.The enduring effects on the long-term space environment have been investigated by various researchers using debris environment models.This paper is focused on the imminent short-term effects of LEO mega constellations on the space operation environment concerned by satellite owners and operators.The effects are measured in terms of the Close Approaches(CAs)and overall collision probability.Instead of using debris environment models,the CAs are determined from integrated orbit positions,and the collision probability is computed for each CA considering the sizes and position covariance of the involving objects.The obtained results thus present a clearer picture of the space operation safety environment when LEO mega constellations are deployed.Many mega constellations are simulated,including a Starlink-like constellation of 1584 satellites,four possible generic constellations at altitudes between 1110 km and 1325 km,and three constellations of 1584 satellites each at the altitudes of 650 km,800 km,and 950 km,respectively,where the Resident Space Object(RSO)spatial density is the highest.The increases in the number of CAs and overall collision probability caused by them are really alarming.The results suggest that highly frequent orbital maneuvers are required to avoid collisions between existing RSOs and constellation satellites,and between satellites from two constellations at a close altitude,as such the constellation operation burden would be very heavy.The study is not only useful for satellite operators but a powerful signal for various stakeholders to pay serious attention to the development of LEO mega constellations.
基金the National 863 Program of China (No. 2006AA12Z325)
文摘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.
基金Founded by the National Natural Science Foundation of China （No.40637034, No.40574004）, the National 863 Program of China （No. 2006AA12Z211） and the Fund of Key Lab of Geodynamic Geodesy of Chinese Academy （No. L06-02）.
文摘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.
基金Supported by the Knowledge Innovation Program of the Chinese Academy of Sciences (Nos.KZCX1-YW-12-04,KZCX2-YW-201)the Instrument Developing Project of the Chinese Academy of Sciences (No.YZ200724)
文摘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.
基金supported jointly by the National Natural Science Foundation of China (Nos. U1939204, 41774056, and 41704046)the Fundamental Research Funds for the Central Universities (No. 2042020kf0010)
文摘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.
基金supported by the National Natural Science Foundation of China (Nos.41504015,41429401)the National 973 Project of China (No.2013CB733302)+2 种基金 China Postdoctoral Science Foundation (No.2015M572146)the National High Technology Research and Development Program of China (No.2011AA060503)the Surveying and Mapping Basic Research Program of National Administration of Surveying,Mapping and Geoinformation (No.15-01-08)
文摘Urban faults in Shenzhen are potential threats to city security and sustainable development. In consideration of the importance of the Shenzhen fault zone, the author provide a detailed interpretation on gravity data model. Bouguer gravity covering the whole Shenzhen City was calculated with a 1-km resolution. Wavelet multi-scale analysis（MSA） was applied to the Bouguer gravity data to obtain the multilayer residual anomalies corresponding to different depths. In addition, 2D gravity models were constructed along three profiles. The Bouguer gravity anomaly shows an NE-striking high-low-high pattern from northwest to southeast, strongly related to the main faults. According to the results of MSA, the correlation between gravity anomaly and faults is particularly significant from 4 to 12 km depth. The residual gravity with small amplitude in each layer indicates weak tectonic activity in the crust. In the upper layers, positive anomalies along most of faults reveal the upwelling of high-density materials during the past tectonic movements. The multilayer residual anomalies also yield important information about the faults, such as the vertical extension and the dip direction. The maximum depth of the faults is about 20 km. In general, NE-striking faults extend deeper than NW-striking faults and have a larger dip angle.
基金This study was funded by the National Natural Science Foundation of China(Grant Nos.41974015,42061134007 and 41474019).
文摘Both the Global Positioning System(GPS)and Gravity Recovery and Climate Experiment(GRACE)/GRACE Follow-On(GFO)provide effective tools to infer surface mass changes.In this paper,we combined GPS,GRACE/GFO spherical harmonic(SH)solutions and GRACE/GFO mascon solutions to analyze the total surface mass changes and terrestrial water storage(TWS)changes in the Shaan-Gan-Ning Region(SGNR)over the period from December 2010 to February 2021.To improve the reliability of GPS inversion results,an improved regularization Laplace matrix and monthly optimal regularization parameter estimation strategy were employed to solve the ill-posed problem.The results show that the improved Laplace matrix can suppress the edge effects better than that of the traditional Laplace matrix,and the corre-lation coefficient and standard deviation(STD)between the original signal and inversion results from the traditional and improved Laplace matrix are 0.84 and 0.88,and 17.49 mm and 15.16 mm,respectively.The spatial distributions of annual amplitudes and time series changes for total surface mass changes derived from GPS agree well with GRACE/GFO SH solutions and mascon solutions,and the correlation coefficients of total surface mass change time series between GPS and GRACE/GFO SH solutions,GPS and GRACE/GFO mascon solutions are 0.80 and 0.77.However,the obvious differences still exist in local regions.In addition,the seasonal characteristics,increasing and decreasing rate of TWS change time series derived from GPS,GRACE/GFO SH and mascon solutions agree well with the Global Land Data Assimilation System(GLDAS)hydrological model in the studied area,and generally consistent with the precipitation data.Meanwhile,TWS changes derived from GPS and GRACE mascon solutions in the SGNR are more reliable than those of GRACE SH solutions over the period from January 2016 to June 2017(the final operation phase of the GRACE mission).
基金supported by the NSFC project (41474046)the DQJB project (DQJB16B05) of the Institute of Geophysics, CEA
文摘On November 13, 2016, an MW7.8 earthquake struck Kaikoura in South Island of New Zealand. By means of back-projection of array recordings, ASTFs-analysis of global seismic recordings, and joint inversion of global seismic data and co-seismic In SAR data, we investigated complexity of the earthquake source. The result shows that the 2016 MW7.8 Kaikoura earthquake ruptured about 100 s unilaterally from south to northeast(~N28°–33°E), producing a rupture area about 160 km long and about 50 km wide and releasing scalar moment 1.01×1021 Nm. In particular, the rupture area consisted of two slip asperities, with one close to the initial rupture point having a maximal slip value ~6.9 m while the other far away in the northeast having a maximal slip value ~9.3 m. The first asperity slipped for about 65 s and the second one started 40 s after the first one had initiated. The two slipped simultaneously for about 25 s.Furthermore, the first had a nearly thrust slip while the second had both thrust and strike slip. It is interesting that the rupture velocity was not constant, and the whole process may be divided into 5 stages in which the velocities were estimated to be 1.4 km/s, 0 km/s, 2.1 km/s, 0 km/s and 1.1 km/s, respectively. The high-frequency sources distributed nearly along the lower edge of the rupture area, the highfrequency radiating mainly occurred at launching of the asperities, and it seemed that no high-frequency energy was radiated when the rupturing was going to stop.
文摘A GIS for ocean applications called ＂the Xiamen Atmospheric and Oceanographic Data Management and Display System （AODMDS）＂ has been designed and developed. The system is based on ArcObjects （AO）, a component-based GIS de- velopment tool. The paper discusses in detail the storage and organization of the atmospheric and oceanographic data, the strategy and methods for the visualization and mapping of oceanographic and atmospheric data, and the implementation of the methods in AODMDS. It also discusses some advanced display control techniques that expand the functions of ArcObjects One of the techniques is ＂gradient-fill-style color-map control,＂ which provides a feasible color-rich display control for all types of raster maps. As a stand-alone desktop GIS system built on AO, AODMDS provides effective data management and powerful mapping and visualization functions for atmospheric and oceanographic data.
