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 transfo...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.展开更多
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 j...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.展开更多
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...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.展开更多
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 me...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.展开更多
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 grav...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.展开更多
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 significan...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.展开更多
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...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.展开更多
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).Compa...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).展开更多
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 ...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.展开更多
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 d...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.展开更多
In ground-based GPS meteorology, Tm is a key parameter to calculate the conversion factor that can convert the zenith wet delay(ZWD) to precipitable water vapor(PWV). It is generally acknowledged that Tm is in an appr...In ground-based GPS meteorology, Tm is a key parameter to calculate the conversion factor that can convert the zenith wet delay(ZWD) to precipitable water vapor(PWV). It is generally acknowledged that Tm is in an approximate linear relationship with surface temperature Ts, and the relationship presents regional variation. This paper employed sliding average method to calculate correlation coefficients and linear regression coefficients between Tm and Ts at every 2°× 2.5° grid point using Ts data from European Centre for Medium-Range Weather Forecasts(ECMWF) and Tm data from "GGOS Atmosphere", yielding the grid and bilinear interpolation-based Tm Grid model. Tested by Tm and Ts grid data, Constellation Observation System of Meteorology, Ionosphere, and Climate(COSMIC) data and radiosonde data, the Tm Grid model shows a higher accuracy relative to the Bevis Tm-Ts relationship which is widely used nowadays. The Tm Grid model will be of certain practical value in high-precision PWV calculation.展开更多
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 t...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/s2, which is around 23% of the present observed value. Studies show in detail that both θ and ω 3 experience complex fluctua-tions 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.展开更多
We compile the GOCE-only satellite model GOSG01 S complete to spherical harmonic degree of 220 using Satellite Gravity Gradiometry(SGG) data and the Satellite-to-Satellite Tracking(SST) observations along the GOCE orb...We compile the GOCE-only satellite model GOSG01 S complete to spherical harmonic degree of 220 using Satellite Gravity Gradiometry(SGG) data and the Satellite-to-Satellite Tracking(SST) observations along the GOCE orbit based on applying a least-squares analysis.The diagonal components(V_(xx),V_(yy),V_(zz)) of the gravitational gradient tensor are used to form the system of observation equations with the band-pass ARMA filter.The point-wise acceleration observations(a_x,a_y,a_z) along the orbit are used to form the system of observation equations up to the maximum spherical harmonic degree/order 130.The analysis of spectral accuracy characteristics of the newly derived gravitational model GOSG01 S and the existing models GOTIM04 S,GODIR04S,GOSPW04 S and JYY_GOCE02S based on their comparison with the ultrahigh degree model EIGEN-6C2 reveals a significant consistency at the spectral window approximately between 80 and 190 due to the same period SGG data used to compile these models.The GOCE related satellite gravity models GOSG01 S,GOTIM05S,GODIR05 S,GOTIM04S,GODIR04 S,GOSPW04S,JYY_-GOCE02 S,EIGEN-6C2 and EGM2008 are also validated by using GPS-leveling data in China and USA.According to the truncation at degree 200,the statistic results show that all GGMs have very similar differences at GPS-leveling points in USA,and all GOCE related gravity models have better performance than EGM2008 in China.This suggests that all these models provide much more information on the gravity field than EGM2008 in areas with low terrestrial gravity coverage.And STDs of height anomaly differences in China for the selected truncation degrees show that GOCE has improved the accuracy of the global models beyond degree 90 and the accuracies of the models improve from 24 cm to 16 cm.STDs of geoid height differences in USA show that GOSGOIS model has best consistency comparing with GPSleveling data for the frequency band of the degree between 20 and 160.展开更多
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.展开更多
Many studies revealed that the Earth medium's lateral heterogeneity can cause considerable effects on the co- and post-seismic deformation field. In this study, the threedimensional finite element numerical method...Many studies revealed that the Earth medium's lateral heterogeneity can cause considerable effects on the co- and post-seismic deformation field. In this study, the threedimensional finite element numerical method are adopted to quantify the effects of lateral heterogeneity caused by material parameters and fault dip angle on the co- and postseismic deformation in the near- and far-field. Our results show that: 1) the medium's lateral heterogeneity does affect the co-seismic deformation, with the effects increasing with the medium's lateral heterogeneity caused by material parameters; 2) the Lame parameters play a more dominant role than density in the effects caused by lateral heterogeneity; 3) when a fault's dip angle is smaller than 90, the effects of the medium's lateral heterogeneity on the hanging wall are greater than on the footwall; 4) the impact of lateral heterogeneity caused by the viscosity coefficient on the post-seismic deformation can affect a large area, including the near- and far-field.展开更多
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.展开更多
According to the quality control needs of filling construction of the face rockfill dam, by means of the global satellite positioning technology, the wireless data communication technology, the computer technology and...According to the quality control needs of filling construction of the face rockfill dam, by means of the global satellite positioning technology, the wireless data communication technology, the computer technology and the data processing and analysis technology, and integrating with the roller compaction machine, the GPS real time supervisory system is developed in this paper. It can be used to real timely supervise the construction quality of the roller compaction for filling engineering. The composition and applied characteristics of GPS system, and the key technique problem and solution of the design are discussed. The height accuracy of GPS system is analyzed and the preliminary application is introduced.展开更多
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).展开更多
According to the space-geodetic data recorded at globally distributed stations over solid land spanning a period of more than 20-years under the International Terrestrial Reference Frame 2008,our previous estimate of ...According to the space-geodetic data recorded at globally distributed stations over solid land spanning a period of more than 20-years under the International Terrestrial Reference Frame 2008,our previous estimate of the average-weighted vertical variation of the Earth's solid surface suggests that the Earth's solid part is expanding at a rate of 0.24 ± 0.05 mm/a in recent two decades.In another aspect,the satellite altimetry observations spanning recent two decades demonstrate the sea level rise(SLR) rate 3.2 ± 0.4 mm/a,of which1.8 ± 0.5 mm/a is contributed by the ice melting over land.This study shows that the oceanic thermal expansion is 1.0 ± 0.1 mm/a due to the temperature increase in recent half century,which coincides with the estimate provided by previous authors.The SLR observation by altimetry is not balanced by the ice melting and thermal expansion,which is an open problem before this study.However,in this study we infer that the oceanic part of the Earth is expanding at a rate about 0.4 mm/a.Combining the expansion rates of land part and oceanic part,we conclude that the Earth is expanding at a rate of 0.35 ± 0.47 mm/a in recent two decades.If the Earth expands at this rate,then the altimetry-observed SLR can be well explained.展开更多
基金supported by the National Natural Science Foundation of China(41874011,42074007)the Fundamental Research Funds for the Central Universities(2042023kfyq01)。
文摘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.
基金supported by the Natural Science Basic Research Program of Shaanxi(No.2023-JC-QN-0306)the Special Fund of the Institute of Geophysics,China Earthquake Administration(No.DQJB21B32)the National Natural Science Foundation of China(No.42174069).
文摘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.
基金the Natural Science Innovation Group Foundation of China under Grants NO.41721003the Science and Technology Support Project of Department of Natural Resources of Hubei Province under Grants NO.ZRZY2022KJ29+1 种基金the Special Fund of Hubei Luojia Laboratory under Grants NO.220100020the National Natural Science Foundation of China under Grants NO.42174030.
文摘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.
基金funded by Key R&D project of Ningxia Hui Autonomous Region(2021BDE931027)Science and technology project of State Grid Ningxia Electric Power Co.Ltd.(229DK2004P).
文摘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.
基金National Key R&D Program of China(Nos.2021YFB3900200,2021YFB3900203)。
文摘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.
基金National Key R&D Program of China(No.2021YFB3901301)National Natural Science Foundation of China(Nos.42074043,42122026,42174038)Youth Innovation Promotion Association of the Chinese Academy of Sciences(No.Y9E006033D)。
文摘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.
基金supported by grants from the National Natural Science Foundation of China(Grant No.sU2139205,41774011,41874011)the National Key Research and Development Program of China(Grant No.2018YFC1503605)。
文摘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.
基金supported by National Natural Science Foundation of China under Grants 42192531 and 42192534the Special Fund of Hubei Luojia Laboratory(China)under Grant 220100001the Natural Science Foundation of Hubei Province for Distinguished Young Scholars(China)under Grant 2022CFA090。
文摘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).
基金supported by the National Natural Science Foundation of China(Grant Nos.42264004,42274033,and 41904012)the Open Fund of Hubei Luojia Laboratory(Grant Nos.2201000049 and 230100018)+2 种基金the Guangxi Universities’1,000 Young and Middle-aged Backbone Teachers Training Program,the Fundamental Research Funds for Central Universities(Grant No.2042022kf1197)the Natural Science Foundation of Hubei(Grant No.2020CFB282)the China Postdoctoral Science Foundation(Grant Nos.2020T130482,2018M630879)。
文摘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.
基金supported by the National Natural Science Foundation of China(Grants:41974022 and 42192531)the Educational Commission of Hubei Province of China(Grant:2020CFA109)+1 种基金the Special Fund of Hubei Luojia Laboratory(grant#220100002)Key Laboratory of Geospace Environment and Geodesy,Ministry of Education,Wuhan University(grant#210204)。
文摘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.
基金supported by National Natural Science Foundation of China(41301377)by the Fundamental Research Funds for the Central Universities(2014214020202)by Surveying and Mapping Basic Research Program of National Administration of Surveying,Mapping and Geoinformation(13-02-09)
文摘In ground-based GPS meteorology, Tm is a key parameter to calculate the conversion factor that can convert the zenith wet delay(ZWD) to precipitable water vapor(PWV). It is generally acknowledged that Tm is in an approximate linear relationship with surface temperature Ts, and the relationship presents regional variation. This paper employed sliding average method to calculate correlation coefficients and linear regression coefficients between Tm and Ts at every 2°× 2.5° grid point using Ts data from European Centre for Medium-Range Weather Forecasts(ECMWF) and Tm data from "GGOS Atmosphere", yielding the grid and bilinear interpolation-based Tm Grid model. Tested by Tm and Ts grid data, Constellation Observation System of Meteorology, Ionosphere, and Climate(COSMIC) data and radiosonde data, the Tm Grid model shows a higher accuracy relative to the Bevis Tm-Ts relationship which is widely used nowadays. The Tm Grid model will be of certain practical value in high-precision PWV calculation.
基金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/s2, which is around 23% of the present observed value. Studies show in detail that both θ and ω 3 experience complex fluctua-tions 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.
基金financially supported by the National Key Basic Research Program of China(973 program,grant no.:2013CB733302,2013CB733301)the Major International(Regional) Joint Research Project(grant no.:41210006)+1 种基金DAAD Thematic Network Project(grant no.:57173947)the National Natural Science Foundation of China(grant No.41374022)
文摘We compile the GOCE-only satellite model GOSG01 S complete to spherical harmonic degree of 220 using Satellite Gravity Gradiometry(SGG) data and the Satellite-to-Satellite Tracking(SST) observations along the GOCE orbit based on applying a least-squares analysis.The diagonal components(V_(xx),V_(yy),V_(zz)) of the gravitational gradient tensor are used to form the system of observation equations with the band-pass ARMA filter.The point-wise acceleration observations(a_x,a_y,a_z) along the orbit are used to form the system of observation equations up to the maximum spherical harmonic degree/order 130.The analysis of spectral accuracy characteristics of the newly derived gravitational model GOSG01 S and the existing models GOTIM04 S,GODIR04S,GOSPW04 S and JYY_GOCE02S based on their comparison with the ultrahigh degree model EIGEN-6C2 reveals a significant consistency at the spectral window approximately between 80 and 190 due to the same period SGG data used to compile these models.The GOCE related satellite gravity models GOSG01 S,GOTIM05S,GODIR05 S,GOTIM04S,GODIR04 S,GOSPW04S,JYY_-GOCE02 S,EIGEN-6C2 and EGM2008 are also validated by using GPS-leveling data in China and USA.According to the truncation at degree 200,the statistic results show that all GGMs have very similar differences at GPS-leveling points in USA,and all GOCE related gravity models have better performance than EGM2008 in China.This suggests that all these models provide much more information on the gravity field than EGM2008 in areas with low terrestrial gravity coverage.And STDs of height anomaly differences in China for the selected truncation degrees show that GOCE has improved the accuracy of the global models beyond degree 90 and the accuracies of the models improve from 24 cm to 16 cm.STDs of geoid height differences in USA show that GOSGOIS model has best consistency comparing with GPSleveling data for the frequency band of the degree between 20 and 160.
基金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.
基金co-supported by the National Natural Science Foundation of China (41431069)the State Key Development Program for Basic Research of China (2013CB733304, 2013CB733303)+1 种基金the Doctoral Fund of Ministry of Education of China (20110141130010)China Postdoctoral Science Foundation funded project (2013M542062)
文摘Many studies revealed that the Earth medium's lateral heterogeneity can cause considerable effects on the co- and post-seismic deformation field. In this study, the threedimensional finite element numerical method are adopted to quantify the effects of lateral heterogeneity caused by material parameters and fault dip angle on the co- and postseismic deformation in the near- and far-field. Our results show that: 1) the medium's lateral heterogeneity does affect the co-seismic deformation, with the effects increasing with the medium's lateral heterogeneity caused by material parameters; 2) the Lame parameters play a more dominant role than density in the effects caused by lateral heterogeneity; 3) when a fault's dip angle is smaller than 90, the effects of the medium's lateral heterogeneity on the hanging wall are greater than on the footwall; 4) the impact of lateral heterogeneity caused by the viscosity coefficient on the post-seismic deformation can affect a large area, including the near- and far-field.
基金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.
文摘According to the quality control needs of filling construction of the face rockfill dam, by means of the global satellite positioning technology, the wireless data communication technology, the computer technology and the data processing and analysis technology, and integrating with the roller compaction machine, the GPS real time supervisory system is developed in this paper. It can be used to real timely supervise the construction quality of the roller compaction for filling engineering. The composition and applied characteristics of GPS system, and the key technique problem and solution of the design are discussed. The height accuracy of GPS system is analyzed and the preliminary application is introduced.
基金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 National 973 Project China(2013CB733305,2013CB733301)National Natural Science Foundation of China(41174011,41429401,41210006,41128003,41021061)
文摘According to the space-geodetic data recorded at globally distributed stations over solid land spanning a period of more than 20-years under the International Terrestrial Reference Frame 2008,our previous estimate of the average-weighted vertical variation of the Earth's solid surface suggests that the Earth's solid part is expanding at a rate of 0.24 ± 0.05 mm/a in recent two decades.In another aspect,the satellite altimetry observations spanning recent two decades demonstrate the sea level rise(SLR) rate 3.2 ± 0.4 mm/a,of which1.8 ± 0.5 mm/a is contributed by the ice melting over land.This study shows that the oceanic thermal expansion is 1.0 ± 0.1 mm/a due to the temperature increase in recent half century,which coincides with the estimate provided by previous authors.The SLR observation by altimetry is not balanced by the ice melting and thermal expansion,which is an open problem before this study.However,in this study we infer that the oceanic part of the Earth is expanding at a rate about 0.4 mm/a.Combining the expansion rates of land part and oceanic part,we conclude that the Earth is expanding at a rate of 0.35 ± 0.47 mm/a in recent two decades.If the Earth expands at this rate,then the altimetry-observed SLR can be well explained.