用单部多普勒雷达径向速度资料反演的风场资料及其应用

利用安装在南京气象学院的天气多普勒雷达收集到的雷达资料,分析了在江淮流域大暴雨后期静止锋附近风场的中尺度结构特征。结果表明:静止锋附近的近地面层为冷平流,往上为暖平流。低层吹东北风,随着高度增大,风速减小,在2km高度附近风速最小,以后转为西南风,并反映出水平风场中存在垂直切变。     

福建省风廓线雷达资料在一次强对流天气过程中的应用

利用福建省永安站的CFL-03边界层风廓线雷达提供的资料,分析了福建省永安市2012年4月1l～12日的一次强对流天气过程,结果表明：风廓线雷达水平风资料可以相对连续地反映测站上空风场垂直结构和变化特点,直观而精细地反映出天气过程的演变特征。风廓线实况中的风向风速变化,可指示高空槽和气旋配合过境的情况,判断中低空急流的强弱;强降水出现前风随高度的变化存在着明显不连续现象,风向风速切变明显。这些特征可作为对降水性质、落区、持续时间等作出短时或临近预报的依据。垂直风资料可反映出降水的开始、结束和降水的强度,其波动发展的高度成为判断对流发展强弱的一个重要指标。降水期间功率谱密度出现双峰谱甚至多峰谱,而降雹时间段波束图出现了速度模糊,证明强对流系统通过测站。风廓线雷达对天气发展趋势提供可靠的指示作用。     

冰雹天气过程的综合分析

选取发生在邯郸辖区内的5次冰雹天气过程,应用常规观测资料并结合雷达回波资料,从天气形势、物理量场和雷达回波演变特征等几方面进行了综合分析。结果表明:冰雹天气出现在对流性不稳定层结条件下,中高层干冷低层暖湿,地面有中尺度辐合切变线配合,0℃层、-20℃层和强垂直风切变的高度适宜。多普勒雷达能很好地监测中尺度天气系统的发展演变过程,回波强度和回波顶高度的变化、速度对和逆风区的出现、1 h降水累积值和垂直积分液态水含量的急剧增加等都对冰雹的出现具有指示意义。     

林芝机场低空风的垂直切变特征及其对飞行安全的影响

为了解林芝机场低空风的垂直切变特征及其对飞行安全的影响,利用林芝机场风廓线雷达资料,计算林芝机场低空风的垂直风切变,讨论了中度及以上风切变的时空分布特征。分析发现,林芝机场低空风切变频繁出现,严重影响林芝机场的飞行安全。林芝机场强顺风垂直切变出现频率为1.0%,严重顺风垂直切变为0.2%,最大顺风垂直切变13.4m/s/30m。强侧风垂直切变出现频率为1.9%,严重侧风垂直切变为0.6%,最大侧风垂直切变16.0m/s/30m。林芝机场低空强垂直切变对起降航班的威胁具有明显的日变化与季节变化特征,强垂直切变主要出现在各个季节的下午及傍晚,8月较弱。     

大连机场典型低空风切变的成因分析

风切变是指空间两点之间的风的矢量差,风切变分为水平风切变和垂直风切变。航空气象学中,低空风切变通常是指近地面600 m高度以下的垂直风切变,低空风切变严重危害航空活动安全,尤其是在飞机起飞和着陆阶段。为了提高对低空风切变的预报准确率和提前预警时间,利用大连机场话音方式航空器报告统计大连机场出现的低空风切变特点。大连机场低空风切变出现最多的季节是秋末和春季,其中以西北大风和偏北大风造成的低空风切变最为典型,这种切变常常出现在地面风向为360°或350°、平均风速9 m·s^-1、阵风风速16 m·s^-1的观测时次附近。选取3个大连机场典型低空风切变的天气个例,利用NCEP/NCAR FNL(1°×1°)资料分析大连机场出现典型低空风切变时的天气形势特点,利用WRF模式输出资料,分析了大连机场出现典型低空风切变时低层的风场特征。分析结果表明,WRF模式输出的风场变化较小,不能体现低空风切变的突变特征,利用模式输出资料计算得到的反映偏差风相对大小的I值与风切变的出现有明显关系,当I大于0.5且逐渐增加时,容易出现低空风切变,当I小于0.5且逐渐减小时,不易出现低空风切变,可以利用I值及其变化情况进行风切变的预报预警。     

基于雷达资料四维变分同化技术对北京地区一次下山突发性增强风暴热动力机制的模拟分析

利用三维数值云模式和雷达资料四维变分(4DVar)同化技术,通过对4部新一代多普勒天气雷达探测资料进行快速更新同化和云尺度模拟,初步分析了2009年8月1日发生在北京地区的一次短生命史、突发性增强风暴的低层动力和热力影响机制。此次风暴过程处于弱天气尺度背景和弱层结背景下,冷池和低层环境风场相互作用是造成山上对流风暴增强传播下山的关键机制,而风暴的短生命史和平原地区上空弱风垂直切变环境有关:在对流风暴产生的初期,由于平原地区局地热力、动力场分布的差异,在平原地区西部近地面形成冷池结构,而冷池的“障碍物”作用进而阻碍环境风场的传播。在此机制下,导致在冷池东南边缘附近形成辐合中心、较强的低层水平风垂直切变和全螺旋度大值中心,有助于风暴传播下山。在风暴临近山边阶段,平原地区原有冷池的“绕流”等机制仍然有助于形成有利于主体风暴传播下山增强的近地面辐合中心、强低层水平风垂直切变和全螺旋度大值中心等环境。此外,随着山上风暴降水产生若干冷池,由于风暴形成的阵风锋抬升作用以及新生冷池与老冷池的逐渐发展并相互靠近,使冷池之间暖空气不断抬升,在冷池之间低层形成较强的辐合中心、全螺旋度大值中心。并且,由于冷池边缘的热力场分布不均匀,同样在冷池边缘形成较大扰动气压和扰动温度,增大了垂直加速度,在冷池之间中高层形成上升气流区,这些机制使北部风暴重新增强和新生风暴产生的同时,最终也导致这些风暴互相靠近,合并组织成带状回波。风暴在平原传播阶段,带状回波产生的冷池进一步增强,并明显扩展。低层风场指示冷池出流(阵风锋)更加强烈且存在明显的“前冲”特征,显现出部分飑线系统特征。但是,由于此时平原地区处于弱风垂直切变环境,此时冷池强于低层风垂直切变,即冷池产生的负涡度大于低层风垂直切变产生的正涡度,因此,冷池前沿的上升气流向后倾斜并导致阵风锋逐渐离开主体风暴,不利于沿着出流边界形成新的对流单体,从而不利于维持对流风暴系统的发展传播。随后,阵风锋和前方东南气流交汇,形成新的孤立单体。并且,基于模拟结果计算了与对流系统发展密切相关的全螺旋度、风垂直切变。结果显示,风垂直切变(尤其是0—3 km)和全螺旋度与风暴发生和传播位置及强度相关性较高,反映出模拟量对带状回波风暴过程具有较好的指示意义。     

一次强风切变过程的分析及成因讨论

2006年6月21日09:04南京禄口机场监测到较强风切变。龙王山(南京)多普勒天气雷达观测图上显示,此时恰有一直径为50 km的对流云团自西向东经过禄口机场。南京市地面自动气象站的观测也表明,此时在南京地区近地层气流出现急剧变化。利用多普勒天气雷达资料对风暴进行立体结构分析,并结合地面自动气象站资料分析得知,这次过程在3 km以下的低层有水平方向的风切变,在1 km高度处有垂直方向的风切变。从多普勒天气雷达的垂直积分液态含水量产品能够更直接地看到对流云团的发展过程,这有助于预报风切变发生的位置。这次雷暴系统中宏下击暴流引起的高空气流下沉,至地面后向两侧辐散,造成水平方向及垂直方向都有强风切变存在。     

伊犁一次区域性暴雨分析

利用伊宁新一代天气雷达资料,结合高空和地面观测资料,分析了2010年7月19日伊犁地区一次局地暴雨天气过程。结果表明：此次强对流天气过程的主要影响系统为500 hPa中亚低槽、200hpa高空急流、低层风速辐合和地面雷暴高压。较强的层结不稳定和低层垂直风切变有利于对流的产生;云图和雷达资料分析表明,此次局地暴雨是由中尺度强对流云团产生,具有典型的对流单体形成、发展成回波短带合并形成带状回波,该带状回波最后演变成一个尺度较大的弓形回波。     

一次弱对流天气降雹成灾的雷达回波特征分析

通过对新疆玛河流域一次冰雹天气雷达资料及地面实况资料的分析,研究弱对流云降雹的雷达反射率因子与径向速度场的回波特征,进一步探讨弱对流云降雹与0℃层高度较低及垂直风切变的关系：弱对流天气当日零度层高度值可作为分析弱冰雹云的依据;高低空风速的差异形成的风切变造成的动压力产生垂直加速度,使对流运动得到加强和维持,也是此次降雹的动力机制。分析此次天气,揭示了新疆玛河流域弱对流云降雹天气过程生消的一些重要特征和变化,对雹云的识别和分析具有一定的指导意义。     

利用激光雷达判别机场晴空风切变事件成因

为了研究高原机场的风切变特征,本文利用丽江机场自动观测站资料、NCEP再分析资料,引入三维激光雷达测风资料,对丽江机场2020年3月29日晴空低空风切变特征进行分析。结果表明,此次低空风切变过程主要是受地面变压风辐合、地面增温湍流加强以及动量下传等因素共同影响而形成的。风切变发生期间低空环境风场与地面风场变化基本吻合。5次低空风切变事件中,出现在11:00—12:20风向转变期间的2次风切变为侧风风向切变,出现在12:20—14:20风速逐渐增大期间的3次风切变为风速切变。风切变指数均大于0.1(1/s),切变强度为中度及以上,严重影响飞行安全。研究结果为机场风切变的激光测风雷达监测预警提供了参考依据。     

Summer dryness in a warmer climate: a process study with a regional climate model

Earlier GCM studies have expressed the concern that an enhancement of greenhouse warming might increase the occurrence of summer droughts in mid-latitudes, especially in southern Europe and central North America. This could represent a severe threat for agriculture in the regions concerned, where summer is the main growing season. These predictions must however be considered as uncertain, since most studies featuring enhanced summer dryness in mid-latitudes use very simple representations of the land-surface processes ("bucket" models), despite their key importance for the issue considered. The current study uses a regional climate model including a land-surface scheme of intermediate complexity to investigate the sensitivity of the summer climate to enhanced greenhouse warming over the American Midwest. A surrogate climate change scenario is used for the simulation of a warmer climate. The control runs are driven at the lateral boundaries and the sea surface by reanalysis data and observations, respectively. The warmer climate experiments are forced by a modified set of initial and lateral boundary conditions. The modifications consist of a uniform 3 K temperature increase and an attendant increase of specific humidity (unchanged relative humidity). This strategy maintains a similar dynamical forcing in the warmer climate experiments, thus allowing to investigate thermodynamical impacts of climate change in comparative isolation. The atmospheric CO ₂ concentration of the sensitivity experiments is set to four times its pre-industrial value. The simulations are conducted from March 15 to October 1st, for 4 years corresponding to drought (1988), normal (1986, 1990) and flood (1993) conditions. The numerical experiments do not present any great enhancement of summer drying under warmer climatic conditions. First, the overall changes in the hydrological cycle (especially evapotranspiration) are of small magnitude despite the strong forcing applied. Second, precipitation increases in spring lead to higher soil water recharge during this season, compensating for the enhanced soil moisture depletion occurring later in the year. Additional simulations replacing the plant control on transpiration with a bucket-type formulation presented increased soil drying in 1988, the drought year. This suggests that vegetation control on transpiration might play an important part in counteracting an enhancement of summer drying when soil water gets limited. Though further aspects of this issue would need investigating, our results underline the importance of land-surface processes in climate integrations and suggest that the risk of enhanced summer dryness in the region studied might be less acute than previously assumed, provided the North American general circulation does not change markedly with global warming.     

Continental vegetation as a dynamic component of a global climate model: a preliminary assessment

A simplified vegetation distribution prediction scheme is used in combination with the Biosphere-Atmosphere Transfer Scheme (BATS) and coupled to a version of the NCAR Community Climate Model (CCM1) which includes a mixed-layer ocean. Employed in an off-line mode as a diagnostic tool, the scheme predicts a slightly darker and slightly rougher continental surface than when BATS' prescribed vegetation classes are used. The impact of tropical deforestation on regional climates, and hence on diagnosed vegetation, differs between South America and S.E. Asia. In the Amazon, the climatic effects of removing all the tropical forest are so marked that in only one of the 18 deforested grid elements could the new climate sustain tropical forest vegetation whereas in S.E. Asia in seven of the 9 deforested elements the climate could continue to support tropical forest. Following these off-line tests, the simple vegetation scheme has been coupled to the GCM as an interactive (or two-way) submodel for a test integration lasting 5.6 yr. It is found to be a stable component of the global climate system, producing only ~ 3% (absolute) interannual changes in the predicted percentages of continental vegetation, together with globally-averaged continental temperature increases of up to + 1.5 °C and evaporation increases of 0 to 5 W m^–2 and no discernible trends over the 67 months of integration. On the other hand, this interactive land biosphere causes regional-scale temperature differences of ± 10 °C and commensurate disturbances in other climatic parameters. Tuning, similar to the q-flux schemes used for ocean models, could improve the simulation of the present-day surface climate but, in the longer term, it will be important to focus on predicting the characteristics of the continental surface rather than simple vegetation classes. The coupling scheme will also have to allow for vegetation responses occurring over longer timescales so that the coupled system is buffered from sudden shocks.     

Numerical analysis of a thermotopographically-induced mesoscale circulation in a mountain basin using a non-hydrostatic model

Summary The boundary-layer wind field during weak synoptic conditions is largely controlled by the nature of the landscape. Mesoscale (sub-synoptic) circulations result from horizontal gradients of sensible heat flux due to variation in local topography, variation in surface-cover, and discontinuities such as land-sea contrasts. Such flows are usually referred to as thermally-driven circulations, and are diurnal in nature and often predictable. In this paper we use a state-of-the-art non-hydrostatic computer model to shed light on the physical mechanisms that drive a persistent easterly wind that develops in the afternoon in the Mackenzie Basin, New Zealand. The easterly – Canterbury Plains Breeze (CPB) – is observed early in the afternoon and is often intense, with mean wind speeds reaching up to 12 m s⁻¹. Although computer modelling in mountainous terrain is extremely challenging, the model is able to simulate this circulation satisfactorily. To further investigate the mechanisms that generate the Canterbury Plains Breeze, two additional idealized model experiments are performed. With each experiment, the effects of the synoptic scale wind and the ocean around the South Island, New Zealand were successively removed. The results show that contrary to previous suggestions, the Canterbury Plains Breeze is not an intrusion of the coastal sea breeze or the Canterbury north-easterly, but can be generated by heating of the basin alone. This conclusion highlights the importance of mountain basins and saddles in controlling near-surface wind regimes in complex terrain.     

A cold air lake formation in a basin — a simulation with a mesoscale numerical model

Summary A formation of a cold air lake in a basin is studied with a mesometeorological model.A dynamic Boussinesq hydrostatic mesoscale numerical model is developed in a staggered orthogonal grid with a horizontal resolution of 1 km and with a varying vertical grid. The topography is presented in a block shape so that computation levels are horizontal.The mesometeorological model is tested in three idealized topography cases (a valley, a single mountain, a basin) and test results are discussed.In an alpine basin surrounded by mountains and plateaus the air is supposed to be stagnant at the beginning of the night. Due to differences in radiation cooling an inversion layer is formed in the basin and local wind circulation is studied by model simulations.With 14 Figures     

Effects of a Fence on Pollutant Dispersion in a Boundary Layer Exposed to a Rural-to-Urban Transition

Simultaneous particle-image velocimetry and laser-induced fluorescence combined with large-eddy simulations are used to investigate the flow and pollutant dispersion behaviour in a rural-to-urban roughness transition. The urban roughness is characterized by an array of cubical obstacles in an aligned arrangement. A plane fence is added one obstacle height h upstream of the urban roughness elements, with three different fence heights considered. A smooth-wall turbulent boundary layer with a depth of 10h is used as the approaching flow, and a passive tracer is released from a uniform line source 1h upstream of the fence. A shear layer is formed at the top of the fence, which increases in strength for the higher fence cases, resulting in a deeper internal boundary layer (IBL). It is found that the mean flow for the rural-to-urban transition can be described by means of a mixing-length model provided that the transitional effects are accounted for. The mixing-length formulation for sparse urban canopies, as found in the literature, is extended to take into account the blockage effect in dense canopies. Additionally, the average mean concentration field is found to scale with the IBL depth and the bulk velocity in the IBL.     

Evaporation from a blanket bog in a foggy coastal environment

Frequent fog severely restricts evaporation from blanket bogs in Newfoundland because it more than halves the radiant energy input, and it eliminates the vapor pressure deficit, resulting in evaporation at the equilibrium rate (average = 0.99 during fog). During these periods, there is no surface resistance to evaporation because the bog has been wetted by fog drip, and although the latent heat flux dominates over sensible heat (average = 0.8), both are small. In contrast, the surface dries during clear periods, increasing the surface resistance to evaporation so that sensible heat becomes more important ( = 1.05). When the mosses are dry, evaporation is below the equilibrium rate ( = 0.87), although the higher available energy ensures that actual evaporation is higher. During clear periods, daily evaporation averaged 2.5 mm, compared to 1.1 and 0.7 mm for fog and rain, respectively. The suppressed evaporation at this site is important in maintaining appropriate hydrological conditions for blanket bog development.     

Parameterization of a momentum source in a tropical cumulus ensemble

Summary An eddy effect of tropical deep convection on the large-scale momentum, resp vorticity budget is investigated. The process is specified by a simple parameterization approach which is based on a concept of rotating clouds exerting a momentum on the large-scale flow. The cloud rotation is associated with the thermal properties of a cloud ensemble by the principle of conservation of potential vorticity. A decomposition of cloud classes is applied in consistency with the thermodynamical parameterization scheme of Arakawa and Schubert (1974).The parameterization is tested with observations of GATE74, Phase III. The data are processed on a B/C-scale grid (55km) in a region within 9N and 16N, and between 21W and 27W, and with a vertical resolution of 41 levels. The parameterization results correspond to the observed patterns, especially in situations with strong large-scale wind shear. The findings suggest that certain large-scalle flow regimes provoke convective scale momentum generation rather than redistributing large-scale momentum by convective circulations.With 10 Figures     

Spatial structure of a jet flow at a river mouth

The present work concentrates on the latest data measured in the Jordan river flow in lake Kinneret. Spectral characteristics of fluctuating velocity components have been obtained and processed. The three-dimensional structure of turbulence along the flow has been described. The main features of the jet flow turbulence in the river mouth are: a) The supply of turbulent energy changes due to different mechanisms along the flow. b) The structure of turbulence formed in the river decays rapidly along the flow, and c) In the sand area and beyond it, a significant generation of turbulent energy occurs. Quantitative estimations of the above effects were carried out.     

Atmospheric dispersion of a heavier-than-air gas near a two-dimensional obstacle

Flow over a two-dimensional obstacle and dispersion of a heavier-than-air gas near the obstacle were studied. Two species, one representing air and the other representing the heavier-than-air gas were treated. Equations for mass and momentum were cast in mass-averaged form, with turbulent Reynolds stresses and mass fluxes modeled using eddy-viscosity and diffusivity hypotheses. A two-equation k- turbulence model was used to determine the effective turbulent viscosity. Streamline curvature and buoyancy corrections were added to the basic turbulence formulation. The model equations were solved using finite difference techniques. An alternating-direction-implicit (ADI) technique was used to solve the parabolic transport equations and a direct matrix solver was used to solve the elliptic pressure equation.Mesh sensitivities were investigated to determine the optimum mesh requirements for the final calculations. It was concluded that at least 10 grid spaces were required across the obstacle width and 15 across the obstacle height to obtain valid solutions. A non-uniform mesh was used to concentrate the grid points at the top of the obstacle.Experimental measurements were made with air flow over a 7.6 by 7.6 cm obstacle in a boundary-layer wind tunnel. Smoke visualization revealed a low-frequency oscillation of the bubble downstream of the obstacle. Hot-wire anemometer data are presented for the mean velocity and turbulent kinetic energy at the mid-plane of the obstacle and the mid-plane of the downstream recirculation bubble. A single hot-wire probe was found to be suitable for determining mean streamwise velocities with an accuracy of 11 %. The downstream recirculation bubble was unsteady and had a length range from 3 to 8 obstacle lengths.The experimental results for flow over the obstacle were compared with numerical calculations to validate the numerical solution procedure. A sensitivity study on the effect of curvature correction and variation of turbulence model constants on the numerical solution was conducted. Calculations that included the curvature correction model gave a downstream recirculation bubble length of 5.9 obstacle lengths while excluding the correction reduced this length to 4.4.In the second part of the study, numerical calculations were performed for the dispersion of a heavier-than-air gas in the vicinity of the two-dimensional obstacle. Characteristics of an adiabatic boundary layer were used in these calculations. The densities of the contaminant gases were 0, 25 and 50% greater than the air density. Calculations were performed with the contaminant injection source upstream and downstream of the obstacle.Use of the pressure gradient model reduced the size of the dense gas cloud by as much as 12%. The curvature correction model also affected the cloud expanse by reducing the effective turbulent viscosity in the downstream recirculation bubble. The location of the injection source had the largest impact on the cloud size. The area of the cloud within the 5 % contour was three times larger for downstream injection than for upstream injection.