Monitoring of ducting by using a ground-based GPS receiver

In this paper,we describe the estimation of low-altitude refractivity structure from simulation and real ground-based GPS delays.The vertical structure of the refractive environment is modeled using three parameters,i.e.,duct height,duct thickness,and duct slope.The refractivity model is implemented with a priori constraints on the duct height,thickness,and strength,which might be derived from soundings or numerical weather-prediction models.A ray propagation model maps the refractivity structure into a replica field.Replica fields are compared with the simulation observed data using a squarederror objective function.A global search for the three environmental parameters is performed using a genetic algorithm.The inversion is assessed by comparing the refractivity profiles from the radiosondes to those estimated.This technique could provide near-real-time estimation of the ducting effect.The results suggest that ground-based GPS provides significant atmospheric refractivity information,despite certain fundamental limitations of ground-based measurements.Radiosondes are typically launched just a few times daily.Consequently,estimates of temporally and spatially varying refractivity that assimilate GPS delays could substantially improve over-estimates caused by using radiosonde data alone.     

Inversion for atmosphere duct parameters using real radar sea clutter

This paper addresses the problem of estimating the lower atmospheric refractivity (M profile) under nonstandard propagation conditions frequently encountered in low altitude maritime radar applications. The vertical structure of the refractive environment is modeled using five parameters and the horizontal structure is modeled using five parameters. The refractivity model is implemented with and without a priori constraint on the duct strength as might be derived from soundings or numerical weather-prediction models. An electromagnetic propagation model maps the refractivity structure into a replica field. Replica fields are compared with the observed clutter using a squared-error objective function. A global search for the 10 environmental parameters is performed using genetic algorithms. The inversion algorithm is implemented on the basis of S-band radar sea-clutter data from Wallops Island, Virginia (SPANDAR). Reference data are from range-dependent refractivity profiles obtained with a helicopter. The inversion is assessed (ⅰ) by comparing the propagation predicted from the radar-inferred refractivity profiles with that from the helicopter profiles, (ⅱ) by comparing the refractivity parameters from the helicopter soundings with those estimated. This technique could provide near-real-time estimation of ducting effects.     

Propagation modeling of ocean-scattered low-elevation GPS signals for maritime tropospheric duct inversion

The maritime tropospheric duct is a low-altitude anomalous refractivity structure over the ocean surface,and it can significantly affect the performance of many shore-based/shipboard radar and communication systems.We propose the idea that maritime tropospheric ducts can be retrieved from ocean forward-scattered low-elevation global positioning system(GPS) signals.Retrieval is accomplished by matching the measured power patterns of the signals to those predicted by the forward propagation model as a function of the modified refractivity profile.On the basis of a parabolic equation method and bistatic radar equation,we develop such a forward model for computing the trapped propagation characteristics of an ocean forward-scattered GPS signal within a tropospheric duct.A new GPS scattering initial field is defined for this model to start the propagation modeling.A preliminary test on the performance of this model is conducted using measured data obtained from a 2009-experiment in the South China Sea.Results demonstrate that this model can predict GPS propagation characteristics within maritime tropospheric ducts and serve as a forward model for duct inversion.     

Sensing of the structure of the radio wave refractivity in the troposphere by a network of satellite navigation system receivers in the city of Kazan

We present the results of an experimental study of the structure of the radio-wave refractivity and its dynamics by a network of seven ground-based GPS and GLONASS receivers in Kazan, Russia. It is shown that the remote sensing results agree well with the data of weather stations, radio sounding, and reanalysis. The standard deviation of the refractivity value, which are obtained using the sensing results from the radio sounding data amounts to 2% of the average value at altitudes of up to 500 m. It is found that the refractivity structure has overnight variations, as well as and mesoscale spatial and temporal variability.     

The estimation of lower refractivity uncertainty from radar sea clutter using the Bayesian-MCMC method

The estimation of lower atmospheric refractivity from radar sea clutter(RFC) is a complicated nonlinear optimization problem.This paper deals with the RFC problem in a Bayesian framework.It uses the unbiased Markov Chain Monte Carlo(MCMC) sampling technique,which can provide accurate posterior probability distributions of the estimated refractivity parameters by using an electromagnetic split-step fast Fourier transform terrain parabolic equation propagation model within a Bayesian inversion framework.In contrast to the global optimization algorithm,the Bayesian-MCMC can obtain not only the approximate solutions,but also the probability distributions of the solutions,that is,uncertainty analyses of solutions.The Bayesian-MCMC algorithm is implemented on the simulation radar sea-clutter data and the real radar seaclutter data.Reference data are assumed to be simulation data and refractivity profiles are obtained using a helicopter.The inversion algorithm is assessed(i) by comparing the estimated refractivity profiles from the assumed simulation and the helicopter sounding data;(ii) the one-dimensional(1D) and two-dimensional(2D) posterior probability distribution of solutions.     

基于雷达海杂波的区域性非均匀蒸发波导反演方法

由于不同海域上空气象条件的不同, 海上蒸发波导在大尺度海面上空发生时通常是区域性非均匀的, 这一特性使得该环境中的电波传播特性相对于水平均匀的蒸发波导环境情况而明显不同, 因此, 进行区域性非均匀的蒸发波导探测反演对正确预测电波传播特性及提高雷达系统的工作性能具有重要的意义. 考虑到实际应用中蒸发波导信息获取手段的多样性, 将中尺度数值气象模式MM5预报的区域性蒸发波导修正折射率剖面作为先验信息, 提出了一种含该先验信息的区域性非均匀蒸发波导的雷达海杂波后验概率估计模型. 该模型使用主分量分析法对蒸发波导的水平非均匀性进行参数化建模, 然后通过贝叶斯理论将修正折射率剖面参数的先验概率分布、后验概率分布和似然函数联系起来, 利用雷达海杂波实现蒸发波导剖面参数的最大后验概率估计反演. 通过我国东海海域的实际区域性非均匀蒸发波导反演测试, 表明该模型能够以更高的精度实现区域性非均匀蒸发波导的反演.     

Refractivity estimations from an angle-of-arrival spectrum

This paper addresses the probability of atmospheric refractivity estimation by using field measurements at an array of radio receivers in terms of angle-of-arrival spectrum.Angle-of-arrival spectrum information is simulated by the ray optics model and refractivity is expressed in the presence of an ideal tri-linear profile.The estimation of the refractivity is organized as an optimization problem and a genetic algorithm is used to search for the optimal solution from various trial refractivity profiles.Theoretical analysis demonstrates the feasibility of this method to retrieve the refractivity parameters.Simulation results indicate that this approach has a fair anti-noise ability and its accuracy performance is mainly dependent on the antenna aperture size and its positions.     

Estimation of lower refractivity uncertainty from radar sea clutter using Bayesian-MCMC method

Estimation of lower atmospheric refractivity from radar sea clutter (RFC) is a complicated nonlinear optimization problem. This paper deals with the RFC problem in a Bayesian framework. It uses unbiased Markov Chain Monte Carlo (MCMC) sampling technique, which can provide accurate posterior probability distributions of the estimated refractivity parameters by using an electromagnetic split-step fast Fourier transform terrain parabolic equation propagation model within a Bayesian inversion framework. In contrast to the global optimization algorithm, the Bayesian-MCMC can obtain not only the approximate solutions, but also the probability distributions of the solutions, that is, uncertainty analyses of solutions. The Bayesian-MCMC algorithm is implemented on the simulation radar sea-clutter data and the real radar sea-clutter data. Reference data are assumed to be simulation data and refractivity profiles obtained with a helicopter. Inversion algorithm is assessed (i) by comparing the estimated refractivity profiles from the assumed simulation and the helicopter sounding data; (ii) the one-dimensional (1D) and two-dimensional (2D) posterior probability distribution of solutions.     

基于高光谱数据和模型反演植被叶面积指数的进展

植被叶面积指数(Leaf Area Index , LAI)是陆面过程中影响陆-气交换的重要参数,也是表征植被冠层结构最基本的参量之一。准确而快速地获取LAI是植被-气候相互作用、植被生态和农作物估产研究不可缺少的工作。本文首先针对LAI和高光谱遥感进行概述,然后从不同平台高光谱传感器数据和不同反演方法两个角度总结了国内外近些年来高光谱遥感LAI的研究进展,最后分析了高光谱遥感反演LAI的未来发展方向。     

大气波导条件下雷达海杂波功率仿真

考虑近海面大气折射率垂直梯度的变化, 采用曲面波频谱估计方法(CWSE)计算掠射角, 并结合修正的GIT海表反射率模型和雷达距离方程对大气波导条件下雷达海杂波功率进行仿真, 最后通过数值实验及与实测数据比较, 验证了CWSE方法的可行性. 关键词： 大气波导 雷达海杂波 曲面波频谱估计     

Remote sensing of atmospheric duct parameters using simulated annealing

Simulated annealing is one of the robust optimization schemes. Simulated annealing mimics the annealing process of the slow cooling of a heated metal to reach a stable minimum energy state. In this paper, we adopt simulated annealing to study the problem of the remote sensing of atmospheric duct parameters for two different geometries of propagation measurement. One is from a single emitter to an array of radio receivers (vertical measurements), and the other is from the radar clutter returns (horizontal measurements). Basic principles of simulated annealing and its applications to refractivity estimation are introduced. The performance of this method is validated using numerical experiments and field measurements collected at the East China Sea. The retrieved results demonstrate the feasibility of simulated annealing for near real-time atmospheric refractivity estimation. For comparison, the retrievals of the genetic algorithm are also presented. The comparisons indicate that the convergence speed of simulated annealing is faster than that of the genetic algorithm, while the anti-noise ability of the genetic algorithm is better than that of simulated annealing.     

Refractivity estimation from radar sea clutter

This paper addresses the problem of estimating lower atmospheric refractivity under the nonstandard propagation conditions frequently encountered in low altitude maritime radar applications. The vertical structure of the refractive environment is modeled by using a five-parameter model, and the horizontal structure is modeled as range-independent. The electromagnetic propagation in the troposphere is simulated by using a split-step fast Fourier transform based on parabolic approximation to the wave equation. A global search marked as a modified genetic algorithm (MGA) for the 5 environmental parameters is performed by using a genetic algorithm (GA) integrated with a simulated annealing technique. The retrieved results from simulated runs demonstrate the ability of this method to make atmospheric refractivity estimations. A comparison with the classical GA and the Bayesian Markov Chain Monte Carlo (Bayesian-MCMC) technique shows that the MGA can not only shorten the inverse time but also improve the inverse precision. For real data cases, the inversion values do not match the reference data very well. The inverted profile, however, can be used to synoptically describe the real refractive structure.     

GPS/LEO无线电掩星开环反演技术

在低对流层,传统的闭环模式容易造成较大的信号追踪误差. 开环模式通过内部或外部方式消除导航资料调制(navigation data modulation, NDM)对大气反演的影响,能够准确地恢复信号的相位. 内部和外部方式的区别是,前者通过信号相邻采样的内在联系移除NDM, 后者利用卫星导航位元数据进行处理. 通过内部和外部两种方式处理COSMIC (constellation observing system for meteorology, ionosphere and climate)掩星数据,利用全谱反演方法获得弯曲角,通过Abel积分变换计算折射率. COSMIC掩星个案分析表明,在低对流层,内部方式可能会导致信号的半周跳现象, 从而造成反演的折射率出现误差. 分别用上述两种方式对2007年第71天至73天共约3130个COSMIC掩星开环数据进行处理. 将折射率反演结果与ECMWF(european centre for medium-range weather forecasts)分析场资料进行统计比较,结果显示:外部处理方式的探测深度比内部方式低100 —200 m. 在热带的3 km以下,内部方式具有比外部方式稍大的折射率负偏差. 因此,外部方式能够得到比内部方式略优的反演结果.     

AIRS红外高光谱卫星数据反演卷云光学厚度和云顶高度

基于大气红外探测器L1B红外高光谱辐射观测资料,结合中分辨率成像光谱仪(moderate resolution imaging spectroradiometer, MODIS)云产品数据,利用通用大气辐射传输模式(combined atmospheric radiative transfer model,CART),根据模式模拟和AIRS实际观测亮温的亮温差,研究从AIRS红外波段1 070～1 135 cm-1高光谱数据反演卷云的光学厚度和云顶高度。将反演的卷云光学厚度与云顶高度作为输入参数模拟计算650～1 150 cm-1波段卷云大气顶的辐射亮温谱,并将模拟值与AIRS观测亮温谱进行了对比分析。将反演的卷云光学厚度和云顶高度和AIRS的760通道(900.56 cm-1,11.1 μm)的亮温以及MODIS卷云反射率进行了对比分析。最后将反演的卷云云顶高度和MODIS云顶高度进行了对比分析。研究结果表明：反演所使用的650～1 150 cm-1波段模式模拟和观测亮温谱吻合得很好,说明CART可以较好的模拟AIRS亮温谱。反演的卷云参数与AIRS在大气窗口区的760通道(900.56 cm-1,11.1 μm)的亮温的分布满足低亮温对应较大的卷云光学厚度和高云顶高度。反演的卷云参数和MODIS卷云的反射率分布满足高卷云光学厚度和云顶高度对应高卷云反射率。反演的卷云云顶高度和MODIS的卷云云顶高度之间线性相关系数相对较高,且都在8.5～11.5 km的概率较高,两者的概率分布趋势一致。说明CART可以用于反演卷云的性质,反演结果具有一定的可靠性。     

GPS掩星一维弯曲角算子在四维变分资料同化系统中的实现方法研究

GPS无线电掩星资料的良好特性使其作为十分重要的观测资料在欧洲 及美国数值天气预报中效果明显. 本文设计了一维弯曲角算子在四维变分资料同化系统中的实现方法和质量控制策略. 详细阐述了一维弯曲角正演算子及其切线性/伴随算子在全球四维变分资料同化系统中的实现步骤和过程, 其中折射率计算中采用Rüeger的折射率最佳平均因子的三项表达式, 计算出更加精确的折射率值. 对2012年8月份的GPS无线电掩星资料进行了观测误差统计, 引入了更加符合四维变分资料同化系统实际情况的一维弯曲角观测误差模型和质量控制策略. 同时, 对比了观测误差统计试验中同化GPS无线电掩星资料的模式背景场和未同化GPS无线电掩星资料的背景场的观测误差统计结果, 发现同化GPS无线电掩星资料之后的观测误差统计结果在全球区域正作用十分明显, 尤其是在中低纬度地区, 高度在10 km 以上的区域观测误差减小十分明显. 关键词： GPS无线电掩星 一维弯曲角算子 四维变分同化系统 质量控制     

基于Kramers-Kronig关系建立金属太赫兹色散模型

提出了一种基于测量反射率谱、使用Kramers-Kronig(KK)关系建立金属太赫兹色散模型的方法.结合合金铝和合金铜4—40 THz的测量反射率谱,通过反射系数振幅和相位的KK关系,采用高频端指数外推,低频端常数外推的方法,反演金属复折射率.以KK反演的复折射率作为实验值,以拟合复折射率和实验值误差最小为准则,使用遗传优化算法,拟合了合金铝和合金铜的Drude色散参数(等离子频率和碰撞频率).基于优化的Drude模型计算了0.1—40 THz材料的复折射率,与椭偏仪的实测结果符合,验证了模型的准确性.该方法理论与实验相互验证,以测量的复折射率作为实验定标,将远红外频段的色散信息拓展到太赫兹频域,确定了太赫兹频段金属的微观物理参数,提供了太赫兹频段色散和散射机理的研究依据.     

含突发E层的电离层模型建立及其在测高中的应用

针对高频射线测高模型未考虑电离层突发E层(sporadic-E,Es层)的问题,从电离层物理结构特性入手,结合实测的电离层垂测数据,在多层准抛物模型的基础上研究了含Es层的电离层模型及其在目标高度测量中的应用.首先,利用反转抛物线特性模拟了含Es的电离层模型,并得到了等离子体频率与高度的关系;然后,利用该模型分析了射线的电离层传输路径与发射仰角/频率的关系以及高频射线微多径特征与目标高度的关系;最后,结合含Es的电离层模型与射线微多径特征,提出了基于分段爬山搜索的快速匹配域测高方法,该方法能大大减少搜索时间.研究结果表明:含Es的电离层模型和提出的测高方法能准确估计出目标高度,并具有较强的实时性.     

AR模型辅助SINS／GPS组合导航的实验研究

激光陀螺的随机漂移噪声类似于白噪声,可以采用AR（2）模型进行数字滤波来减小其对激光陀螺精度的影响。首先介绍了AR模型,然后利用IMU单元三组陀螺信号的AR（2）模型给出了车载SINS／GPS组合导航系统误差模型,并在GPS的辅助下,对SINS／GPS组合导航进行了实验研究。最后利用Kalman滤波器对车载SINs／GPS组合导航实际测量数据进行了离线的半实物仿真。实验结果表明,采用AR（2）模型辅助的SINS／GPS组合导航系统能够得到较高的姿态精度,并能有效抑制速度误差的发散。     

双声道波导中低频环境噪声分布特性

针对北极冰下双声道波导的噪声分布特性研究,采用垂直阵对双声道波导内的噪声场进行了观测,北极冰下噪声物理模型使用混合加性噪声模型描述,统计分布模型使用α稳定分布描述.并将噪声划分为3个频段,按噪声来源分为背景噪声、冲击噪声以及航船噪声3种场景,分析其统计特性,给出了直方图统计、核函数估计与正态分布的比较结果,α稳定分布参数估计结果,峰度与偏度的计算结果以及累积分布函数检验结果。经分析讨论,给出结论:双声道波导中的噪声统计分布可用α稳定分布来描述,但上下表面声道的分布参数具有明显差别,上表面波导背景噪声明显偏离高斯分布而下表面声道的背景噪声接近于高斯分布;上下表面声道中不同频带的噪声分布差异也很明显,其中2~30 Hz差异性最大,30~100 Hz差异性较大,而100~500 Hz差异性相对较小,并且可以认为是对称α稳定分布.此外,上下表面声道中不同频带噪声α稳定特征参数变化也存在差异,上表面波导受冲击噪声影响α值变大,拖尾变薄;而下表面波导受冲击噪声影响α值则会变小,拖尾变重。      

“Refractivity-from-clutter” based on local empirical refractivity model

Constructing sophisticated refractivity models is one of the key problems for the RFC(refractivity from clutter)technology. If prior knowledge of the local refractivity environment is available, more accurate parameterized model can be constructed from the statistical information, which in turn can be used to improve the quality of the local refractivity retrievals. The validity of this proposal was demonstrated by range-dependent refractivity profile inversions using the adjoint parabolic equation method to the Wallops' 98 experimental data.