Impact of a novel hybrid accelerometer on satellite gravimetry performance

The state-of-the-art electrostatic accelerometers (EA) used for the retrieval of non-gravitational forces acting on a satellite constitute a core component of every dedicated gravity field mission. However, due to their difficult-to-control thermal drift in the low observation frequencies, they are also one of the most limiting factors of the achievable performance of gravity recovery. Recently, a hybrid accelerometer consisting of a regular EA and a novel cold atom interferometer (CAI) that features a time-invariant observation stability and constantly recalibrates the EA has been developed in order to remedy this major drawback. In this paper we aim to assess the value of the hybrid accelerometer for gravity field retrieval in the context of GRACE-type and Bender-type missions by means of numerical closed-loop simulations where possible noise specifications of the novel instrument are considered and different components of the Earth’s gravity field signal are added subsequently. It is shown that the quality of the gravity field solutions is mainly dependent on the CAI’s measurement accuracy. While a low CAI performance (10^?8 to 10^?9?m/s²/Hz^1/2) does not lead to any gains compared to a stand-alone EA, a sufficiently high one (10^?11?m/s²/Hz^1/2) may improve the retrieval performance by over one order of magnitude. We also show that improvements which are limited to low-frequency observations may even propagate into high spherical harmonic degrees. Further, the accelerometer performance seems to play a less prominent role if the overall observation geometry is improved as it is the case for a Bender-type mission. The impact of the accelerometer measurements diminishes further when temporal variations of the gravity field are introduced, pointing out the need for proper de-aliasing techniques. An additional study reveals that the hybrid accelerometer is – contrary to a stand-alone EA – widely unaffected by scale factor instabilities.     

Evaluation of uncertainty in gravity wave potential energy calculations through GPS radio occultation measurements

The application of the Global Positioning System (GPS) radio occultation (RO) method to the atmosphere enables the determination of height profiles of temperature, among other variables. From these measurements, gravity wave activity is usually quantified by calculating the potential energy through the integration of the ratio of perturbation and background temperatures between two given altitudes in each profile. The uncertainty in the estimation of wave activity depends on the systematic biases and random errors of the measured temperature, but also on additional factors like the selected vertical integration layer and the separation method between background and perturbation temperatures. In this study, the contributions of different parameters and variables to the uncertainty in the calculation of gravity wave potential energy in the lower stratosphere are investigated and quantified. In particular, a Monte Carlo method is used to evaluate the uncertainty that results from different GPS RO temperature error distributions. In addition, our analysis shows that RO data above 30 km height becomes dubious for gravity waves potential energy calculations.     

Modeling and analysis of Earth’s gravity field measurement performance by inner-formation flying system

The Earth’s gravity field can be measured with high precision by constructing the purely gravitational orbit of the inner-satellite in Inner-formation Flying System (IFS), which is independently proposed by Chinese scholars and offers a new way to carry out gravity field measurement by satellite without accelerometers. In IFS, for the purpose of quickly evaluating the highest degree of recovered gravity field model and geoid error as well as analyzing the influence of system parameters on gravity field measurement, an analytical formula was established by spectral analysis method. The formula can reflect the analytical relationship between gravity field measurement performance and system parameters such as orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and total measurement time. This analytical formula was then corrected by four factors introduced from numerical simulation of IFS gravity field measurement. By comparing computation results from corrected analytical formula and the actual gravity field measurement performance by CHAMP, the correctness and rationality of this analytical formula were verified. Based on this analytical formula, the influences of system parameters on IFS gravity field measurement were analyzed. It is known that gravity field measurement performance is a monotone decreasing function of orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and the reciprocal of total measurement time. There is a match relationship between the inner-satellite orbit determination error and residual disturbances, in other words, the change rate of gravity field measurement performance with one of them is seriously restricted by their relative size. The analytical formula can be used to quantitatively evaluate gravity field measurement performance fast and design IFS parameters optimally. It is noted that the analytical formula and corresponding conclusions are applied to any gravity satellite which measures gravity field by satellite perturbation orbit.     

Estimating low resolution gravity fields at short time intervals to reduce temporal aliasing errors

The Gravity Recovery and Climate Experiment (GRACE) satellite mission has been estimating temporal changes in the Earth’s gravitational field since its launch in 2002. While it is not yet fully resolved what the limiting source of error is for GRACE, studies on future missions have shown that temporal aliasing errors due to undersampling signals of interest (such as hydrological variations) and errors in atmospheric, ocean, and tide models will be a limiting source of error for missions taking advantage of improved technologies (flying drag-free with a laser interferometer). This paper explores the option of reducing the effects of temporal aliasing errors by directly estimating low degree and order gravity fields at short time intervals, ultimately resulting in data products with improved spatial resolution. Three potential architectures are considered: a single pair of polar orbiting satellites, two pairs of polar orbiting satellites, and a polar orbiting pair of satellites coupled with a lower inclined pair of satellites. Results show that improvements in spatial resolution are obtained when one estimates a low resolution gravity field every two days for the case of a single pair of satellites, and every day for the case of two polar pairs of satellites. However, the spatial resolution for these cases is still lower than that provided by simply destriping and smoothing the solutions via standard GRACE post-processing techniques. Alternately, estimating daily gravity fields for the case of a polar pair of satellites coupled with a lower inclined pair results in solutions with superior spatial resolution than that offered by simply destriping and smoothing the solutions.     

Application of HF Doppler measurements for the investigation of internal atmospheric waves in the ionosphere

We present results of the spectral analysis of data series of Doppler frequency shifted signals reflected from the ionosphere, using experimental data received at Kazan University, Russia. Spectra of variations with periods from 1 min to 60 days have been calculated and analyzed for different scales of periods. The power spectral density for spring and winter differs by a factor of 3–4. Local maxima of variation amplitude are detected, which are statistically significant. The periods of these amplitude increases range from 6 to 12 min for winter, and from 24 to 48 min for autumn. Properties of spectra for variations with the periods of 1–72 h have been analyzed. The maximum of variation intensity for all seasons and frequencies corresponds to the period of 24 h. Spectra of variations with periods from 3 to 60 days have been calculated. The maxima periods of power spectral density have been detected by the MUSIC method for the high spectral resolution. The detected periods correspond to planetary wave periods. Analysis of spectra for days with different level of geomagnetic activity shows that the intensity of variations for days with a high level of geomagnetic activity is higher.     

The response of sporadic E-layer to the total solar eclipse of July 22, 2009 over the equatorial ionization anomaly region of the Indian zone

The digital ionosonde located in Bhopal (23.2°N, 77.2°E), India has been used to investigate the responses of the Es layer in the equatorial ionization anomaly (EIA) crest to the total solar eclipse (TSE) of July 22, 2009. Results show the presence of intense Es layer during and after the eclipse period. The gravity waves induced by the solar eclipse propagated upward in the Es layer and produced the periodic disturbance. The results of the wavelet analysis display the presence of dominant oscillation of about 24–32, 16–20 and 8 min. The appearance of intense sporadic-E concomitantly with the signatures of gravity wave suggests that the wind shear introduced by the solar eclipse induced gravity wave might be the plausible mechanism behind the intensification of Es-layer ionization.     

现代小卫星的重力梯度姿态稳定系统

本文首先论述现代小卫星重力梯度稳定系统组成，其次建立数学模型，最后进行计算机仿真，研究初始姿态角速度、惯量比、飞轮角动量和伸杆时间等因素对姿态动态特性的影响。     

Contribution of satellite laser ranging to combined gravity field models

In the framework of satellite-only gravity field modeling, satellite laser ranging (SLR) data is typically exploited to recover long-wavelength features. This contribution provides a detailed discussion of the SLR component of GOCO02S, the latest release of combined models within the GOCO series. Over a period of five years (January 2006 to December 2010), observations to LAGEOS-1, LAGEOS-2, Ajisai, Stella, and Starlette were analyzed. We conducted a series of closed-loop simulations and found that estimating monthly sets of spherical harmonic coefficients beyond degree five leads to exceedingly ill-posed normal equation systems. Therefore, we adopted degree five as the spectral resolution for real data analysis. We compared our monthly coefficient estimates of degree two with SLR and Gravity Recovery and Climate Experiment (GRACE) time series provided by the Center for Space Research (CSR) at Austin, Texas. Significant deviations in C₂₀ were noted between SLR and GRACE; the agreement is better for the non-zonal coefficients. Fitting sinusoids together with a linear trend to our C₂₀ time series yielded a rate of (−1.75 ± 0.6) × 10⁻¹¹/yr; this drift is equivalent to a geoid change from pole to equator of 0.35 ± 0.12 mm/yr or an apparent Greenland mass loss of 178.5 ± 61.2 km³/yr. The mean of all monthly solutions, averaged over the five-year period, served as input for the satellite-only model GOCO02S. The contribution of SLR to the combined gravity field model is highest for C₂₀, and hence is essential for the determination of the Earth’s oblateness.     

Atmosphere-ionosphere coupling from convectively generated gravity waves

Ionospheric variability impacts operational performances of a variety of technological systems, such as HF communication, Global Positioning System (GPS) navigation, and radar surveillance. The ionosphere is not only perturbed by geomagnetic inputs but is also influenced by atmospheric tides and other wave disturbances propagating from the troposphere to high altitudes. Atmospheric Gravity Waves (AGWs) excited by meteorological sources are one of the largest sources of mesoscale variability in the ionosphere. In this paper, Total Electron Content (TEC) data from networks of GPS receivers in the United States are analyzed to investigate AGWs in the ionosphere generated by convective thunderstorms. Two case studies of convectively generated gravity waves are presented. On April 4, 2014 two distinct large convective systems in Texas and Arkansas generated two sets of concentric AGWs that were observed in the ionosphere as Traveling Ionospheric Disturbances (TIDs). The period of the observed TIDs was 20.8 min, the horizontal wavelength was 182.4 km, and the horizontal phase speed was 146.4 m/s. The second case study shows TIDs generated from an extended squall line on December 23, 2015 stretching from the Gulf of Mexico to the Great Lakes in North America. Unlike the concentric wave features seen in the first case study, the extended squall line generated TIDs, which exhibited almost plane-parallel phase fronts. The TID period was 20.1 min, its horizontal wavelength was 209.6 km, and the horizontal phase speed was 180.1 m/s. The AGWs generated by both of these meteorological events have large vertical wavelength (>100 km), which are larger than the F2 layer thickness, thus allowing them to be discernible in the TEC dataset.     

空间太阳能电站姿态运动-结构振动耦合建模与分析

针对多旋转关节空间太阳能电站构型，利用基于能量等效原理的连续体等效方法将其等效为柔性梁模型，并考虑重力梯度影响，建立了姿态运动与结构振动的耦合动力学模型；结合Runge-Kutta 法和Newmark法的优点，提出了适用于求解姿态运动与结构振动耦合动力学方程的改进算法，相比于经典Runge-Kutta 法大幅提高了效率；利用改进算法得到了不同参数下的动力学响应。在此基础上，推导了结构振动量级随结构尺寸的六次方量级增加的规律，仿真结果表明尺寸过大引发不稳定现象；分析了姿态运动和重力梯度对结构振动频率和振幅的影响；发现了姿态运动周期受结构柔性影响而增大的现象，这种现象在低轨以及大初始姿态角下影响更为明显。