查看更多>>摘要:Variations in ice mass deform the Earth and modify its gravity field,a process known as Glacial Isostatic Adjustment(GIA).GIA in Antarctica remains poorly constrained due to the cumulative effect of past and present ice-mass changes,the unknown history of the past ice-mass change,and the uncertainties on the mechanical properties of the Earth.This paper investigates the effect of using Andrade and Burgers viscoelastic rheologies,rather than the commonly used Maxwell rheology,to model GIA-induced deformation in Antarctica.The Love number and Green's function formalism are used to compute the radial surface displacements and the gravity changes induced by the past and present ice-mass changes.We consider an Earth model whose elastic properties and radial structure are averaged from the Preliminary Reference Earth Model and two viscosity profiles to account for the recently published results on the present ice-mass changes.Using the three rheological laws affects the temporal response of the Earth differently,leading to smaller discrepancies than those induced by the two viscosity structures.The differences are the largest between Maxwell and Burgers rheologies during the 100-1000 years following the beginning of the surface-mass change.Results show that using the Andrade and Burgers rheologies allows the Earth to respond on decennial to centennial time scales,up to 10 m more than Maxwell.Considering only the recent ice-mass changes,the deformation rates derived from Burgers and Andrade rheologies are several times larger than those estimated by Maxwell rheology.
查看更多>>摘要:Modeling the earth's fluid and elastic response to the melting of the glaciers of the last ice age is the most direct way to infer the earth's radial viscosity profile.Here,we compare two methods for calculating the viscoelastic response to surface loading.In one,the elastic equation of motion is converted to a visco-elastic equation using the Correspondence Principle.In the other,elastic deformation is added to the viscous flow as isostatic adjustment proceeds.The two modeling methods predict adjustment histories that are different enough to potentially impact the interpretation of the observed glacial isostatic adjustment(GIA).The differences arise from buoyancy and whether fluid displacements are subjected to hydrostatic pre-stress.The methods agree if they use the same equations and boundary conditions.The origin of the differences is determined by varying the boundary conditions and pre-stress application.
查看更多>>摘要:While the geodetic excitation x(t)of polar motion p(t)is essential to improve our understanding of global mass redistributions and relative motions with respect to the terrestrial frame,the widely adopted method to derive x(t)from p(t)has biases in both amplitude and phase responses.This study has developed a new simple but more accurate method based on the combination of the frequency-and time-domain Liouville's equation(FTLE).The FTLE method has been validated not only with 6-h sampled synthetic excitation series but also with daily and 6-h sampled polar motion measurements as well as x(t)produced by the interactive webpage tool of the International Earth Rotation and Reference Systems Service(IERS).Numerical comparisons demonstrate that x(t)derived from the FTLE method has superior performances in both the time and frequency domains with respect to that obtained from the widely adopted method or the IERS webpage tool,provided that the input p(t)series has a length around or more than 25 years,which presents no practical limitations since the necessary polar motion data are readily available.The FTLE code is provided in the form of MatLab function.
查看更多>>摘要:Due to environmental noise and human factors,magnetic data collected in the field often contain various noises and interferences that significantly affect the subsequent data processing and inter-pretation.Empirical Mode Decomposition(EMD),an adaptive multiscale analysis method for nonlinear and non-stationary signals,is widely used in geophysical and geodetic data processing.Compared with traditional EMD,Improved Complete Ensemble EMD with Adaptive Noise(ICEEMDAN)is more effective in addressing the problem of mode mixing.Based on the principles of ID ICEEMDAN,this paper presents an alternative algorithm for 2D ICEEMDAN,extending its application to two-dimensional scenarios.The effectiveness of the proposed approach is demonstrated through syn-thetic signal experiments,which show that the 2D ICEEMDAN exhibits a weaker mode mixing effect compared to the traditional bidimensional EMD(BEMD)method.Furthermore,to improve the per-formance of the denoising method based on 2D ICEEMDAN and preserve useful signals in high-frequency components,an improved soft thresholding technique is introduced.Synthetic magnetic anomaly data testing indicates that our denoising method effectively preserves signal continuity and outperforms traditional soft thresholding methods.To validate the practical application of this improved threshold denoising method based on 2D ICEEMDAN,it is applied to ground magnetic survey data in the Yandun area of Xinjiang.The results demonstrate the effectiveness of the method in removing noise while retaining essential information from practical magnetic anomaly data.In particular,practical applications suggest that 2D ICEEMDAN can extract trend signals more accurately than the BEMD.In conclusion,as a potential tool for multi-scale decomposition,the 2D ICEEMDAN is versatile in processing and analyzing 2D geophysical and geodetic data.
查看更多>>摘要:The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occur-rences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Mul-tiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.
查看更多>>摘要:The precise correction of atmospheric zenith tropospheric delay(ZTD)is significant for the Global Navigation Satellite System(GNSS)performance regarding positioning accuracy and convergence time.In the past decades,many empirical ZTD models based on whether the gridded or scattered ZTD products have been proposed and widely used in the GNSS positioning applications.But there is no comprehensive evaluation of these models for the whole China region,which features complicated topography and climate.In this study,we completely assess the typical empirical models,the IGGtropSH model(gridded,non-meteorology),the SHAtropE model(scattered,non-meteorology),and the GPT3 model(gridded,meteorology)using the Crustal Movement Observation Network of China(CMONOC)network.In general,the results show that the three models share consistent per-formance with RMSE/bias of 37.45/1.63,37.13/2.20,and 38.27/1.34 mm for the GPT3,SHAtropE and IGGtropSH model,respectively.However,the models had a distinct performance regarding geographical distribution,elevation,seasonal variations,and daily variation.In the southeastern re-gion of China,RMSE values are around 50 mm,which are much higher than that in the western region,approximately 20 mm.The SHAtropE model exhibits better performance for areas with large variations in elevation.The GPT3 model and the IGGtropSH model are more stable across different months,and the SHAtropE model based on the GNSS data exhibits superior performance across various UTC epochs.
查看更多>>摘要:Utilizing the adopted average topographic density of 2670 kg/m3 in the reduction of gravity anomalies introduces errors attributed to topographic density variations,which consequently affect geoid modeling accuracy.Furthermore,the mean gravity along the plumbline within the topography in the definition of Helmert orthometric heights is computed approximately by applying the Poincaré-Prey gravity reduction where the topographic density variations are disregarded.The Helmert orthometric heights of bench-marks are then affected by errors.These errors could be random or systematic depending on the specific geological setting of the region where the leveling network is physically established and/or the geoid model is determined.An example of systematic errors in orthometric heights can be given for large regions characterized by sediment or volcanic deposits,the density of which is substantially lower than the adopted topographic density used in Helmert's definition of heights.The same applies to geoid modeling errors.In this study,we investigate these errors in the Hong Kong territory,where topographic density is about 20%lower than the density of 2670 kg/m3.We use the digital rock density model to estimate the effect of topographic density variations on the geoid and orthometric heights.Our results show that this effect on the geoid and Helmert orthometric heights reach maxima of about 2.1 and 0.5 cm,respectively.Both results provide clear evidence that rock density models are essential in physical geodesy applications involving gravimetric geoid modeling and orthometric height determination despite some criticism that could be raised regarding the reliability of these density models.However,in regions dominated by sedimentary and igneous rocks,the geological information is essential in these applications because topographic densities are substantially lower than the average density of 2670 kg/m3,thus introducing large systematic errors in geoid and orthometric heights.
查看更多>>摘要:Frequency modulated continuous wave(FMCW)radar is an advantageous sensor scheme for target estimation and environmental perception.However,existing algorithms based on discrete Fourier transform(DFT),multiple signal classification(MUSIC)and compressed sensing,etc.,cannot achieve both low complexity and high resolution simultaneously.This paper proposes an efficient 2-D MUSIC algo-rithm for super-resolution target estimation/tracking based on FMCW radar.Firstly,we enhance the efficiency of 2-D MUSIC azimuth-range spectrum estimation by incorporating 2-D DFT and multi-level resolution searching strategy.Secondly,we apply the gradient descent method to tightly integrate the spatial continuity of object motion into spectrum estimation when processing multi-epoch radar data,which improves the efficiency of continuous target tracking.These two approaches have improved the algorithm efficiency by nearly 2-4 orders of magnitude without losing accuracy and resolution.Simu-lation experiments are conducted to validate the effectiveness of the algorithm in both single-epoch estimation and multi-epoch tracking scenarios.
查看更多>>摘要:Extreme hydrological events such as droughts and floods have been increasingly influenced by abnormal atmospheric disturbances caused by human activity and global warming.The Dongting Lake Basin(DLB)has experienced challenging circumstances over the past 20 years due to complex climatic variations,leading to extreme flooding and drought.This study aims to investigate the spatiotemporal variation in terrestrial water storage anomalies(TWSA)over the DLB using data from the GRACE/GRACE-FO and GLDAS spanning the latest two decades.A significant decline trend in TWSA is unveiled from July 2019 to May 2023,with the rate of change determined as-1.94 cm/year and-1.99 cm/year based on the GRACE/GRACE-FO and GLDAS,respectively.The GRACE-Drought Severity Index(DSI)is employed to identify and evaluate the severity and spatiotemporal evolution of the 2022 drought event in the DLB.The results accurately capture the drought event,which began in July 2022 and continued until March 2023,with the most severe conditions occurring in October 2022,when the GRACE-DSI value stood at-2.06 and the TWSA decreased by 15.24 cm and 33.51 cm relative to the same month in 2021 and 2020,respectively.Additionally,the daily water level variation at the Chenglingji hydrological gauging station in 2022 broke previous records,reaching a minimum of only 19 m.Comparing the 2022 drought event with the drought events in 2006 and 2011,the impact of drought on vegetation growth conditions was relatively small,but there was still significant vegetation degradation across the DLB.
查看更多>>摘要:The longitudinal dependence of the behavior of ionospheric parameters has been the subject of a number of works where significant variations are discovered.This also applies to the prediction of the iono-spheric total electron content(TEC),which neural network methods have recently been widely used.However,the results are mainly presented for a limited set of meridians.This paper examines the lon-gitudinal dependence of the TEC forecast accuracy in the equatorial zone.In this case,the methods are used that provided the best accuracy on three meridians:European(30° E),Southeastern(110° E)and American(75° W).Results for the stations considered are analyzed as a function of longitude using the Jet Propulsion Laboratory Global Ionosphere Map(JPL GIM)for 2015.These results are for 2 h ahead and 24 h ahead forecast.It was found that in this case,based on the metric values,three groups of archi-tectures can be distinguished.The first group included long short-term memory(LSTM),gated recurrent unit(GRU),and temporal convolutional networks(TCN)models as a part of unidirectional deep learning models;the second group is based on the recurrent models from the first group,which were supple-mented with a bidirectional algorithm,increasing the TEC forecasting accuracy by 2-3 times.The third group,which includes the bidirectional TCN architecture(BiTCN),provided the highest accuracy.For this architecture,according to data obtained for 9 equatorial stations,practical independence of the TEC prediction accuracy from longitude was observed under the following metrics(Mean Absolute Error MAE,Root Mean Square Error RMSE,Mean Absolute Percentage Error MAPE):MAE(2 h)is 0.2 TECU approximately;MAE(24 h)is 0.4 TECU approximately;RMSE(2 h)is less than 0.5 TECU except Niue station(RMSE(2 h)is 1 TECU approximately);RMSE(24 h)is in the range of 1.0-1.7 TECU;MAPE(2 h)<1%except Darwin station,MAPE(24 h)<2%.This result was confirmed by data from additional 5 stations that formed latitudinal chains in the equatorial part of the three meridians.The complete cor-respondence of the observational and predicted TEC values is illustrated using several stations for disturbed conditions on December 19-22,2015,which included the strongest magnetic storm in the second half of the year(min Dst=-155 nT).