差分GPS定向校准方法的研究

提出了一种差分GPS定向校准方法,即利用差分GPS测得固定站和运动站的坐标(WGS-84),应用计算模型求出运动站相对固定站的大地方位角,再与由一等天文点计算给出的大地方位角进行比对。给出了计算模型,试验设备,试验方法和试验数据,验证了思路的正确性和可行性。     

基于小波变换的GPS载波相位周跳检测与修复

载波相位周跳检测与修复是实现GPS高精度定位的难点之一。小波变换具有空间局部性好的特点,本文利用小波变换把载波相位观测量分解为高频和低频分量,提出了利用小波系数局部极值点检测相位周跳的方法。基于载波相位和多普勒频率共同约束的埃尔米特多项式,外推周跳时刻载波相位,与观测值比较确定周跳的大小,从而修复周跳。利用飞机实际观测数据和仿真的高动态数据进行了计算,给出了计算结果,表明该方法可以对高动态情况下载波相位周跳进行准确探测和修复。     

Analysis of ionospheric TEC response to solar and geomagnetic activities at different solar activity stages

Studying the relationship of total electron content (TEC) to solar or geomagnetic activities at different solar activity stages can provide a reference for ionospheric modeling and prediction. On the basis of solar activity indices, geomagnetic activity parameters, and ionospheric TEC data at different solar activity stages, this study analyzes the overall variation relationships of solar and geomagnetic activities with ionospheric TEC, the characteristics of the quasi-27-day periodic oscillations of the three variables at different stages, and the delayed TEC response of solar activity by conducting correlation analysis, Butterworth band-pass filtering, Fourier transform, and time lag analysis. The following results are obtained. (1) TEC exhibits a significant linear relationship with solar activity at different solar activity stages. The correlation coefficients |R| are arranged as follows: |R|_EUV > |R|_F10.7 > |R|sunspot number. No significant linear relationship exists between TEC and geomagnetic activity parameters (|R| < 0.35). (2) TEC, solar activity indices, and geomagnetic activity parameters have a period of 10.5 years. The maximum amplitudes of the Fourier spectrum for TEC and solar activity indices are nearly 27 days and those of geomagnetic activity parameters are nearly 27 and 13.5 days. (3) The deviations of the quasi-27-day significant periodic oscillation of TEC and solar activity indices are consistent. (4) No evident relationship exists between the quasi-27-day periodic oscillation of TEC and geomagnetic activity parameters. (5) The delay time of TEC for the 10.7 cm solar radio flux and extreme ultraviolet is always consistent, whereas that for sunspot number varies at each stage.     

Ionospheric TEC prediction using hybrid method based on ensemble empirical mode decomposition (EEMD) and long short-term memory (LSTM) deep learning model over India

Total electron data (TEC) from GPS nowadays can be used as a tool for understanding the space weather phenomena. The development of prediction model for TEC is quiet crucial and challenging due to the dynamic behavior of the ionosphere, since it depends on different factors such as seasonal, diurnal and spatial variations, solar geomagnetic conditions etc. In this paper, an attempt is made for predicting the GPS derived TEC values for different GNSS stations over India using a hybrid method based on Ensemble empirical mode decomposition (EEMD) and Long Short-Term Memory (LSTM) deep learning method. The daily TEC time series data from the IISc Bangalore (Latitude 13.021, Longitude 77.570), Lucknow (Latitude 26.912, Longitude 80.956) and Hyderabad (Latitude 17.417, Longitude 78.551) stations over India during the period 2008 to 2015 of solar cycle 23 and 24 is used for analysis. The assessment of model performance for testing predicted output compared with LSTM and EMD-LSTM models, and their comparison results show that the hybrid EEMD-LSTM model presents better than the other models.     

飞机载体的杆臂效应对GPS测速精度的影响

杆臂效应是由于测量体的安装位置与运动载体质心不重合而引起的。从理论上分析了杆臂效应产生的原理,推导得出杆臂效应引起的速度误差公式。通过仿真试验和飞行试验数据分析了飞机载体不同运动状态下,杆臂效应对GPS速度测量精度的影响。通过对杆臂误差补偿后,GPS测速精度得到很大提高。在惯导系统飞行试验中引入该方法,可以科学合理地给出系统试飞评定结果。     

Precise orbit determination and precision comparison for FY3C and FY3D with receiver antenna GPS and BDS PCV using spaceborne BDS and GPS observation data

The FY3C and FY3D satellites were equipped with global navigation satellite occultation detector (GNOS) receivers that received both GPS and BDS-2 signals. For further improving precise orbit determination (POD) precisions, we estimated receiver GPS and BDS signal phase center variations (PCV) models with 2° and 5° resolutions and set the different weights for GPS and BDS-2 observations in the combined POD. The BDS-based POD precision using BDS-2 satellite antenna phase center offset (PCO) values from the China Satellite Navigation Office (CSNO) are not as accurate as those obtained from the International GNSS Service (IGS) Multi-GNSS experiments project (MGEX). The estimated receiver GPS and BDS PCV models with 2° and 5° resolutions were estimated from the GPS phase residuals of GPS-based POD and BDS phase residuals of combined POD, respectively. In most cases, the POD precisions using the estimated PCVs with 2° resolution are superior to those with 5° resolution. The precisions of the BDS-based POD and combined POD were both improved by introducing the receiver BDS PCV models. The weighting for GPS and BDS-2 observations can further improve the precision of the combined POD. The tested results of selected weights are better than those with equal weight in the combined POD. The experiment results show that orbital precisions of FY3C are worse than those of FY3D.     

Negative VTEC with spherical harmonic expansion model under different solar activity conditions

The Earth’s ionosphere can be described by a spherical harmonic (SH) expansion up to a specific degree. However, there exist negative vertical total electron content (VTEC) values in the global ionosphere map (GIM) with the SH expansion model. In this contribution, we specifically investigated the negative VTEC values that are induced by the SH expansion model and validated the performance of the inequality-constrained least squares (ICLS) method in eliminating the negative VTEC values. The GPS data from 2004 to 2017 was selected to cover one solar cycle and the experiments under different solar activity conditions were analyzed. The results in our work show that the occurrence of the negative VTEC values is attributed to the deficiency of the SH expansion model when the VTEC itself is small instead of the unevenly distribution of the GNSS stations. The negative VTEC values appear periodically in the temporal domain, showing apparently one year and half year periods. During one year, two peaks in June and December can be observed in the time series of the negative VTEC values. The number of negative VTEC values in June is obvious larger than that in December. During one solar cycle, the number of negative VTEC values under quiet solar activity condition is obvious larger than that under strong solar activity condition. In the spatial domain, the appearance of the negative VTEC values is strongly related with the movement of the subsolar point. In the latitude of the subsolar point has the largest magnitude, the negative values will appear on the opposite hemisphere and the further from the subsolar point the more negative values. The maximum number of the negative VTEC values in the southern hemisphere appears in June, while the peak value in the northern hemisphere appears in December. The maximum number of negative VTEC values in the southern hemisphere is generally larger than that in the northern hemisphere. In addition, the negative VTEC values are distributed both at middle latitude and high latitude in the southern hemisphere, while they are mainly distributed at high latitude in the northern hemisphere. When the ICLS method is used, the negative VTEC values can be eliminated efficiently and it has nearly no influence on the positive VTEC values. The ICLS method can also improve the receiver’s differential code bias (DCB) and significantly decrease the unreasonable negative slant TEC (STEC) values along the lines of sight. Using the final GIM product of the Jet Propulsion Laboratory (JPLG) as a reference, the root mean square (RMS) of the ICLS solution shows maximum 25%, 20% and 45% improvement relative to the least squares (LS) solution at northern high latitude, southern middle latitude and southern high latitude, respectively.     

捷联式重力仪在海洋测量中的应用与数据处理

针对海洋重力测量对重力测量系统的要求,给出了一种捷联式海洋重力仪SAG-2M,论述了系统组成、工作原理、特点和数据处理流程。利用近期获得的某海域重力测量数据,评估了捷联式重力仪精度。测量结果表明,重复测线的重力异常曲线吻合度较高,交叉点不符值精度约为1.2mGal,满足海洋重力测量的指标要求。     

SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

Aspects related to variability of radiative cooling by NO in lower thermosphere,TEC and O/N2 correlation,and diffusion of NO into mesosphere during the Halloween storms

Nitric Oxide is a very important trace species which plays a significant role acting as a natural thermostat in Earth’s thermosphere during strong geomagnetic activity. In this paper, we present various aspects related to the variation in the NO Infrared radiative flux (IRF) exiting the thermosphere by utilizing the TIMED/SABER (Thermosphere Ionosphere Mesosphere Energetics and Dynamics/ Sounding of the Atmosphere using Broadband Emission Radiometry) observational data during the Halloween storm which occurred in late October 2003. The Halloween storm comprised of three intense-geomagnetic storms. The variability of NO infrared flux during these storm events and its connection to the strength of the geomagnetic storms were found to be different in contrast to similar super storms. The connection between the quantum of energy outflux from the upper atmosphere into space in terms of NO IRF and the duration of storms is established. The NO radiative cooling, and the closely correlated depletion in O/N₂ ratio are controlled by the Joule heating intensity (proxied by AE-index). The collisional excitation rate of NO, calculated using the modelled datasets of WACCM-X (Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension), correlates well with the observed pattern of radiative emission by NO. Observational datasets from TIMED/GUVI (Global Ultra-Violet Imager) and MIT Haystack observatory madrigal GNSS (Global navigation satellite system) total electron content (TEC) database shows that the TEC and O/N₂ enhancement in low-mid northern hemispheric latitudes are mainly controlled by the z-component of Interplanetary magnetic field (IMF-B_z). The penetration of eastward electric field during the storm events is found to be responsible for the overall enhancement of TEC. The contribution of enhanced day-side TEC in observed variation of O/N₂ ratio by GUVI is also reported. It is also seen that during substorms related events the night-time polar region experiences more cooling due to NO than the daytime polar region. The connections between the mid- and low-latitude enhancement in NO IRF with the propagation of LSTIDs (Large-scale traveling ionospheric disturbances) in combination with the O/N₂ variability, and the altitudinal variation in NO flux with the progression of the storm is also investigated. This study presents the evidence on the role of diffusion processes in the large scale enhancement of NO in the mesospheric altitudes.