山基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无线电掩星技术反演大气参数方法对比

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

地基GPS低高度角观测反演大气折射率廓线的模拟仿真

用价值函数的最优估计方法,发展了地基GPS弯曲角和相位观测联合反演大气折射率廓线的算法. 在价值函数中提出了一个比较合理的权系统;通过四参数形式的折射率模型,采用变分同化技术,反演出大气折射率廓线. 模拟结果显示,地基GPS观测反演得到的低层大气折射率廓线与真值吻合较好,并可以应用到大气波导的探测.     

COSMIC数据验证AMSU平流层低层观测的初步分析结果

基于Global Positioning System (GPS)掩星数据在平流层具有较高准确性、稳定性的优势,本文尝试用新一代GPS掩星观测——the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC)资料验证不同卫星平台上先进的微波探测仪（AMSU）的平流层观测结果.通过COSMIC大气温度廓线与AMSU辐射传输模式结合,得到模拟亮温,然后与AMSU平流层观测进行匹配比较.分析表明GPS掩星数据能够作为一个相对独立的参量检验NOAA15、16、18卫星平台内部的偏差.通过一年数据的比较验证,初步显示不同卫星平台的AMSU观测亮温在平流层低层都偏低,并且NOAA18平台的亮温偏低程度明显大于NOAA15、16.AMSU亮温偏差在极地冬季较为显著,尤其南极地区NOAA18的偏差幅度达到1.8 K.结合24小时内AMSU观测亮温偏差变化及其样本分布特征,可以看到明显的太阳辐射差异可能是导致AMSU观测亮温在极地偏差显著的主要原因.     

GPS测量误差对大气掩星反演精度影响分析

在无线电GPS大气掩星观测中,GPS和LEO卫星的位置误差、速度误差、卫星钟的漂移、载波相位观测误差等会给其反演的大气参量带来一定的误差.为定量分析上述GPS测量误差对大气掩星反演精度的影响,本文采用模拟分析的方法,利用CHAMP实测轨道数据和大气与电离层经验模式,采用三维射线追踪方法模拟得到大气附加相位观测数据;然后根据各种GPS测量误差的不同性质,在模拟掩星观测大气附加相位数据中加入GPS测量误差并进行反演,将其反演结果与未加入误差时的反演结果进行比较,分析各测量误差所带来的掩星反演误差的特点;利用一天的151个分布全球的掩星事件的误差结果进行统计分析.结果表明:(1)±10 m的GPS径向位置误差引起的温度反演误差的平均误差在30 km的高度为(-+)0.8K,标准偏差为0.05K;±0.4 m的LEO径向位置误差引起的温度反演误差的平均误差在30 km的高度为(-+)1.2K,标准偏差为0.07K.(2)仅考虑卫星速度误差对掩星几何关系的影响,则±30 m/s的GPS速度误差引起的温度反演误差的平均误差最大为±0.02K,标准偏差为0.3K;±10 m/s的LEO速度误差引起的温度反演误差的平均误差最大为±0.4K,标准偏差为0.06K.(3)当卫星钟漂为1×10-9时,反演温度误差的平均误差在30 km的高度为263K,标准偏差为16K,当卫星钟漂为1×10-14时,反演温度误差的平均误差在30km的高度为0.019K,标准偏差为0.004K.(4)当L1和L2载波相位观测误差分别为1 mm和2 mm时,引起的温度反演的平均误差在30 km的高度为0.01K,标准偏差为0.4K.(5)利用较接近真实观测的三维大气和电离层背景下的模拟附加相位数据,同时加入国际最新指标的各种GPS测量误差进行反演,误差分析结果为:在30 km以下,折射率和密度的相对误差小于0.3%(1σ),压强的相对误差小于0.5%(1σ),温度误差小于1K(1σ).     

基于COSMIC掩星探测资料的云底高反演研究

本文基于相对湿度廓线进入云层时的突变实现云底高反演的思想,采用2008年11月至2009年1月的COSMIC掩星湿空气数据反演全球云底高度,并与探空资料反演结果进行对比分析,得出以下重要结论:(1)当温度－40 ℃相似文献   

地基GNSS VTEC约束的大气掩星电离层误差修正方法

介绍了一种考虑电离层沿GNSS掩星射线路径分布不对称且兼顾一阶和二阶项的大气掩星电离层误差修正新方法,它综合地基GNSS VTEC的水平变化信息和电离层模式的垂直变化信息,在沿"入射线"与"出射线"双边局部球对称假设下,估算GNSS掩星射线路径上的电子密度;进而计算一阶和二阶项弯曲角电离层误差廓线.采用太阳活动低年的2008年7月15日和太阳活动较高年的2013年7月15日两天的MetOp-A和GRACE掩星观测资料和IGS的GNSS VTEC数据,计算了GNSS大气掩星弯曲角一阶和二阶项电离层误差廓线.对比分析表明：新方法经验模型和理论模型的二阶项弯曲角电离层残差,以及经验模型与实测数据的一阶和二阶项弯曲角电离层误差均具有良好的一致性,因此该方法可用于L2信号数据质量较差或L2信号中断的掩星事件,同时修正大气掩星一阶和二阶电离层误差,从而提高弯曲角精度和掩星观测资料的利用率.     

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

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

电离层水平不均匀性对无线电掩星资料反演的影响

电离层水平不均匀结构是引起掩星资料反演的主要误差源.本文模拟不均匀结构对无线电掩星资料反演的影响.首先利用3-D全球射线追踪技术进行了模拟.结果指出:电离层电子密度20%的变化可能引起GPS的时延增加约12m.其次,将90%的掩星时延对应的路径长度,定义为电离层掩星测量的水平分辨率,其值为1200-2000 km.最后,计算了不同纬度电离层当地正午、子夜和日出时的电离层水平不对称指数.它表明:低纬夏季日出时,不对称指数达-78%,反演结果可能有较大误差.     

Characteristics of inertia-gravity waves around jet stream from radiosonde observations in Wuhan (30.5°N, 114.4°E)

The 5 years’ radiosonde data obtained from January 2000 to December 2004 in Wuhan (30.5°N, 114.4°E) have been used for studying the behaviors of inertia-gravity waves in the vicinity of the jet stream. It is observed that the wave intensity has a similar seasonal variation with the jet stream intensity with a strong winter maximum and a summer minimum. Moreover, a similar inter-annual trend for both the wave intensity and jet stream intensity is also found. These results suggest that the jet stream may be the predominant source of the inertia-gravity waves in the troposphere and lower stratosphere over Wuhan in the period of the 5 years. It is noticed from 28 radiosonde profiles during wintertime that the energy of inertia-gravity waves exhibits upward and downward propagation respectively above and below the jet stream. This indicates that the source of the inertia-gravity waves is within the jet stream. In these cases, the twin waves below and above the jet stream usually hold similar amplitudes. The horizontal propagation of the twin waves also shows some interesting relationship.     

Atmospheric temperature profiling by joint Raman,Rayleigh and Fe Boltzmann lidar measurements

Simultaneous and complete temperature profiles from near ground to about 100 km are essential for studying the dynamical coupling between different atmospheric layers. They are acquired by combining three different lidar techniques at Wuhan, China (30.5°N, 114.4°E). The atmospheric temperatures from about 3 to 25 km are calculated from the nitrogen molecule density profiles obtained from the N₂ vibrational Raman backscatter, while the atmospheric temperatures between 30 and ∼75 km are calculated by the standard Rayleigh scattering method. The temperatures in the 80–100 km altitude region are derived from the Fe Boltzmann technique. The temperature profiles measured by our lidar systems exhibit good agreement when compared with the radiosonde and satellite data, as well as the model. A Lomb–Scargle spectral analysis of the normalized temperature perturbations in the altitude range from 4 to 60 km shows that the spectral slopes of the vertical wave number spectra tended to −3 for large vertical wave numbers. This is consistent with the model predictions of saturated gravity wave spectra.     

Argo data assimilation in ocean general circulation model of Northwest Pacific Ocean

The Argo temperature and salinity profiles in 2005–2009 are assimilated into a coastal ocean general circulation model of the Northwest Pacific Ocean using the ensemble adjustment Kalman filter (EAKF). Three numerical tests, including the control run (CTL) (without data assimilation, which serves as the reference experiment), ensemble free run (EnFR) (without data assimilation), and EAKF experiment (with Argo data assimilation using EAKF), are carried out to examine the performance of this system. Using the restarts of different years as the initial conditions of the ensemble integrations, the ensemble spreads from EnFR and EAKF are all kept at a finite value after a sharp decreasing in the first few months because of the sensitive of the model to the initial conditions, and the reducing of the ensemble spread due to Argo data assimilation is not much. The ensemble samples obtained in this way can well represent the probabilities of the real ocean states, and no ensemble inflation is necessary for this EAKF experiment. Different experiment results are compared with satellite sea surface temperature (SST) data and the Global Temperature-Salinity Profile Program (GTSPP) data. The comparison of SST shows that modeled SST errors are reduced after data assimilation; the error reduction percentage after assimilating the Argo profiles is about 10?% on average. The comparison against the GTSPP profiles, which are independent of the Argo profiles, shows improvements in both temperature and salinity. The comparison results indicated a great error reduction in all vertical layers relative to CTL and the ensemble mean of EnFR; the maximum value for temperature and salinity reaches to 85?% and 80?%, respectively. The standard deviations of sea surface height are employed to examine the simulation ability, and it is shown that the mesoscale variability is improved after Argo data assimilation, especially in the Kuroshio extension area and along the section of 10°N. All these results suggest that this system is potentially useful for improving the simulation ability of oceanic numerical models.     

Rayleigh lidar temperature measurements in the upper troposphere and lower stratosphere

A single-wavelength Rayleigh lidar system has been used to measure the temperatures in the upper troposphere and lower stratosphere in the night in the altitude range from about 8 to $\text{30}\text{km}$. The temperature derivation is based on an inversion algorithm of the pure Rayleigh backscatter. Calculations include the derivation of the air molecular concentration by an iterative method and the backscattered signals corrected by the background aerosol, which is now found to be low and stable. The uncertainties in estimating the temperature using this method are discussed in detail.The temperature profiles and the tropopause characteristics derived by using the lidar measurements are compared with the radiosonde data. Good agreement is found between these two measurements revealing the potential of this method. The comparison with radiosonde data shows that the lidar measured tropopause temperature is lower by $\text{0.8±1.5}\text{K}$ and the tropopause height is higher by $\text{0.45±0.8}\text{km}$ than the radiosonde measurements. The climatology of local tropopause (24.57°N,121.13°E) is briefly discussed in terms of a double tropopause formation and seasonal variations of the tropopause height and temperature.     

Changes of the 0°C isotherm and the equilibrium line altitude in central Chile during the last quarter of the 20th century / Changements de l'isotherme 0°C et de la ligne d'équilibre des neiges dans le Chili central durant le dernier quart du 20ème siècle

Abstract

An overall retreat of glaciers has been observed in the Andes of central Chile during the last ~100 years. Precipitation is mainly of frontal origin and concentrates in winter months. Analysis of precipitation data shows a decrease until 1976, an increase thereafter north of 34°S and a decrease south of 34°S, but overall no significant trends during the last quarter of the 20th century. Analysis of radiosonde data of central Chile shows mid-tropospheric warming with an elevation increase of the 0°C isotherm of 122 ± 8 m and 200 ± 6 m in winter and summer, respectively, during the 27-year period between 1975 and 2001. The results point to a snowline elevation increase in the region during the last quarter of the 20th century and a concurrent rise of the equilibrium line altitude (ELA) and suggest that mid-troposphere warming is the main cause for glacier retreat in central Chile.     

A new hydrostatic mapping function for tropospheric delay estimation

A new hydrostatic mapping function is proposed by modifying the analytical solution of the Chapman grazing incident (Chi) function. Applicability of this mapping function is evaluated over the Indian region by comparing it with the true value of mapping function estimated through ray-tracing the refractivity profile of neutral atmosphere as well as by comparing it with the Niell and Global Mapping Functions. This comparison shows that the new mapping function, agrees very well with the ray-traces values and it is in par with if not better than the other global mapping functions for very high zenith angles up to ～88°.     

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     

Quality Assessment of Atmospheric Motion Vectors Over the Indian Ocean

Because conventional observations over the oceans are not available, especially during tropical cyclones, multi-spectral atmospheric motion vectors (AMVs) estimated from geostationary satellites are routinely assimilated in the numerical weather prediction models at different operational centres across the globe. The derived AMVs are generally validated with radiosonde observations available over land at synoptic hours; however, over the ocean there is a limited scope to assess the quality of AMVs. Over ocean, AMVs can be validated with radiosonde data available from opportunistic ships or using dropsonde data available from aircrafts. In this study, the accuracy of the AMVs derived from the geostationary satellites Kalpana-1 and Meteosat-7 is evaluated over the oceanic region. Radiosonde data available from a ship cruise held in the Bay of Bengal during the period 09 July–08 August 2012 and from the Cal/Val site situated at Kavaratti Island (72.62°E, 10.57°N) in the southern Indian Ocean are used to assess the AMV accuracy. In this study, 83 radiosonde profiles are used to validate the Kalpana-1 AMVs, to allow a better understanding of AMV errors over the Indian Ocean. The RMSVD of Kalpana-1 AMVs for the high-, mid- and low-levels are found to be 7.9, 9.4 and 5.3 m s^?1, respectively, while the corresponding RMSVD for Meteosat-7 AMVs are 9.1, 5.5 and 3.7 m s^?1. A similar accuracy is observed when the AMVs are validated against the NCEP analyses collocated with the nearest radiosonde locations. The high RMSVD and bias for Kalpana-1 AMVs at the mid-level and Meteosat-7 AMVs at the high-level are associated with the limitation of satellite winds to resolve the upper-level easterly jet in conjunction with errors in the height assignment. This study could help the numerical modellers to assign appropriate observation error over this region during the assimilation of AMVs into the NWP models.     

Inertial gravity wave parameters for the lower stratosphere from radiosonde data over China

Characterization of gravity wave(GW)parameters for the stratosphere is critical for global atmospheric circulation models.These parameters are mainly determined from measurements.Here,we investigate variation in inertial GW activity with season and latitude in the lower stratosphere(18-25 km)over China,using radiosonde data with a high vertical resolution over a 2-year period.Eight radiosonde stations were selected across China,with a latitudinal range of 22°-49°N.Analyses show that the GW energy in the lower stratosphere over China has obvious seasonal variation and a meridional distribution,similar to other regions of the globe.The GW energy is highest in winter,and lowest in summer;it decreases with increasing latitude.Velocity perturbations with longitude and latitude are almost the same,indicating that GW energy is horizontally isotropic.Typically,85%of the vertical wavelength distribution is concentrated between elevations of 1 and 3 km,with a mean value of 2 km;it is almost constant with latitude.Over 80%of all the horizontal wavelengths occur in the range 100-800 km,with a mean value of 450km;they show a weak decrease with increasing latitude,yielding a difference of about 40 km over the 22°-49°N range.The ratio of horizontal wavelength over vertical wavelength is about 200:1,which implies that inertial GWs in the lower stratosphere propagate along nearly horizontal planes.Ratios of their intrinsic frequency to the Coriolis parameter decrease with increasing latitude;most values are between 1 and 2,with a mean value of 1.5.Study of the propagation directions of GW energy shows that upward fractions account for over 60%at all stations.In contrast,the horizontal propagation direction is significantly anisotropic,and is mainly along prevailing wind directions;this anisotropy weakens with increasing latitude.     

Using ensemble adjustment Kalman filter to assimilate Argo profiles in a global OGCM

An ensemble adjustment Kalman filter (EAKF) is used to assimilate Argo profiles of 2008 in a global version of the Modular Ocean Model version 4. Four assimilation experiments are carried out to compare with the simulation without data assimilation, which serves as the control experiment. All experiment results are compared with dataset of Global Temperature–Salinity Profile Program and satellite sea surface temperature (SST). The first experiment (Exp 1) is implemented by perturbing temperature of upper layers in the initial conditions (ICs) with an amplitude of 1.0°C and no ensemble inflation. The results from Exp 1 show that the simulated temperature (salinity) deviation in the upper 400 m (500 m) is reduced through Argo data assimilation; however, these deviations are increased in deeper layers. The error reduction in SST is much greater during January to June than during the rest of the year. Three more experiments are designed to understand the responses in different layers and months. Two of them test model sensitivities to ICs by perturbing them vertically: one over the vertical extent of the whole water column (Exp 2) and the other employs smaller perturbation amplitude of 0.1°C (Exp 3). Exp 2 shows that the simulated temperature and salinity deviations are systematically improved in the whole water column. Comparison between Exps 2 and 3 suggests that perturbation amplitude is important. Exp 4 tests the influence of the optimal inflation factor of 5%, which is determined by other set of numerical tests. Exp 4 improves assimilation performance much more than the other three experiments without inflation. Therefore, we conclude that the perturbation should be introduced to all model layers, proper perturbation amplitude is important for Ocean data assimilation using EAKF, and the ensemble inflation by an optimal inflation is critical to improve the skill of the EAKF analysis.