后向散射辐射探测大气臭氧总量反演方法的改进

以SBDART模式作为辐射传输计算模式,在各种大气气溶胶存在的情况下利用现有的大气后向散射反演大气臭氧总量的Version7方法,模拟计算得到各种气溶胶存在时的大气臭氧总量,并分析了反演结果误差产生的原因.进一步提出根据斜气柱臭氧总量sΩ0,选择不同连续波长区间测量得到的后向散射强度来对臭氧总量的初估值Ω0进行修正的方法.将得到的结果与Version7的进行比较,证明了新方法的有效性.      

紫外后向散射探测臭氧总量的反演误差研究

通过探测紫外后向散射获取臭氧总量目前已经成为臭氧总量探测的主要方法之一. 本文以模拟计算为主要方法, 从模式计算误差, 反演算法误差和仪器误差三方面着手, 对影响臭氧总量反演精度的各个参数进行分析研究, 并对由其引起的臭氧总量反演误差进行定量计算, 研究结果对臭氧总量反演精度的提高和实际观测资料的处理有着重要的指导意义.      

大气臭氧地基遥感反演的退卷积方法

本文提出了一种大气臭氧地基遥感信息处理的新方法,它利用非线性退卷积方法和等效仪器函数提高光谱分辨率,利用非线性最小二乘法进行大气臭氧反演计算,可以显着提高大气臭氧遥感的可探测高度和精度。      

中国不同地区剩余臭氧及臭氧探空订正因子

根据中国不同地点臭氧探空数据，研究气球炸点臭氧浓度定值（CMR）法、卫星（SBUV和MLS）纬向平均法确定的剩余臭氧Ω_res及其对订正因子C_ref的影响，同时检验臭氧垂直分布对C_ref的贡献.结果显示: CMR法对气球炸点高度依赖性明显，且易高估Ω_res使C_ref整体低于100%；卫星纬向平均Ω_res对气球炸点高度不敏感，但在中国东部的臭氧总量高值区或青藏高原及低纬度臭氧低值地区，Ω_res呈现近10DU以上低值，这是经向臭氧总量及其垂直分布差异在卫星遥感数据上的反映.地面到100hPa的对流层臭氧（Ω_tro），100~10hPa的平流层臭氧（Ω_str）以及10hPa以上的Ω_res对C_ref贡献平均分别为（16±3.4）%，（65±2.3）%，（19±3.3）%.表明基于C_ref评估或订正探空仪平流层臭氧测值时，需考虑对流层臭氧及确定Ω_res方法的影响.卫星纬向平均法，特别是近似实测的SBUV臭氧廓线的值适用于确定Ω_res.      

太赫兹火星大气临边探测仿真研究

对火星大气进行连续高分辨率观测是研究火星大气物理和化学过程的重要手段.太赫兹临边探测技术通过测量火星大气中的风和光化学循环中的重要气体（CO,O₃,H₂O,H₂O₂等）提高对火星的认知.针对火星大气遥感的探测需求,分析了300～1000GHz频段的频谱特征.针对探测卫星对于载荷质量、功耗等参数的要求,提出一个560GHz频段的火星大气太赫兹临边探测仪设计方案,并利用辐射传输模型ARTS中的行星工具箱进行仿真.仿真结果显示:火星大气温度的反演精度优于4K,其中45km高度以下优于2K;H₂O丰度的反演精度在90km以下优于50%,30km以下优于2%;H₂O₂的反演精度在40km以下优于50%;O₃的反演精度在50km以下优于60%;大气风速度的反演精度在65km以上优于5m·s-1,最高可以达到2m·s-1.研究结果表明,利用太赫兹波段的吸收谱线可以很好地探测火星大气中各成分的丰度、变化趋势以及中高层大气的风,可为后续火星表面及大气探测提供参考.      

小尺度重力波引起夜光云反照率变化的统计特征

利用AIM卫星搭载的CIPS云成像探测器获得的云图数据,提取2008-2009年南北半球共6664个小尺度重力波（波长10～150km）个例,通过重力波区域与背景云层反照率变化值的对比分析,研究重力波引起云层反照率的变化特征.结果表明,重力波引起的反照率变化值以正值为主,最大平均值4.48×10-6sr-1出现在南半球降交轨道.反照率变化值与IWC变化值正相关,相关系数均在0.85以上.重力波引起的反照率变化呈现出很强的纬度和时间依赖性,且几乎均为正值.反照率变化值在中期阶段（冬/夏至日之后的50天）的高纬地区（>70°）更大,但在中期以外始末阶段的低纬地区（<70°）逐渐变小,甚至开始出现负值.随着背景云层的增强,反照率平均值呈线性增大,小尺度重力波能够引起背景云层反照率约14.6%～28.8%的变化量.当重力波引起的反照率周期性变化的振幅逐渐增大时,反照率变化值也线性增大,变化率约为0.909%～1.194%.南半球的变化率整体比北半球稍小,这与背景大气条件的差异有关      

基于FY-3C/MWHTS观测资料反演中低纬度海面气压

利用风云三号C星（FY-3C）微波温湿探测仪（MWHTS）的实测亮温数据，开展了中低纬度（40°S-40°N）区域海面气压反演研究.MWHTS 118.75GHz氧气通道的辐射亮温测量值与氧气气柱总量密切相关，可用于反演海面气压.根据辐射传输方程分析了MWHTS 8个氧气通道对海面气压的敏感性.结果表明，与位于氧气吸收带中心的通道相比，位于吸收带翼区的探测通道对海面气压的变化更敏感.基于神经网络方法建立了中低纬度海面气压反演算法，通过将反演结果与ERA-Interim再分析数据以及原位观测数据进行对比分析，发现建立的反演算法在中低纬度晴空、云天、雨天条件下，对海面气压的估计精度分别为2.0，3.0和3.5hPa.最后，开展了生成初期热带气旋的反演试验，结果表明反演的海面气压资料对热带低压的判别有一定帮助.      

地基Fabry-Perot中高层大气风速反演及误差分析

基于子午工程地基法布里-帕罗干涉仪（Fabry-Perot Interferometer,FPI）的气辉观测数据,结合地基独特的观测模式（天顶角为0° 的天顶方向和天顶角为45°的东西北南四个方向）对地基中高层大气风速进行反演,包括数据预处理、干涉环圆心确定、干涉环半径计算和风速反演. 将2010年5月6-13日8天十个环（十个干涉环同时参与反演）的反演结果与地基FPI风速实测数据进行比较,得到557.7nm,630.0nm,892.0nm三种谱线气辉的反演平均偏差分别为2.7m·s-1,5.5m·s-1,7.7m·s-1. 此外,基于反演算法对上述反演精度影响因素进行了分析. 研究发现,气辉辐射强度对风速的反演精度影响较大,气辉辐射越强,外环的半径计算精度越高,可参与的反演环数越多,则最终的风速反演精度越高. 而圆心偏差± 2pixel（五个环）和± 1pixel（十个环）及焦距变化（±10mm）对风速反演精度的影响相对较小,但当超出这一偏差范围,风速反演偏差会迅速增大.      

临近空间大气环境落球探测中的科氏力影响

气象火箭落球探测技术是临近空间大气环境探测的重要方法.在落球探测数据处理过程中,通常忽略科氏力项的影响.本文利用经验预报模式构建落球探测正演仿真模型,并根据落球探测原理建立参数反演模型,在此基础上仿真模拟了落球探测数据处理过程中忽略科氏力项对大气参数反演精度的影响.在95～100 km高度范围内,忽略科氏力项将引起温度...     

风云-3A卫星TOU臭氧总量地基对比验证分析

通过地基臭氧总量观测,对风云三号(FY-3A)气象卫星TOU(2009年7月至2013年12月)臭氧总量开展检验分析.结果显示,TOU臭氧总量与地基测值之间的线性相关系数达到0.95以上,相对差别(RD)为(-0.16+4.3)%.赤道地区RD随地基测值的上升明显,太阳天顶角(SZA)为70°以上时TOU臭氧总量明显下降,75°~85°有一4%~-10%不等的降幅.云量变化(0~100%)或云高变化(1~15 km)对RD的影响在±2%以内.RD随着吸收性气溶胶(AI为正值)有-2%~-5%的变化,在南半球中纬度地区明显.TOU与地基Dobson,Brewer及SAOZ光谱仪测值的RD分别为(0.6±0.17)%,(0.05±0.2)%,(0.04±0.01)%.当地基臭氧总量在200 DU以下时,TOU值低于Dobson或SAOZ测值;而在430~500 DU时,TOU高于SAOZ值2%~6%.TOU分别与Dobson,Brewer,SAOZ的RD对SZA(70°以下)、云量变化均不敏感,但TOU分别与Brewer,SAOZ的RD在云高5 km以上时或随AI的增加均出现一定的趋势特征.     

First results of ozone profiles between 35 and 65 km retrieved from SCIAMACHY limb spectra and observations of ozone depletion during the solar proton events in October/November 2003

Ozone density profiles between 35 and 65 km altitude are derived from scattered sunlight limb radiance spectra measured by the SCIAMACHY instrument on the Envisat satellite. The method is based on the inversion of normalized limb radiance profiles in the Hartley absorption bands of ozone at selected wavelengths between 250 and 310 nm. It employs a non-linear Newtonian iteration version of Optimal Estimation (OE) coupled with the radiative transfer model SCIARAYS. The limb scatter technique combined with a classical OE retrieval in the short-wave UV-B and long-wave UV-C delivers reliable results as shown by a first comparison with MIPAS V4.61 profiles yielding agreement within 10% between 38 and 55 km. An overview of the methodology and an initial error analysis are presented. Furthermore the effect of the solar proton storm between 28 October and 6 November 2003 on the ozone concentration profiles is shown. They indicate large depletion of ozone of about 60% at 50 km in the Northern hemisphere, a weaker depletion in the Southern hemisphere and a dependence of the depletion on the Earth’s magnetic field.     

平流层臭氧和辐射场的季节分布特征

利用美国NCAR化学气候耦合模式WACCM3对平流层温度场、风场、臭氧及辐射场进行了模拟.结果表明,在适宜飞艇长期驻留的准零风层高度20～22km(对应大气压强范围为50～30hPa,以下均采用气压值表征对应大气高度),7-8月风速小于5m·s-1的风带可长期稳定在40°N以北.臭氧空间分布显示,在30hPa气压高度处中国地区臭氧浓度出现了带状分布,30hPa高度以下低纬度地区臭氧浓度低于中纬度地区.平流层太阳加热率的时空变化表明,在平流层上层,太阳加热率可达100×10-6K·s-1,而在平流层下层,只有10×10-6K·s-1.6-8月中国区域的太阳加热率大于9月;在100～30hPa高度内,中纬度地区太阳加热率高于低纬度地区,在30hPa高度以上,低纬度地区太阳加热率高于中纬度地区;8-9月30～40hPa高度处,太阳加热率的空间变化较小.在30hPa高度上,太阳加热率在40°N昼夜变化最大;50hPa高度处,太阳加热率的昼夜变化小于30hPa高度处,而且白天太阳加热率出现极大值的纬度明显靠北.平流层低纬度地区的长波加热率小于中纬度地区.青藏高原由于地形特殊,其6-7月的臭氧浓度、太阳加热率和长波加热率均小于同纬度其他地区.      

基于蒙特卡罗频率法的葡萄籽总酚含量高光谱测量变量选择

在利用高光谱建立葡萄籽总酚含量的预测模型中,为解决变量过多、模型复杂度高等问题,需依据光谱特点进行有效地数据降维。提出了一种蒙特卡罗频率法（MCF）对高光谱数据进行波长选择,并建立了葡萄籽总酚的支持向量回归（SVR）预测模型。该方法首先采用蒙特卡罗采样（MCS）选择波长子集;然后建立大量SVR子模型,并选出均方根误差（RMSE）较小的子模型,统计每个波长出现的频次;最后根据指数递减函数确定波长个数,选取频次最高的波长子集作为特征波长。结果表明,采用MCF可以在降维的同时提高模型的预测性能,波长数目由原始的196个减少到9个,波长范围均在950~1 400 nm,RMSE值从0.42减少到0.37,预测精度优于SPA等其他波长选择方法。因此,提出的基于MCF在高光谱数据处理中能有效选择特征波长,为准确建立预测模型提供了一种有效的方法。      

Synergetic retrieval of terrestrial AOD from MODIS images of twin satellites Terra and Aqua

Aerosol optical depth (AOD) is one of the most important indicators of atmospheric pollution. It can be retrieved from satellite imagery using several established methods, such as the dark dense vegetation method and the deep blue algorithm. All of these methods require estimation of surface reflectance prior to retrieval, and are applicable to a certain pre-designated type of surface cover. Such limitations can be overcome by using a synergetic method of retrieval proposed in this study. This innovative method is based on the fact that the ratio K of surface reflectance at different angles/geometries is independent of wavelength as reported by Flowerdew and Haigh (1995). An atmospheric radiative transfer model was then established and resolved with the assistance of the ratio K obtained from two Moderate Resolution Imaging Spectroradiometer (MODIS) spectral bands acquired from the twin satellites of Terra and Aqua whose overpass is separated by three hours. This synergetic method of retrieval was tested with 20 pairs of MODIS images. The retrieved AOD was validated against the ground observed AOD at the Taihu station of the AErosol RObotic NETwork (AERONET). It is found that they are correlated with the observations at a coefficient of 0.828 at 0.47 μm and 0.921 at 0.66 μm wavelengths. The retrieved AOD has a mean relative error of 25.47% at 0.47 μm and 24.3% at 0.66 μm. Of the 20 samples, 15 and 17 fall within two standard error of the line based observed AOD data on the ground at the 0.47 μm and 0.66 μm, respectively. These results indicate that this synergetic method can be used to reliably retrieve AOD from the twin satellites MODIS images, namely Terra and Aqua. It is not necessary to determine surface reflectance first.     

Vertical wavenumber spectra of atmospheric ozone measured from ozonesonde observations

Vertical profiles of ozone have been measured at balloon altitudes. Our purpose is to examine the character of vertical wavenumber spectra of ozone fluctuations, to assess the possible roles of gravity wave field in ozone fluctuations, and to determine dominant vertical wavelengths of ozone spectra. Vertical wavenumber spectra of 12 ozone fluctuations obtained during June–August 2003 are presented. Results indicate that mean spectral slopes in the wavenumber range from 4.69 × 10⁻⁴ to 2.50 × 10⁻³ cyc/m are about −2.91 in the troposphere and −2.87 in the lower stratosphere, which is close to the slope of −3 predicted by current gravity wave saturation models. The consistency of the observed spectral slopes with the value of −3 predicted by current gravity wave saturation models suggests that the observed ozone fluctuations are due primarily to atmospheric gravity waves. At m = 1/(1000 m) the mean spectral amplitude is over 30 times larger in the lower stratosphere than in the troposphere. Mean vertical wavenumber spectra in area-preserving form reveal dominant vertical wavelengths of ∼2.6 km in the troposphere and ∼2.7 km in the lower stratosphere, which is consistent with the values varying between 1.5 and 3.0 km estimated from the velocity field and temperature field at these heights.     

Observations of the ozone and nitrogen dioxide profiles in the TROICA-4 experiment

To investigate the vast area of Russia, a mobile scientific facility based in a railway carriage was developed. It is capable to perform continues measurements being coupled in a passenger train traveling along railroads. It was first equipped with a spectrometer for remote sensing of ozone and nitrogen dioxide in the atmosphere for the transcontinental observations into the chemistry of the atmosphere-4 expedition performed from February 18 to March 5, 1998. A twilight DOAS method, which was applied for retrieval of the nitrogen dioxide profiles basing on spectral measurements at the visible wavelengths (434–451 nm), is described in the paper. Main features of a new algorithm for retrieval of the ozone profile and total content using the differential structure of the UV spectrum (310–335 nm) are presented. The ozone and nitrogen dioxide contents are obtained and shortly validated against available alternative data.     

Spaceborne lidar measurement accuracy: Simulation of aerosol,cloud, molecular density,and temperature retrievals

Spaceborne lidar measurements and retrievals are simulated using realistic errors in signal, conventional density information, atmospheric transmission, and lidar calibration. We find that by day, independent analysis of returns at wavelengths of 0.53 and 1.06 μm yields vertical profiles (0.1- to 1-km resolution) of tenuous clouds and boundary-layer, Saharan, and strong volcanic stratospheric aerosols to accuracies of 30% or better, provided particulate optical depth does not exceed ?0.3. By night all these constituents are retrieved, plus noctilucent clouds, mesospheric aerosols, and upper tropospheric/nonvolcanic stratospheric (UT/NVS) areosols. Molecular-density uncertainties are a dominant source of error for UT/NVS retrievals.To reduce these errors and also to provide density and temperature profiles, we developed a procedure that combines returns at 0.35 and 1.06 μm. This technique significantly improves UT/NVS aerosol retrieval accuracy and also yields useful density and temperature profiles there. Strong particulate contamination limits the technique to the cloud-free upper troposphere and above.