COSMIC计划及掩星数据误差分析

掩星观测资料的误差特性是GPS气象学研究的热点之一,COSMIC（the Constellation Observing System for Meteorology,Ionosphere and Climate）计划在技术上做了很多改进,包括采用“开环”跟踪技术,以改进对低对流层掩星信号的跟踪能力,从而减小掩星数据的反演误差。本文采用2006年7月29日至12月31号的COSMIC掩星数据与全球探空资料进行对比,统计了掩星数据在不同地区、不同高度层的分布以及其折射率的误差特性,结果表明：（1）与探空资料匹配后的掩星事件主要集中在中纬度地区陆地上;在垂直方向5km以上,掩星探测趋于稳定。（2）“开环”技术的应用修正了以前掩星计划中出现在低对流层的折射率负偏差,但是却引进了折射率正偏差。在低纬地区,这种正偏差最大,在地表到1km范围内达到0．6％;随着纬度增加,正偏差减小,在中纬度地区10km以下最大值为0．3％;到高纬度地区,正偏差减小到0．2％以内。     

COSMIC反演精度和有关特性的检验

COSMIC是一个由6颗低轨卫星组成的用于天气、气候和电离层观测的空基GPS星座观测系统, 从2006年9月开始每天可提供覆盖全球的约2000~3000个掩星点, 掩星过程可提供从40 km高空到近地面的大气温、压、湿的廓线资料。为了有效利用这些资料, 以探空资料为基准, 对2007年1-10月我国及邻近区域的COSMIC掩星资料进行精度、可用性和全天候性的检验。结果表明: COSMIC反演的温度和折射率的精度很高, 水汽压的精度相对较差。与NCEP/NCAR再分析资料相比, 折射率和湿度廓线有更高精度。     

GNSS remote sensing of the Australian tropopause

Radio occultation (RO) techniques that use signals transmitted by Global Navigation Satellite Systems (GNSS) have emerged over the past decade as an important tool for measuring global changes in tropopause temperature and height, a valuable capacity given the tropopause??s sensitivity to temperature variations. This study uses 45,091 RO data from the CHAMP (CHAllenging Minisatellite Payload, 80 months), GRACE (Gravity Recovery And Climate Experiment, 23 months) and COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate, 20 months) satellites to analyse the variability of the tropopause??s height and temperature over Australia. GNSS RO temperature profiles from CHAMP, GRACE, and COSMIC are first validated using radiosonde observations provided by the Bureau of Meteorology (Australia). These are compared to RO soundings from between 2001 and 2007 that occurred within 3 h and 100 km of a radiosonde. The results indicate that RO soundings provide data of a comparable quality to radiosonde observations in the tropopause region, with temperature deviations of less than 0.5 ± 1.5 K. An analysis of tropopause height and temperature anomalies indicates a height increase over Australia as a whole of ca. 4.8 ± 1.3 m between September 2001 and April 2008, with a corresponding temperature decrease of ?0.019 ± 0.007 K. A similar pattern of increasing height/decreasing temperature was generally observed when determining the spatial distribution of the tropopause height and temperature rate of change over Australia. Although only a short period has been considered in this study, a function of the operating time of these satellites, the results nonetheless show an increase in the height of the tropopause over Australia during this period and thus may indicate regional warming. Several mechanisms could be responsible for these changes, such as an increase in the concentration of greenhouse gases in the atmosphere, and lower stratospheric cooling due to ozone loss, both of which have been observed during the last decades.     

利用COSMIC掩星资料分析青藏高原地区对流层顶季节变化特征

根据WMO(1957)对流层高度的判别方法,利用2007—2013年COSMIC掩星资料计算了对流层高度,并用无线电探空资料对结果进行检验,分析了青藏高原地区对流层顶季节变化特征.结果表明,COSMIC掩星资料和无线电探空资料判定的对流层高度具有很高的线性相关关系,相关系数高达0.976,平均值偏差为0.448 km,...     

The Impact of Assimilating FY-3C GNOS GPS Radio Occultation Observations on GRAPES Forecasts

In the present study, a gross quality control (QC) procedure is proposed for the Global Navigation Satellite System Occultation Sounder (GNOS) Global Positioning System radio occultation (GPS RO) refractivity data to remove abnormal data before they are assimilated. It consists of a climate extreme check removing data outside the range of the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) climate maxima and minima over approximately five years, and a vertical gradient check that rejects profiles containing super-refractions. These two QC steps were applied sequentially to identify outliers in GNOS GPS RO refractivity data during boreal winter 2013/2014. All of the abnormal refractivity profiles and the outliers at each level of the GNOS GPS RO observations were effectively removed by the proposed QC procedure. The post-QC GNOS GPS RO refractivity observations were then assimilated in the Global/ Regional Analysis and PrEdiction System (GRAPES) using the three-dimensional variational (3D-Var) system. The impacts of the GNOS refractivity observation on GRAPES analysis and forecasting were evaluated and analyzed using an observation system experiment run over one whole winter season of 2013 / 2014. The experiment results demonstrated a positive impact of GNOS GPS RO data on analysis and forecast quality. The root mean squared error of GRAPES analysis temperature was reduced by 1% in the Southern Hemisphere (SH) extratropics and in the tropics, and the anomaly correlation scores of the forecasted 500-hPa geopotential height over the SH increased significantly during days 1 to 5. Overall, the benefits of using GNOS GPS RO data are significant in the SH and tropics.     

COSMIC掩星数据与L波段探空数据的对比分析

COSMIC(Constellation Observation System for Meteorology,Ionosphere and Climate)每天可以提供全球2000～3000条从40 km高空到近地面的大气温、压、湿的廓线资料,有效地弥补了常规探空资料在时间和空间上分辨率的不足。通过对2008年5月20日至2008年11月26日COSMIC资料与L波段探空秒数据进行比对,结果表明,在10 km高度以下,COSMIC反演的湿廓线资料与L波段探空数据偏差较小,温度偏差为-0.5℃,均方根误差为1.5℃;折射率偏差为1.4N,均方根误差为5.9N;气压偏差为2.0 hPa,均方根误差为4.7 hPa;水汽压偏差为0.1 hPa,均方根误差为1.1 hPa。COSMIC干廓线资料与L波段探空相比,在10～30 km高度内,温度偏差为-0.3℃,均方根误差为1.9℃;折射率偏差为0.4N,均方根误差为0.9N;气压偏差为1.4 hPa,均方根误差为2.6 hPa。表明COSMIC资料既具有较高的时空分辨率,又具有较好的精度,在数值模式中具有重要的潜在应用。     

COSMIC反演大气温湿廓线与模式客观分析场的比较

利用2009年不同季节COSMIC湿反演的大气温度和相对湿度廓线数据,分别与时、空相匹配的ECMWF(European Centre for Medium-Range Weather Forecasts,欧洲中尺度天气预报中心)、NCEP(National Centers for Environmental Prediction,美国环境预报中心)模式客观分析场和无线电探空观测数据,进行全球范围的比较分析.初步研究表明,无论夏季还是冬季,各种资料源之间相互比较的偏差和标准差分布相似,与季节无关.就温度而言,三种资料源的温度水平、垂直分布都很接近,ECMWF模式数据比NCEP不论是温度廓线还是湿度廓线都更接近COSMIC反演值.模式的水汽客观分析场在对流层基本上都比无线电探空观测值偏湿,对流层中高层在大部分海洋地区也比COSMIC反演场偏湿.COSMIC反演的相对湿度相对于无线电探空整层偏大,具有明显正偏差,在300 hPa偏差达最大值(约30%).     

Planetary-Scale Wave Structures of the Earth's Atmosphere Revealed from the COSMIC Observations

GPS radio occultation(GPS RO) method,an active satellite-to-satellite remote sensing technique,is capable of producing accurate,all-weather,round the clock,global refractive index,density,pressure,and temperature profiles of the troposphere and stratosphere.This study presents planetary-scale equatorially trapped Kelvin waves in temperature profiles retrieved using COSMIC(Constellation Observing System for Meteorology,Ionosphere,and Climate) satellites during 2006-2009 and their interactions with background atmospheric conditions.It is found that the Kelvin waves are not only associated with wave periods of higher than 10 days(slow Kelvin waves) with higher zonal wave numbers(either 1 or 2),but also possessing downward phase progression,giving evidence that the source regions of them are located at lower altitudes.A thorough verification of outgoing longwave radiation(OLR) reveals that deep convection activity has developed regularly over the Indonesian region,suggesting that the Kelvin waves are driven by the convective activity.The derived Kelvin waves show enhanced(diminished) tendencies during westward(eastward) phase of the quasi-biennial oscillation(QBO) in zonal winds,implying a mutual relation between both of them.The El Nino and Southern Oscillation(ENSO) below 18 km and the QBO features between 18 and 27km in temperature profiles are observed during May 2006-May 2010 with the help of an adaptive data analysis technique known as Hilbert Huang Transform(HHT).Further,temperature anomalies computed using COSMIC retrieved temperatures are critically evaluated during different phases of ENSO,which has revealed interesting results and are discussed in light of available literature.     

通过影响高度多路径质量控制改进热带对流层低层GPS掩星资料同化效果

与观测气象、电离层和气候的GPS卫星计划COSMIC-1相比,COSMIC-2有更大功率的GPS接收天线,高2倍的采样率(100 Hz),小3倍的轨道倾角(24°)。因此,COSMIC-2和未来小卫星将为热带地区提供前所未有的大量无线电掩星观测资料。因为热带对流层低层水汽含量高,折射率有局地大梯度,违背大气球对称性假定,造成影响多路径现象,所以热带地区的掩星资料同化具有较大挑战性。本文首先利用之前研究中提出的影响参数多路径质量控制方法,去除热带对流层低层弯角模拟中出现影响参数多路径现象的资料,然后使用有限射线长度二维(2D)路径追踪观测算子和一维(1D)Abel变换观测算子,进行GPS弯角资料同化。结果表明,如果射线路径限制在离近地面点±300~400 km的水平距离以内,有限射线长度2D路径追踪观测算子得到的弯角模拟与2D射线追踪观测算子得到的弯角廓线之差不仅无偏,而且其标准差远比1D Abel的弯角模拟精度高。此外,无论采用1D Abel变换还是有限射线长度2D路径追踪观测算子来模拟弯角,影响参数多路径质量控制方案都能进一步改进GPS掩星弯角资料同化和预报精度。     

Planetary-Scale Wave Structures of the Earth’s Atmosphere Revealed from the COSMIC Observations

GPS radio occultation（GPS RO） method,an active satellite-to-satellite remote sensing technique,is capable of producing accurate,all-weather,round the clock,global refractive index,density,pressure,and temperature profiles of the troposphere and stratosphere.This study presents planetary-scale equatorially trapped Kelvin waves in temperature profiles retrieved using COSMIC（Constellation Observing System for Meteorology,Ionosphere,and Climate） satellites during 2006-2009 and their interactions with background atmospheric conditions.It is found that the Kelvin waves are not only associated with wave periods of higher than 10 days（slow Kelvin waves） with higher zonal wave numbers（either 1 or 2）,but also possessing downward phase progression,giving evidence that the source regions of them are located at lower altitudes.A thorough verification of outgoing longwave radiation（OLR） reveals that deep convection activity has developed regularly over the Indonesian region,suggesting that the Kelvin waves are driven by the convective activity.The derived Kelvin waves show enhanced（diminished） tendencies during westward（eastward） phase of the quasi-biennial oscillation（QBO） in zonal winds,implying a mutual relation between both of them.The El Nino and Southern Oscillation（ENSO） below 18 km and the QBO features between 18 and 27km in temperature profiles are observed during May 2006-May 2010 with the help of an adaptive data analysis technique known as Hilbert Huang Transform（HHT）.Further,temperature anomalies computed using COSMIC retrieved temperatures are critically evaluated during different phases of ENSO,which has revealed interesting results and are discussed in light of available literature.     

利用COSMIC掩星资料研究青藏高原地区大气边界层高度

以往关于青藏高原边界层的研究都是基于个别站点的常规观测,对青藏高原边界层的整体性认识受限。GPS掩星资料具有测量精度高和垂直分辨率高的特性,其廓线中含有大量有价值的边界层信息。利用2007—2013年COSMIC掩星资料,通过计算大气折射率最小梯度来确定边界层高度,并用无线电探空资料对结果进行了检验。在此基础上,对青藏高原地区边界层高度的特征及其形成机制展开了研究,比较了COSMIC掩星确定的边界层高度和ERA-Int的差别,讨论了最小梯度法用于边界层研究的不确定性。结果表明:青藏高原上COSMIC掩星和无线电探空数据检测的边界层高度相关系数为0.786,平均值偏差为0.049 km,均方根误差为0.363 km,COSMIC掩星数据检测的边界层高度和无线电探空的结果非常接近。青藏高原上边界层高度呈现西高东低的分布特征,高原中西部边界层高度主要为1.8—2.3 km,而高原东部边界层为1.4—1.8 km,最大值在高原西南部。青藏高原地区边界层有明显的季节差异,冬季高原上大部分地区边界层高度超过2.0 km;春季大部分地区高度降低,但在受印度季风影响的高原南部有明显的抬升,最大值可超过3.0 km;夏季高原上边界层高度开始升高,大部分地区超过1.8 km;秋季又开始回落。青藏高原以北塔克拉玛干沙漠和高原以南印度季风活动区是两个高值区,北部的沙漠地区边界层高度在夏季最高,南部印度季风活动区在季风爆发前（4月）达到全年最大值。青藏高原中西部地区有水平风辐合以及广泛的上升运动,为边界层的发展提供了动力条件,而东部的下沉运动对边界层的发展有抑制作用。青藏高原边界层各个季节的空间分布与地表感热通量分布一致。COSMIC掩星资料确定的边界层高度和ERA-Int相比,空间分布基本一致但ERA-Int边界层高度明显偏低。当有系统性强逆温存在的时候,或者云中液态水或冰水含量较大时,用最小梯度法检测的边界层高度不确定性增加。      

一次降水天气过程的GPS掩星资料在GSI同化系统中的应用

为了探索全球定位卫星观测系统无线电掩星技术(Global Positioning System/Radio Occulta-tion,简称GPSRO)在业务数值预报模式中的应用效果,以美国NCEP(National Centers for Environ-mental Prediction)业务同化系统为平台,对有/无COSMIC(Constellation Observing System for Meteor-ology,Ionosphere and Climate)GPS无线电掩星观测资料进行了一次降水过程的详细的同化对比试验。试验中GPSRO观测算子主要采用了相对简单但计算省时的局地(Local)同化技术。该个例试验表明:COSMICGPS掩星观测对分析影响较大的地区是常规观测资料较稀缺地区(如洋面),因此在这些地区GPSRO资料可以被看作是对常规观测很好的补充。同时还注意到COSMIC GPS RO资料在对流层下层对湿度场分析和预报的影响显著,在对流层中上层对风场分析也产生正面影响。这证实GPS RO资料在业务数值天气预报中应用的有效性,为我国数值预报未来开展业务GPS掩星观测资料同化提供了有益的科学参考依据。     

一种新的COSMIC大气折射率资料观测误差估计方法及在GRAPES全球三维变分同化中的应用

观测误差协方差是变分同化系统中决定分析及预报效果的关键参数之一,观测误差的估计精度直接影响变分同化分析和预报效果。分析了新息增量法(H-L法)估计全球定位系统无线电掩星这类观测点不固定资料的观测误差的适用条件,并利用1年的气象、电离层及气候星座观测系统(COSMIC)折射率资料,针对局地观测算子,估计了COSMIC折射率在南、北半球高、中、低6个纬度带四季的观测误差,分析了COSMIC折射率观测误差的纬度、高度和季节变化的特点,并将估计的折射率观测误差应用于GRAPES(Global/Regional Assimilation and Prediction Enhanced System)三维变分同化系统。结果表明,折射率观测误差随高度和纬度有明显变化;在中、高纬度带,折射率观测误差有显著的季节变化:夏季折射率的观测误差约为冬季2倍,春、秋两季折射率误差具有较好的南北对称性,冬、夏两季折射率观测误差南、北半球差异较大。与GRAPES原来使用的全球平均单一的折射率观测误差相比,在GRAPES全球三维变分同化系统中使用本研究估计的较高精度的随纬度变化的COSMIC折射率观测误差能够提高GRAPES全球变分同化系统的预报水平。     

COSMIC-2掩星反演数据质量分析

为评价COSMIC-2掩星反演数据精度,利用探空站点资料,对比验证了基于COSMIC-2附加相位数据由无线电掩星数据处理软件ROPP计算得到的大气折射率及温湿廓线。实验结果表明:COSMIC-2数据经反演后的大气廓线质量较高,大气折射率相对偏差标准差为1.5%~4.3%,大气温度和相对湿度均方根误差分别为1.3~2.2 K、10%~15%;低纬地区样本精度总体优于中纬地区; 季节性差异不明显,但冬春季样本精度总体略优于夏秋季。反演数据与COSMIC-2二级产品相比,折射率与温度在部分高度层精度更优,相对湿度在总体上精度更优。      

The impact of Global Positioning System data on the prediction of an extratropical cyclone: an observing system simulation experiment

In this paper, we report a series of observing system simulation experiments that we conducted to assess the potential impact of Global Positioning System/meteorology (GPS/MET) refractivity data on short-range numerical weather prediction. We first conducted a control experiment using the Penn State/NCAR mesoscale model MM5 at 90-km resolution on an extratropical cyclone known as the ERICA (Experiment on Rapidly Intensifying Cyclones over the Atlantic) IOP 4 storm. The results from the control experiment were then used to simulate GPS/MET refractivity observations with different spatial resolution and measurement characteristics. The simulated refractivity observations were assimilated into an 180-km model during a 6-h period, which was followed by a 48-h forecast integration. Key findings can be summarized as follows:

• The assimilation of refractivity data at the 180-km resolution can recover important atmospheric structures in temperature and moisture fields both in the upper and lower troposphere, and, through the internal model dynamical processes, also the wind fields. The assimilation of refractivity data led to a considerably more accurate prediction of the cyclone.

• Distributing the refractivity randomly in space and applying a line averaging did not alter the results significantly, while reducing the spatial resolution from 180 km to 360 km produced a moderately degraded result. Even at the 360-km resolution, the GPS-type refractivity data still have a notable positive impact on cyclone prediction.

• Restricting the refractivity data to altitude 3 km and above considerably degraded its impact on cyclone prediction. This degradation was greater than the combined effects of distributing the refractivity data randomly, performing line averaging, and reducing the resolution to 360 km.

These results showed that the GPS/MET refractivity data is likely to have a significant impact on short-range operational numerical weather prediction. The random distribution and line averaging associated with the inherent GPS occultation do not pose a problem for effective assimilation. On the other hand, these results also argue that we need to improve the GPS/MET retrieval algorithm in order to recover useful data in the lower troposphere, and to increase the number of low-earth-orbiting satellites carrying GPS receivers in order to increase the density of GPS soundings, so that the potential impact of GPS/MET refractivity data on numerical weather prediction can be fully realized.     

全谱反演技术在山基掩星数据处理中的应用

山基掩星观测是一种探测特定区域气象数据的新技术。该文从山基掩星探测技术的概念出发,阐述了山基掩星探测技术的特点和优势,介绍了山基掩星观测数据几何光学反演方法; 在详细介绍全谱反演技术的基础上,提出了采用全谱反演技术来处理山基掩星观测数据的新方法。采用该方法对2005年8月在河北兴隆雾灵山获取的山基掩星观测试验探测数据进行了处理,成功获得了大气折射率剖面, 将全谱反演结果与几何光学反演结果和同时获取的时空匹配的探空数据进行了比对分析。结果表明:全谱反演结果与几何光学反演结果平均相对偏差小于2%,全谱反演结果偏小,标准偏差低于3%;全谱反演结果与常规探空结果的平均相对偏差为8.15%,全谱反演结果偏小,标准偏差为1.4%。     

Retrieval of Vertical Distribution of Tropospheric Refractivity through Ground-Based GPS Observation简

In the present reported study, the vertical distributions of local atmospheric refractivity were retrieved from ground- based GPS observations at low elevation angles. An improved optimization method was implemented at altitudes of 0-10 km to search for a best-fit refractivity profile that resulted in atmospheric delays most similar to the delays calculated from the observations. A ray-tracing model was used to simulate neutral atmospheric delays corresponding to a given refractivity profile. We initially performed a ＂theoretical retrieval＂, in which no observation data were involved, to verify the optimization method. A statistical relative error of this ＂theoretical retrieval＂ （-2% to 2%） indicated that such a retrieval is effective. In a practical retrieval, observations were obtained using a dual-frequency GPS receiver, and its initial value was provided by CIRA86aQ_UoG data. The statistical relative errors of the practical retrieval range from -3% to 5% were compared with co-located radiosonde measurements, Results clearly revealed diurnal variations in local refractivity prc,files, The results also suggest that the general vertical distribution of refractivity can be derived with a high temporal resolution. However, further study is needed to describe the vertical refractivity gradient clearly.     

COSMIC掩星反演大气温湿资料的质量特征分析

针对COSMIC(Constellation Observation System of Meteorology Ionosphere and Climate)掩星反演的大气温度和水汽压二级资料,利用常规探空观测和NCEP(National Centers for Environmental Prediction)分析资料分别进行质量检验分析,以揭示反演资料质量的海陆差异、随纬度和高度变化等三维空间特征。结果显示:我国区域反演资料的温度略低于探空观测与NCEP分析资料,温度与水汽压相对于检验资料的均方根误差较小;全球范围内,反演大气温湿资料的质量随高度和纬度的不同而存在明显差异,而且海洋和陆地上质量的水平和垂直分布特征也存在显著不同。总体上看,COSMIC反演大气温湿资料具有良好的可靠性与精确度,可作为我国数值预报资料同化的新资料,反演温湿资料的质量特征也可为COSMIC资料同化的质量控制和垂直稀疏化方案的设计提供科学依据。     

4D tomographic reconstruction of the tropospheric wet refractivity using the concept of virtual reference station,case study: northwest of Iran

Iran enjoys a variety of climatological conditions. Moreover, numerical weather prediction (NWP) models are not assimilated with the meteorological data in Iran, the country suffering from poor spatial and temporal resolution of radiosonde measurements. These facts make modeling of troposphere impossible using the measurements and NWP. On the other hand, the global positioning system (GPS) has been emerged as a valuable tool for modeling and remote sensing of Earth’s atmosphere. This research is the first attempt to address the tropospheric wet refractivity modeling by GPS measurements in Iran. Changes of topography in the study area are taken into account. As a leading work, virtual reference stations (VRS) are used to fix the rank deficiency of the problem. The model space resolution matrix is used to achieve the optimum spatial resolution of the tomographic model and the optimum number of VRS stations. The accuracy of the developed model (KNTU1) is investigated by deploying radiosonde measurements.     

基于无线电掩星数据,探空资料以及再分析资料估算大气边界层高度的比较

用11年的全球无线电掩星数据(COSMIC),无线电探空数据(IGRA)以及欧洲中心再分析资料(ERA-Interim)对全球大气边界层高度(PBLH)进行估算比较.结果表明:(1)在1200 UTC和0000 UTC,由ERA-Interim和IGRA数据估算得到的全球PBLH空间分布较为一致,相关性较好,在白天正午时候太阳辐射能力较强,对流活动频繁,估算得到的大气边界层高度较高.(2)由COSMIC掩星数据估算得到的边界层高度比探空数据和再分析数据估算结果整体偏大.(3)COSMIC掩星数据,IGRA探空数据以及ERA-Interim再分析资料估算结果都表明边界层高度在低纬度地区偏大,高纬度地区偏小.(4)分析不同数据估算边界层高度纬向季节性差异表明,IGRA探空数据和COSMIC数据间差异为-1700m至-500m,IGRA与ERA-Interim之间的差异为-500m至250m.此外,对于大多数纬度而言,三个数据集之间的差异在冬季较大,在夏季较小.