NeQuick模型算法研究及性能比较

研究了NeQuick2算法改进及其实现方法,从不同角度分析了NeQuick2模型在全球区域和中国区域内的性能优势。一个太阳活动周期内,中国区域NeQuick2模型计算的电子总含量（total electron content,TEC）比NeQuick1模型精度有显著提升,改正精度与太阳活动水平具有较强的相关性,低年比高年的改善效果更为显著。以全球电离层数据（global ionosphere maps,GIM）为参考标准,中国中高纬区域太阳活动低年NeQuick2模型TEC的系统年平均偏差减少了76%,年平均均方根（root mean square,RMS）值减少了约72%。太阳活动高年NeQuick2模型TEC的系统年平均偏差减少了38%,平均RMS减少了13%左右,且中高纬区域改正精度优于低纬区域11%~13%。全球区域太阳活动峰值期间NeQuick2模型TEC比NeQuick1模型日平均偏差改善了25%,日平均RMS改善了30%左右。分别用NeQuick1和NeQuick2模型得出F₂层顶部区域在太阳活动峰值期电子密度随高度剖面分布,顶部电子密度剖面精度改善近40%。最后分别得出了两个模型中国区域中高纬地区E和F₁层区域在100 km、150 km和200 km高度的电子密度分布图,结果显示NeQuick2模型改善了电子密度分布状况,有效避免了NeQuick1在底部区域电子密度梯度不连续以及电离层异常结构的情况。     

国际参考电离层2016模型在陆海区域的精度分析

为了研究IRI 2016模型在陆地、海洋及全球整体的预报精度和差异性问题,该文利用IGS组织提供的20002019年全球电离层TEC数据和12个GNSS跟踪站(陆地区域跟踪站6个,海洋区域跟踪站6个)实测TEC数据,基于陆地和海洋独立研究的方法,借助数理统计、相关系数及时间序列,分析了IRI 2016模型在陆地与海洋区域的精度特征.结果 表明:IRI 2016模型精度与研究区域内跟踪站的数量、纬度有密切关系,跟踪站密集区域、低纬度地区模型精度较高;太阳活动强度与IRI 2016模型精度高度相关,2008年和2019年为太阳活动低年,模型的精度较高.IRI 2016模型在全球范围内,相较于海洋区域,陆地区域模型的精度较高;与春秋冬三季相比,夏季TEC预测值与CODE GIM统计差值最小,模型的精度最高.     

LEO卫星单频精密定轨电离层模型改进算法

电离层延迟的有效改正是LEO卫星单频精密定轨的关键所在。目前主要采用电离层比例因子法进行LEO卫星电离层延迟改正,但该方法在电子密度峰值高度确定时未考虑太阳活动、经纬度、昼夜变化、季节等因素的影响。IRI2012模型虽然考虑了上述因素对电子密度峰值高度的影响,但因其与电离层薄层高度选择的标准不一致,通常它们之间存在系统性偏差而无法直接使用。为此本文提出将电离层薄层高度作为约束条件对IRI2012模型确定的电子密度峰值高度的均值进行参数约束估计,得到一种改进的电离层模型算法,并利用Swarm卫星GPS观测数据对其进行验证。结果表明:改进的电离层模型对Swarm卫星径向、切向和法向定轨精度均有不同程度的提高,尤其对轨道径向和法向精度改善最为明显,分别提高了31.6%和32.0%;同时较大幅度地降低了轨道的系统性偏差,尤其是在径向和法向,分别平均降低了65.0%和54.7%。     

基于深度学习构建的全球电离层NmF2模型

采用递归神经网络对空基COSMIC和地基垂测站数据建立了全球电离层峰值电子密度模型,模型均方根误差达到1.3×10~5 el/cm~3。在春夏秋冬4个季节内,人工神经网络ANN模型预测精度比IRI模型分别提高了25.7%、19.7%、33.3%和21.8%。另外,ANN模型不仅能够有效地模拟全球电离层时空变化特征,也能够成功地模拟电离层的诸多区域物理变化特性,如赤道电离异常、威德尔海异常、中纬度夜间异常和冬季异常。ANN模型可以为改正单频接收机的电离层延迟发挥一定的作用。     

重力固体潮信号的独立成分分析与谱相关解调

针对提取重力固体潮潮汐谐波分量,分析各谐波分量间调制关系的问题,该文结合重力固体潮的产生机制建立了一种重力固体潮信号正交分解模型。利用独立分量分析算法和谱相关方法,分离出与模型相对应的独立分量信号,并从独立成分中解调出载波分量和被调制分量。从而完整地描述了地球自转产生的潮汐谐波与月球、太阳相对于地球轨道变化产生的潮汐谐波的正交分解关系。实验表明,此方法能有效地提取出重力固体潮信号中所包含的潮汐谐波信息,并揭示出该重力固体潮信号中存在3个线性叠加的独立分量,在月波分量和年波分量都以乘性调制的方式调制到半日波分量和日波分量上。     

利用卫星重力梯度数据确定地球重力场的单层位模型

本文提出了处理卫星重力梯度数据以确定高分辩力重力场模型的单层位法并对其中的独立估计法进行了误差分析，数字结果显示：当卫星高度为200km，卫星数据网格宽度为15′，卫星重力梯度数据的精度为2×10~(-3)E时，利用独立估计法可得到分辩力为1°×1°（100km）的全球重力场模型，其重力异常精度小于1（mgal）；若卫星高度降至160km，卫星重力梯度数据的精度达到3×10~(-4)E，则获得的重力场模型的分辩力可提高到0.5°×0.5°（50km），其重力异常精度仍小于1（mgal）。     

利用掩星技术研究南极地区顶部电离层特性

掩星观测能够提供地面到低轨卫星轨道高度处的整个电离层电子密度剖面,对于顶部电离层的研究有重要的作用。本文利用COSMIC(constellation observing system for meteorology ionosphere and climate)掩星数据反演了电子密度剖面,提取了F₂层峰值高度(hmF2)、F₂层峰值密度(NmF2)、垂直标尺高(vertical scale height,VSH)等电离层参数,研究了南极地区的F₂层在太阳活动周期内的变化、年际变化、周日变化等,并且重点分析了南极地区的顶部电离层的垂直结构特征,尤其是威德尔海异常在垂直方向上的变化。结果表明,整个南极的hmF2每日均值在250~300 km左右,NmF2每日均值在1~8×1011 el/m3之间,VSH每日均值在100~250 km,威德尔海异常主要表现在顶部电子密度的增强和底部电子密度的减少。     

联合多源数据确定中国海及周边海底地形模型!

联合船测海深、重力异常和垂直重力梯度异常数据,构建了中国海及周边地区1′×1′海底地形模型。船测海深数据主要用于构建波长大于200km波段的海底地形,垂直重力梯度异常数据用于反演100~200km波段内的海底地形,而重力异常数据用于反演波长小于100km波段的海底地形。以船测海深为检核参考,考察了本文模型在南中国海和菲律宾海地区的精度。结果表明,本文模型精度优于ETOPO1、GEBCO和DTU10模型,在南中国海地区与斯克里普斯海洋研究所(Scripps Institute of Oceanography,SIO)的海底地形模型V15.1精度一致,在菲律宾海地区较V15.1模型精度提高了约8.2%。     

AR模型谱在超导重力数据信号检测中的分析研究

在介绍AR（auto-regression）模型谱分析原理的基础上,分别采用AR模型谱和周期图法对法国Strasbourg、澳大利亚Mt Stromlo和日本Matsushiro三个站的超导重力数据进行信号检测,以半日波的理论值为依据,运用两种方法进行半日波信号检测、分析与比较。结果表明,在超导重力数据信号检测分析中,AR模型谱比周期图法更准确、稳定,且受数据量的影响较小。     

南极中山站的重力和区域海洋潮汐特征

本文利用中山站弹簧重力仪记录的重力潮汐时间序列、验潮站数据、CATS2008区域和Eot11a全球海潮模型研究重力和海洋潮汐特征。结果表明,在周日频段,潮波O1的海潮振幅达到28 cm,4个主要潮波（Q1、O1、P1和K1）的全球模型与验潮站潮高差之和为4.2 cm,区域模型与验潮站潮高差之和为4.4 cm;在半日频段,潮波M2的海潮振幅达到20 cm,4个主要潮波（N2、M2、S2和K2）的潮高差之和分别为7.7 cm和5.1 cm,说明利用区域模型修正全球模型的重要性。经区域模型修正的全球海潮负荷改正后,重力主波K1、M2和S2的最终残差振幅分别下降了9.84%、56.14%和37.08%,说明区域海潮模型更能反映海洋潮汐的真实特征,用区域模型修正全球海潮模型的有效性得到验证。     

The international reference ionosphere today and in the future

The international reference ionosphere (IRI) is the internationally recognized and recommended standard for the specification of plasma parameters in Earth’s ionosphere. It describes monthly averages of electron density, electron temperature, ion temperature, ion composition, and several additional parameters in the altitude range from 60 to 1,500 km. A joint working group of the Committee on Space Research (COSPAR) and the International Union of Radio Science (URSI) is in charge of developing and improving the IRI model. As requested by COSPAR and URSI, IRI is an empirical model being based on most of the available and reliable data sources for the ionospheric plasma. The paper describes the latest version of the model and reviews efforts towards future improvements, including the development of new global models for the F2 peak density and height, and a new approach to describe the electron density in the topside and plasmasphere. Our emphasis will be on the electron density because it is the IRI parameter most relevant to geodetic techniques and studies. Annual IRI meetings are the main venue for the discussion of IRI activities, future improvements, and additions to the model. A new special IRI task force activity is focusing on the development of a real-time IRI (RT-IRI) by combining data assimilation techniques with the IRI model. A first RT-IRI task force meeting was held in 2009 in Colorado Springs. We will review the outcome of this meeting and the plans for the future. The IRI homepage is at .     

A preliminary evaluation of the performance of multiple ionospheric models in low- and mid-latitude regions of China in 2010–2011

Ionospheric delay is a dominant error source in Global Navigation Satellite System (GNSS). Single-frequency GNSS applications require ionospheric correction of signal delay caused by the charged particles in the earth’s ionosphere. The Chinese Beidou system is developing its own ionospheric model for single-frequency users. The number of single-frequency GNSS users and applications is expected to grow fast in the next years in China. Thus, developing an appropriate ionospheric model is crucially important for the Chinese Beidou system and worldwide single-frequency Beidou users. We study the performance of five globally accessible ionospheric models Global Ionospheric Map (GIM), International Reference Ionosphere (IRI), Parameterized Ionospheric Model (PIM), Klobuchar and NeQuick in low- and mid-latitude regions of China under mid-solar activity condition. Generally, all ionospheric models can reproduce the trend of diurnal ionosphere variations. It is found that all the models have better performances in mid-latitude than in low-latitude regions. When all the models are compared to the observed total electron content (TEC) data derived from GIM model, the IRI model (2012 version) has the best agreement with GIM model and the NeQuick has the poorest agreement. The RMS errors of the IRI model using the GIM TEC as reference truth are about 3.0–10.0 TECU in low-latitude regions and 3.0–8.0 TECU in mid-latitude regions, as observed during a period of 1 year with medium level of solar activity. When all the ionospheric models are ingested into single-frequency precise point positioning (PPP) to correct the ionospheric delays in GPS observations, the PIM model performs the best in both low and mid-latitudes in China. In mid-latitude, the daily single-frequency PPP accuracy using PIM model is ~10 cm in horizontal and ~20 cm in up direction. At low-latitude regions, the PPP error using PIM model is 10–20 cm in north, 30–40 cm in east and ~60 cm in up component. The single-frequency PPP solutions indicate that NeQuick model has the lowest accuracy among all the models in both low- and mid-latitude regions of China. This study suggests that the PIM model may be considered for single-frequency GNSS users in China to achieve a good positioning accuracy in both low- and mid-latitude regions.     

Simulation study of a follow-on gravity mission to GRACE

The gravity recovery and climate experiment (GRACE) has been providing monthly estimates of the Earth’s time-variable gravity field since its launch in March 2002. The GRACE gravity estimates are used to study temporal mass variations on global and regional scales, which are largely caused by a redistribution of water mass in the Earth system. The accuracy of the GRACE gravity fields are primarily limited by the satellite-to-satellite range-rate measurement noise, accelerometer errors, attitude errors, orbit errors, and temporal aliasing caused by un-modeled high-frequency variations in the gravity signal. Recent work by Ball Aerospace & Technologies Corp., Boulder, CO has resulted in the successful development of an interferometric laser ranging system to specifically address the limitations of the K-band microwave ranging system that provides the satellite-to-satellite measurements for the GRACE mission. Full numerical simulations are performed for several possible configurations of a GRACE Follow-On (GFO) mission to determine if a future satellite gravity recovery mission equipped with a laser ranging system will provide better estimates of time-variable gravity, thus benefiting many areas of Earth systems research. The laser ranging system improves the range-rate measurement precision to ~0.6 nm/s as compared to ~0.2 μm/s for the GRACE K-band microwave ranging instrument. Four different mission scenarios are simulated to investigate the effect of the better instrument at two different altitudes. The first pair of simulated missions is flown at GRACE altitude (~480 km) assuming on-board accelerometers with the same noise characteristics as those currently used for GRACE. The second pair of missions is flown at an altitude of ~250 km which requires a drag-free system to prevent satellite re-entry. In addition to allowing a lower satellite altitude, the drag-free system also reduces the errors associated with the accelerometer. All simulated mission scenarios assume a two satellite co-orbiting pair similar to GRACE in a near-polar, near-circular orbit. A method for local time variable gravity recovery through mass concentration blocks (mascons) is used to form simulated gravity estimates for Greenland and the Amazon region for three GFO configurations and GRACE. Simulation results show that the increased precision of the laser does not improve gravity estimation when flown with on-board accelerometers at the same altitude and spacecraft separation as GRACE, even when time-varying background models are not included. This study also shows that only modest improvement is realized for the best-case scenario (laser, low-altitude, drag-free) as compared to GRACE due to temporal aliasing errors. These errors are caused by high-frequency variations in the hydrology signal and imperfections in the atmospheric, oceanographic, and tidal models which are used to remove unwanted signal. This work concludes that applying the updated technologies alone will not immediately advance the accuracy of the gravity estimates. If the scientific objectives of a GFO mission require more accurate gravity estimates, then future work should focus on improvements in the geophysical models, and ways in which the mission design or data processing could reduce the effects of temporal aliasing.     

Validation of international reference ionosphere models using in situ measurements from GRACE K-band ranging system and CHAMP planar Langmuir probe

The in situ measurements of electron contents from GRACE K-band (dual-frequency) ranging system and CHAMP planar Langmuir probe were used to validate the international reference ionosphere (IRI) models. The comparison using measurements from year 2003 to 2007 shows a general agreement between data and the model outputs. The improvement in the newer IRI model (IRI-2007) is evident with the measurements from the GRACE satellites orbiting at the higher altitude. We present the comparison between the models and data comprehensively for various cases in solar activity, local time, season, and latitude. The IRI models do not well predict the electron density in the years 2006 and later, when the solar activity is extremely low. The IRI models generally overestimate the electron density during local winter while they underestimate during local summer. In the equatorial region, the large difference at local sunrise lasts for all years and all seasons. The IRI models do not perform well in predicting the anomaly in the polar region such as the Weddell Sea Anomaly. These discrepancies are likely due to smoothed (12-month averaged) solar activity indices used in the IRI models and due to insufficient spherical harmonic representation not able to capture small spatial scales. In near future, further improvement on the IRI models is expected by assimilating those in situ satellite data by implementing higher resolution (spatial and temporal) parameterizations.     

利用卫星重力梯度与地面数据确定地球重力场的分析

本文利用简捷的球谐分析方法讨论了重力场元在地面和空间的谱分布特征和向下延拓问题,分析了各类测量数据求定重力场的最高分辨率及精度。结果表明,在一个低轨道卫星上以适当的精度(优于10~(-2)E)的重力梯度测量可以获得空间分辨率为100公里、精度高于5mgal和10cm的重力场和大地水准面。     

临近空间大气参数误差特性分析

为定量分析临近空间大气环境参数的准确性,利用10年以上的TIMED和ENVISAT卫星探测数据,以及根据理想气体状态方程和地转风公式计算出的密度和风场结果,统计得到大气温度、密度、纬向风、经向风和合成风的分布,并与中国参考大气开展对比验证,分析中国区域内临近空间大气参数误差随高度、纬度和经度的季节变化规律,对于临近空间卫星数据应用、环境特性分析和气象保障具有重要意义。研究表明:温度偏差在55 km高度以上的春季和秋季沿纬向逐渐减小,密度偏差随高度降低逐渐增大,沿经向在30 km高度以下存在偏差较大的经度带。纬向风偏差在夏季随高度的变化特征明显不同,沿纬向在40—70 km高度逐渐减小;经向风偏差在40 km高度以下沿纬向均匀分布,且季节性差异较小;合成风偏差随高度的震荡特征明显,沿纬向在秋季的40—60 km高度先减小后增大,沿经向在春季的30—45 km高度呈不断减小趋势。     

月球重力场模型研究进展和我国将来月球卫星重力梯度计划实施

本文介绍了基于国际探月观测数据建立的月球重力场模型:8×4、15×8、13×13、5×5、7×7、16×16-1/2/3、Lun60d、GLGM-1/2、LP75D/G、LP100K/J、LP165P、LP150Q和SGM90d;通过对比SST-HL/LL-Doppler-VLBI和SST-HL/SGG-Doppler-VLBI跟踪观测模式的优缺点,建议我国将来首期月球卫星重力测量计划采用SST-HL/SGG-Doppler-VLBI较优;其次,通过对比静电悬浮、超导和量子卫星重力梯度仪的优缺点,建议我国将来首期月球卫星重力梯度计划采用静电悬浮重力梯度仪;并建议我国将来首颗月球重力梯度卫星的轨道高度(50～100 km)选择在已有月球探测卫星的测量盲区,轨道倾角(90°±3°)设计为有利于月球卫星观测数据全球覆盖的近极轨模式。     

不同GPS掩星电离层剖面产品相关性分析

将在一定时空限定范围内的不同低轨卫星COSMIC、GRACE、CHAMP、FY3C的电离层掩星电子密度剖面定义为一个掩星对来对比分析不同类型掩星电离层产品。结果表明：COSMIC掩星对之间的电子密度剖面整体轮廓符合得很好，电子密度剖面主要在250 km以下和500 km以上存在较大的偏差，250~500 km的电子密度整体偏差较小，统计得到的COSMIC掩星对的电子密度参量NmF2和hmF2的相关系数能分别达到0.99和0.97，具有高度相关性，不同COSMIC卫星之间没有明显的系统误差；COSMIC、GRACE、CHAMP和FY3C不同低轨卫星间的电子密度剖面略有差异，通过统计电子密度参量NmF2和hmF2之间的相关系数，COSMIC和CHAMP的相关系数分别为0.95和0.86，COSMIC和GRACE的相关系数分别为0.98和0.94，COSMIC和FY3C的相关系数分别为0.96和0.92，不同掩星类型之间的电子密度参量之间也具有高度相关性，验证了不同卫星任务GPS掩星电离层剖面的一致性。     

利用COSMIC掩星资料分析大气重力波参数的全球分布特征

本文介绍了从单一温度廓线提取重力波的水平波长、垂直波长、势能和动量通量的方法,利用2007~2008年COSMIC干温廓线数据计算得到了可以表征中尺度重力波活动特征的相关参数值在各年夏季和冬季的全球分布。结果与国内外已有的相关结果在形态分布上符合较好。结果表明,在20~30 km大气层,重力波在北半球处于冬季时的活动明显强于夏季,重力波的活动呈现纬向分布特征,夏半球高纬地区的重力波活动有明显的季节性变化,重力波的活动主要受地形和对流作用影响。COSMIC数据计算的动量通量分布大致与势能分布趋于一致。