山基GPS掩星与自动气象站观测结果比较

山基GPS掩星技术可以精确获得低层大气折射率信息,其观测原理是:在高山的山顶安装GPS接收机,跟踪GPS卫星的低仰角和负仰角信号,通过载波相位的变化求出大气折射造成的弯曲角剖面,从而利用科学反演方法获得观测点高度以下的大气折射率剖面.为了对山基掩星测量数据进行验证,2005年8月1日～29日,在河北雾灵山(40.60°N,117.48°E,2118m)开展了为期一个月的山基GPS掩星观测与其他探测的时比实验.实验所用的JAVAD双频GPS接收机共工作576小时,接收机朝南观测,记录到山基掩星事件1136次,其中上升掩星621次,下降掩星515次.实验期间,在雾灵山顶利用自动气象观测仪每天24小时观测当地的温度、湿度和气压,由此可计算大气折射率.本文利用雾灵山GPS掩星实验所获得的观测数据,成功反演得到大气折射率剖面,并将所得接收机高度的折射率与雾灵山顶自动气象站观测结果进行比对,分析结果表明:山基GPS掩星和自动气象站观测结果是一致的.山基GPS掩星观测可为低层大气环境监测提供大量数据,具有潜在的应用前景.     

大气GPS掩星观测反演方法

大气GPS掩星观测可获得全球的大气气象参量剖面信息.本文阐述了地球大气GPS掩星观测反演原理,详细介绍了其几何光学近似反演方法和全谱反演方法,提出了将几何光学反演方法和全谱反演方法结合以形成可以处理多路径掩星数据的反演新方案.该方案和几何光学反演方法应用于GPS/MET和CHAMP大气掩星数据反演,成功地获得了大气参量剖面.结果表明,新反演方案是可行的、有效的GPS大气掩星反演方案.     

COSMIC大气掩星开环数据反演方法

COSMIC星座GPS无线电掩星探测利用GPS开环接收技术提高低层大气观测数据的质量和对上升掩星事件的跟踪能力.开环掩星观测数据受到GPS导航数据调制的影响,在其数据后处理中必须消除该影响以获得高质量的科学反演结果.利用GPS导航数据调制码数据和利用开环数据本身内在的关系等两种方法可以消除该影响.将上述方法应用于COSMIC的掩星事件个例反演,获得了修复的大气附加相位数据;并利用几何光学近似反演方法和全谱反演方法,获得了射线弯曲角. 全谱反演方法获得的弯曲角及其温度反演结果与COSMIC数据中心的结果一致,说明我们的算法是有效的.     

大气掩星反演误差特性初步分析

GPS大气掩星探测技术可以获得全球大气折射率、气压、密度、温度和水汽压等气象参数,该技术基本原理是基于几何光学近似的Abel积分反演.地球扁率、电离层传播时间延迟、大气大尺度水平梯度、多路径传播现象等因素在某些高度范围影响大气反演的精度.本文采用模拟的方法,分析其中地球扁率及电离层对反演结果的影响,并讨论局部圆弧修正及电离层修正的效果.利用CHAMP掩星实测轨道数据和有关电离层和大气经验模式、采用三维射线追踪方法模拟计算几种情形下的GPS掩星观测附加相位数据,对模拟数据进行反演,将反演气象参量剖面与模拟时给定模式剖面进行比较,得到了0～60 km高度范围内的反演误差.误差统计分析结果表明,局部圆弧中心的修正以及电离层修正,对于高精度的GPS掩星反演是非常重要的;电离层修正残差仍是制约30～60 km高度范围内反演精度的重要因素,进一步完善和优化大气掩星反演需要发展新的电离层修正算法.     

电离层GPS掩星观测反演技术

在电离层局部地区球对称假设下,推导了利用双频和单频无线电掩星观测数据,反演电离层电子密度剖面的两种方法. 双频反演的误差来自于载波相位的观测误差,单频反演误差则主要由伪距的观测精度决定. 由于载波相位测量精度比伪距测量精度高两个量级,因此双频反演的精度一般比单频反演的高些. 不过,两载波信号L1和L2之间的传播路径差异会给双频方法带来误差. 利用三维射线追踪的程序模拟的无线电掩星数据来评估这些方法,结果表明,反演出的电离层剖面与给定的模式电离层非常吻合,验证了两种方法的可靠性和准确性. 将这两种反演方法应用于处理实测的GPS/MET掩星观测数据,均能获取合理的电离层剖面信息. 且单频方法得到的反演剖面与双频方法相当一致, 这为利用LEO星载单频GPS接收机进行电离层掩星观测提供了理论基础.     

GNSS掩星中大气水汽的非线性反演

研究利用无线电掩星探测数据反演大气水汽参数,文中使用一维非线性方法,反演对流层水汽压廓线.首先,定义了一维目标函数,并以几何光学假设为前提,反演出电波弯曲角;将电波弯曲角作为掩星探测量,大气模式输出的温、湿参量作为初始场,代入目标函数,对目标函数求最优反演出水汽压廓线.文中给出了部分反演结果,经讨论分析认为:非线性方法反演水汽压对初始场的精度依赖较小,能够同时反演出大气温度\水汽压廓线及相应的误差分析,对数值天气预报的应用及数据同化研究具有重要意义.     

山基GPS掩星折射率与探空折射率比较

本文简要阐述了山基GPS掩星数据处理过程,介绍了雾灵山(40.60.N,117.48.E,～2118 m)山基GPS掩星观测实验的基本情况,将雾灵山山基GPS掩星折射率与探空折射率进行了比较.按照山基掩星和无线电探空在空间(经纬度)上相差在2.以内、时间上相差在1.5 h以内的条件,共选出了206个山基掩星折射率廓线样本用于比较分析.山基掩星折射率与探空折射率的个例比较,表明山基掩星折射率廓线合理,且和探空折射率基本一致.统计结果显示,在2100 m 以下,山基掩星折射率与探空折射率相对偏差的平均值为-2.0%,标准偏差为7.3%.     

电离层GPS掩星观测改正TEC反演方法

电离层掩星观测中,当低轨卫星（LEO）轨道高度较低时,轨道以上的电离层电子总含量（TEC）对掩星反演的影响不能忽略.目前,一般采用指数函数等外推方法来处理该问题,对反演结果可能引起较大误差.为提高电离层掩星反演精度,本文研究利用LEO处于非掩星一侧GPS观测数据的改正TEC新反演方法.用三维射线追踪程序计算出电离层掩星观测模拟数据,分别应用改正TEC方法和外推方法进行反演,将反演结果与所用模式值进行比较.结果表明：对于轨道高度约800km的GPS/MET掩星模拟数据,外推方法和改正TEC方法反演结果都与模式值基本一致;对于轨道高度约400km的CHAMP掩星模拟数据,外推方法误差较大,改正TEC方法反演结果与模式值相符得较好.将改正TEC方法应用于GPS/MET实测数据的反演,取得了合理的结果.这些说明,改正TEC算法是一种有效的电离层掩星反演方法,尤其是对于轨道较低的LEO的电离层掩星观测反演特别有用.     

GPS地面台网和掩星观测结合的时变三维电离层层析

本文给出GPS地面台网和掩星观测结合的时变三维电离层层析的原理、算法和基于实测数据的反演结果.反演结果的比较表明,联合地基GPS与掩星观测数据进行重建,电子密度整体图像的重建质量特别是其垂直结构的重建质量得到了明显改善.在平静日和磁暴期间两种条件下利用实测数据的重建结果表明,GPS地面台网和掩星观测结合的电离层层析可以获得电离层电子密度在高度-纬度-经度-时间四维空间中的变化.重建结果清晰地显示了磁暴期间电离层负相暴效应,表明结合GPS地面台网和掩星观测的时变三维电离层层析可以有效地监测扰动条件下的大尺度电离层结构.     

GPS无线电掩星技术反演大气参数方法对比

无线电掩星技术是一种探测地球大气的崭新途径.菲涅耳衍射反演是GPS掩星技术中基于物理光学原理的新的反演方法.文章研究了菲涅耳及几何光学方法的反演原理,利用GPS/MET数据分别反演得到了大气折射率及温度廓线.对比分析了菲涅耳衍射方法在10 Hz、25 Hz和50 Hz三个采样频率下不同的反演结果,通过与几何光学反演结果及相应探空数据的对比分析,表明菲涅耳衍射反演得到的结果在对流层,特别是水汽含量丰富的区域中更为准确可靠.     

Inversion, error analysis, and validation of GPS/MET occultation data

The global positioning system meteorology (GPS/MET) experiment was the first practical demonstration of global navigation satellite system (GNSS)-based active limb sounding employing the radio occultation technique. This method measures, as principal observable and with millimetric accuracy, the excess phase path (relative to propagation in vacuum) of GNSS-transmitted radio waves caused by refraction during passage through the Earth’s neutral atmosphere and ionosphere in limb geometry. It shows great potential utility for weather and climate system studies in providing an unique combination of global coverage, high vertical resolution and accuracy, long-term stability, and all-weather capability. We first describe our GPS/MET data processing scheme from excess phases via bending angles to the neutral atmospheric parameters refractivity, density, pressure and temperature. Special emphasis is given to ionospheric correction methodology and the inversion of bending angles to refractivities, where we introduce a matrix inversion technique (instead of the usual integral inversion). The matrix technique is shown to lead to identical results as integral inversion but is more directly extendable to inversion by optimal estimation. The quality of GPS/MET-derived profiles is analyzed with an error estimation analysis employing a Monte Carlo technique. We consider statistical errors together with systematic errors due to upper-boundary initialization of the retrieval by a priori bending angles. Perfect initialization and properly smoothed statistical errors allow for better than 1 K temperature retrieval accuracy up to the stratopause. No initialization and statistical errors yield better than 1 K accuracy up to 30 km but less than 3 K accuracy above 40 km. Given imperfect initialization, biases ≫ 2 K propagate down to below 30 km height in unfavorable realistic cases. Furthermore, results of a statistical validation of GPS/MET profiles through comparison with atmospheric analyses of the European Centre for Medium-range Weather Forecasts (ECMWF) are presented. The comparisons indicate the high utility of the occultation data in that very good agreement of upper troposphere/lower stratosphere temperature (better than 1.5 K rms, ≪ 0.5 K bias) is found for a region (Europe+USA) where the ECMWF analyses are known to be good, but poorer agreement for a region (Southern Pacific) where the analyses are known to be degraded.     

Inversion of GPS meteorology data

The GPS meteorology (GPS/MET) experiment, led by the Universities Corporation for Atmospheric Research (UCAR), consists of a GPS receiver aboard a low earth orbit (LEO) satellite which was launched on 3 April 1995. During a radio occultation the LEO satellite rises or sets relative to one of the 24 GPS satellites at the Earths horizon. Thereby the atmospheric layers are successively sounded by radio waves which propagate from the GPS satellite to the LEO satellite. From the observed phase path increases, which are due to refraction of the radio waves by the ionosphere and the neutral atmosphere, the atmospheric parameter refractivity, density, pressure and temperature are calculated with high accuracy and resolution (0.5–1.5 km). In the present study, practical aspects of the GPS/MET data analysis are discussed. The retrieval is based on the Abelian integral inversion of the atmospheric bending angle profile into the refractivity index profile. The problem of the upper boundary condition of the Abelian integral is described by examples. The statistical optimization approach which is applied to the data above 40 km and the use of topside bending angle profiles from model atmospheres stabilize the inversion. The retrieved temperature profiles are compared with corresponding profiles which have already been calculated by scientists of UCAR and Jet Propulsion Laboratory (JPL), using Abelian integral inversion too. The comparison shows that in some cases large differences occur (5 K and more). This is probably due to different treatment of the upper boundary condition, data runaways and noise. Several temperature profiles with wavelike structures at tropospheric and stratospheric heights are shown. While the periodic structures at upper stratospheric heights could be caused by residual errors of the ionospheric correction method, the periodic temperature fluctuations at heights below 30 km are most likely caused by atmospheric waves (vertically propagating large-scale gravity waves and equatorial waves).Present address: Communications Research Laboratory, Upper Atmosphere Section, 4-2-1 Nukui- Kita, Koganei- shi, Tokyo 184, Japan     

Radio occultation data analysis by the radioholographic method

The radioholographic method is briefly described and tested by using data of 4 radio occultation events observed by the GPS/MET experiment on 9 February 1997. The central point of the radioholographic method (Pavelyev, 1998) is the generation of a radiohologram along the LEO satellite trajectory which allows the calculation of angular spectra of the received GPS radio wave field at the LEO satellite. These spectra are promising in view of detection, analysis and reduction of multipath/diffraction effects, study of atmospheric irregularities and estimation of bending angle error. Initial analysis of angular spectra calculated by the multiple signal classification (MUSIC) method gives evidence that considerable multibeam propagation occurs at ray perigee heights below 20 km and at heights around 80–120 km for the 4 GPS/MET occultation events. Temperature profiles obtained by our analysis (radioholographic method, Abel inversion) are compared with those of the traditional retrieval by the UCAR GPS/MET team (bending angle from slope of phase front, Abel inversion). In 3 of 4 cases we found good agreement (standard deviation σ_T∼1.5°K between both retrievals at heights 0–30 km).     

The Assimilation of GPS Radio Occultation Data and its Impact on Rainfall Prediction along the West Coast of India during Monsoon 2002

In this study, the Weather Research and Forecasting (WRF-2.0.3.1) model with three-dimensional variational data assimilation (3DVAR) was utilized to study a heavy rainfall event along the west coast of India with and without the assimilation of GPS occultation refractivity soundings in the monsoon period of 2002. The WRF model is a next-generation mesoscale numerical weather prediction system designed to serve both operational forecasting and atmospheric research communities. The Global Positioning System (GPS) radio occultation (RO) refractivity data, processed by UCAR, were obtained from the CHAMP and SAC-C missions. This study investigates the impact of thirteen GPS occultation refractivity soundings only, as assimilated into the WRF model with 3DVAR, on the rainfall prediction over the western coastal mountain of India. The model simulation, with the finest resolution of 10 km, was in good agreement with rainfall observations, up to 72-h forecast. There are some subtle but important differences in predicted rainfalls between the control run CN (without the assimilation of refractivity soundings) and G13 (with the assimilation of thirteen GPS RO soundings). In general, the assimilation run G13 gives a better prediction in terms of both rainfall locations and amounts at later times. The moisture increments were analyzed at the initial and forecast times to assess the impact of GPS RO data assimilation. The results indicate that remote soundings in the forcing region could have significant impacts on distant downstream regions. It is anticipated, based on this study, that considerably occultation soundings available from the six-satellite constellation of FORMOSAT-3/COSMIC would have even more significant impacts on weather prediction in this region.     

GPS无线电掩星后向传播方法

在大气多路径传播条件下,讨论了两种无线电信号的反演方法:几何光学方法和后向传播方法.当大气存在多路径效应时,采用几何光学方法反演大气参数会引入较大误差.后向传播方法将无线电信号从多路径区域反推至单路径区域,减少了大气多路径效应的影响.为了比较不同的反演方法,利用多相位屏模型,数值模拟了大气多路径条件下无线电信号在大气中...