机载微波大气温度探测仪多高度飞行观测试验结果分析

机载微波大气温度探测仪可以机动灵活地获取大气温度廓线信息。针对一次机载微波大气温度探测仪的多高度飞行观测试验,基于逐线积分模式和大气参数廓线库,建立用于不同飞行高度的快速辐射传输模式,分析了仪器观测亮温的质量并对仪器观测进行了订正;建立了基于神经网络的微波大气温度廓线反演算式,分析了不同高度、不同通道选择对于大气温度廓线反演性能的影响。研究结果表明:（1）较低飞行高度计算得到的各地表敏感通道地表比辐射率之间具有较好的一致性;（2）采用订正算式订正后,不同飞行高度的模拟亮温与观测亮温具有较好的一致性;（3）机载微波大气温度反演最优通道组合依赖于平台飞行高度;（4）采用最优的通道组合,4 200 m、3 200 m和2 500 m高度层温度反演均方根误差范围分别为0.5~1.8 K、0.5~1.3 K和0.4~1.0 K。      

利用地基微波辐射计资料反演0-10km大气温湿廓线试验研究

实测与模拟的微波辐射计亮温存在偏差,导致基于BP神经网络模型的大气温湿廓线反演精度的降低。研究了一种基于资料订正后的BP神经网络反演大气温湿廓线的方法。首先,基于2014年6月南京江宁探空资料,利用MonoRTM模式,模拟中心频率在22.2GHz~58.8GHz范围内22通道亮温;对比模拟和实测南京站微波辐射计资料,建立实测微波辐射计资料订正模型。然后,以南京地区2011-2013年探空资料为输入,模拟22通道亮温数据,并基于模拟的22通道亮温数据和当地探空资料,利用BP神经网络算法,建立大气温度、水汽密度、相对湿度廓线反演模型。最后,利用构建的订正模型,对2014年7月试验获取的微波辐射计资料进行订正,并将订正后的微波辐射计资料输入BP神经网络反演模型,反演0-10km高度58层的大气温度、水汽密度和相对湿度,对比实际探空资料以及微波辐射计二级产品,评估分析反演效果。实验结果表明：所建的反演模型提高了大气温湿廓线反演精度,大气温度、水汽密度和相对湿度均方根误差范围分别为1.0~2.0K、0.20 ~1.93g/m3和2.5%~18.6%。     

地基微波辐射计遥感大气廓线的BP神经网络反演方法研究

讨论了12通道地基微波辐射计遥感反演温度、相对湿度和云液态水廓线的BP神经网络反演方法,利用探空资料,对北京春、夏、秋、冬四个季节的大气廓线进行神经网络训练,并对训练好的网络的反演能力进行数值检验,分析了反演精度;对北京南郊观象台12通道微波辐射计的观测亮温资料进行实际反演,结果表明,神经网络(BPNN)反演的廓线与微波辐射计自带RadiomeNN的相比更加接近实际。     

基于一维变分算法的地基微波辐射计遥感大气温湿廓线研究

为研究地基微波辐射计遥感温、湿度廓线的一维变分算法的反演能力,用北京地区2010—2011年00和12时(世界时)的多通道地基微波辐射计亮温资料进行试验。首先,利用同时次的地面观测资料、红外亮温(由地基微波辐射计自带红外传感器测得)及探空观测数据,给出提取无云样本的方案,得到432个无云样本;再以辐射传输模式计算得到的模拟亮温为参考,对无云条件下的观测亮温进行质量控制;然后利用探空数据进行模拟试验,结果发现,一维变分算法对3 km以下的温度廓线有较大调整。使反演结果更加接近探空,而对湿度廓线在0—10 km都有不同程度的优化;最后利用一维变分算法对地基微波辐射计观测亮温进行大气温湿廓线反演,将结果与探空对比可以看出,温度廓线的均方根误差小于2.9 K,绝对湿度的均方根误差小于0.47 g/m~3;进一步与地基微波辐射计自带神经网络的反演结果比较表明,一维变分的反演结果更接近实际大气。     

从118.75 GHz附近六通道亮温反演大气温度廓线的数值模拟研究

为了能在静止气象卫星上实现微波被动遥感探测大气温度廓线,并保持一定的地面空间分辨率(如视场小于60 km),就需要使用高频微波及大天线.欧洲和美国下一代静止气象卫星上都已考虑采用118.75 GHz附近通道.为了充分了解118.75 GHz附近通道遥感反演温度廓线的能力,为仪器研制及今后资料的解释反演提供必要的基础数据,作者开展了采用118.75 GHz附近六个通道遥感反演大气温度廓线的数值模拟研究.统计反演的数值试验表明,118.75 GHz附近六通道对温度垂直分布有一定的遥感反演能力;温度反演较好的层次对应于权重函数峰值所在的位置.     

预估临边探测大气温度和水汽廓线的反演精度

针对在研仪器——大气辐射超高光谱探测仪的临边探测模式,模拟计算了大气温度和水汽的权重函数。以此为基础,利用信息量和权重函数线性化方法,结合仪器的可探测亮温阈值0.3 K,计算并分析6种大气状态下,大气温度和水汽混合比廓线在不同反演精度条件下可获得的光谱通道数,在满足最佳光谱通道数200的要求下,理论上预估其反演精度。温度廓线整体反演精度为0.6 K,水汽混合比廓线反演精度可达到5%,但热带大气在16~20 km高度的水汽廓线反演精度仅为10%。反演精度预估,仅提供了一种全面认识仪器性能的方法,精度的确定还有赖于真实探测数据的获取和反演方法。     

辐射传输模式对地基微波辐射计观测亮温的模拟能力分析

提供高时间分辨率大气温度湿度廓线的地基微波辐射计近年来广泛使用,多通道观测亮温的数据质量是大气廓线产品合理性的基本保障。一般定期液氮绝对定标可以更好维护亮温数据质量,但实际操作颇为不易。辐射传输模式作为一种辅助工具,可以检验和认识地基微波辐射计观测亮温的数据质量。本文针对三个辐射传输模式:MonoRTM、ARTS和MWRT,结合北京探空观测资料、北京观象台和河北香河站同类型的德国RPG地基微波辐射计观测资料,分析比较了三个模式的模拟与观测亮温差异,评估不同辐射传输模式对地基微波辐射计观测的模拟能力。地基微波辐射计14个通道观测亮温与模式模拟的差异统计比较发现:三个模式的模拟结果与地基微波辐射计大部分通道的观测亮温都很接近,与观测结果具有很好一致性（如相关系数高达0.99）,而对温度通道ch8（51.26 GHz）和ch9（52.28 GHz）,三个模式模拟与观测相关系数明显较低（<0.80）,并且存在显著的绝对偏差（4～5 K）,表明模式在这两个通道的模拟能力有待提高。三个模式中,MonoRTM模式在温度通道ch8、ch9和ch10（53.86 GHz）存在明显的系统性偏差,尤其是ch8高达5 K;ARTS模式对水汽通道ch1（22.24 GHz）的模拟能力相对较弱;MWRT模拟与观测亮温在多个通道上相对更为接近和稳定,尤其系统性偏差最小。此外,探空廓线与地基观测站的空间位置不一致,对地基微波辐射计水汽通道的模拟结果影响较为显著,而对水汽不敏感的温度通道影响甚微。两地观测亮温与模式模拟的比对,初步表明北京观象台地基微波辐射水汽通道的观测质量有待改进。     

云区卫星微波通道亮度温度的数值模拟

目前多数资料同化系统中对卫星的观测值都是采用晴空模拟,然而用晴空辐射传输模式模拟云区卫星微波通道的辐射值会造成与观测较大的偏差,导致大量云区卫星资料被直接抛弃无法进入同化系统,因而有必要改进云区卫星辐射亮温的模拟能力,进而提高同化系统中云区卫星资料利用率。以2010年台风“圆规”、“凡比亚”和“鲇鱼”为例,基于先进的微波扫描辐射计AMSR-E观测应用一维变分算法反演台风区域的云宏观参数,包括云液水含量廓线、云冰水含量廓线和雨水含量廓线;然后,以大气温度、湿度廓线及这些反演的云参数作为快速辐射传输模式CRTM的输入参数,模拟AMSR-E各通道的辐射亮温。通过对比晴空、有云两种情况下模拟亮度温度与实际观测亮度温度间的偏差,发现增加云参数作为辅助参数、启动辐射传输的散射模块,可以有效地改进台风外围云区卫星辐射亮温的模拟效果,大幅减少模拟亮温与观测亮温间的偏差,增加了同化进数值预报系统的卫星观测数据量。      

利用神经网络从118.75 GHz附近通道亮温反演大气温度

为了准确快速地从118.75 GHz附近六通道亮温计算大气温度,作者开展了利用人工神经网络技术反演大气温度的数值模拟研究。与线性统计反演算法比较,海面上大气温度反演的总体均方根误差减小17%,陆面上大气温度反演的总体均方根误差减小15%。两种下垫面条件下的温度反演结果表明,近陆面的温度反演结果优于近海面的温度反演结果。另外,对温度廓线垂直结构反演性能的分析结果表明,对于具有较厚逆温层结构的温度廓线,神经网络反演对廓线的复现能力优于线性统计反演。     

重庆不同天气条件下地基微波辐射计探测特征

利用2012年1月至2014年8月重庆沙坪坝站的微波辐射计和探空数据,通过数值模拟检验微波辐射计的亮温精度,并统计分析晴空、有云和降水天气条件下微波辐射计反演产品的变化特征。结果表明:(1)有云时微波辐射计氧气通道53.85、54.00 GHz亮温与探空观测温度相关性较好;晴空和有云时MonoRTM模拟亮温与微波辐射计观测亮温相关性较好。(2)不同天气条件下,微波辐射计反演温度与探空观测值的相关性都较高,降水时4.0 km以下微波辐射计反演温度明显偏高,有云和晴空时3.8 km以下的温度平均绝对误差小于2℃。微波辐射计反演的相对湿度与探空观测值的相关性较同高度层温度的相关性差,有云时1.0~2.6 km高度反演的相对湿度平均误差很小,降水时4.5 km以下平均误差也较小且稳定。降水时4.0 km以下微波辐射计反演的水汽密度平均误差明显偏大,有云时多数高度层平均误差较小。(3)4.2 km以下降水时08:00微波辐射计反演温度的平均误差较大,有云时08:00微波辐射计反演温度和水汽密度的平均误差均较小。说明微波辐射计反演的大气廓线具有可用性,且在稳定大气环境中反演效果更好。     

Retrieving atmospheric temperature profiles from AMSU-A data with neural networks

Back propagation neural networks are used to retrieve atmospheric temperature profiles from NOAA-16 Advanced Microwave Sounding Unit-A (AMSU-A) measurements over East Asia. The collocated radiosonde observation and AMSU-A data over land in 2002-2003 are used to train the network, and the data over land in 2004 are used to test the network. A comparison with the multi-linear regression method shows that the neural network retrieval method can significantly improve the results in all weather conditions. When an offset of 0.5 K or a noise level of ±0.2 K is added to all channels simultaneously, the increase in the overall root mean square (RMS) error is less than 0.1 K. Furthermore, an experiment is conducted to investigate the effects of the window channels on the retrieval. The results indicate that the brightness temperatures of window channels can provide significantly useful information on the temperature retrieval near the surface. Additionally, the RMS errors of the profiles retrieved with the trained neural network are compared with the errors from the International Advanced TOVS (ATOVS) Processing Package (IAPP). It is shown that the network-based algorithm can provide much better results in the experiment region and comparable results in other regions. It is also noted that the network can yield remarkably better results than IAPP at the low levels and at about the 250-hPa level in summer skies over ocean. Finally, the network-based retrieval algorithm developed herein is applied in retrieving the temperature anomalies of Typhoon Rananim from AMSU-A data.     

A COMPARISON OF THE RETRIEVAL OF ATMOSPHERIC TEMPERATURE PROFILES USING OBSERVATIONS OF THE 60 GHZ AND 118.75 GHZ ABSORPTION LINES

The Microwave Temperature Sounder-Ⅱ(MWTS-Ⅱ) and Microwave Humidity and Temperature Sounder(MWHTS) onboard the Fengyun-3 C(FY-3 C) satellite can be used to detect atmospheric temperature profiles. The MWTS-II has 13 temperature sounding channels around the 60 GHz oxygen absorption band and the MWHTS has 8 temperature sounding channels around the 118.75 GHz oxygen absorption line. The data quality of the observed brightness temperatures can be evaluated using atmospheric temperature retrievals from the MWTS-Ⅱ and MWHTS observations. Here, the bias characteristics and corrections of the observed brightness temperatures are described. The information contents of observations are calculated, and the retrieved atmospheric temperature profiles are compared using a neural network(NN) retrieval algorithm and a one-dimensional variational inversion(1 D-var) retrieval algorithm. The retrieval results from the NN algorithm show that the accuracy of the MWTS-Ⅱ retrieval is higher than that of the MWHTS retrieval, which is consistent with the results of the radiometric information analysis. The retrieval results from the 1 D-var algorithm show that the accuracy of MWTS-Ⅱ retrieval is similar to that of the MWHTS retrieval at the levels from 850-1,000 h Pa, is lower than that of the MWHTS retrieval at the levels from 650-850 h Pa and 125-300 h Pa, and is higher than that of MWHTS at the other levels. A comparison of the retrieved atmospheric temperature using these satellite observations provides a reference value for assessing the accuracy of atmospheric temperature detection at the 60 GHz oxygen band and 118.75 GHz oxygen line. In addition, based on the comparison of the retrieval results, an optimized combination method is proposed using a branch and bound algorithm for the NN retrieval algorithm, which combines the observations from both the MWTS-Ⅱand MWHTS instruments to retrieve the atmospheric temperature profiles. The results show that the optimal combination can further improve the accuracy of MWTS-Ⅱ retrieval and enhance the detection accuracy of atmospheric temperatures near the surface.     

GMS-5 估计可降水量的研究

文章证明了由静止气象卫星GMS-5的分裂窗通道和水汽通道亮温反演可降水量的可行性,探讨了GMS-5红外通道亮温与可降水量的关系,建立了由3个通道亮温反演可降水量的经验公式。用60组大气平均廓线,对公式模拟检验误差为0.18 g/cm2,而用实际124组探空和对应的GMS-5亮温资料进行检验,误差0.40 g/cm2。用得到的经验公式可反演大范围的晴空可降水量分布。     

地基微波辐射计探测大气边界层高度方法

采用2013年中国科学院大气物理研究所香河大气综合观测试验站的地基微波辐射计和激光雷达观测数据,以激光雷达探测的大气边界层高度为参考,分别利用非线性神经网络和多元线性回归方法建立微波亮温直接反演大气边界层高度的算法,并对比两种方法的反演能力, 同时分析非线性神经网络算法在不同时段及不同天气状况下反演结果的差异。结果表明:非线性神经网络算法的反演能力优于多元线性回归算法,其反演结果与激光雷达探测的大气边界层高度有较好一致性,冬、春季的相关系数达到0.83,反演精度比线性回归算法约高26%;对于不同时段和不同天气条件,春季的反演结果最好,晴空的反演结果好于云天; 四季和不同天气状况的划分也有利于提高反演精度。     

Retrieval of Sea Surface Wind Speed by One-Dimensional Synthetic Aperture Microwave Radiometer

One-dimensional synthetic aperture microwave radiometers have higher spatial resolution and record measurements at multiple incidence angles. In this paper, we propose a multiple linear regression method to retrieve sea surface wind speed at an incidence angle between 0° ~65°. We assume that a one-dimensional synthetic aperture microwave radiometer operates at frequencies of 6.9, 10.65, 18.7, 23.8 and 36.5 GHz. Then, the microwave radiative transfer forward model is used to simulate the measured brightness temperatures. The sensitivity of the brightness temperatures at 0°~65° to the sea surface wind speed is calculated. Then, vertical polarization channels(VR), horizontal polarization channels(HR) and all channels(AR) are used to retrieve the sea surface wind speed via a multiple linear regression algorithm at 0° ~65°, and the relationship between the retrieval error and incidence angle is obtained. The results are as follows:(1) The sensitivity of the vertical polarization brightness temperature to the sea surface wind speed is smaller than that of the horizontal polarization.(2) The retrieval error increases with Gaussian noise. The retrieval error of VR first increases and then decreases with increasing incidence angle, the retrieval error of HR gradually decreases with increasing incidence angle, and the retrieval error of AR first decreases and then increases with increasing incidence angle.(3) The retrieval error of AR is the lowest and it is necessary to retrieve the sea surface wind speed at a larger incidence angle for AR.     

云对地基微波辐射计反演湿度廓线的影响

利用中国气象局大气探测试验基地的L波段探空数据和微波辐射计观测数据,采用MonoRTM辐射传输模型作为正演亮温模型,BP (back propagation) 神经网络作为反演工具,在由亮温反演大气湿度廓线的过程中,添加与样本匹配的云底高度和云厚度信息,建立新的反演模型,使新反演模型得到的反演湿度廓线和未添加云信息的反演湿度廓线分别与探空数据进行对比,获取两种反演方法各高度层的均方根误差,分析云信息对反演大气湿度廓线的影响。对比结果表明:未添加云信息时,测试样本的反演湿度廓线与探空廓线的相关系数平均值为0.685,而添加云信息后,相关系数平均值为0.805。相比未添加云信息的反演廓线,添加云信息之后多数高度层的均方根误差均有不同程度减小,而在有云以上高度层表现尤为明显。     

微波辐射计温湿廓线反演方法改进试验

为提升地基微波辐射计在不同天气条件下, 特别是云天条件下温湿廓线的反演精度, 利用2011年1月—2016年12月中国气象局北京国家综合气象观测试验基地探空数据, 在微波辐射计反演温湿度廓线的过程中通过区分晴天和云天条件并引入全固态Ka波段测云仪云高及云厚信息, 对反演输入亮温进行质量控制和偏差订正, 建立BP神经网络模型, 采用2017年1月—2018年3月微波辐射计探测数据评估检验, 结果表明:在亮温订正前提下, 晴天温度模型、云天温度模型、晴天相对湿度模型和云天相对湿度模型反演结果与探空的相关系数分别为0.99, 0.99, 0.80和0.78, 均方根误差为2.3℃, 2.3℃, 9%和16%, 较微波辐射计自带产品（LV2产品）减小约0.4℃, 0.3℃, 11%和9%, 准确性提升约30%, 28%, 64%和45%;温度模型偏差在±2℃以内、湿度模型偏差在±20%以内的占比分别为68%, 70%和95%, 78%, 较LV2产品分别提高了7%, 5%和27%, 23%, 其中相对湿度改善明显。可见亮温订正、区分天气类型训练反演模型有利于改善地基微波辐射温湿廓线反演精度。     

EFFECTS OF A CLOUD FILTERING METHOD FOR FENGYUN-3C MICROWAVE HUMIDITY AND TEMPERATURE SOUNDER MEASUREMENTS OVER OCEAN ON RETRIEVALS OF TEMPERATURE AND HUMIDITY

For Microwave Humidity and Temperature sounder (MWHTS) measurements over the ocean, a cloud filtering method is presented to filter out cloud- and precipitation-affected observations by analyzing the sensitivity of the simulated brightness temperatures of MWHTS to cloud liquid water, and using the root mean square error (RMSE) between observation and simulation in clear sky as a reference standard. The atmospheric temperature and humidity profiles are retrieved using MWHTS measurements with and without filtering by multiple linear regression (MLR), artificial neural networks (ANN) and one- dimensional variational (1DVAR) retrieval methods, respectively, and the effects of the filtering method on the retrieval accuracies are analyzed. The numerical results show that the filtering method can improve the retrieval accuracies of the MLR and the 1DVAR retrieval methods, but have little influence on that of the ANN. In addition, the dependencies of the retrieval methods upon the testing samples of brightness temperature are studied, and the results show that the 1DVAR retrieval method has great stability due to that the testing samples have great impact on the retrieval accuracies of the MLR and the ANN, but have little impact on that of the 1DVAR.     

Physical retrievals of over-ocean rain rate from multichannel microwave imagery. Part I: Theoretical characteristics of normalized polarization and scattering indices

Summary Microwave rain rate retrieval algorithms have most often been formulated in terms of the raw brightness temperatures observed by one or more channels of a satellite radiometer. Taken individually, single-channel brightness temperatures generally represent a near-arbitrary combination of positive contributions due to liquid water emission and negative contributions due to scattering by ice and/or visibility of the radiometrically cold ocean surface. Unfortunately, for a given rain rate, emission by liquid water below the freezing level and scattering by ice particles above the freezing level are rather loosely coupled in both a physical and statistical sense. Furthermore, microwave brightness temperatures may vary significantly (30–70 K) in response to geophysical parameters other than liquid water and precipitation. Because of these complications, physical algorithms which attempt to directly invert observed brightness temperatures have typically relied on the iterative adjustment of detailed microphysical profiles or cloud models, guided by explicit forward microwave radiative transfer calculations.In support of an effort to develop a significantly simpler and more efficient inversion-type rain rate algorithm, the physical information content of two linear transformations of single-frequency, dual-polarization brightness temperatures is studied: thenormalized polarization difference P of Petty and Katsaros (1990, 1992), which is intended as a measure of footprint-averaged rain cloud transmittance for a given frequency; and ascattering index S (similar to the polarization corrected temperature of Spencer et al., 1989) which is sensitive almost exclusively to ice. A reverse Monte Carlo radiative transfer model is used to elucidate the qualitative response of these physically distinct single-frequency indices to idealized 3-dimensional rain clouds and to demonstrate their advantages over raw brightness temperatures both as stand-alone indices of precipitation activity and as primary variables in physical, multichannel rain rate retrieval schemes.As a byproduct of the present analysis, it is shown that conventional plane-parallel analyses of the well-known footprint-filling problem for emission-based algorithms may in some cases give seriously misleading results.With 11 Figures     

Latitudinal and Scan-dependent Biases of Microwave Humidity Sounder Measurements and Their Dependences on Cloud Ice Water Path

The relationship between differences in microwave humidity sounder(MHS)–channel biases which represent measured brightness temperatures and model-simulated brightness temperatures, and cloud ice water path(IWP) as well as the influence of the cloud liquid water path(LWP) on the relationship is examined. Seven years(2011–17) of NOAA-18 MHS-derived measured brightness temperatures and IWP/LWP data generated by the NOAA Comprehensive Large Array-data Stewardship System Microwave Surface and Precipitation Products System are used. The Community Radiative Transfer Model, version2.2.4, is used to simulate model-simulated brightness temperatures using European Center for Medium-Range Weather Forecasts reanalysis data as background fields. Scan-angle deviations of the MHS window channel biases range from-1.7 K to1.0 K. The relationships between channels 2, 4, and 5 biases and scan angle are symmetrical about the nadir. The latitudedependent deviations of MHS window channel biases are positive and range from 0–7 K. For MHS non-window channels,the latitudinal deviations between measured brightness temperatures and model-simulated brightness temperatures are larger when the detection height is higher. No systematic warm or cold deviations are found in the global spatial distribution of difference between measured brightness temperatures and model-simulated brightness temperatures over oceans after removing scan-angle and latitudinal deviations. The corrected biases of five different MHS channels decrease differently with respect to the increase in IWP. This decrease is stronger when LWP values are higher.