The elasto-gravitational deformation response of the Earths solid parts to the perturbations of the pressure and gravity on the core-mantle boundary (CMB) and the solid inner core boundary (ICB), due to the dynamical ...The elasto-gravitational deformation response of the Earths solid parts to the perturbations of the pressure and gravity on the core-mantle boundary (CMB) and the solid inner core boundary (ICB), due to the dynamical behaviors of the fluid outer core (FOC), is discussed. The internal load Love numbers, which are formulized in a general form in this study, are employed to describe the Earths deformation. The preliminary reference Earth model (PREM) is used as an example to calculate the internal load Love numbers on the Earths surface, CMB and ICB, respectively. The characteristics of the Earths deformation variation with the depth and the perturbation periods on the boundaries of the FOC are also investigated. The numerical results indicate that the internal load Love numbers decrease quickly with the increasing degree of the spherical harmonics of the displacement and depend strongly on the perturbation frequencies, especially on the high frequencies. The results, obtained in this work, can be used to construct the boundary conditions for the core dynamics of the long-period oscillations of the Earths fluid outer core.展开更多
We apply ambient noise tomography to significant seismic data resources in a region including the northeastern Tibetan plateau, the Ordos block and the Sichuan basin. The seismic data come from about 160 stations of t...We apply ambient noise tomography to significant seismic data resources in a region including the northeastern Tibetan plateau, the Ordos block and the Sichuan basin. The seismic data come from about 160 stations of the provincial broadband digital seismograph networks of China. Ambient noise cross-correlations are performed on the data recorded between 2007 and 2009 and high quality inter-station Rayleigh phase velocity dispersion curves are obtained between periods of 6 s to 35 s. Resulting Rayleigh wave phase velocity maps possess a lateral resolution between 100 km and 200 kin. The phase velocities at short periods （〈20 s） are lower in the Sichuan basin, the northwest segment of the Ordos block and the Weihe graben, and outline sedimentary deposits. At intermediate and long periods （〉25 s）, strong high velocity anomalies are observed within the Ordos block and the Sichuan basin and low phase velocities are imaged in the northeastern Tibetan plateau, reflecting the variation of crustal thickness from the Tibetan plateau to the neighboring regions in the east. Crustal and uppermost mantle shear wave velocities vary strongly between the Tibetan plateau, the Sichuan basin and the Ordos block. The Ordos block and the Sichuan basin are dominated by high shear wave velocities in the crust and uppermost mantle. There is a triangle-shaped low velocity zone located in the northeastern Tibetan plateau, whose width narrows towards the eastern margin of the plateau. No low velocity zone is apparent beneath the Qinling orogen, suggesting that mass may not be able to flow eastward through the boundary between the Ordos block and the Sichuan basin in the crust and uppermost mantle.展开更多
Based on the elastic dislocation theory, multilayered crustal model, and rupture model obtained by seismic waveform inversion, we calculated the coand post-seismic surface deformation and gravity changes caused by the...Based on the elastic dislocation theory, multilayered crustal model, and rupture model obtained by seismic waveform inversion, we calculated the coand post-seismic surface deformation and gravity changes caused by the Yushu M W 6.9 earthquake occurred on April 14, 2010. The observed GPS velocity field and gravity field in Yushu areas are disturbed by the coand post-seismic effects induced by Yushu earthquake, thus the theoretical coand post-seismic deformation and gravity changes will provide important modification for the background tectonic movement of Yushu and surrounding regions. The time relaxation results show that the influences of Yushu earthquake on Yushu and surrounding areas will last as long as 30 to 50 years. The maximum horizontal displacement, vertical uplift and settlement are about 1.96, 0.27 and 0.16 m, respectively, the maximal positive and negative value of gravity changes are 8.892×10-7 m·s-2 and -4.861×10-7 m·s-2 , respectively. Significant spatial variations can be found on the coand post-seismic effects： The co-seismic effect mainly concentrates in the region near the rupture fault, while viscoelastic relaxation mostly acts on the far field. Therefore, when using the geodetic data to research tectonic motion, we should not only consider the effect of co-seismic caused by earthquake, but also pay attention to the effect of viscoelastic relaxation.展开更多
The buildings and structures of mines were monitored automatically using modern surveying technology.Through the analysis of the monitoring data,the deformation characteristics were found out from three aspects contai...The buildings and structures of mines were monitored automatically using modern surveying technology.Through the analysis of the monitoring data,the deformation characteristics were found out from three aspects containing points,lines and regions,which play an important role in understanding the stable state of buildings and structures.The stability and deformation of monitoring points were analysed,and time-series data of monitoring points were denoised with wavelet analysis and Kalman filtering,and exponent function and periodic function were used to get the ideal deformation trend model of monitoring points.Through calculating the monitoring data obtained,analyzing the deformation trend,and cognizing the deformation regularity,it can better service mine safety production and decision-making.展开更多
A three-dimensional crustal density model beneath North China is determined using P-wave traveltimes and gravity datum by sequential inversion method in this paper. To improve the uniqueness of the solution, we used a...A three-dimensional crustal density model beneath North China is determined using P-wave traveltimes and gravity datum by sequential inversion method in this paper. To improve the uniqueness of the solution, we used a linear relationship between velocity and density to achieve mutual conversions and constraints between velocity difference and density difference. Algebraic reconstruction technique （ART） was used in density inversion, which highly improved the calculation speed comparing with common least squares method. The inversion results indicate that the crustal density beneath North China is extremely inhomogeneous with its distribution coherent with regional tectonics. The lower crust of Taihang mountain tectonic zone shows an obvious low density characteristic. We proposed that it may be an indicator of upwelling of mantle materials or partial melting of lowermost crust.展开更多
On April 20, 2013, an earthquake with mag- nitude 7.0 occurred in the southwest of the Longmenshan fault system in and around Lushan County, Sichuan Province, China. This devastating earthquake killed hun- dreds of pe...On April 20, 2013, an earthquake with mag- nitude 7.0 occurred in the southwest of the Longmenshan fault system in and around Lushan County, Sichuan Province, China. This devastating earthquake killed hun- dreds of people, injured 10 thousand others, and collapsed countless buildings. In order to analyze the potential risk of this big earthquake, we calculate the co- and post-seismic surface deformation and gravity changes of this event. In this work, a multilayered crustal model is designed, and the elastic dislocation theory is utilized to calculate the co- and post-seismic deformations and gravity changes. During the process, a rupture model obtained by seismic waveform inversion （Liu et al. Sci China Earth Sci 56（7）： 1187-1192, 2013） is applied. The time-dependent relaxation results show that the influences on Lushan and its surrounding areas caused by the Ms7.0 Lushan earthquake will last as long as 10 years. The maximum horizontal displacement, vertical uplift, and settlement are about 5 cm, 21.24 cm, and 0.16 m, respectively; the maximal positive and nega- tive values of gravity changes are 45 and -0.47 μGal, respectively. These results may be applied to evaluate the long-term potential risk caused by this earthquake and to provide necessary information for post-earthquake reconstruction.展开更多
We compare the factors which affect the movement of Tibetan Plateau by building three types of finite element models： an elastic materials （M-EC）, a continuous model composed by non-linear materials （M-PC）, and a...We compare the factors which affect the movement of Tibetan Plateau by building three types of finite element models： an elastic materials （M-EC）, a continuous model composed by non-linear materials （M-PC）, and an elastic model with discontinuous fault movements （M-ET）. Both in M-ET and M-EC, the materials in Qiangtang and Lhasa block are elastic, and in M-ET, discontinuous movement of faults is considered for evaluating the effects of strike-slip faults. In model M-PC Druker-Prager plastic materials are used in Qiangtang and Lhasa block. Comparisons of the numerical simulation and the GPS observations show following characteristics： （1） Under present tectonic environment, short-term deformation of Tibetan Plateau can be simulated well by elastic models; （2） Discontinuous fault activities increase the lateral extrusion of the eastern part of Tibetan Plateau, reduce the stress field level in Qiangtang, Tarim and Qaidam blocks and strengthen the E-W extensional force in the east and the west parts of Qiangtang block; （3） Properties of plastic materials reduce the total stress field and the E-W extensional force, thus, the normal fault earthquakes in southern Tibet is mainly owed to the effect of active fault movement. Based on the numerical simulations we speculate that faults movement may play a more important role on the kinematic pattern of Tibetan Plateau than bulk properties.展开更多
The crustal structure of Xiachayu-Gonghe geophysical profile in eastern Tibetan plateau is simulated with Bouguer anomaly corrected for sediments and lithosphere. The forward simulation shows that the thickness of upp...The crustal structure of Xiachayu-Gonghe geophysical profile in eastern Tibetan plateau is simulated with Bouguer anomaly corrected for sediments and lithosphere. The forward simulation shows that the thickness of upper crust in eastem Tibetan plateau is about 20 km, and the density is 2.78 × 10^3 kg/m^3. The bottom interface of middle crust changes from 30 km to 40 km, the density of middle crust is 2.89 × 10^3 kg/m^3. The materials with low density of 2.78 × 10^3 kg/m^3 exist in middle crust, and those with high density of 3.33 × 10^3 kg/m^3 exist at the bottom of middle crust between Wenquan and Tanggemu. The density is 3.10× 10^3 kg/m^3 in lower crust. The shallowest depth of Moho interface is about 56 km, and the deepest one is about 74 km, the undulation of interface is large, the deep Moho is located in Xiachayu, Chayu, Nujiang, and Wenquan. The crustal density of eastern Tibetan plateau is larger than that of central section; the low velocity layers are located in middle crust and bottom in eastern Tibetan plateau and at the bottom of the upper crust in the central plateau.展开更多
We present a three-dimensional geometry model for tortuosity of streamlines in porous media with randomly placed cylindrical particles. The proposed model is expressed as functions of porosity and geometrical paramete...We present a three-dimensional geometry model for tortuosity of streamlines in porous media with randomly placed cylindrical particles. The proposed model is expressed as functions of porosity and geometrical parameters with no empirical constant. This might be helpful for understanding the physical mechanism for tortuosity of streamlines in three-dimensional porous media. The model predictions are found to be in good agreement with the experimental data available.展开更多
Firstly, the new combined error model of cumulative geoid height influenced by four error sources, including the inter-satellite range-rate of an interferometric laser （K-band） ranging system, the orbital position a...Firstly, the new combined error model of cumulative geoid height influenced by four error sources, including the inter-satellite range-rate of an interferometric laser （K-band） ranging system, the orbital position and velocity of a global positioning system （GPS） receiver and non-conservative force of an accelerometer, is established from the perspectives of the power spectrum principle in physics using the semi-analytical approach. Secondly, the accuracy of the global gravitational field is accurately and rapidly estimated based on the combined error model; the cumulative geoid height error is 1.985× 10^-1 m at degree 120 based on GRACE Level 1B measured observation errors of the year 2007 published by the US Jet Propulsion Laboratory （JPL）, and the cumulative geoid height error is 5.825 × 10^-2 m at degree 360 using GRACE Follow-On orbital altitude 250 km and inter-satellite range 50 km. The matching relationship of accuracy indexes from GRACE Follow-On key payloads is brought forward, and the dependability of the combined error model is validated. Finally, the feasibility of high-accuracy and high-resolution global gravitational field estimation from GRACE Follow-On is demonstrated based on different satellite orbital altitudes.展开更多
Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate...Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate regularized matrix, and a singular decomposition method was used to generate regularization parameters. Numerical test results suggest that the regularized ambiguity float solution is more stable and reliable than the least-squares float solution. The mean square error matrix of the new method possesses a lower correlation than the variancecovariance matrix of the least-squares estimation. The size of the ambiguity search space is reduced and the search efficiency is improved. The success rate of the integer ambiguity searching process is improved significantly when the ambiguity resolution by using constraint equation method is used to determine the correct ambiguity integervector. The ambiguity resolution by using constraint equation method requires an initial input of the ambiguity float solution candidates which are obtained from the LAMBDA method in the new method. In addition, the observation time required to fix reliable integer ambiguities can he significantly reduced.展开更多
Based on the nearly diurnal resonance in the tidal gravity observations, the temporal variations in period of the Earth＇s free core nutation （FCN） are investigated by using the tidal gravity observations of 18-year...Based on the nearly diurnal resonance in the tidal gravity observations, the temporal variations in period of the Earth＇s free core nutation （FCN） are investigated by using the tidal gravity observations of 18-year duration recorded continu- ously with a superconducting gravimeter （SG） at Brussels. The effects of the global oceanic tide loading and local barometric pressure on the SG observations have been removed by using eleven high-precision global digital models of oceanic tides and barometric pressure measurements recorded simultaneously at the same site. The results indicate that there exist decade-scale variations in the FCN period. The results should be further confirmed by the measurements using other space-based geodetic techniques （such as the very long baseline interferometry） and the SG observations from globally distributed stations.展开更多
The method of extracting Green＇s function between stations from cross correlation has proven to be effective theoretically and experimentally. It has been widely applied to surface wave tomography of the crust and up...The method of extracting Green＇s function between stations from cross correlation has proven to be effective theoretically and experimentally. It has been widely applied to surface wave tomography of the crust and upmost mantle. However, there are still controversies about why this method works. Snieder employed stationary phase approximation in evaluating contribution to cross correlation function from scatterers in the whole space, and concluded that it is the constructive interference of waves emitted by the scatterers near the receiver line that leads to the emergence of Green＇s function. His derivation demonstrates that cross correlation function is just the convolution of noise power spectrum and the Green＇s function. However, his derivation ignores influence from the two stationary points at infinities, therefore it may fail when attenuation is absent. In order to obtain accurate noise-correlation function due to scatters over the whole space, we compute the total contribution with numerical integration in polar coordinates. Our numerical computation of cross correlation function indicates that the incomplete stationary phase approximation introduces remarkable errors to the cross correlation function, in both amplitude and phase, when the frequency is low with reasonable quality factor Q. Our results argue that the dis- tance between stations has to be beyond several wavelengths in order to reduce the influence of this inaccuracy on the applications of ambient noise method, and only the station pairs whose distances are above several （〉5） wavelengths can be used.展开更多
We analyzed the seismic waveforms from the December 26, 2004 Sumatra-Andaman earthquake recorded at broadband seismic stations in western Europe. Previous studies involving of the beam-forming technique and high frequ...We analyzed the seismic waveforms from the December 26, 2004 Sumatra-Andaman earthquake recorded at broadband seismic stations in western Europe. Previous studies involving of the beam-forming technique and high frequency analysis suggest that the earthquake ruptured with a duration of around 500 s. This very long duration makes P wave overlap with later arrivals such as PP wave, which follows P in about 200 s. Since P waves are crucial for modeling earthquake processes, we propose an iterative method to separate P and PP waveforms. The separated P waveform confirms a second large energy release around 300 s after the initial rupture. The iterative signal separation technique is particularly useful for mixed signals that are not independent and the number of recording stations far exceeds number of mixed signal sources.展开更多
In this paper,we develop an approach to study the effect of second-order ionospheric delay on GPS positioning based on the ionosphere-free combination (abbreviated to Lc) of GPS dual-frequency carrier phase observable...In this paper,we develop an approach to study the effect of second-order ionospheric delay on GPS positioning based on the ionosphere-free combination (abbreviated to Lc) of GPS dual-frequency carrier phase observables,in which the first-order ionospheric delay has been eliminated.GPS data from IGS WUHN tracking station during April 9 23,2003 is used to perform the above approach,and results show that the second-order ionospheric delay in GPS so-called ionosphere-free observables will result in the regular southward shift of GPS positioning solution.Additionally,the influence of the Earth's magnetic field on second-order ionospheric delay over the Wuhan area is discussed,and it indicates that the Earth's magnetic field is one of the main factors why second-order ionospheric delay results in the regular southward shift of GPS positioning solution over Wuhan area.展开更多
Sumatra-Andaman Large Earthquake on Dec. 26,2004 generated not only the Indian Ocean Tsunami but also the Earth's free oscillations (EFO). The signals of Earth's free oscillations were perfectly re-corded by t...Sumatra-Andaman Large Earthquake on Dec. 26,2004 generated not only the Indian Ocean Tsunami but also the Earth's free oscillations (EFO). The signals of Earth's free oscillations were perfectly re-corded by the superconducting gravimeter C0-32 at Wuhan station in China. After the pre-treatment and spectral analysis on the observational data from Wuhan station,we obtained more than ninety EFO modes including 42 fundamental modes,2 radial modes and 49 harmonic modes. On the basis of the discussions on some observed harmonic modes and abnormal splitting phenomena,we considered that the real rigidity might be lower than the theoretical prediction of PREM model in the inner core and however the anisotropy of compressive wave was brightly higher than the present estimations in the inner core. This suggested that the anisotropy of the inner core could be much more complicated than our present understanding,and there might be some new geophysical phenomena in the formation process of the inner core.展开更多
During the period when a GPS satellite,the Earth and the Sun are approximately collinear,the phenomenon of eclipsing the satellite occurs,when the satellite yaw attitude deviates from its nominal case,i.e. the body X-...During the period when a GPS satellite,the Earth and the Sun are approximately collinear,the phenomenon of eclipsing the satellite occurs,when the satellite yaw attitude deviates from its nominal case,i.e. the body X-axis cannot point towards the Sun for Block II&IIA or away from it for Block IIR satellites. The yaw attitude of the eclipsing satellites has a significant influence on both the satellite clock estimation at each International GNSS Service (IGS) Analysis Center (AC) and for users of the precise point positioning (PPP) implementations. It is known that,during the eclipsing periods,inconsistent yaw attitude models among the ACs will contribute to the errors of the IGS combined clock products. As for the PPP user,the influence of the eclipsing satellite is two-fold. First,as the satellite clocks are always kept fixed during PPP implementation,the above-mentioned problematic IGS clocks will inevitably be passed on to the PPP estimates. Second,the improper yaw attitude modeling of the eclipsing satellite will cause a correction bias exceeding 1 dm for the two kinds of attitude-related systematic errors,namely the phase wind-up and satellite antenna phase center offset,which will further deteriorate the accuracy of the PPP solutions. A yaw attitude model is introduced in this paper with the aim of improving the reliability of PPP solutions during the satellite eclipsing period.展开更多
Ocean tide loading(OTL)displacements of eight principal constituents at 12 sites in Hong Kong were determined using more than eight years of continuous GPS observations.Parameters of the OTL displacements were estimat...Ocean tide loading(OTL)displacements of eight principal constituents at 12 sites in Hong Kong were determined using more than eight years of continuous GPS observations.Parameters of the OTL displacements were estimated using daily GPS solutions obtained with precise point positioning(PPP)technique.The results were compared with predictions from seven latest global ocean tide models.Gravity measurements of OTL in Hong Kong were also used to validate the GPS results.The study shows that the results from both the GPS and gravity measurements agree best with the GOT4.7 and NAO99b models,when the K1 and K2 constituents are excluded.The agreements between the GPS and the model estimates are generally at sub-millimeter level in both the horizontal and vertical directions,except for S2,K2 and K1 constituents that have relatively larger errors.After removing the systematic biases between the GPS and the model estimates,the misfits of M2,S2,N2,O1,P1 and Q1 at all sites are within 0.5 and 1.0 mm in the horizontal and the vertical directions,respectively,whereas K1 and K2 show relatively larger misfits of up to 2.3 mm.Both the GPS and the gravity S2 estimates have large biases with unknown reasons when compared with the modeled values,which needs to be further investigated.The study demonstrates that GPS is capable of estimating the OTL displacements with the same accuracy as the model predictions,especially for coastal areas.展开更多
We consider the characteristics of long-term changes in non-tidal gravity and their implication to the local perturbations in barometric pressure and water storage and to the local vertical crustal movement using the ...We consider the characteristics of long-term changes in non-tidal gravity and their implication to the local perturbations in barometric pressure and water storage and to the local vertical crustal movement using the long-term continuous gravity observations from a superconducting gravimeter (SG) at Wu-han station,together with the co-site measurements from a Global Positioning System (GPS) receiver and an absolute gravimeter FG5. The observation results indicate that there are obvious seasonal variations in the long-term gravity changes measured with the SG. About 70 percent of the whole sea-sonal changes come from the contribution of the local disturbances in air pressure and water storage,while over 95 percent of the annual changes are attributed to the loading effects of these environmental perturbations. Due to the absence of direct measurements of the local water storage,especially those of the underground water,the global assimilating models of land water LaD (Land Dynamics) and GLDAS (Global Land Data Assimilation System) cannot virtually describe the real hydrologic distur-bances around the station. The resulting gravity changes,which are simulated theoretically by means of convolution integration of the loading Green’s functions and water models LaD and GLDAS,show significantly time delay of about 55 days from those measured with the SG. Compared with the meas-urements of the absolute gravity with the FG5,the long-term drift rate of the SG is determined as about 17.13 nms-2/a. From the co-site GPS measurements,it is found that the local crust is slowly subsiding at a rate of 3.71±0.16 mm/a,and the related gravity variation is estimated as 13.88±0.22 nms-2/a. In other words,the ratio of the changes in gravity and altitude related to the local vertical crustal movement is about -37.41 nms-2/cm. It implies that a considerable mass adjustment may be associated with the local vertical crustal movement,and its dynamic mechanism should be investigated further.展开更多
文摘The elasto-gravitational deformation response of the Earths solid parts to the perturbations of the pressure and gravity on the core-mantle boundary (CMB) and the solid inner core boundary (ICB), due to the dynamical behaviors of the fluid outer core (FOC), is discussed. The internal load Love numbers, which are formulized in a general form in this study, are employed to describe the Earths deformation. The preliminary reference Earth model (PREM) is used as an example to calculate the internal load Love numbers on the Earths surface, CMB and ICB, respectively. The characteristics of the Earths deformation variation with the depth and the perturbation periods on the boundaries of the FOC are also investigated. The numerical results indicate that the internal load Love numbers decrease quickly with the increasing degree of the spherical harmonics of the displacement and depend strongly on the perturbation frequencies, especially on the high frequencies. The results, obtained in this work, can be used to construct the boundary conditions for the core dynamics of the long-period oscillations of the Earths fluid outer core.
基金supported by Chinese Academy of Sciences grant kzcx2-yw-142 and Y009021002National Natural Science Foundation of China under grant No. 40974034+1 种基金NSF-EAR award 0944022NSF-OISE 0730154 at the University of Colorado at Boulder
文摘We apply ambient noise tomography to significant seismic data resources in a region including the northeastern Tibetan plateau, the Ordos block and the Sichuan basin. The seismic data come from about 160 stations of the provincial broadband digital seismograph networks of China. Ambient noise cross-correlations are performed on the data recorded between 2007 and 2009 and high quality inter-station Rayleigh phase velocity dispersion curves are obtained between periods of 6 s to 35 s. Resulting Rayleigh wave phase velocity maps possess a lateral resolution between 100 km and 200 kin. The phase velocities at short periods （〈20 s） are lower in the Sichuan basin, the northwest segment of the Ordos block and the Weihe graben, and outline sedimentary deposits. At intermediate and long periods （〉25 s）, strong high velocity anomalies are observed within the Ordos block and the Sichuan basin and low phase velocities are imaged in the northeastern Tibetan plateau, reflecting the variation of crustal thickness from the Tibetan plateau to the neighboring regions in the east. Crustal and uppermost mantle shear wave velocities vary strongly between the Tibetan plateau, the Sichuan basin and the Ordos block. The Ordos block and the Sichuan basin are dominated by high shear wave velocities in the crust and uppermost mantle. There is a triangle-shaped low velocity zone located in the northeastern Tibetan plateau, whose width narrows towards the eastern margin of the plateau. No low velocity zone is apparent beneath the Qinling orogen, suggesting that mass may not be able to flow eastward through the boundary between the Ordos block and the Sichuan basin in the crust and uppermost mantle.
基金supported by Chinese Academy of Sciences (Nos.KZCX2-YW-116 and KZCX2-YW-142)National Natural Science Foundation of China (No. 40974034)
文摘Based on the elastic dislocation theory, multilayered crustal model, and rupture model obtained by seismic waveform inversion, we calculated the coand post-seismic surface deformation and gravity changes caused by the Yushu M W 6.9 earthquake occurred on April 14, 2010. The observed GPS velocity field and gravity field in Yushu areas are disturbed by the coand post-seismic effects induced by Yushu earthquake, thus the theoretical coand post-seismic deformation and gravity changes will provide important modification for the background tectonic movement of Yushu and surrounding regions. The time relaxation results show that the influences of Yushu earthquake on Yushu and surrounding areas will last as long as 30 to 50 years. The maximum horizontal displacement, vertical uplift and settlement are about 1.96, 0.27 and 0.16 m, respectively, the maximal positive and negative value of gravity changes are 8.892×10-7 m·s-2 and -4.861×10-7 m·s-2 , respectively. Significant spatial variations can be found on the coand post-seismic effects： The co-seismic effect mainly concentrates in the region near the rupture fault, while viscoelastic relaxation mostly acts on the far field. Therefore, when using the geodetic data to research tectonic motion, we should not only consider the effect of co-seismic caused by earthquake, but also pay attention to the effect of viscoelastic relaxation.
基金Project(40771175)supported by the National Nature Science Foundation of China
文摘The buildings and structures of mines were monitored automatically using modern surveying technology.Through the analysis of the monitoring data,the deformation characteristics were found out from three aspects containing points,lines and regions,which play an important role in understanding the stable state of buildings and structures.The stability and deformation of monitoring points were analysed,and time-series data of monitoring points were denoised with wavelet analysis and Kalman filtering,and exponent function and periodic function were used to get the ideal deformation trend model of monitoring points.Through calculating the monitoring data obtained,analyzing the deformation trend,and cognizing the deformation regularity,it can better service mine safety production and decision-making.
基金financially supported by the National Natural Science Foundation of China(No.90814009)
文摘A three-dimensional crustal density model beneath North China is determined using P-wave traveltimes and gravity datum by sequential inversion method in this paper. To improve the uniqueness of the solution, we used a linear relationship between velocity and density to achieve mutual conversions and constraints between velocity difference and density difference. Algebraic reconstruction technique （ART） was used in density inversion, which highly improved the calculation speed comparing with common least squares method. The inversion results indicate that the crustal density beneath North China is extremely inhomogeneous with its distribution coherent with regional tectonics. The lower crust of Taihang mountain tectonic zone shows an obvious low density characteristic. We proposed that it may be an indicator of upwelling of mantle materials or partial melting of lowermost crust.
基金supported by NSFC Grant Nos 41021003,41174086,41074052key Program from Chinese Academy of Sciences Grant No KZZD-EW-TZ-05
文摘On April 20, 2013, an earthquake with mag- nitude 7.0 occurred in the southwest of the Longmenshan fault system in and around Lushan County, Sichuan Province, China. This devastating earthquake killed hun- dreds of people, injured 10 thousand others, and collapsed countless buildings. In order to analyze the potential risk of this big earthquake, we calculate the co- and post-seismic surface deformation and gravity changes of this event. In this work, a multilayered crustal model is designed, and the elastic dislocation theory is utilized to calculate the co- and post-seismic deformations and gravity changes. During the process, a rupture model obtained by seismic waveform inversion （Liu et al. Sci China Earth Sci 56（7）： 1187-1192, 2013） is applied. The time-dependent relaxation results show that the influences on Lushan and its surrounding areas caused by the Ms7.0 Lushan earthquake will last as long as 10 years. The maximum horizontal displacement, vertical uplift, and settlement are about 5 cm, 21.24 cm, and 0.16 m, respectively; the maximal positive and nega- tive values of gravity changes are 45 and -0.47 μGal, respectively. These results may be applied to evaluate the long-term potential risk caused by this earthquake and to provide necessary information for post-earthquake reconstruction.
基金jointly supported by Chinese Academy of Sciences(Nos.KZCX2-YW-116 and KZCX2-YW-142)National Natural Science Foundation of China (Nos.40974034and 40064004)
文摘We compare the factors which affect the movement of Tibetan Plateau by building three types of finite element models： an elastic materials （M-EC）, a continuous model composed by non-linear materials （M-PC）, and an elastic model with discontinuous fault movements （M-ET）. Both in M-ET and M-EC, the materials in Qiangtang and Lhasa block are elastic, and in M-ET, discontinuous movement of faults is considered for evaluating the effects of strike-slip faults. In model M-PC Druker-Prager plastic materials are used in Qiangtang and Lhasa block. Comparisons of the numerical simulation and the GPS observations show following characteristics： （1） Under present tectonic environment, short-term deformation of Tibetan Plateau can be simulated well by elastic models; （2） Discontinuous fault activities increase the lateral extrusion of the eastern part of Tibetan Plateau, reduce the stress field level in Qiangtang, Tarim and Qaidam blocks and strengthen the E-W extensional force in the east and the west parts of Qiangtang block; （3） Properties of plastic materials reduce the total stress field and the E-W extensional force, thus, the normal fault earthquakes in southern Tibet is mainly owed to the effect of active fault movement. Based on the numerical simulations we speculate that faults movement may play a more important role on the kinematic pattern of Tibetan Plateau than bulk properties.
基金Supported by the National Natural Science Foundation of China (No. 40874036,No. 90814009,No. 40774011)the National Science & Technology Major Project of China (No. 2008ZX05008-006-062)
文摘The crustal structure of Xiachayu-Gonghe geophysical profile in eastern Tibetan plateau is simulated with Bouguer anomaly corrected for sediments and lithosphere. The forward simulation shows that the thickness of upper crust in eastem Tibetan plateau is about 20 km, and the density is 2.78 × 10^3 kg/m^3. The bottom interface of middle crust changes from 30 km to 40 km, the density of middle crust is 2.89 × 10^3 kg/m^3. The materials with low density of 2.78 × 10^3 kg/m^3 exist in middle crust, and those with high density of 3.33 × 10^3 kg/m^3 exist at the bottom of middle crust between Wenquan and Tanggemu. The density is 3.10× 10^3 kg/m^3 in lower crust. The shallowest depth of Moho interface is about 56 km, and the deepest one is about 74 km, the undulation of interface is large, the deep Moho is located in Xiachayu, Chayu, Nujiang, and Wenquan. The crustal density of eastern Tibetan plateau is larger than that of central section; the low velocity layers are located in middle crust and bottom in eastern Tibetan plateau and at the bottom of the upper crust in the central plateau.
基金Supported by the Natural Science Foundation of Hubei Province of China under Grant No 2009CDB187, the Hubei Province Key Laboratory of Refractories and Ceramics Ministry-Province Jointly-Constructed Cultivation Base for State key Laboratory （No G201009）, the Open Research Fund Program of Hubei Province Key Laboratory of Systems Science in Metallurgical Process （No C201019）, and the National Natural Science Foundation of China under Grant No 10932010.
文摘We present a three-dimensional geometry model for tortuosity of streamlines in porous media with randomly placed cylindrical particles. The proposed model is expressed as functions of porosity and geometrical parameters with no empirical constant. This might be helpful for understanding the physical mechanism for tortuosity of streamlines in three-dimensional porous media. The model predictions are found to be in good agreement with the experimental data available.
基金supported by the National Natural Science Foundation of China (Grant No 40674038)the Funds of the Chinese Academy of Sciences for Key Topics in Innovation Engineering (Grant Nos KZCX2-YW-143 and KZCX2-YW-202)+1 种基金the National High Technology Research and Development Program of China (863) (Grant Nos 2009AA12Z138 and 2006AA09Z153)the Grant-in-Aid for Scientific Research of Japan (Grant No B19340129)
文摘Firstly, the new combined error model of cumulative geoid height influenced by four error sources, including the inter-satellite range-rate of an interferometric laser （K-band） ranging system, the orbital position and velocity of a global positioning system （GPS） receiver and non-conservative force of an accelerometer, is established from the perspectives of the power spectrum principle in physics using the semi-analytical approach. Secondly, the accuracy of the global gravitational field is accurately and rapidly estimated based on the combined error model; the cumulative geoid height error is 1.985× 10^-1 m at degree 120 based on GRACE Level 1B measured observation errors of the year 2007 published by the US Jet Propulsion Laboratory （JPL）, and the cumulative geoid height error is 5.825 × 10^-2 m at degree 360 using GRACE Follow-On orbital altitude 250 km and inter-satellite range 50 km. The matching relationship of accuracy indexes from GRACE Follow-On key payloads is brought forward, and the dependability of the combined error model is validated. Finally, the feasibility of high-accuracy and high-resolution global gravitational field estimation from GRACE Follow-On is demonstrated based on different satellite orbital altitudes.
文摘Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate regularized matrix, and a singular decomposition method was used to generate regularization parameters. Numerical test results suggest that the regularized ambiguity float solution is more stable and reliable than the least-squares float solution. The mean square error matrix of the new method possesses a lower correlation than the variancecovariance matrix of the least-squares estimation. The size of the ambiguity search space is reduced and the search efficiency is improved. The success rate of the integer ambiguity searching process is improved significantly when the ambiguity resolution by using constraint equation method is used to determine the correct ambiguity integervector. The ambiguity resolution by using constraint equation method requires an initial input of the ambiguity float solution candidates which are obtained from the LAMBDA method in the new method. In addition, the observation time required to fix reliable integer ambiguities can he significantly reduced.
基金supported by the Key Project of the Knowledge Innovation of Chinese Academy of Sciences (Grant No. KZCX2-YW-133, KZCX2-YW-Q08-2)National Natural Science Foun-dation of China (Grant No. 40874038 and 40730316)
文摘Based on the nearly diurnal resonance in the tidal gravity observations, the temporal variations in period of the Earth＇s free core nutation （FCN） are investigated by using the tidal gravity observations of 18-year duration recorded continu- ously with a superconducting gravimeter （SG） at Brussels. The effects of the global oceanic tide loading and local barometric pressure on the SG observations have been removed by using eleven high-precision global digital models of oceanic tides and barometric pressure measurements recorded simultaneously at the same site. The results indicate that there exist decade-scale variations in the FCN period. The results should be further confirmed by the measurements using other space-based geodetic techniques （such as the very long baseline interferometry） and the SG observations from globally distributed stations.
基金supported by the National Natural Science Foundation of China (No. 40674027)CAS outstanding 100 research program,MOST program 2007FY220100
文摘The method of extracting Green＇s function between stations from cross correlation has proven to be effective theoretically and experimentally. It has been widely applied to surface wave tomography of the crust and upmost mantle. However, there are still controversies about why this method works. Snieder employed stationary phase approximation in evaluating contribution to cross correlation function from scatterers in the whole space, and concluded that it is the constructive interference of waves emitted by the scatterers near the receiver line that leads to the emergence of Green＇s function. His derivation demonstrates that cross correlation function is just the convolution of noise power spectrum and the Green＇s function. However, his derivation ignores influence from the two stationary points at infinities, therefore it may fail when attenuation is absent. In order to obtain accurate noise-correlation function due to scatters over the whole space, we compute the total contribution with numerical integration in polar coordinates. Our numerical computation of cross correlation function indicates that the incomplete stationary phase approximation introduces remarkable errors to the cross correlation function, in both amplitude and phase, when the frequency is low with reasonable quality factor Q. Our results argue that the dis- tance between stations has to be beyond several wavelengths in order to reduce the influence of this inaccuracy on the applications of ambient noise method, and only the station pairs whose distances are above several （〉5） wavelengths can be used.
基金supported by CAS fund(KZCX2-YW-116-1)National Natural Science Foundation of China(40821160549 and 41074032)China Earthquake Administration fund(200808078)
文摘We analyzed the seismic waveforms from the December 26, 2004 Sumatra-Andaman earthquake recorded at broadband seismic stations in western Europe. Previous studies involving of the beam-forming technique and high frequency analysis suggest that the earthquake ruptured with a duration of around 500 s. This very long duration makes P wave overlap with later arrivals such as PP wave, which follows P in about 200 s. Since P waves are crucial for modeling earthquake processes, we propose an iterative method to separate P and PP waveforms. The separated P waveform confirms a second large energy release around 300 s after the initial rupture. The iterative signal separation technique is particularly useful for mixed signals that are not independent and the number of recording stations far exceeds number of mixed signal sources.
基金Supported by Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KZCX2-YW-202)National Natural Science Foundation of China (Grant No. 40674038)Program of Key Laboratory of Space Environment and Geodesy of Ministry of Education of China
基金supported by the National Science Foundation for Distinguished Young Scholars (Grant No. 40625013)the National Natural Science Foundation of China (Grant No. 40890160)+1 种基金the National High-Tech Research and Development Program of China (Grant No. 2007AA12Z311)the Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No.0909141014)
文摘In this paper,we develop an approach to study the effect of second-order ionospheric delay on GPS positioning based on the ionosphere-free combination (abbreviated to Lc) of GPS dual-frequency carrier phase observables,in which the first-order ionospheric delay has been eliminated.GPS data from IGS WUHN tracking station during April 9 23,2003 is used to perform the above approach,and results show that the second-order ionospheric delay in GPS so-called ionosphere-free observables will result in the regular southward shift of GPS positioning solution.Additionally,the influence of the Earth's magnetic field on second-order ionospheric delay over the Wuhan area is discussed,and it indicates that the Earth's magnetic field is one of the main factors why second-order ionospheric delay results in the regular southward shift of GPS positioning solution over Wuhan area.
基金the National Natural Science Foundation of China (Grant Nos. 40404005 and 40374029)the Excellent Prize of President ScholarshipHundred Talents Program of the Chinese Academy of Sciences
文摘Sumatra-Andaman Large Earthquake on Dec. 26,2004 generated not only the Indian Ocean Tsunami but also the Earth's free oscillations (EFO). The signals of Earth's free oscillations were perfectly re-corded by the superconducting gravimeter C0-32 at Wuhan station in China. After the pre-treatment and spectral analysis on the observational data from Wuhan station,we obtained more than ninety EFO modes including 42 fundamental modes,2 radial modes and 49 harmonic modes. On the basis of the discussions on some observed harmonic modes and abnormal splitting phenomena,we considered that the real rigidity might be lower than the theoretical prediction of PREM model in the inner core and however the anisotropy of compressive wave was brightly higher than the present estimations in the inner core. This suggested that the anisotropy of the inner core could be much more complicated than our present understanding,and there might be some new geophysical phenomena in the formation process of the inner core.
基金supported by the National Science Foundation for Distinguished Young Scholars (40625013)the National Natural Science Foundation of China (40890160 and 40874009)the National High-Tech Research and Development Program of China (2007AA12Z311,2007A-A12-Z305)
文摘During the period when a GPS satellite,the Earth and the Sun are approximately collinear,the phenomenon of eclipsing the satellite occurs,when the satellite yaw attitude deviates from its nominal case,i.e. the body X-axis cannot point towards the Sun for Block II&IIA or away from it for Block IIR satellites. The yaw attitude of the eclipsing satellites has a significant influence on both the satellite clock estimation at each International GNSS Service (IGS) Analysis Center (AC) and for users of the precise point positioning (PPP) implementations. It is known that,during the eclipsing periods,inconsistent yaw attitude models among the ACs will contribute to the errors of the IGS combined clock products. As for the PPP user,the influence of the eclipsing satellite is two-fold. First,as the satellite clocks are always kept fixed during PPP implementation,the above-mentioned problematic IGS clocks will inevitably be passed on to the PPP estimates. Second,the improper yaw attitude modeling of the eclipsing satellite will cause a correction bias exceeding 1 dm for the two kinds of attitude-related systematic errors,namely the phase wind-up and satellite antenna phase center offset,which will further deteriorate the accuracy of the PPP solutions. A yaw attitude model is introduced in this paper with the aim of improving the reliability of PPP solutions during the satellite eclipsing period.
基金supported by the Research Grants Council(RGC)of the Hong Kong Special Administrative Region(Grant Nos.PolyU5157/05E,PolyU5161/06E)the Scientific Research Foundation of The Hong Kong Polytechnic University(Grant No.GY-F81)the Fundamental Research Funds for the Central Universi-ties(Grant Nos.SWJTU09BR029,SWJTU09BR031)
文摘Ocean tide loading(OTL)displacements of eight principal constituents at 12 sites in Hong Kong were determined using more than eight years of continuous GPS observations.Parameters of the OTL displacements were estimated using daily GPS solutions obtained with precise point positioning(PPP)technique.The results were compared with predictions from seven latest global ocean tide models.Gravity measurements of OTL in Hong Kong were also used to validate the GPS results.The study shows that the results from both the GPS and gravity measurements agree best with the GOT4.7 and NAO99b models,when the K1 and K2 constituents are excluded.The agreements between the GPS and the model estimates are generally at sub-millimeter level in both the horizontal and vertical directions,except for S2,K2 and K1 constituents that have relatively larger errors.After removing the systematic biases between the GPS and the model estimates,the misfits of M2,S2,N2,O1,P1 and Q1 at all sites are within 0.5 and 1.0 mm in the horizontal and the vertical directions,respectively,whereas K1 and K2 show relatively larger misfits of up to 2.3 mm.Both the GPS and the gravity S2 estimates have large biases with unknown reasons when compared with the modeled values,which needs to be further investigated.The study demonstrates that GPS is capable of estimating the OTL displacements with the same accuracy as the model predictions,especially for coastal areas.
基金the Key Project of the Knowledge Innovation of Chinese Academy of Sciences (Grant No.KZCX2-YW-133)the National Natural Science Foundation of China (Grant Nos.40574034 and 40730316)
文摘We consider the characteristics of long-term changes in non-tidal gravity and their implication to the local perturbations in barometric pressure and water storage and to the local vertical crustal movement using the long-term continuous gravity observations from a superconducting gravimeter (SG) at Wu-han station,together with the co-site measurements from a Global Positioning System (GPS) receiver and an absolute gravimeter FG5. The observation results indicate that there are obvious seasonal variations in the long-term gravity changes measured with the SG. About 70 percent of the whole sea-sonal changes come from the contribution of the local disturbances in air pressure and water storage,while over 95 percent of the annual changes are attributed to the loading effects of these environmental perturbations. Due to the absence of direct measurements of the local water storage,especially those of the underground water,the global assimilating models of land water LaD (Land Dynamics) and GLDAS (Global Land Data Assimilation System) cannot virtually describe the real hydrologic distur-bances around the station. The resulting gravity changes,which are simulated theoretically by means of convolution integration of the loading Green’s functions and water models LaD and GLDAS,show significantly time delay of about 55 days from those measured with the SG. Compared with the meas-urements of the absolute gravity with the FG5,the long-term drift rate of the SG is determined as about 17.13 nms-2/a. From the co-site GPS measurements,it is found that the local crust is slowly subsiding at a rate of 3.71±0.16 mm/a,and the related gravity variation is estimated as 13.88±0.22 nms-2/a. In other words,the ratio of the changes in gravity and altitude related to the local vertical crustal movement is about -37.41 nms-2/cm. It implies that a considerable mass adjustment may be associated with the local vertical crustal movement,and its dynamic mechanism should be investigated further.
