星载雷达在订正地基天气雷达标定误差中的应用

在雷达反射率数据定量应用中,标定误差是导致结果产生偏差的重要原因之一。星载雷达（TRMM PR）长期工作的稳定性和连续性已被验证,本文将星载雷达数据转换到S波段,通过对比星载雷达和南京雷达同时段不同高度（相对于融化层的位置）不同降水类型（层云降水或对流降水）的数据,得到两部雷达对融化层以下层云降水的观测相关性高,差异稳定。通过对比分析星载雷达和南京雷达同时段零度层以下层状云降水的观测数据,得出两部雷达反射率因子值的回归关系式。以星载雷达观测数据为基准,使用该关系式对南京雷达反射率数据进行订正,并通过雨量计数据对订正结果进行验证,结果显示使用本文订正关系式订正后雷达估计的降雨量更接近雨量计的观测值。     

94/220 GHz星载雷达双波长比对非球形冰晶云参数敏感性分析

利用不同形状冰晶的散射特性,获得了非球形冰晶云的94/220 GHz测云雷达双波长比,探讨了非球形冰晶云的双波长比与云内微物理参数的关系,分析了衰减前后的星载雷达反射率因子及双波长比的垂直廓线。结果表明:（1）双波长比可以反映小到0.1 mm中值尺度的冰粒子,对粒子总数、谱的形状参数不敏感,对粒子大小、形状、云衰减较敏感。（2）雷达灵敏度一定时,星载雷达可测云厚与雷达波长、冰含水量（IWC）的垂直分布、云厚及衰减有关;没有进行衰减订正时,双波长比和衰减有关,冰含水量越大,波长越短,衰减越大,双波长比最大值与可探测云厚有关。两部雷达可探测冰含水量为0.001—0.1 g/m3、厚2 km的冰云;当云厚5 km、冰含水量垂直分布在0.001—0.2 g/m3时,云厚的94%基本可以被220 GHz云雷达探测到。（3）如果两部雷达气象方程中用水的介电因子,测量回波强度应进行介电因子的订正后再计算双波长比。      

A 35-GHz Polarimetric Doppler Radar and Its Application for Observing Clouds Associated with Typhoon Nuri

Millimeter-wavelength radar has proved to be an effective instrument for cloud observation and research. In this study, 8-mm-wavelength cloud radar (MMCR) with Doppler and polarization capabilities was used to investigate cloud dynamics in China for the first time. Its design, system specifications, calibration, and application in measuring clouds associated with typhoon are discussed in this article. The cloud radar measurements of radar reflectivity (Z), Doppler velocity (Vr), velocity spectrum width (Sw) and the depolar-ization ratio (LDR) at vertical incidence were used to analyze the microphysical and dynamic processes of the cloud system and precipitation associated with Typhoon Nuri, which occurred in southern China in August 2008. The results show the reflectivity observed using MMCR to be consistent with the echo height and the melting-layer location data obtained by the nearby China S-band new-generation weather radar (SA), but the Ka-band MMCR provided more detailed structural information about clouds and weak precipitation data than did the SA radar. The variation of radar reflectivity and LDR in vertical structure reveals the transformation of particle phase from ice to water. The vertical velocity and velocity spectrum width of MMCR observations indicate an updraft and strong turbulence in the stratiform cloud layer. MMCR provides a valuable new technology for meteorological research in China.     

Application of Radar Reflectivity Factor in Initializing Cloud-Resolving Mesoscale Model. Part Ⅰ: Retrieval of Microphysical Elements and Vertical Velocity

Assuming that cloud reaches static state in the warm microphysical processes, water vapor mixing ratio(qv), cloud water mixing ratio (qc), and vertical velocity (w) can be calculated from rain water mixing ratio (qr)- Through relation of Z-qr, qr can be retrieved by radar reflectivity factor (Z). Retrieval results indicate that the distributions of mixing ratios of vapor, cloud, rain, and vertical velocity are consistent with radar images, and the three-dimensional spatial structure of the convective cloud is presented. Treating q,v saturated at the echo area, the retrieved qr is about 0.1 g kg-1, qc is always less than 0.3 g kg-1, w is usually below 0.5 m s-1, and rain droplet terminal velocity (vr) is around 5.0 m s-1 in the place where radar reflectivity factor is about 25 dBz; in the place where echo is 45 dBz, the retrieved qr and qc are always about 3.0 g kg-1, w is greater than 5.0 m s-1, and vr is around 7.0 m s-1. In the vertical, the maximum updraft velocity is greater than 3.0 m s-1 at the height of around 5.0 kin, the maximum cloud water content is about 3.0 g kg-1 above 5 km and the maximum rain water content is about 3.0 g kg-1 below 6 kin. Due to the assumption that the cloud is in static state, there will be some errors in the retrieved variables within the clouds which axe rapidly growing or dying-out, and in such cases, more sophisticated radar data control technique will help to improve the retrieval results.     

NUMERICAL STUDY OF REMOTE SENSING OF CLOUD LIQUID WATER CONTENT BY SPACE-BORNE COMBINED RADAR-RADIOMETER

A retrieval method of microwave(MW)space-borne remote sensing of cloud liquid watercontent by the combined radar-radiometer is suggested.A three-layer cloud model is chosen torepresent the typical stratified precipitating cloud.The retrieval method mainly follows ourpreviously suggested scheme with some modifications.Numerical comparative study shows that inspace-borne remote sensing of cloud liquid water content by the combined method is much betterthan by radar only;Also the retrieval accuracy of cloud liquid water content may be improvedwhen using the three-layer cloud model in the combined method.     

NUMERICAL STUDY OF REMOTE SENSING OF CLOUD LIQUID WATER CONTENT BY SPACE-BORNE COMBINED RADAR-RADIOMETER*

A retrieval method of microwave(MW) space-borne remote sensing of cloud liquid water content by the combined radar-radiometer is suggested.A three-layer cloud model is chosen to represent the typical stratified precipitating cloud.The retrieval method mainly follows our previously suggested scheme with some modifications.Numerical comparative study shows that in space-borne remote sensing of cloud liquid water content by the combined method is much better than by radar only;Also the retrieval accuracy of cloud liquid water content may be improved when using the three-layer cloud model in the combined method.     

SACOL站冰云粒子下降末速度的反演及其时空分布特征研究

冰云是影响气候变化最为重要的因子之一,其生命周期的变化在很大程度上决定了冰云的气候辐射效应。冰云粒子下降末速度是影响冰云生命周期的关键参数。为了开展对冰云粒子下降末速度的研究,利用兰州大学半干旱气候与环境监测站Ka波段毫米波云雷达2013年8月至2015年7月连续观测数据,反演了冰云粒子的下降末速度（V_t）,并根据雷达反射率因子（Z）与V_t的关系计算了拟合因子a、b的值;在此基础上应用聚类分析方法,对比分析了4种不同特性冰云Z、V_t和拟合因子a、b的时、空分布特征,进而尝试通过参数垂直分布特征识别研究云中不同位置上云微物理过程的变化。结果表明:冰云粒子下降末速度的分布与雷达反射率因子有很好的对应,最大频率都出现在距离地面约7 km高度处,且具有显著的季节变化,粒子下降末速度在暖季较冷季可增大25%,峰值出现在6月和9月;云层较厚且持续时间长的第一、三类冰云,其雷达反射率因子、粒子下降末速度及拟合因子a和b的平均值都显著大于云层较薄且持续时间短的第二、四类云。垂直方向上,Z、V_t和拟合因子b从云顶到云底随着高度的降低呈现先增大后减小的趋势,体现了云粒子在云顶区域成核和水汽凝华效应,随着粒子在下落过程中碰并增长,云滴粒子逐渐增大,水汽的凝华和粒子的聚合起主要作用,最后在云底部分,云粒子蒸发、升华减小消亡的过程。由此表明中纬度干旱半干旱地区冰云是从云顶到云底自上而下的形成过程。      

基于CloudSat卫星对地基云雷达反射率因子的校准

地基云雷达是云的重要探测手段，但随着运行时间的增加，雷达发射机、接收器等参数的变化，会使观测数据产生漂移偏差，从而对云物理特性的反演产生显著影响，因此云雷达数据的校准是一个重要的基础问题。针对KAZR（Ka Band Zenith Radar，K波段云雷达）云雷达特征，本文在Pavlos等提出的雷达数据校准方法基础上进行改进，优化了对弱云和降水的信号识别，利用CloudSat星载雷达观测的反射率因子，气体衰减校正等数据，对兰州大学半干旱气候环境监测站（Semi Arid Climate and Environment Observatory of Lanzhou University，SACOL）KAZR云雷达2013年8月至2017年5月反射率因子数据进行了校准，建立了KAZR雷达反射率因子46个月的历史资料校准数据库，并对校准周期的变化进行了对比分析。校准数据库的建立对SACOL站云的长期观测研究具有重要意义，同时为不同波段地基雷达的对比增加了可行性。     

Relationship between cloud characteristics and radar reflectivity based on aircraft and cloud radar co-observations

Cloud properties were investigated based on aircraft and cloud radar co-observation conducted at Yitong, Jilin, Northeast China. The aircraft provided in situ measurements of cloud droplet size distribution, while the millimeter-wavelength cloud radar vertically scanned the same cloud that the aircraft penetrated. The reflectivity factor calculated from aircraft measurements was compared in detail with simultaneous radar observations. The results showed that the two reflectivities were comparable in warm clouds, but in ice cloud there were more differences, which were probably associated with the occurrence of liquid water. The acceptable agreement between reflectivities obtained in water cloud confirmed that it is feasible to derive cloud properties by using aircraft data, and hence for cloud radar to remotely sense cloud properties. Based on the dataset collected in warm clouds, the threshold of reflectivity to diagnose drizzle and cloud particles was studied by analyses of the probability distribution function of reflectivity from cloud particles and drizzle drops. The relationship between reflectivity factor (Z) and cloud liquid water content (LWC) was also derived from data on both cloud particles and drizzle. In comparison with cloud droplets, the relationship for drizzle was blurred by many scatter points and thus was less evident. However, these scatters could be partly removed by filtering out the drop size distribution with a large ratio of reflectivity and large extinction coefficient but small effective radius. Empirical relationships of Z-LWC for both cloud particles and drizzle could then be derived.     

雷达反射率因子在中尺度云分辨模式初始化中的应用Ⅰ：云微物理量和垂直速度的反演

对于水平网格距小于10 km的高分辨率非静力平衡的显式云分辨模式,云微物理量的初始化以及物理量之间的相互协调是十分重要的,并且一直是云分析领域的一个难题。考虑在云雨处于定常状态的前提下,根据暖云过程,可得到云中水成物之间相互转化以及与垂直速度的约束关系,而雨水含量跟雷达反射率因子(Z)有关。因此,采用合肥新一代天气雷达2003年7月5日02时(北京时)的观测资料,针对雷达探测的特点对雷达原始数据进行了坐标转换,并进行了插值处理和简单的质量控制,然后依据Z-qr关系和定常暖云方案,反演了雨水混合比(qr)、云水混合比(qc)、水汽混合比(qv)和垂直速度(w)。结果表明,由此得到的雷达回波主要特征及量值分布与雷达站的CAPPI图像基本一致;水汽、云水、雨水和垂直速度量值大小的分布与雷达回波强度的分布也是吻合的,体现出了云中水成物和垂直速度的三维分布结构;各物理量的量值也符合梅雨锋暴雨的特点,梅雨锋积层混合云系中层状云和对流云的差别十分明显。在较强的回波区,雨水在6 km以下含量较大,最大值位于4 km左右,可超过3.0 g/kg;上升速度在5 km左右最大,最大值超过5 m/s;云水含量大值区位于5 km以上,在7 km高度上达到最大,超过3.0 g/kg;雨滴末速度虽然上下比较一致,一般几米每秒,但在5 km左右为大值区。     

多普勒天气雷达资料对中尺度模式短时预报的影响

利用中尺度模式ARPS(The Advanced Regional Prediction System)及其资料分析系统ADAS(ARPS Data Analysis System),将国内新一代多普勒雷达(CINRAD)反射率及径向风资料直接用于中尺度数值模拟,通过一次华北地区暴雨过程的模拟对比试验,分析了雷达资料对初始场的改进效果及其对模拟结果的影响,结果表明:(1)利用雷达径向风资料对初始风场进行调整后,自近地面到对流层顶的u,v,w都发生了变化,调整后的初始风场在对流层中层变化最大.(2)利用雷达反射率进行微物理调整和云分析能调整初始场中的云水信息,使得雷达回波附近3 km以下的水汽混合比(qv)增加,4 km以下的雨水混合比(qr)增加,对流层(约10 km以下)的云水混合比(qc)增加,4～9 km的对流层上部云冰混合比(qi)和雪混合比(qs)增加.ADAS通过非绝热初始化调整温度场,从而得到了一个动力和热力上平衡的初始场.(3)模拟的1 h雨量与实况的对比表明,同时利用雷达反射率和径向风改进过的初始场能明显增强3 h内的降水强度和落区预报,改善中尺度数值模式短时定量降水预报.模拟的1 h流场对比分析表明,经雷达径向风调整后,能够在初始场中增加气旋性涡旋等中小尺度风信息,明显减少模式的spin-up时间.(4)通过对雷达径向风和反射率对模式初始场和模拟结果影响的对比分析发现,雷达径向风主要是改进初始风场,而雷达反射率主要是改进初始场中的湿度参数,增加初始场中云水等的含量,调整温度场.通过模拟的6 h降水对比发现,利用雷达径向风调整初始场后,对降水模拟有一定的改进,但效果不甚明显,而雷达反射率资料对定量降水预报改进效果明显,同时使用雷达径向风和反射率资料改进初始场后对降水的模拟效果最明显.     

The First Observed Cloud Echoes and Microphysical Parameter Retrievals by China's 94-GHz Cloud Radar

By using the cloud echoes first successfully observed by China's indigenous 94-GHz SKY cloud radar, the macrostructure and microphysical properties of drizzling stratocumulus clouds in Anhui Province on 8 June2013 are analyzed, and the detection capability of this cloud radar is discussed. The results are as follows.(1) The cloud radar is able to observe the time-varying macroscopic and microphysical parameters of clouds,and it can reveal the microscopic structure and small-scale changes of clouds.(2) The velocity spectral width of cloud droplets is small, but the spectral width of the cloud containing both cloud droplets and drizzle is large. When the spectral width is more than 0.4 m s-1, the radar reflectivity factor is larger(over –10 dBZ).(3) The radar's sensitivity is comparatively higher because the minimum radar reflectivity factor is about–35 dBZ in this experiment, which exceeds the threshold for detecting the linear depolarized ratio(LDR) of stratocumulus(commonly –11 to –14 dBZ; decreases with increasing turbulence).(4) After distinguishing of cloud droplets from drizzle, cloud liquid water content and particle effective radius are retrieved. The liquid water content of drizzle is lower than that of cloud droplets at the same radar reflectivity factor.     

A Ka-band Solid-state Transmitter Cloud Radar and Data Merging Algorithm for Its Measurements

This study concerns a Ka-band solid-state transmitter cloud radar, made in China, which can operate in three different work modes, with different pulse widths, and coherent and incoherent integration numbers, to meet the requirements for cloud remote sensing over the Tibetan Plateau. Specifically, the design of the three operational modes of the radar(i.e., boundary mode M1, cirrus mode M2, and precipitation mode M3) is introduced. Also, a cloud radar data merging algorithm for the three modes is proposed. Using one month's continuous measurements during summertime at Naqu on the Tibetan Plateau,we analyzed the consistency between the cloud radar measurements of the three modes. The number of occurrences of radar detections of hydrometeors and the percentage contributions of the different modes' data to the merged data were estimated.The performance of the merging algorithm was evaluated. The results indicated that the minimum detectable reflectivity for each mode was consistent with theoretical results. Merged data provided measurements with a minimum reflectivity of -35 dBZ at the height of 5 km, and obtained information above the height of 0.2 km. Measurements of radial velocity by the three operational modes agreed very well, and systematic errors in measurements of reflectivity were less than 2 dB. However,large discrepancies existed in the measurements of the linear depolarization ratio taken from the different operational modes.The percentage of radar detections of hydrometeors in mid- and high-level clouds increased by 60% through application of pulse compression techniques. In conclusion, the merged data are appropriate for cloud and precipitation studies over the Tibetan Plateau.     

闪电定位和雷达观测资料在云分析中的应用及数值试验

本文针对2016年6月山西、山东地区一次强对流活动,通过GRAPES（Global/Regional Assimilation andPrediction Enhanced System）云分析系统融合了多普勒雷达反射率三维组网拼图资料和ADTD（Active DirectoryTopology Diagrammer）闪电定位资料,对模式初始场的云水、云冰等云微物理变量进行了调整和分析,并对这次强对流过程设计了3组数值试验。结果表明:（1）云分析系统在加入雷达反射率资料以后,能在模式初始场中较准确的计算模式的初始水物质;（2）闪电定位资料转换成的替代雷达反射率对原有的雷达反射率有一定的补充作用,使模式的初始场与真实状况更接近;（3）从模式降水的角度看,在加入雷达资料以后,GRAPES云分析可以显著提高6小时以内降水的模拟效果,而闪电定位资料的加入,对中雨以上量级的降水模拟有进一步的改进作用。     

Analysis and very short range forecast of cyclone “AILA” with radar data assimilation with rapid intermittent cycle using ARPS 3DVAR and cloud analysis techniques

In this study, both reflectivity and radial velocity are assimilated into the Weather Research and Forecasting (WRF) model using ARPS 3DVAR technique and cloud analysis procedure for analysis and very short range forecast of cyclone ÁILA. Doppler weather radar (DWR) data from Kolkata radar are assimilated for numerical simulation of landfalling tropical cyclone. Results show that the structure of cyclone AILA has significantly improved when radar data is assimilated. Radar reflectivity data assimilation has strong influence on hydrometeor structures of the initial vortex and precipitation pattern and relatively less influence is observed on the wind fields. Divergence/convergence conditions over cyclone inner-core area in the low-to-middle troposphere (600–900 hPa) are significantly improved when wind data are assimilated. However, less impact is observed on the moisture field. Analysed minimum sea level pressure (SLP) is improved significantly when both reflectivity and wind data assimilated simultaneously (RAD-ZVr experiment), using ARPS 3DVAR technique. In this experiment, the centre of cyclone is relocated very close to the observed position and the system maintains its intensity for longer duration. As compared to other experiments track errors are much reduced and predicted track is very much closer to the best track in RAD-ZVr experiment. Rainfall pattern and amount of rainfall are better captured in this experiment. The study also reveals that cyclone structure, intensification, direction of movement, speed and location of cyclone are significantly improved and different stages of system are best captured when both radar reflectivity and wind data are assimilated using ARPS 3DVAR technique and cloud analysis procedure. Thus optimal impact of radar data is realized in RAD-ZVr experiment. The impact of DWR data reduces after 12 h forecast and it is due to the dominance of the flow from large-scale global forecast system model. Successful coupling of data assimilation package ARPS 3DVAR with WRF model for Indian DWR data is also demonstrated.     

毫米波云雷达功率谱密度数据的检验和 在弱降水滴谱反演中的应用研究

本文首先利用数值模拟的方法,分析了利用毫米波云雷达功率谱密度反演雨滴谱时,降水粒子米散射效应、空气湍流、空气上升速度等对雨滴谱和液态水含量等参数反演的影响;建立了功率谱密度处理及其直接反演雨滴谱、液态水含量、降水强度和空气上升速度的方法;并利用2012年7月在云南腾冲观测的二次弱降水数据,采用毫米波雷达和Ku波段微降水雷达观测的回波强度、径向速度垂直廓线以及780 m高度上的功率谱密度对比的方法,以及毫米波云雷达观测的780 m高度上功率谱密度、回波强度与地面雨滴谱计算得到的这些量的对比方法,分析了毫米波雷达数据的可靠性;并将780 m高度上毫米波雷达反演的雨滴谱与地面雨滴谱数据进行了对比,分析了毫米波雷达反演的雨滴谱的准确性;分析了毫米波雷达回波强度偏弱的原因,讨论了该高度以下降水对毫米波雷达衰减的影响。结果表明:空气湍流对弱降水微物理参数反演影响不大,而空气上升速度和米散射效应均对反演结果有一定影响;毫米波雷达观测到的径向速度和功率谱密度与微降水雷达比较一致,回波强度的垂直廓线的形状与微降水雷达也比较一致,但毫米波雷达观测的回波强度偏弱;与雨滴谱计算值相比,毫米波雷达观测的低层的回波强度也偏弱,天线上的积水是造成毫米波雷达回波强度变弱的主要原因。毫米波雷达观测的低层的功率谱密度与地面雨滴谱观测的数据形状比较一致,但有一定的位移。毫米波雷达反演的雨滴谱与地面观测的谱型和粒子大小也比较一致。这些结果初步验证了毫米波雷达观测的功率谱密度及其反演方法的可靠性。     

太赫兹频段雷达探测的冰云微物理参数反演算法模拟研究

在假设冰云粒子呈球形以及粒子谱服从对数正态分布的条件下,利用离散偶极子近似法(DDA),计算出太赫兹频段(220 GHz)冰云粒子的雷达反射率因子,及其与瑞利假设下雷达反射率因子的比值。忽略衰减和多次散射的影响,根据太赫兹波段冰云的雷达反射率因子,基于最优估计理论反演冰云的微物理参数,并验证该算法的可靠性。反演结果表明,当冰云粒子大小在设定的尺度范围内时,有效粒子半径(r_e)的反演误差小于4%,粒子谱宽(σ)的误差小于2.5%、粒子数密度(N_T)的误差小于1%,冰水含量(IWC)的误差小于5%。还分析了当N_T和σ为定值时,反演结果随粒子尺寸的变化情况,当冰云粒子尺寸在模拟计算设定的范围内时,r_e的反演误差小于0.04%,σ的反演误差小于0.02%,N_T的反演误差小于0.50%,IWC的反演误差小于0.08%,如果冰云粒子大小超出模拟计算设置的范围,反演误差随着r_e增加而增大。该结果证明了基于最优估计理论反演得到的冰云微物理参数与模拟设定值有良好的一致性,说明该方法可应用于太赫兹频段云雷达的冰云观测及云微物理参数的反演和研究。      

WRF模式云参数化方案对一次深对流系统模拟的验证和改进

使用中尺度数值模式WRF中的双参数云微物理方案WDM6针对2008年台风“凤凰”登陆过程中造成的强降水进行数值模拟,通过卫星模拟器利用MTSAT-1R和TRMM卫星观测的红外云顶黑体亮温TBB、PR雷达反射率资料使用统计方法验证模拟结果。通过修改云水向雨水自动转化过程、冰晶核化过程、雪和霰的下落末速度、雪和霰的截距进行敏感性试验,减小模拟结果和卫星观测结果的差异。研究结果表明:WDM6方案模拟的台风“凤凰”登陆后的降水,强对流云系及对流柱状雷达回波基本符合实况,但模拟结果局部偏强。WDM6方案模拟产生了较多的浅对流云,低估了对流云系的出现频率。不同云类型模拟的雷达回波均偏强,对流云系雷达回波垂直分布接近观测。敏感性试验结果说明修改WDM6方案中云水向雨水自动转化率有效地改善了模拟效果。同时发现云滴初始数浓度影响云水向雨水自动转化率并最终影响云系结构和雷达反射率的模拟结果,过高的云滴初始数浓度会使模拟结果变差。     

双线偏振雷达资料在数值模式中的应用: 模拟器的构建

基于Rayleigh-Gans散射原理构造一个S波段双线偏振雷达模拟器。模拟器考虑了云冰、雪、雨和雹4种水凝物,以水凝物的混合比和数浓度以及水凝物粒子的轴比、相对介电常数、下落倾角为输入量,计算得到水平/垂直偏振反射率因子、差分反射率、比差分相位等偏振量。通过二维理想飑线模拟的试验结果表明,模拟器合理地再现了二维理想飑线系统成熟期的主要偏振特征:雹的反射率较高,差分反射率较低(又称Zdr洞);对流云区的Zdr柱;层云区的反射率和比差分相位的0℃层亮带特征;雨滴反射率与差分反射率因子的“雨线”统计特征。利用该模拟器建立了模式变量和偏振雷达观测的联系,有助于未来将偏振雷达观测资料应用于模式预报效果评估及对流尺度资料同化等方面的研究和应用。