降雨强度对激光雷达测风精度的影响研究

多普勒激光雷达在大气、环境以及风能等领域得到越来越广泛的应用,但对于复杂天气下多普勒激光雷达的适用性仍然有待深入研究。为此本研究采用2020年8—10月期间位于福建三沙的地基多普勒激光雷达与边界层高塔所搭载的超声风温仪观测的风场数据进行对比,发现多普勒激光雷达在水平风速、风向方面具有稳定的高精度探测性能,与超声风温仪之间相关系数达到了0.948和0.984。相比之下,激光雷达垂直风速的探测误差较大,与超声风温仪之间相关系数仅有0.353。研究发现,降雨强度与激光雷达垂直风速误差呈正相关关系,强降雨下垂直风速偏差最大可达到9 m/s。     

测风激光雷达与风廓线雷达的探测性能评估及数据融合

测风激光雷达和风廓线雷达作为L波段探空测风的有效补充,均可以提供高时空分辨率的大气风场信息,然而由于工作原理和适用条件存在明显差异,在探测性能上各有优缺点,单一设备的探测数据已不能满足精细化预报的要求。本研究使用2020年1—5月北京南郊观象台的L波段探空资料对同址观测的测风激光雷达和风廓线雷达进行了数据质量评估,结果表明测风激光雷达与探空的一致性较高,U、V分量的相关系数分别为0.97和0.98,均方根误差分别为1.1和0.95 m·s-1,然而在2 km以上数据获取率较低且偏差较大;风廓线雷达与探空相比,U、V分量的相关系数分别为0.94和0.93,均方根误差分别为2.94和2.91 m·s-1,风廓线雷达的探测距离虽然更远,但在0.5 km以下和6 km以上的测量偏差较大。考虑到两种测风雷达在不同探测高度上的性能优缺点,提出分段曲面拟合法对两者的水平风资料进行融合处理,并选取个例对融合效果进行验证,结果表明,融合后的风廓线与融合前相比,风向和风速的一致性均得到明显提升。     

单多普勒雷达二维风场反演——Extended VPP方法

文中提出的EVPP(Extended Velocity Plan Processing)方法是根据风场连续及局部分片光滑的性质,假定在一很小的分析窗内风矢量保持线性变化关系,使得在分析窗内反演得到的风场在矢径上的投影与实际观测的多普勒雷达径向速度场之间的误差达到最小,并利用零速度线这一特点对问题的求解做出约束,使得所得解尽可能达到全局最优,即:(1)利用变分的思想求出水平风速的大小及方向;(2)利用风场局部分片光滑,对风的方向做出矫正。虽然由于风场受到线性假设这一限制,使得其使用的范围得到了局限,但是在本文中利用零速度线这一特征使得反演的结果得到了进一步的提高。文中利用模拟的理想个例和实际观测到的个例分别进行了研究分析。为了证明其理论的可行性,首先我们利用三维冰雹云模式模拟了辽宁2005年7月8日的一次冰雹云过程,考虑到个例的代表性,我们选择了辐合、辐散、旋转以及均匀风场等4种类型的流场结构进行了模拟分析,同时为了验证该种方法对于资料质量的依赖性以及噪声的敏感性,我们在理想质量中加入了随机噪声。从对模拟风场的反演,证明了该方法理论上的可行性,而且噪声对其结构的影响不大。然后针对台风、飑线、强对流等实际天气过程作了进一步分析。最后对反演结果进行了误差分析,结果表明这一方法能较好地反演出中尺度系统的水平流场结构。     

单部多普勒天气雷达风场反演及其在数值预报中的应用试验

针对1998年7月5日下午到夜间北京地区出现的一次强降雨过程, 采用“局部VAD法”, 对中国气象科学研究院多普勒雷达测得的1998年7月5日12:00UTC的径向风反演成水平风, 并将反演风应用到数值天气预报 (NWP) 同化系统中, 进行了同化预报试验。结果表明, 反演的雷达风基本符合事实, 将其应用到NWP中, 有助于分析和预报时空尺度较小的中尺度系统。     

多普勒激光雷达风场反演方法研究

采用三维变分同化反演(3DVAR)、 四维变分同化反演(4DVAR) 对多普勒激光雷达资料反演风场的方法进行了研究, 利用车载多普勒激光雷达在2008年残奥会测试赛期间外场试验取得的数据, 反演了海面10 m高度处的风场, 并将风场反演结果与浮标资料进行了对比分析, 结果表明: 3DVAR、4DVAR风场反演方法均能实现近海面风的精细化风场反演, 并能反映出风向的变化, 反演风场与浮标数据基本一致, 在风速较大的天气情况, 3DVAR与4DVAR反演风场的一致性要好于风速较小的天气情况; 4DVAR反演方法中以浮标资料作为背景场, 使得其与浮标的符合程度要好于3DVAR方法反演风场; 反演风场的风向与浮标风向具有很好的相关关系, 反演风场的风速与浮标风速具有一定的相关关系, 反演风场的风向、风速与浮标的风向、风速之间平均均方根误差和平均绝对误差表明, 这两序列之间具有一定差别, 在风速较小的天气情况下使用时需要注意。     

Effects of a change of terrain height and roughness on a wind profile

We have measured profiles of an onshore wind 200 m upwind and 200 m downwind from an abruptly rising shore using a remote-sensing Doppler lidar anemometer. Data were taken at heights between 4.7 and 66.5 m above sea level. Results show that the onshore velocity vector slopes upward 16.6 to 9.6°, the amount depending on the height of measurement, due to the combined effects of a 1.7-m high bluff shoreline and the frictionally decelerated flow over land. The profile 200 m inland has the expected deceleration at lower levels because of increased surface roughness and implies a velocity vector at 66.5 m height with an upward slope of approximately 18° (2.6 m s^-1 upward component, 8.4 m s^-1 vector magnitude), an acceleration to 0.3 m s^-1 greater than the upwind value, or a combination of both effects. All three options are consistent with mass continuity. The experiment exhibits the usefulness and limitations of a backscatter Doppler lidar for boundary-layer profile measurements in a horizontally inhomogeneous environment.     

从单Doppler速度场反演风矢量场的VAP方法

从单Doppler速度场反演风矢量场的VAP(Velocity Azimuth Processing)方法是在假定相邻方位角的风矢量相等的前提下,从单部Doppler雷达所测得的各个距离圈上的径向速度随方位角分布的廓线推算风向和风速的方法。这个假定可以保证反演所得的风矢量场中仍保留次天气尺度的风场信息。 对VAP方法反演公式的误差分析表明,反演出的风向误差为20—30度,与常规测风的误差相近。风速误差一般为10~0m/s,但在某些特定的情况下可达到与风速相同的量级(10~1m/s)。 用VAP方法反演风矢量场前必须首先消除原始Doppler速度资料中的脉动。它可以通过简单的一维线性平滑实现。 文中给出了一个反演实例,从中可以分析出次天气尺度的切变线和辐合线。常规天气图和UHF雷达的风廓线观测表明,这些次天气尺度系统是存在的。它们与反射率因子分布的配置还表明这些次天气尺度系统是和降水相联系的。     

The Application of Sea-Surface Wind Detection with Doppler Lidar in Olympic Sailing

The mobile incoherent Doppler lidar(MIDL),which was jointly developed by State Key Laboratory of Severe Weather(LaSW) of the Chinese Academy of Meteorological Sciences(CAMS) and Ocean University of China,provided meteorological services during the Olympic sailing events in Qingdao in 2008.In this study,two experiments were performed based on these measurements.First,the capabilities of MIDL detection of sea-surface winds were investigated by comparing its radial velocities with those from a sea buoy.MIDL radial velocity was almost consistent with sea-buoy data;both reflected the changes in wind with time.However,the MIDL data was 0.5 m s 1 lower on average than the sea-buoy data due to differences in detection principle,sample volume,sample interval,spatial and temporal resolution.Second,the wind fields during the Olympic sailing events were calculated using a four-dimensional variation data assimilation(4DVAR) algorithm and were evaluated by comparing them with data from a sea buoy.The results show that the calculations made with the 4DVAR wind retrieval method are able to simulate the fine retrieval of sea-surface wind data-the retrieved wind fields were consistent with those of sea-buoy data.Overall,the correlation coefficient of wind direction was 0.93,and the correlation coefficient of wind speed was 0.70.The distribution of retrieval wind fields was consistent with that of MIDL radial velocity;the root-mean-square error between them had an average of only 1.52 m s 1.     

激光测风雷达PPI扫描监测风切变的数值模拟

低空风切变严重威胁飞机起降安全,而激光测风雷达在晴空风切变监测方面具有显著优势,但其监测性能受到风场特点和监测算法的影响。由于无法对局地空间风场进行时空同步精细化探测,且风场具有阵性,导致风切变监测会出现虚警率、漏报等现象,只能通过数值模拟仿真的方式,对不同风场的监测性能进行分析评估。本文首先仿真构建了四种具有代表性的风场类型,然后模拟激光测风雷达平面—位置显示（Plane Position Indication,PPI）扫描方式获取不同仿真风场下的径向风速数据,最后利用八邻域系统算法对径向风速数据进行风切变监测,验证不同类型风场下激光测风雷达PPI扫描方式监测风切变算法的有效性。结果表明,激光测风雷达PPI扫描结果能够真实地反映风场情况;八邻域系统算法能够有效识别风切变;在进行风切变检测时,应尽可能地考虑到多个方位上的切变值,避免漏报。     

Retrieval of atmospheric and oceanic parameters and the relevant numerical calculation

It is well known that retrieval of parameters is usually ill-posed and highly nonlinear, so parameter retrieval problems are very difficult. There are still many important theoretical issues under research, although great success has been achieved in data assimilation in meteorology and oceanography. This paper reviews the recent research on parameter retrieval, especially that of the authors. First, some concepts and issues of parameter retrieval are introduced and the state-of-the-art parameter retrieval technology in meteorology and oceanography is reviewed briefly, and then atmospheric and oceanic parameters are retrieved using the variational data assimilation method combined with the regularization techniques in four examples： retrieval of the vertical eddy diffusion coefficient; of the turbulivity of the atmospheric boundary layer; of wind from Doppler radar data, and of the physical process parameters. Model parameter retrieval with global and local observations is also introduced.     

Vertical Structure of the Urban Boundary Layer over Marseille Under Sea-Breeze Conditions

During the UBL-ESCOMPTE program (June–July 2001), intensive observations were performed in Marseille (France). In particular, a Doppler lidar, located in the north of the city, provided radial velocity measurements on a 6-km radius area in the lowest 3 km of the troposphere. Thus, it is well adapted to document the vertical structure of the atmosphere above complex terrain, notably in Marseille, which is bordered by the Mediterranean sea and framed by numerous massifs. The present study focuses on the last day of the intensive observation period 2 (26 June 2001), which is characterized by a weak synoptic pressure gradient favouring the development of thermal circulations. Under such conditions, a complex stratification of the atmosphere is observed. Three-dimensional numerical simulations, with the Méso-NH atmospheric model including the town energy balance (TEB) urban parameterization, are conducted over south-eastern France. A complete evaluation of the model outputs was already performed at both regional and city scales. Here, the 250-m resolution outputs describing the vertical structure of the atmosphere above the Marseille area are compared to the Doppler lidar data, for which the spatial resolution is comparable. This joint analysis underscores the consistency between the atmospheric boundary layer (ABL) observed by the Doppler lidar and that modelled by Méso-NH. The observations and simulations reveal the presence of a shallow sea breeze (SSB) superimposed on a deep sea breeze (DSB) above Marseille during daytime. Because of the step-like shape of the Marseille coastline, the SSB is organized in two branches of different directions, which converge above the city centre. The analysis of the 250-m wind fields shows evidence of the role of the local topography on the local dynamics. Indeed, the topography tends to reinforce the SSB while it weakens the DSB. The ABL is directly affected by the different sea-breeze circulations, while the urban effects appear to be negligible.     

多普勒雷达风场反演VAP方法的资料预处理

在分析了三类单多普勒雷达反演风场方法的基础上, 认为VAP方法是近期内可望在业务中应用的方法, 并指出VAP方法的关键是能否将原始资料中的中-U及以上尺度的风场信息过滤出来, 以消除中-V尺度以下风场信息对反演结果的影响.为此提出了二维区域滑动低通滤波方法对资料进行预处理.此外还对VAP方法的风速计算公式进行了改进以减小计算误差.     

An Optimal Inverse Method Using Doppler Lidar Measurements to Estimate the Surface Sensible Heat Flux

Inverse methods are widely used in various fields of atmospheric science. However, such methods are not commonly used within the boundary-layer community, where robust observations of surface fluxes are a particular concern. We present a new technique for deriving surface sensible heat fluxes from boundary-layer turbulence observations using an inverse method. Doppler lidar observations of vertical velocity variance are combined with two well-known mixed-layer scaling forward models for a convective boundary layer (CBL). The inverse method is validated using large-eddy simulations of a CBL with increasing wind speed. The majority of the estimated heat fluxes agree within error with the proscribed heat flux, across all wind speeds tested. The method is then applied to Doppler lidar data from the Chilbolton Observatory, UK. Heat fluxes are compared with those from a mast-mounted sonic anemometer. Errors in estimated heat fluxes are on average 18 %, an improvement on previous techniques. However, a significant negative bias is observed (on average $-63\,\%$ ) that is more pronounced in the morning. Results are improved for the fully-developed CBL later in the day, which suggests that the bias is largely related to the choice of forward model, which is kept deliberately simple for this study. Overall, the inverse method provided reasonable flux estimates for the simple case of a CBL. Results shown here demonstrate that this method has promise in utilizing ground-based remote sensing to derive surface fluxes. Extension of the method is relatively straight-forward, and could include more complex forward models, or other measurements.     

基于多普勒测风激光雷达的低空急流观测研究

基于相干多普勒测风激光雷达于2018年8月在山东德州获取的为期一个月的风廓线观测数据,进行了低空急流的判定、识别与统计分析。参考BONNER对低空急流的判定标准,对1 500 m高度以下的每10 min平均风廓线数据进行低空急流识别与统计,急流发生频率仅为3.6%。参考张世丰对低空急流的判定标准,统计了350 m高度以下10 min平均风廓线的低空急流风速、高度、风向及风切变等结构特征。急流发生频率为24.9%,急流速度主要介于6～10 m·s-1之间,急流高度出现3个峰值,分别位于110 m、160 m和220 m左右,急流风向主要为偏东风和偏南风。结果表明,多普勒激光雷达可以获取高时空分辨率的风廓线数据,进而可以有效检测低空急流结构的存在及其特征。     

Remote Sensing Of Three-Dimensional Winds With Elastic Lidar: Explanation Of Maximum Cross-Correlation Method

Maximum cross-correlation provides a method to remotely determine highly resolved three-dimensional fields of horizontal winds with elastic lidar throughout large volumes of the planetary boundary layer (PBL). This paper details the technique and shows comparisons between elastic lidar winds, remotely sensed laser Doppler velocimeter (LDV) wind profiles, and radiosonde winds. Radiosonde wind data were acquired at Barcelona, Spain, during the Barcelona Air-Quality Initiative (1992), and the LDV wind data were acquired at Sunland Park, New Mexico (N.M.), during the Border Area Air-Quality Study (1994). Comparisons show good agreement between the different instruments, and demonstrate the method useful for air pollution management at the local/regional scale. Elastic lidar winds could thus offer insight into aerosol and pollution transport within the PBL. Lidar wind fields might also be used to nudge or improve initialization and evaluation of atmospheric meteorological models.     

Verification of a Synthesized Method for the Calculation of Low-Level Climatological Wind Fields Using a Mesoscale Boundary-Layer Model

Low-level climatological wind fields over the La Plata River region of South America are synthesized with a dry, hydrostatic mesoscale boundary-layer numerical model. The model is forced at the upper boundary with the 1200 UTC local radiosonde observations and at the lower boundary with a land-river differential heating function defined from the daily meteorological observations of the region. The climatological wind field is defined as the mean value of a series of individual daily forecasts, employing two methods. The simplified method considers a 192-member ensemble (16 wind directions and 12 wind-speed classes at the upper boundary). Each member has a probability of occurrence that is determined from the 1959–1984 observations; the daily method uses a total of 3,248 days with available data during the same period. In both methods each realization is a daily forecast from which the mean wind distributions at 0300, 0900, 1500 and 2100 local standard time are calculated and compared to the observations of five meteorological stations in the region. The validation of the climatological wind fields for both methods is evaluated by means of the root-mean-square error of the wind-direction frequency distribution and mean wind speed by wind sector. The results obtained with the two methods are similar, and the errors in wind speed are always smaller than those in wind direction. The combined errors of wind direction and wind speed show that the ensemble method is outperformed by the daily method, on average by meteorological station in only one out of five of them, and on average by the time of the day in only one out of 4 h. The conclusion of the study is that the ensemble method is an appropriate methodology for determining high resolution, low-level climatological wind fields, with the boundary-layer model applied to a region with a strong diurnal cycle of surface thermal contrast. The proposed methodology is of particular utility for synthesizing wind fields over regions with limited meteorological observations, since the 192-member matrix can be easily defined with few observing points, as well as in the case of relatively incomplete records.     

新疆地区一次对流性降水的三维中尺度风场研究

利用2004年外场试验获得的乌鲁木齐和五家渠C波段双多普勒雷达资料, 分析了双多普勒雷达风场反演方法和资料的可靠性, 研究了2004年8月8日发生在乌鲁木齐和五家渠的一次强对流性降水的回波和风场中尺度结构及演变过程。结果表明:这两部雷达观测的回波强度相关很好, 雷达基线上的径向速度基本一致, 资料可靠, 适合进行双多普勒雷达观测; Cressman插值的影响半径的变化对风场的中尺度结构基本没有影响, 径向速度误差引起的风场反演误差与该点所处的位置有关, 1 m/s径向速度误差也不会改变风场的中尺度结构。该过程为对流单体发展为对流带状回波的过程, 在对流单体的左侧生成新的对流单体, 逐步发展为长度约90 km范围的带状对流系统, 该系统恰与较强的东北风和较弱的西风形成的辐合相对应, 上升气流与强对流回波相对应, 不同对流单体有各自相独立的风场结构。用双多普勒雷达观测得到对流系统的内部风场有利于了解对流系统的内部动力过程, 从而探讨降水的形成和演变机理。     

机载多普勒激光测风雷达风场反演研究

中国科学院上海光学精密机械研究所研制出机载全固态多普勒激光测风雷达,并成功进行了试验飞行。本文介绍了机载多普勒激光测风雷达的基本参数和探测三维风场的工作原理,并对利用激光回波信号和机载辅助数据反演三维风场的详细算法进行了研究,获得了3 000 m以下高度具有较高精度的三维风场,同时用经纬仪气球测风法进行了对比验证实验,验证了机载多普勒激光测风雷达测量数据的准确性。     

变分方法反演双多普勒雷达低层二维风场

提出用变分方法和正则化方法相结合反演双多普勒雷达资料的低层二维风场的方法,即在目标泛函中加入正则化项（先验信息项）且合理选择正则化参数。数值试验结果表明：正则化项的引入能有效滤除反演风场中因奇点和观测资料引起的误差,提高反演精度,得到最优的全局反演风场。     

Observations of the Relation Between Upslope Flows and the Convective Boundary Layer in Steep Terrain

Slope flow mechanisms are crucial for the transport of air pollutants in complex terrain. Previous observations in sloping terrain showed upslope flows filling the entire convective boundary layer (CBL) and reducing air pollution concentrations by venting air pollutants out of the CBL into the free atmosphere. During the Pacific 2001 Air Quality Field Study in the Lower Fraser Valley, British Columbia, Canada, we observed slope flows during weak synoptic winds, clear skies, and strong daytime solar heating. With a Doppler sodar we measured the three wind components at the foot of a slope having an average angle of 19° and a ridge height of 780 m. We operated a scanning lidar system and a tethersonde at a nearby site on the adjacent plain to measure backscatter of particulate matter, temperature, wind speed, wind direction, and specific humidity. Strong daytime upslope flows of up to 6 m s⁻¹ through a depth of up to 500 m occurred in the lower CBL, but with often equally strong and deep return flows in the upper part of the CBL. The mass transport of upslope flow and return flow approximately balanced over a 4-h morning period, suggesting a closed slope-flow circulation within the CBL. These observations showed that air pollutants can remain trapped within a CBL rather than being vented from the CBL into the free atmosphere.