Ensemble retrieval of atmospheric temperature profiles from AIRS

ABSTRACT Satellite-based observations provide great opportunities for improving weather forecasting. Physical retrieval of atmo spheric profiles from satellite observations is sensitive to the uncertainty of the first guess and other factors. In order to improve the accuracy of the physical retrieval, an ensemble methodology was developed with an emphasis on perturbing the first guess. In the methodology, a normal probability density function （PDF） is used to select the optimal profile from the ensemble retrievals. The ensemble retrieval algorithm contains four steps： （1） regression retrieval for original first guess; （2） perturbation of the original first guess to generate new first guesses （ensemble first guesses）; （3） using the ensemble first guesses and nonlinear iterative physical retrieval to generate ensemble physical results; and （4） the final optimal profile is selected from the ensemble physical results by using PDE Temperature eigenvectors （EVs） were used to generate the pertur- bation and generate the ensemble first guess. Compared with the regular temperature profile retrievals from the Atmospheric InfraRed Sounder （AIRS）, the ensemble retrievals RMSE of temperature profiles selected by the PDF was reduced between 150 and 320 hPa and below 400 hPa, with a maximum improvement of 0.3 K at 400 hPa. The bias was also reduced in many layers, with a maximum improvement of 0.69 K at 460 hPa. The combined optimal （CombOpt） profile and a mean optimal （MeanOpt） profile of all ensemble physical results were improved below 150 hPa. The MeanOpt profile was better than the CombOpt profile, and was regarded as the final optimal （FinOpt） profile. This study lays the foundation for improving temperature retrievals from hyper-spectral infrared radiance measurements.     

卫星热红外遥感技术在火山区地热探测中的应用——以内蒙古锡林郭勒火山区为例

能源的全球性短缺以及全球环境问题的日益加剧,使得绿色新能源的研究与开发成为当今世界许多国家的一项重要任务。地热资源作为其中一员其重要性日益凸显。在地热资源的勘查中,传统的地质和地球物理勘探投入高、难度大,而遥感技术具有经济、高效、覆盖范围广等优势,已逐渐成为地热资源探测与调查的一种新技术手段。地热能开发的首选地区往往是新生代火山活动区,因其在相对较小的深度内就能获得很高的温度。本文综述了国内外利用遥感技术探测火山地热的发展历史和研究现状,介绍了这项技术在美国黄石公园、日本阿苏火山以及我国长白山火山和腾冲火山的应用。同时选取内蒙古锡林郭勒盟新生代火山活动构造区作为实例,利用美国NASA陆地卫星Landsat-7 ETM+第6波段的红外遥感数据反演了研究区的地表温度,结合ETM+7-4-1波段组合影像的地质构造解译结果以及其它相关资料对研究区地表温度分布进行了分析。结果表明,区内红外遥感地表温度异常区多位于北东东与北西向构造的交汇处或与环形构造密切相关。本文圈定了7个地表温度异常区作为这一地区下一步地热资源勘查的远景区。     

Retrieval of wheat leaf area index from AWiFS multispectral data using canopy radiative transfer simulation

Accurate representation of leaf area index (LAI) from high resolution satellite observations is obligatory for various modelling exercises and predicting the precise farm productivity. Present study compared the two retrieval approach based on canopy radiative transfer (CRT) method and empirical method using four vegetation indices (VI) (e.g. NDVI, NDWI, RVI and GNDVI) to estimate the wheat LAI. Reflectance observations available at very high (56 m) spatial resolution from Advanced Wide-Field Sensor (AWiFS) sensor onboard Indian Remote Sensing (IRS) P6, Resourcesat-1 satellite was used in this study. This study was performed over two different wheat growing regions, situated in different agro-climatic settings/environments: Trans-Gangetic Plain Region (TGPR) and Central Plateau and Hill Region (CPHR). Forward simulation of canopy reflectances in four AWiFS bands viz. green (0.52–0.59 μm), red (0.62–0.68 μm), NIR (0.77–0.86 μm) and SWIR (1.55–1.70 μm) were carried out to generate the look up table (LUT) using CRT model PROSAIL from all combinations of canopy intrinsic variables. An inversion technique based on minimization of cost function was used to retrieve LAI from LUT and observed AWiFS surface reflectances. Two consecutive wheat growing seasons (November 2005–March 2006 and November 2006–March 2007) datasets were used in this study. The empirical models were developed from first season data and second growing season data used for validation. Among all the models, LAI-NDVI empirical model showed the least RMSE (root mean square error) of 0.54 and 0.51 in both agro-climatic regions respectively. The comparison of PROSAIL retrieved LAI with in situ measurements of 2006–2007 over the two agro-climatic regions produced substantially less RMSE of 0.34 and 0.41 having more R² of 0.91 and 0.95 for TGPR and CPHR respectively in comparison to empirical models. Moreover, CRT retrieved LAI had less value of errors in all the LAI classes contrary to empirical estimates. The PROSAIL based retrieval has potential for operational implementation to determine the regional crop LAI and can be extendible to other regions after rigorous validation exercise.     

中性大气掩星标准反演技术

结合JPL和GFZ的中性大气标准反演算法，发展了上海天文台中性大气掩星标准反演算法流程．对CHAMP观测资料进行中性大气反演，并将反演结果与ECMWF(the European Center for Medium Range Weather Forecasting)的预报大气数据进行比较．作为讨论，指出将来的发展方向．     

基于贝叶斯网络的Web访问模式挖掘模型系统

根据用户的访问记录,以及每次检索输入的关键词、用户的动作来动态更新用户的兴趣点.提出了一个能较好吻合这一思想的数据挖掘模型,以用户的停留时间、访问频率和用户输入的检索关键词为数据对象,进行用户兴趣分析,动态跟踪用户需求意向,简化检索过程,提高WEB资源的使用效率.     

基于星载激光雷达ICESat-2的雪深监测方法研究

目前,被动微波数据是积雪深度反演的主要数据源,受其较粗空间分辨率的影响,反演雪深存在较大的不确定性。激光雷达由于其较高的测高精度,在雪深监测方面具有一定的潜力,基于星载激光雷达ICESat-2数据对北疆地区2018年10月至2019年9月积雪季的雪深进行了提取。由于很难获取ICESat-2轨迹点的雪深观测资料,因此首先利用地面雪深观测数据对目前流行的被动微波雪深反演产品进行验证,获取精度可靠的雪深产品并与发展的ICESat-2监测雪深数据产品进行对比。结果表明：AMSR2雪深产品在北疆地区误差较大,整体存在高估现象,中国雪深长时间序列雪深产品精度相对可靠,以作为对比ICESat-2模拟雪深的参考数据;ICESat-2雪深与中国雪深长时间序列雪深产品在空间上以及变化趋势方面吻合度较高,但ICESat-2雪深变化更加连续,说明ICESat-2激光雷达数据不但可以提取区域积雪的深度,对积雪深度的空间变化也比被动微波数据更加敏感,可以获取更加详细的积雪深度空间变化细节,为精细化的积雪深度空间分布提供数据支撑。     

一次1月山东半岛东部极端海效应暴雪的发生机制分析

利用海上浮标站、自动站、多普勒天气雷达、L波段雷达探空、NCEP/NCAR再分析逐6 h和降水等观测资料,结合EVAP雷达风场反演获得的水平风场资料,对2018年1月9—11日一次渤海海效应暴雪过程的产生机制进行了分析。结果表明:此次海效应暴雪过程是一次极端降雪事件,具有强降雪持续时间长、降雪量大、暴雪分布近γ中尺度等特征。暴雪发生前后受两次强冷空气影响,渤海和山东半岛地区持续降温,850 hPa温度降至-18～-16℃,是产生强海效应降雪的有利条件;此次冷空气明显强于12月渤海海效应暴雪,1月产生海效应暴雪的850 hPa温度中位数较12月低约5℃。受强冷空气影响时,海气温差明显增大,海洋向低层大气输送的最大感热通量可达226.8 W·m~(-2),低层大气高湿饱和,导致大气层结不稳定,不稳定局限于850 hPa以下,为浅层对流。雷达反射率因子图上具有明显的"列车效应"。造成窄带回波的原因在于出现了低层切变线,即在山东半岛北部沿海的小范围区域内出现了东北风及西西南风,形成了西北风与东北风、西西南风与东北风的切变线,触发暴雪产生。而东北风达到的高度不超过1.2 km,多在0.6 km以下。通过此次极端暴雪过程的综合观测资料分析,揭示了较少出现的1月海效应暴雪的特征,其形成的环流形势、热力不稳定、动力条件等与常见的12月海效应暴雪基本类似,主要差异在于冷空气强度较12月偏强,这可成为1月海效应暴雪的首要预报着眼点;海上浮标站、雷达风场反演技术是定量揭示海效应暴雪中尺度特征的良好资料和方法。     

Raindrop Size Distribution Parameters Retrieved from Xinfeng C-Band Polarimetric Radar Observations

The objective of this research was to acquire a raindrop size distribution (DSDs) retrieved from C-band polarimetric radar observations scheme for the first time in south China.An observation period of the precipitation process was selected,and the shape-slope (μ-Λ) relationship of this region was statistically analyzed using the raindrop sample observations from the two-dimensional video disdrometer (2DVD) at Xinfeng Station,Guangdong Province.Simulated data of the C-band polarimetric radar reflectivity Z_(HH)and differential reflectivity Z_(DR)were obtained through scattering simulation.The simulation data were combined with DSD fitting to determine the Z_(DR)-Λand log_(10)(Z_(HH)/N_0)-Λrelationships.Using Xinfeng C-band polarimetric radar observations Z_(DR)and Z_(HH),the raindrop Gamma size distribution parametersμ,Λ,and N_0were retrieved.A scheme for using C-band polarimetric radar to retrieve the DSDs was developed.This research revealed that during precipitation process,the DSDs obtained using the C-band polarimetric radar retrieval scheme are similar to the 2DVD observations,the precipitation characteristics of rainfall intensity (R),mass-weighted mean diameter (D_m) and intercept parameter (N_w) with time obtained by radar retrieval are basically consistent with the observational results of the 2DVD.This scheme establishes the relationship between the observations of the C-band polarimetric radar and the physical quantities of the numerical model.This method not only can test the prediction of the model data assimilation system on the convective scale and determine error sources,but also can improve the microphysical precipitation processes analysis and radar quantitative precipitation estimation.The present research will facilitate radar data assimilation in the future.     

Kinematic features of a bow echo in southern China observed with doppler radar

A bow echo is a type of mesoscale convective phenomenon that often induces extreme weather and appears with strong reflectivity on radar images. A strong bow echo that developed from a supercell was observed over Foshan City in southern China on 17 April 2011. The intense gusty winds and showers caused huge losses of property and severely affected human lives. This paper presents an analysis of this strong meso- n-scale convective system based on Doppler radar observations. The isolated bow echo exhibited a horizontal scale of about 80 km in terms of reflectivity above 40 dBZ, and a life span of 8 hours. The system originated from the merging of a couple of weakly organized cells in a shear line, and developed into an arch shape as it moved through the shear zone. Sufficient surface moisture supply ensured the convective instability and development of the bow echo. The low-altitude winds retrieved from single Doppler radar observations showed an obvious rear-inflow jet along the notch area. Different from the conventional definition, no book- end anticyclone was observed throughout the life cycle. Very strong slantwise updrafts and downdrafts were recognizable from the retrieved winds, even though the spatial scale of the bow echo was small. Strong winds and induced damage on the surface are considered to have been caused by the mid-level rear-inflow jet and intense convective downdrafts.     

从NOAA卫星AVHRR资料反演中国区域地表反照率

由NOAA卫星AVHRR短波通道1、2反射率反演地表反照率需要3个反演模式,分别是窄-宽波段反射率转换模式、大气顶双向反射模式、大气订正模式.基于模式和国家卫星气象中心接收处理的NOAA-18 AVHRR1B数据,处理了2006年1月至2010年12月的中国区域地表反照率,由于云的影响,15天合成技术用来形成周期为15天的地表反照率数据文件.2006年、2010年2年的处理结果与MODIS同类产品对比,RMS为0.028～0.074、相关系数为0.76～0.93,误差较大出现在冬季,原因是两者15天合成方法不同;5年的日平均地表反照率与21个中国地面气象一级辐射站的观测测值作对比,结果是：RMS为0.053、相关系数为0.88.反演模式系统误差以及云和气溶胶影响是卫星反演地表反照率的主要误差来源.