GOCE卫星监测的中国区域重力异常

介绍了利用最新一代重力卫星GOCE(Gravity field and steady-state Ocean Circulation Explorer)重力场模型数据计算重力异常的原理和方法,采用最新发布的GOCE重力场模型数据(2009年10月-2010年7月)计算了中国区域(70°～130°E,15°～55°N)的重力异常,为了进一步分析重力异常与区域地质构造及地震活动性的对应关系,将计算结果与中国区域地形、地震活动区域等资料进行了对比分析,结果表明利用GOCE重力场模型数据计算的重力异常能够较好的反映区域地质构造分布特征,而且强震震中通常位于重力异常变化剧烈的高梯度带上.     

Gravity field contribution analysis of GOCE gravitational gradient components

A gravity field model is computed from the four accurate gravitational gradient components of GOCE (Gravity field and steady-state Ocean Circulation Explorer), combined with the analysis of the kinematic orbits, and some moderate constraint (or stabilization) in the polar areas where no observation from GOCE is available due to the orbit geometry. The normal matrix of each component is computed individually in order to study its contribution to the combined solution. The results show that the contribution of V_zz is the largest, with an average value of 32.74% of the total solution; the second and the third largest are V_zz and V_yy, with average values of 28.04% and 26.08%, respectively; the component V_xz contributes 11.81%. Validation with external data shows that each component has its characteristic value and that the information content of the component V_xz is not negligible and should be included for gravity field recovery. The orbit part as derived from high-low satellite-to-satellite tracking (SST-hl) to the GPS contributes mostly to the coefficients below degree and order (d/o) 20, and to non-zonal coefficients from d/o 20 to 80. The mean value of the contribution of the polar stabilization is the smallest with a value of 0.22%, nevertheless it is important. In addition to the contribution analysis in terms of the normal matrices, each individual component of the gradiometer has been combined with SST and polar stabilization, to give a set of single component gravity field models. These partially combined solutions are compared to the fully combined solution in terms of geoid differences. They show that the partially combined solution with V_zz is closest to the complete solution. Even closer is a combination with V_xx and V_yy. In addition to the GOCE-only solution, a GOCE-GRACE (Gravity Recovery And Climate Experiment) combined gravity field model is derived and the information content of GOCE and an available set of normal equations of GRACE are investigated. Results show that, as expected, GRACE dominates the solution below degree 90 and GOCE above degree 140.     

A comparison of GOCE gravitational models with EGM2008

Four new gravity field models from GOCE, two of them combined with GRACE, are compared here with EGM2008. The objectives are to look into the differences in consecutive ranges of the spherical harmonic expansion globally as well as in selected geographical regions and in the regions of the various data sources used for EGM2008. In general, GOCE is able to contribute to improved global gravity models in the spherical harmonic range between 120 and 200 (and above). The agreement between EGM2008 and the GOCE models is very good in well-surveyed regions such as North America, Europe and Australia, with geoid RMS-differences on the order of 4–6 cm. In other regions, where the surface gravity data available for the development of EGM2008 were poor, such as South America, Africa, South-East Asia or China the RMS-differences are on a level of 30 cm. Here GOCE leads to a significant improvement. These findings are confirmed by the analysis of the areas of the various EGM2008 data sources. In the regions of the so-called “fill-in” data of EGM2008 RMS-geoid height differences are high. In Antarctica GOCE also gives important improvements in terms of spatial resolution and accuracy. In general, the agreement between EGM2008 and the GOCE-models up to degree and order (d/o) 200 is good, with a global (excluding the polar gaps of GOCE orbits, throughout) geoid difference RMS of 11 cm, in the ocean areas 8 cm and 20 cm in the continental areas. GOCE models are better suited for ocean circulation studies because no prior ocean information enters into the data reduction process, as it is the case when deducing gravity anomalies from an altimetric mean sea surface. On the other hand, the good consistency between GOCE-models and EGM2008 in ocean areas very likely indicates that the influence of ocean circulation information on EGM2008 is rather small. The four tested GOCE models behave similarly except at the highest latitudes where GOCE lacks data due to its orbit inclination of 96.5° and some form of regularization which has to be applied.     

The earth-moon-sun disturbing function in satellite orbit dynamics

am ¶rt; a uu u m nu uu n ¶rt;uu um nmua u aumau n m -a-. ¶rt;a ¶rt; 4- n¶rt;a um.     

Investigation of the instability of graveyard orbits due to the earth's gravitational perturbations

Summary The Banach theorem is applied to the Lagrange planetary equation for the semimajor axis of a geostationary satellite orbit to estimate the stability of near-geostationary satellite orbits. To achieve a graveyard (disposal) orbit, which will not interfere (=cross) the initial geostationary orbit, the geostationary semi-major axis a_g have to be increased at least by 50 km. Numerical results for a variety of graveyard orbits show that the increase of a_g by about 100 km will yield sufficiently stable orbits (accounting for the Earth's gravitational perturbations only) during the next 150 years.Dedicated to the 75th Birthday of Professor Academician Tibor Kolbenheyer     

基于时空域混合法利用Kaula正则化精确和快速解算GOCE地球重力场

为了研究卫星重力梯度技术对中高频地球重力场反演精度的影响,本文基于时空域混合法,利用Kaula正则化反演了250阶GOCE地球重力场.模拟结果表明:第一,时空域混合法是精确和快速求解高阶地球重力场的有效方法;第二,Kaula正则化是降低正规阵病态性的重要方法;第三,基于改进的预处理共轭梯度迭代法可快速求解大型线性方程组...     

Upward continuation of Dome-C airborne gravity and comparison with GOCE gradients at orbit altitude in east Antarctica

An airborne gravity campaign was carried out at the Dome-C survey area in East Antarctica between the 17th and 22nd of January 2013, in order to provide data for an experiment to validate GOCE satellite gravity gradients. After typical filtering for airborne gravity data, the cross-over error statistics for the few crossing points are 11.3 mGal root mean square (rms) error, corresponding to an rms line error of 8.0 mGal. This number is relatively large due to the rough flight conditions, short lines and field handling procedures used. Comparison of the airborne gravity data with GOCE RL4 spherical harmonic models confirmed the quality of the airborne data and that they contain more high-frequency signal than the global models. First, the airborne gravity data were upward continued to GOCE altitude to predict gravity gradients in the local North-East-Up reference frame. In this step, the least squares collocation using the ITGGRACE2010S field to degree and order 90 as reference field, which is subtracted from both the airborne gravity and GOCE gravity gradients, was applied. Then, the predicted gradients were rotated to the gradiometer reference frame using level 1 attitude quaternion data. The validation with the airborne gravity data was limited to the accurate gradient anomalies (T_XX, T_YY, T_ZZ and T_XZ) where the long-wavelength information of the GOCE gradients has been replaced with GOCO03s signal to avoid contamination with GOCE gradient errors at these wavelengths. The comparison shows standard deviations between the predicted and GOCE gradient anomalies T_XX, T_YY, T_ZZ and T_XZ of 9.9, 11.5, 11.6 and 10.4 mE, respectively. A more precise airborne gravity survey of the southern polar gap which is not observed by GOCE would thus provide gradient predictions at a better accuracy, complementing the GOCE coverage in this region.     

DEMETER卫星观测到的云南普洱地震前的电离层扰动

北京时间2007年6月3日5时34分56秒（2007年6月2日21时34分56秒UT）在中国西南部的云南省普洱市宁洱县(23.0deg;N, 101.1deg;E)发生了MS6.4地震, 本文分析了DEMETER数据中心从此次地震前10天（5月23日）至地震当天共11天的数据, 寻找可能与此次地震有关的电离层前兆. DEMETER卫星在此期间有284个半轨,本文选择研究了飞过震中周围约1 888 km范围上空的29个半轨（在此范围内卫星到震中的最大距离为2 000 km）. 其中,在29个半轨中的7个半轨中发现了7次异常事件,它们的电场频谱、电子和离子的温度及密度与周围空间相比表现了同步的变化. 其中5次事件出现在研究的区域内,即震中周围1 888 km范围的空域中, 而另外两次事件则远离研究的区域. 分析表明这5次事件可能与普洱地震有关.      

利用GOCE和GRACE卫星观测数据确定静态重力场模型

高精度静态卫星重力场模型在全球海洋环流研究、全球/区域数字高程基准面确定等领域有重要应用,本文研究仅利用GOCE卫星和联合GRACE卫星观测数据确定高精度高阶次静态重力场模型.利用GOCE卫星全周期高精度引力梯度分量(V_xx、V_yy、V_zz和V_xz)观测值基于直接最小二乘法构建300阶次的SGG(Satellite Gravity Gradiometry)法方程,并利用卫星跟踪卫星观测值基于点域加速度法构建130阶SST(Satellite-to-Satellite Tracking)法方程,然后利用方差分量估计联合SGG和SST法方程确定300阶次纯GOCE卫星重力场模型GOSG02S.利用全周期GRACE观测数据由动力学方法解算了180阶次的SWPU-GRACE2021S模型,并将其对应法方程与GOCE卫星法方程联合解算了GRACE和GOCE的联合模型WHU-SWPU-GOGR2022S.分别基于XGM2019模型和GPS水准数据对本文解算的三个模型GOSG02S、SWPU-GRACE2021S和WHU-SWPU-GOGR2022S在频域和空域进行了精度分析,结果表明,GOSG02S和WHU-SWPU-GOGR2022S模型与GO_CONS_GCF_2_DIR_R6、GO_CONS_GCF_2_TIM_R6、GO_CONS_GCF_2_SPW_R5、GOCO06s和Tongji-GMMG2021S等使用了GOCE卫星全周期数据的模型精度相当,精度差异基本都在毫米量级;SWPU-GRACE2021S模型在160阶次之前与国际主流GRACE卫星重力场模型ITSG-Grace2018s和Tongji-Grace02s精度相当.

     

基于B-spline和正则化算法的低轨卫星轨道平滑

本文提出了一个利用纯几何轨道和力模型的新算法来计算精确且相对平滑的卫星轨道. 该法将一个纯几何轨道表达为一个B-spline的线性组合,线性组合的系数可以由最小二乘法估计获得. 力模型通过计算加速度来附加约束. 为了平衡几何轨道的点位误差和加速度的不精确,一个基于“广义交互确认（GCV,generalized cross-validation）”的正则化算法运用其中. 由于B-spline的本地控制性,该方法的计算效率相当高. 本文的数值分析表明了该法的有效性. 模拟计算的结论是：带加速度约束较不带加速度约束的平滑效果好. 力模型越精确,平滑的轨道就越精确. 三个月的CHAMP实测轨道数据处理结果表明,平滑后的轨道改进了重力场模型.     

The global stress field in the lithosphere obtained from the satellite gravitational harmonics

In this paper, the lithosphere is considered to be a homogeneous elastic spherical shell for the sake of simplicity and the stress equations for the base of the lithosphere are taken as boundary conditions. Then the stress equations are obtained for use in the computation of the stress field in the lithosphere with the satellite gravitational harmonic coefficients. The 5 × 5° global stress field in the lithosphere is computed from harmonics of 2–30°. The directions of principal stresses of this stress field agree favourably with the directions of principal stresses indicated by mid-plate earthquake mechanisms, in situ stress measurements and sensitive geological features. This result indicates that the drag forces exerted on the base of the lithosphere, due to gravitational mantle convection, may be among the driving forces of plate motion and a major source for the stress field in the lithosphere.     

行星内部地震波及其对卫星轨道的影响

本文指出地震学在天文和行星学科里的重要作用.我们主要介绍最近提出的“潮汐—地震波共振”(tidal-seismic resonance)效应,并且讨论它对卫星轨道演化的作用.当在同步轨道以下周期运动的卫星引起的引潮力的频率和行星内部自由震荡频率吻合时,就会发生潮汐—地震波共振.此时,行星内部的地震波将被激发并引起行星表面的显著位移.升高和下降的地面会对卫星产生一个力矩从而使得卫星轨道下降.因为潮汐共振引起的动态地面位移可以比单纯引潮力引起的位移大两个数量级,所以潮汐共振会显著加速卫星下降速率.我们用我们开发的三维地震波场模拟程序AstroSeis数值计算了潮汐—地震波共振对轨道的影响,进而推测这一共振效应可能对行星早期吸积速度有显著影响.另外,因为行星内部的Q值和S波的波速对潮汐共振影响很大,未来研究微重力环境下的小行星或陨石内部地震波的速度和Q值对研究行星演化和太阳系的形成至关重要.     

基于轨道扰动引力谱分析的方法确定低低观测重力卫星反演重力场的空间分辨率

基于低低卫-卫跟踪重力卫星的轨道特性,从垂直和水平两个方向计算了重力卫星高空扰动引力,并根据其谱特性及星载加速度的测量噪声水平分析了重力卫星能反演重力场的阶数.利用EGM96重力场模型分别计算了400 km、450 km和500 km 轨道高度处重力卫星受到的扰动引力谱及扰动引力谱的平均量级,分析其垂直特性表明:在三个轨道高度处能分别能反演150、140和130阶的重力场模型.利用两颗同轨重力卫星相距220 km的特性,计算了400 km、450 km和500 km 轨道高度处纬度相差2°的两颗卫星纬向扰动引力差,即扰动引力水平分量,分析其谱特性,表明:重力卫星能反演至117阶的地球重力场模型.     

Parameters of traveling ionospheric disturbances estimated from satellite beacon observations in low-earth orbit

Based on beacon sounding the ionosphere using coherent signals of low-orbiting navigation satellites, the following parameters of medium-scale traveling ionospheric disturbances have been determined: the effective thickness of the atmospheric waveguide and the altitude of its longitudinal axis, the horizontal spatial period, the velocity, the maximum amplitudes of disturbances, and the inclination of the disturbance phase front. It has been found that the horizontal spatial period of studied disturbances increases with increasing distance to an initiating impact location and its delay. In this case the minimum value of the horizontal spatial period always exceeds ∼150 km, and the disturbance amplitude increases with increasing horizontal period and distance from the initiating disturbance location and with decreasing altitude of the main ionization maximum. All disturbance parameters are independent of the initiating impact nature. It has been found that disturbances with relative amplitudes of 0.1–0.7 are often observed. Disturbances that are registered more often travel southward, and their absolute velocities are 7–60 m/s.     

Directions of axes of the Earth's ellipsoid of inertia derived from satellite orbit dynamics

Studia Geophysica et Geodaetica - Обсуж¶rt;aеmся 3 возможносmь...     

基于引力位系数相对权重的卫星重力场测量分析

卫星重力场测量已成为最有效的全球重力场测量手段.本文结合典型的重力卫星和重力卫星研究计划,分析了卫星重力测量的三种原理,并基于各阶位系数的相对权重讨论了各种原理的应用优势.分析可知,卫星受摄轨道适用于恢复长波重力场,低轨星间距离变化率适用于恢复中长波重力场,重力梯度适用于恢复中短波重力场.针对中长波高精度重力场测量的需要,设计了综合获取低轨星间距离变化率与受摄轨道的重力卫星方案,该方案由两组内编队组成星星跟踪复合编队,轨道高度为250km,星间距离为50～100km.     

Centimeter-level precise orbit determination for the HY-2A satellite using DORIS and SLR tracking data

The HY-2A satellite is the first ocean dynamic environment monitoring satellite of China. Centimeter-level radial accuracy is a fundamental requirement for its scientific research and applications. To achieve this goal, we designed the strategies of precise orbit determination (POD) in detail. To achieve the relative optimal orbit for HY-2A, we carried out POD using DORIS-only, SLR-only, and DORIS + SLR tracking data, respectively. POD tests demonstrated that the consistency level of DORIS-only and SLR-only orbits with respect to the CNES orbits were about 1.81 cm and 3.34 cm in radial direction in the dynamic sense, respectively. We designed 6 cases of different weight combinations for DORIS and SLR data, and found that the optimal relative weight group was 0.2 mm/s for DORIS and 15.0 cm for SLR, and RMS of orbit differences with respect to the CNES orbits in radial direction and three-dimensional (3D) were 1.37 cm and 5.87 cm, respectively. These tests indicated that the relative radial and 3D accuracies computed using DORIS + SLR data with the optimal relative weight set were obviously higher than those computed using DORIS-only and SLR-only data, and satisfied the requirement of designed precision. The POD for HY-2A will provide the invaluable experience for the following HY-2B, HY-2C, and HY-2D satellites.     

卫星轨道参数对LEO－LEO掩星事件数量及分布影响的模拟研究

文章对LEO－LEO卫星掩星事件的特性进行了深入讨论,通过仿真计算定量讨论了轨道参数对掩星事件次数及分布的影响.结果指出：1)对一年掩星事件数量及其纬度分布影响大的轨道参数依次是轨道倾角、轨道高度、升交点赤经和近地点角距;而这些参数对一年掩星事件的经度分布影响都不大;2)发射卫星和接收卫星倾角之间的相互关系对掩星事件数量和纬度分布影响最大,两卫星倾角接近互补时掩星事件最多, 两颗LEO卫星中与90°相差较大的倾角决定了LEO－LEO掩星事件的纬度分布范围;3)发射卫星和接收卫星处于同一轨道高度时,LEO－LEO掩星事件数取得极大值,发射卫星和接收卫星都采用600 km的轨道高度,一年LEO－LEO掩星事件数最多;4) 发射卫星和接收卫星升交点赤经(RAAN)的相互关系对LEO－LEO掩星事件数量影响很大,发射卫星和接收卫星升交点赤经之差为120°或240°时掩星事件数量最多,RAAN对掩星事件纬度和经度分布影响不大;5)近地点角距对掩星事件数量和分布影响都不大.该研究结果对LEO－LEO掩星探测卫星的轨道设计有一定的参考作用.