On Sea Surface Roughness Parameterization and Its Effect on Tropical Cyclone Structure and Intensity

A new parameterization scheme of sea surface momentum roughness length for all wind regimes, including high winds, under tropical cyclone (TC) conditions is constructed based on measurements from Global Positioning System (GPS) dropsonde. It reproduces the observed regime transition, namely, an increase of the drag coefficient with an increase in wind speed up to 40 m s-1 , followed by a decrease with a further increase in wind speed. The effect of this parameterization on the structure and intensity of TCs...     

不同下垫面空气动力学参数的研究

文中利用中国科学院沙漠研究所与日本国家农业环境技术研究所合作于1990—1994年在中国内蒙古自治区奈曼市半干旱地区沙丘和植被区下垫面观测的微气象数据,根据Monin-Obukhov相似性理论,计算了重度干扰草原、中度干扰草原、轻度干扰草原、无干扰草原、沙丘、沙丘内地、草地、稻田、小麦田、大豆田和玉米田11种下垫面的空气动力学参数粗糙度长度z0,零平面位移d,摩擦速度u*,并分析了它们与水平风速u和Richardson数的关系,比较了不同人为干扰草原生态系统条件下的空气动力学特征。结果表明:地表生物量和覆盖率随着人为干扰强度的增加而减少。不同人为干扰下垫面的粗糙长度与生物量和植被高度以及地表起伏程度有着密切关系;Richardson数也是其影响因子。风速、粗糙度都与摩擦速度成正相关,但对于不同下垫面有所不同,从中可以看到草地对沙漠化有一定的防治作用。同一种下垫面不同时期的空气动力学参数也存在差异。这些结果对建立陆面过程和区域气候模式具有重要的意义。     

上海徐家汇地区台风“灿鸿”与良态风风速特征

上海徐家汇地区建筑分布密集且高低不一,是典型的具有非均一下垫面特征的城市地貌。本文利用该地区的地理信息研究了粗糙度长度的分布规律,并基于80m高度的风速实测数据,对台风"灿鸿"和良态风作用下的平均风速、湍流强度、阵风因子等参数与粗糙度长度之间的关系进行了分析。结果表明:不同风向角对应的计算扇区内建筑物高度、分布密度的差别导致了粗糙度长度值随风向角发生明显变化,但是变化幅度随着计算扇区的增大而减小;台风"灿鸿"作用下的平均风速最大值大于良态风,两者对应的粗糙度长度变化范围差别甚微,但是台风作用下的粗糙度长度中位数较小且分布相对集中;台风"灿鸿"作用下各向湍流强度均随着平均风速的增加呈明显的减小趋势,但不随粗糙度长度变化;良态风作用下,各向湍流强度不随平均风速变化,而随着粗糙度长度的增加而增加;台风"灿鸿"和良态风作用下,各向阵风因子均随湍流强度的增加而增大,但前者作用下的阵风因子略大于后者。     

Effects of different terrain on velocity standard deviations

The standard deviations of wind velocity components are calculated and compared based on the measurements of turbulence fluctuations over three underlying surfaces: uniform, inhomogeneous and urban. Statistical analysis shows that there are the differences between them which prove the influences of surfaces roughness on turbulence.     

Study of Aerodynamic Parameters on Different Underling Surfaces

Aerodynamic parameters including the zero-plane displacement (d), roughness length (20), and friction velocity (u*) on the different underlying surfaces of heavy-grazing site, medium-grazing site, light-grazing site, no-grazing site, dune, inter-dune, grassland, rice paddy site, wheat site, soybean site, and maize site have been computed based on the Monin-Obukhov similarity theory by utilizing the micrometeorologically observed data of dune and vegetation in the semi-arid area at Naiman, Inner Mongolia of China, conducted jointly by the Institute of Desert Research, Chinese Academy of Sciences and the National Institute of Agro-Environmental Sciences of Japan in 1990-1994. And their relationships between wind speed and Richardson number are analyzed. The aerodynamic characteristics of different man-made disturbed grassland ecosystems are also compared. Result shows that the vegetation coverage and the above-ground biomass decrease with the increase in man-made stress of the grassland. The roughness length for different underlying surfaces is closely related to vegetation height, above-ground biomass, and ground surface undulation, and Richardson number Ri is also its influencing factor. The friction velocity varies largely on different underlying surfaces, and it is positively proportional to wind speed and roughness length. The aerodynamic parameters of various times on the same underlying surface are different, too. Above results indicate that grassland and vegetation are of significance in preventing desertification, especially in the arid and semi-arid land ecosystems. And the results of this paper are also important for constructing the land surface physical process as well as regional climate model.     

基于风雨影响评估的深圳市热带气旋数据库系统

风雨影响评估是热带气旋预报预警服务中的一项重要工作,而热带气旋风雨影响的历史气候数据是支撑这一工作的重要依据。以深圳为例,提出了建立热带气旋数据库系统的技术思路,该系统实现了以下几种功能:1)热带气旋数据的后台自动入库;2)关键区和相似离度相结合的相似台风的自动检索和风雨影响分析;3)直观展示广东沿海各市登陆的热带气旋的风雨影响气候概况;4)热带气旋风雨影响等级的自动提示。通过该系统,预报员可以高效、快速地查询历史相似路径或相同登陆点的热带气旋风雨影响情况,从而为分析当前热带气旋可能造成的风雨影响评估提供参考依据。     

Characteristics of Near Surface Winds over Different Underlying Surfaces in China: Implications for Wind Power Development

Accurate wind and turbulence information are essential to wind energy research and utilization, among which wind shear and turbulence intensity/scale have seldom been investigated. In this paper, the observational data from the100-m high wind towers in Xilinhot in Inner Mongolia(2009–10; grassland region), Huanghua in Hebei Province(2009–10; coastal flat region), and Xingzi County in Jiangxi Province(2010–11; mountain–lake region) are used to study the variations in near surface winds and turbulence characteristics related to the development of local wind energy over different underlying surfaces. The results indicate that(1) the percentage of the observed wind shear exponents exceeding 0.3 for the grassland region is 6%, while the percentage is 13% for the coastal flat region and 10%for the mountain–lake region. In other words, if the wind speed at 10 m is 10 ms^–1, the percentage of the wind speed at 100 m exceeding 20 ms^–1 for the grassland region is 6%, while the percentage is 13% for the coastal flat region and 10% for the mountain–lake region.(2) In terms of the turbulent intensity in the zonal, meridional, and vertical directions(I_u, I_v, and I_w, respectively), the frequencies of I_v/I_u < 0.8 in the grassland, coastal flat, and mountain–lake regions are 23%–29%, 32%–38%, and 30%–37%, respectively. Additionally, the frequencies of I_w/I_u < 0.5 in the grassland, coastal flat, and mountain–lake regions are 45%–75%, 52%–70%, and 43%–53%, respectively. The frequencies of I_v/I_u < 0.8 and Iw/I_u < 0.5 in each region mean that I_u is large and the air flow is unstable and fluctuating,which will damage the wind turbines. Therefore, these conditions do not meet the wind turbine design requirements,which must be considered separately.(3) At 50-and 70-m heights, the value of the turbulence scale parameter Λ in the grassland region is greater than that in the coastal flat region, and the latter is greater than that in the mountain–lake region. Therefore, under the same conditions, some parameters, e.g., the extreme directional change and extreme operating gust at the hub height in the grassland region, are greater than those in the coastal flat region,which are greater than those in the mountain–lake region. These results provide a reference for harnessing local wind energy resources and for the selection and design of wind turbines.     

Turbulent Flow Properties Around a Staggered Wind Farm

The fundamental properties of turbulent flow around a perfectly staggered wind farm are investigated in a wind tunnel. The wind farm consisted of a series of 10 rows by 2–3 columns of miniature wind turbines spaced 5 and 4 rotor diameters in the streamwise and spanwise directions respectively. It was placed in a boundary-layer flow developed over a smooth surface under thermally neutral conditions. Cross-wire anemometry was used to obtain high resolution measurements of streamwise and vertical velocity components at various locations within and above the wind farm. The results show that the staggered configuration is more efficient in terms of momentum transfer from the background flow to the turbines compared to the case of an aligned wind turbine array under similar turbine separations in the streamwise and spanwise directions. This leads to improved power output of the overall wind farm. A simplified analysis suggests that the difference in power output between the two configurations is on the order of 10%. The maximum levels of turbulence intensity in the staggered wind farm were found to be very similar to that observed in the wake of a single wind turbine, differing substantially with that observed in an aligned configuration with similar spacing. The dramatic changes in momentum and turbulence characteristics in the two configurations show the importance of turbine layout in engineering design. Lateral homogenization of the turbulence statistics above the wind farm allows for the development of simple parametrizations for the adjustment of flow properties, similar to the case of a surface roughness transition. The development of an internal boundary layer was observed at the upper edge of the wind farm within which the flow statistics are affected by the superposition of the ambient flow and the flow disturbance induced by the wind turbines. The adjustment of the flow in this layer is much slower in the staggered situation (with respect to its aligned counterpart), implying a change in the momentum/power available at turbine locations. Additionally, power spectra of the streamwise and vertical velocity components indicate that the signature of each turbine-tip vortex structure persists to locations deep within the wind farm.     

不同下垫面的大气边界层湍流扩散特征

利用海洋、戈壁、城郊和城市4种不同下垫面上的湍流观测资料,运用多尺度湍流概念和方法计算分析了风速谱和扩散参数等湍流特征量。结果表明:下垫面的粗糙度对大气边界层湍流特征影响很大,随着粗糙度的增大,速度谱的峰值频率向高频方向移动,扩散参数增大。     

A Wind-Tunnel Investigation of Wind-Turbine Wakes: Boundary-Layer Turbulence Effects

Wind-tunnel experiments were performed to study turbulence in the wake of a model wind turbine placed in a boundary layer developed over rough and smooth surfaces. Hot-wire anemometry was used to characterize the cross-sectional distribution of mean velocity, turbulence intensity and kinematic shear stress at different locations downwind of the turbine for both surface roughness cases. Special emphasis was placed on the spatial distribution of the velocity deficit and the turbulence intensity, which are important factors affecting turbine power generation and fatigue loads in wind energy parks. Non-axisymmetric behaviour of the wake is observed over both roughness types in response to the non-uniform incoming boundary-layer flow and the effect of the surface. Nonetheless, the velocity deficit with respect to the incoming velocity profile is nearly axisymmetric, except near the ground in the far wake where the wake interacts with the surface. It is found that the wind turbine induces a large enhancement of turbulence levels (positive added turbulence intensity) in the upper part of the wake. This is due to the effect of relatively large velocity fluctuations associated with helicoidal tip vortices near the wake edge, where the mean shear is strong. In the lower part of the wake, the mean shear and turbulence intensity are reduced with respect to the incoming flow. The non-axisymmetry of the turbulence intensity distribution of the wake is found to be stronger over the rough surface, where the incoming flow is less uniform at the turbine level. In the far wake the added turbulent intensity, its positive and negative contributions and its local maximum decay as a power law of downwind distance (with an exponent ranging from −0.3 to −0.5 for the rough surface, and with a wider variation for the smooth surface). Nevertheless, the effect of the turbine on the velocity defect and added turbulence intensity is not negligible even in the very far wake, at a distance of fifteen times the rotor diameter.     

基于风廓线雷达的广东登陆台风边界层高度特征研究

针对8个登陆广东省的热带气旋,利用经过数据质量控制的风廓线雷达连续、高时空分辨率的风场观测数据,对热带气旋边界层特征进行了分析。研究结果表明:热带气旋边界层中切向风速大值区垂直范围越大、风速越强、持续时间越久,则热带气旋强度越大、登陆后强度维持时间越久。眼区外入流层厚度越大,入流层气流越强,热带气旋登陆后强度维持时间则越久。风廓线雷达信噪比垂直梯度对大气湍流信息有一定的指示作用,对于入流层高度在2000 m以下的热带气旋,其入流层顶所在高度与信噪比梯度最大值所在高度相近,对于入流层较为深厚的热带气旋,用信噪比垂直梯度确定的边界层高度虽接近入流层顶高,但仍有一定差距。不同特点的热带气旋其边界层高度并不相同,对于登陆后强度迅速减弱的热带气旋边界层高度在500~1000 m;登陆后强度持续时间短的热带气旋,其边界层高度约1000~2000 m;登陆后强度持续时间长的热带气旋,其边界层高度在2000 m之上,最高可达5000~7000 m。这些结果加深了对登陆台风边界层高度演变特征的认识。     

沿海风工程设计风速中泊松-耿贝尔法的应用

每年西北太平洋热带气旋(TC)发生的次数、移动路径和强度都是随机的,我国东南沿海各地每年受TC影响的次数便构成了某种离散型分布,而TC影响下的最大风速则可以构成某种连续型分布。该文采用上海台风研究所提供的1961—2006年TC中心风速和TC影响期间各台站大风资料,利用泊松-耿贝尔联合极值风速计算方法,计算了沿海各气象站TC影响大风的多年一遇风工程设计最大风速。结果表明:当观测资料样本序列较短,特别是像TC这样随机性很强的天气事件,泊松-耿贝尔联合极值算法更具优势;我国沿海地区有53.9%的台站50年一遇最大风速在25 m/s以下,最大风速大于42.5 m/s以上的台站分布于浙江的大陈岛、嵊山、石浦,福建的北茭和台山,广东的遮浪、上川岛和海南的西沙岛,在这些地区进行风电开发风险较大,应引起足够重视。     

A Comparative Study of Stationary and Non-stationary Wind Models Using Field Measurements

We present a comparative study of the conventional stationary wind speed model and a newly proposed non-stationary wind speed model using field measurements. The concept of, and the differences between, the two wind models are briefly reviewed. Wind data recorded by a field measurement system for wind turbulence parameters (FMS-WTP) of 1-year duration are analyzed using the two wind models. Comparisons were made between the wind characteristics obtained from the two models, including hourly mean wind speed, turbulence intensity, the wind spectrum, integral length scale, root coherence function and probability density function. The effects of wind types (monsoon or typhoon), statistical properties (stationary or non-stationary), and surface roughness (open-sea fetch or overland fetch) on wind characteristics are discussed. The comparative study demonstrates that the non-stationary wind model appears to be more appropriate than the conventional stationary wind speed model for characterizing turbulent winds of one-hour duration over complex terrain.     

Possible sources of forecast errors generated by the global/regional assimilation and prediction system for landfalling tropical cyclones. Part I: Initial uncertainties

This paper investigates the possible sources of errors associated with tropical cyclone (TC) tracks forecasted using the Global/Regional Assimilation and Prediction System (GRAPES). The GRAPES forecasts were made for 16 landfalling TCs in the western North Pacific basin during the 2008 and 2009 seasons, with a forecast length of 72 hours, and using the default initial conditions (“initials”, hereafter), which are from the NCEP-FNL dataset, as well as ECMWF initials. The forecasts are compared with ECMWF forecasts. The results show that in most TCs, the GRAPES forecasts are improved when using the ECMWF initials compared with the default initials. Compared with the ECMWF initials, the default initials produce lower intensity TCs and a lower intensity subtropical high, but a higher intensity South Asia high and monsoon trough, as well as a higher temperature but lower specific humidity at the TC center. Replacement of the geopotential height and wind fields with the ECMWF initials in and around the TC center at the initial time was found to be the most efficient way to improve the forecasts. In addition, TCs that showed the greatest improvement in forecast accuracy usually had the largest initial uncertainties in TC intensity and were usually in the intensifying phase. The results demonstrate the importance of the initial intensity for TC track forecasts made using GRAPES, and indicate the model is better in describing the intensifying phase than the decaying phase of TCs. Finally, the limit of the improvement indicates that the model error associated with GRAPES forecasts may be the main cause of poor forecasts of landfalling TCs. Thus, further examinations of the model errors are required.     

Possible Sources of Forecast Errors Generated by the Global/Regional Assimilation and Prediction System for Landfalling Tropical Cyclones. Part Ⅰ: Initial Uncertainties

This paper investigates the possible sources of errors associated with tropical cyclone(TC) tracks forecasted using the Global/Regional Assimilation and Prediction System(GRAPES). The GRAPES forecasts were made for 16 landfalling TCs in the western North Pacific basin during the 2008 and 2009 seasons, with a forecast length of 72 hours, and using the default initial conditions("initials", hereafter), which are from the NCEP-FNL dataset, as well as ECMWF initials. The forecasts are compared with ECMWF forecasts. The results show that in most TCs, the GRAPES forecasts are improved when using the ECMWF initials compared with the default initials. Compared with the ECMWF initials, the default initials produce lower intensity TCs and a lower intensity subtropical high, but a higher intensity South Asia high and monsoon trough, as well as a higher temperature but lower specific humidity at the TC center. Replacement of the geopotential height and wind fields with the ECMWF initials in and around the TC center at the initial time was found to be the most efficient way to improve the forecasts. In addition, TCs that showed the greatest improvement in forecast accuracy usually had the largest initial uncertainties in TC intensity and were usually in the intensifying phase. The results demonstrate the importance of the initial intensity for TC track forecasts made using GRAPES, and indicate the model is better in describing the intensifying phase than the decaying phase of TCs. Finally, the limit of the improvement indicates that the model error associated with GRAPES forecasts may be the main cause of poor forecasts of landfalling TCs. Thus, further examinations of the model errors are required.     

基于测风塔实测台风威马逊登陆过程的强风特性分析

利用广东省徐闻县西连镇90 m测风塔在1409号超强台风威马逊登陆期间获取的具备完整的台风代表性的观测数据以及处于台风外围的广东省茂名市博贺镇100 m测风塔的观测数据,对台风威马逊的近地层强风特性进行了分析,西连测风塔结果表明:风速时程曲线呈明显的"M"型分布特征,台风中心经过测风塔前后,风向沿逆时针方向大幅偏转约170°。风速随高度增加而增大,风速廓线较好地符合对数和幂指数律;台风过境前后,各强风区的风速廓线幂指数和粗糙长度呈先减小后增大的特点;粗糙陆地下垫面的风速廓线幂指数和粗糙长度较大。湍流强度和阵风系数在前外围强风区或后外围强风区较大,在前眼壁强风区或后眼壁强风区较小,湍流强度和阵风系数随高度增加而减小,基本符合指数为负值的幂指数律;粗糙下垫面对湍流强度和阵风系数有增大的作用。外围强风区和眼壁强风区的10 min风向变率变化较为平稳,而在眼区变动较为剧烈,在眼区,当风速达到最低值或次低值时,10 min风向变率幅值达到最大值。博贺测风塔结果表明其总体上与西连测风塔台风前外围和前眼壁强风区的情形相似。     

深圳356米气象塔观测数据的质量控制方法和空气动力学参数研究

气象高塔数据资料弥足珍贵,对其进行质量控制将为后续科学研究和业务工作的开展提供便利;此外,利用塔基观测资料计算空气动力学参数有助于校正模式空气动力学参数理论值。对2017—2018年深圳356 m气象梯度观测塔共13层的每10 s风速、风向、相对湿度、温度探测资料进行数据质量控制,基于莫宁-奥布霍夫相似理论和数据质量控制后的气象梯度观测塔近地层(10 m、20 m、40 m、50 m和80 m)1分钟平均的风温资料,利用最小二乘法拟合迭代计算了近中性条件下深圳气象梯度观测塔下垫面空气动力学粗糙度(z₀)和零平面位移(d)。结果表明:深圳气象梯度观测塔的气象探测资料数据质量很高,连续两年平均数据缺失率为1.28%,数据错误率为0.01%。近中性边界层条件下,深圳气象梯度观测塔下垫面空气动力学粗糙度均值为0.35 m,零平面位移均值为5.33 m,结果合理可信。研究表明空气动力学参数受下垫面非均匀性、植株柔软性、气流来向、风速等的共同影响。     

The Effect of Wind-Turbine Wakes on Summertime US Midwest Atmospheric Wind Profiles as Observed with Ground-Based Doppler Lidar

We examine the influence of a modern multi-megawatt wind turbine on wind and turbulence profiles three rotor diameters ( $D$ D ) downwind of the turbine. Light detection and ranging (lidar) wind-profile observations were collected during summer 2011 in an operating wind farm in central Iowa at 20-m vertical intervals from 40 to 220 m above the surface. After a calibration period during which two lidars were operated next to each other, one lidar was located approximately $2D$ 2 D directly south of a wind turbine; the other lidar was moved approximately $3D$ 3 D north of the same wind turbine. Data from the two lidars during southerly flow conditions enabled the simultaneous capture of inflow and wake conditions. The inflow wind and turbulence profiles exhibit strong variability with atmospheric stability: daytime profiles are well-mixed with little shear and strong turbulence, while nighttime profiles exhibit minimal turbulence and considerable shear across the rotor disk region and above. Consistent with the observations available from other studies and with wind-tunnel and large-eddy simulation studies, measurable reductions in wake wind-speeds occur at heights spanning the wind turbine rotor (43–117 m), and turbulent quantities increase in the wake. In generalizing these results as a function of inflow wind speed, we find the wind-speed deficit in the wake is largest at hub height or just above, and the maximum deficit occurs when wind speeds are below the rated speed for the turbine. Similarly, the maximum enhancement of turbulence kinetic energy and turbulence intensity occurs at hub height, although observations at the top of the rotor disk do not allow assessment of turbulence in that region. The wind shear below turbine hub height (quantified here with the power-law coefficient) is found to be a useful parameter to identify whether a downwind lidar observes turbine wake or free-flow conditions. These field observations provide data for validating turbine-wake models and wind-tunnel observations, and for guiding assessments of the impacts of wakes on surface turbulent fluxes or surface temperatures downwind of turbines.     

应用印痕函数研究测风塔资料代表性

影响测风塔代表性的因素很多,为了合理地确定测风塔观测结果能够代表的区域范围,需要通过适当的方法综合考虑这些因素。湍流通量代表性研究中常用的印痕函数是风速、稳定度、下垫面粗糙度等的函数,可以定量地刻画单点观测结果的代表性。利用印痕函数方法对中国风能资源观测网中来自不同下垫面的44个测风塔2010年的观测结果进行了测风塔资料的代表性分析。首先,根据两层(10、70 m)的风速和温度计算稳定度参数L_p并确定中性层结,粗糙度的计算选取了中性层结条件下满足对数风廓线的资料。然后,利用中性层结下70 m高度主导风向和风向频率较大方位的资料分析了两种印痕函数方法(Schuepp和Hsieh方法)的适用性及代表性范围与粗糙度的关系。两种方法得到的印痕函数对应距离(即代表性范围)比较接近。但是,在复杂山地Schuepp方法已不适用。当粗糙度小于10 m时,粗糙度与印痕函数对应距离的关系可用对数-线性函数描述。根据经验推荐了利用印痕函数积分90%对应的距离以及下垫面显著变化位置来确定测风塔代表范围。但是,印痕函数积分阈值仍需要通过示踪实验或数值模拟检验。