Momentum flux in the subcloud layer of a microburst-producing thunderstorm determined from JAWS Dual-Doppler data

Dual-Doppler winds at 1647 MDT for the 14 July, 1982 convective storm collected during the Joint Airport Weather Studies (JAWS) project at Denver's Stapleton International Airport were objectively analyzed to produce a three-dimensional wind field. The domain of interest had dimensions of 10 × 10 × 8.5 km centered on the microburst. Vertical velocities were computed by integrating the anelastic continuity equation downward from the storm's top. A variational approach was then employed to adjust the derived three-dimensional wind field. Subsequently, fields of deviation perturbation pressure and virtual temperature were retrieved from a detailed wind field using the three momentum equations. These retrieved fields were subjected to internal consistency checks to determine the level of confidence before interpetation. The fields were then used to calculate the generation of the vertical transport of horizontal momentum in the subcloud layer of a microburst-producing storm during the quasi-steady mature stage. Results show that the microburst occurrence in the atmospheric boundary layer (ABL) enhances eddy transfer of momentum. Direct calculation of the vertical transport of u- and v-momentum reveals that momentum was being transferred downward from the mid-levels of the storm to the microburst. The dominant processes contributing to the generation/dissipation of horizontal momentum flux were the total buoyancy production, pressure effects, vertical mean wind shear and vertical transport of momentum. The above processes play an important role in maintaining the strength of the microburst outflow in the ABL during the quasi-steady mature stage of the microburst life cycle.     

The subcloud-layer eddy kinetic energy budget of a microburst-producing thunderstorm determined from jaws dual-Doppler measurements

Dual-Doppler data collected from 1646 to 1648 MDT on 14 July, 1982 in Colorado were employed to study the eddy kinetic energy budget in the subcloud layer of a microburst-producing thunderstorm during its mature stage. Each term in the budget equation was computed from the Doppler-derived winds and retrieved thermodynamic fields within the 10 by 10 km horizontal domain. Results show that in the atmospheric boundary layer (ABL) where the microburst dominates, the turbulent flow extracts energy from the mean flow in order for the microburst to maintain its strong diverging outflow at low levels. The vertical transport of eddy kinetic energy is predominantly downward in the low layer due to the organized downdrafts in the microburst area. The horizontal flux convergence (divergence) of eddy kinetic energy by the mean and eddying motions is approximately balanced by that of the vertical flux divergence (convergence). Similarly, the contributions from the pressure and buoyancy production terms are nearly in balance. As a result, a net change of the eddy kinetic energy generation in the subcloud layer is relatively small in comparison with the individual term in the budget equation.     

2017年7月一次微下击暴流背景的多尺度特征

利用中尺度数值模式WRFV3.6模拟了2017年7月12日宁波机场附近的一次微下击暴流天气过程,结合浙江省自动站资料、机场自动观测数据、多普勒天气雷达资料等分析微下击暴流成因,结果表明：此次微下击暴流是由海风锋触发的强雷雨引起的,机场自动观测数据显示的风、温、压等气象要素的变化均呈现出明显的微下击暴流特征;数值模式较好地模拟出微下击暴流的水平风场结构;拖曳作用、下沉过程中冰雹融化、液态水和雨水持续蒸发降温作用是形成此次下击暴流的重要原因;低层位温扰动加强了垂直运动,中低层位涡异常增大区与雷暴强降水区域有较好对应关系。     

Boundary-layer structure near the cold front of a marine cyclone during “ERICA”

The boundary layer in the warm sector of a moderately deepening winter cyclone during the Experiment on Rapidly Intensifying Cyclones over the Atlantic (ERICA) is studied near the cold front. Data from the National Center for Atmospheric Research Electra research aircraft are used to examine mean and turbulence quantities. The aircraft data and supplemental data from ships, drifting buoys and moored buoys reveal an equivalent-barotropic pressure field. The area is found to be dominated by gradients in temperature and in turbulent fluxes, with changes occurring over 100 km horizontally being comparable to changes over 350 m vertically. The horizontal components of the gradients are found to be a maximum in a direction perpendicular to the front. Cross-sections perpendicular to the front are used to illustrate boundary-layer structure. Profiles of wind speed, stress, wind direction and stress direction are estimated from an Ekman model that is modified to take into account the equivalent-barotropic pressure field. Comparison of profiles from the model to the aircraft-measured data show reasonable agreement far from the front (100 km) when the model uses a constant eddy viscosity of approximately 6 kg m^–1 s^–1. Near the front there is less agreement with the model. Profiles of turbulent fluxes of momentum, heat and latent heat are divergent, with along-wind momentum flux negative and decreasing upward, cross-wind momentum flux positive and increasing upward, and heat flux and latent heat flux small, positive and decreasing upward. Far from the front, the turbulent kinetic energy budget shows that dissipation balances shear production. However, near-front behavior has an imbalance at low altitude, with shear production appearing as a TKE sink.     

Role of Downward Momentum Transport in the Formation of Severe Surface Winds

The Weather Research and Forecasting (WRF) model was used to investigate the role of downward momentum transport in the formation of severe surface winds for a squall line on 3-4 June 2009 across regions of the Henan and Shandong Provinces of China. The results show that there was a strong westerly jet belt with a wind speed greater than 30 m s 1 and a thickness of 5 km at an altitude of 11-16 km. The jet belt was accelerated, and it descended while the squall line convective system occurred. It was found that the appearance of strong negative perturbation pressure accompanied by the squall line caused the acceleration of the upper-level westerly jet and increased the horizontal wind speed by a maximum of 18%. Meanwhile, the negative buoyancy due to the loading, melting, and evaporation of cloud hydrometeors induced the downward momentum transport from the upper levels. The downward momentum transport contributed approximately 70% and the surface cold pool 30% to the formation of severe surface winds.     

Temperature stratification of the atmospheric boundary layer over the Deccan Plateau,India, during the summer monsoon

The vertical and horizontal temperature structure of the atmospheric boundary layer (ABL) were studied using aircraft observations made in the lowest 2.4 km above ground level during the summer monsoon.The vertical temperature structure of the ABL in the region may be classified into the following four categories.Category The ABL consisted of two layers of thickness 700–900 m separated by a thin transition layer. The lapse rates in the former two layers were dry adiabatic.Category The lowest layer of the ABL of thickness 400–600 m was adiabatically stratified and the overlying layer was stable with gradients of potential temperature 4–5°C km^–1. The stable layer contained a thin adiabatic stratified layer of 200–300 m thickness at a height of 1.5 km.Category The lowest 200–400 m layer of the ABL was adiabatically stratified and the overlying layer was stable with potential temperature gradients of 5–6 °C km¹.Category The ABL was mainly stable with potential temperature gradients of 6 °C km^–1 or greater. Occasionally thin layers with adiabatic stratification were found embedded in the ABL.The temperature distribution of the horizontal temperature at 900 m was mainly normal. The high-frequency portion of the spectra lying between 0.05 and 0.16 Hz (corresponding to wave length 1 km to 300 m) oscillated around the –\2/3 power law line. The spectral curve showed a significant peak at 0.011 Hz having a wave-length of 5 km.Department of Geoscience, North Carolina State University, Raleigh, NC, 27650, U.S.A.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.     

Regional shear stress of broken forest from radiosonde wind profiles in the unstable surface layer

Mean wind speed profiles were measured by tracking radiosondes in the unstable atmospheric boundary layer (ABL) over the forested Landes region in southwestern France. New Monin-Obukhov stability correction functions, recently proposed following an, analysis by Kader and Yaglom, as well as the Businger-Dyer stability formulation were tested, with wind speeds in the surface sublayer to calculate the regional shear stress. These profile-derived shear stresses were compared with eddy correlation measurements gathered above a mature forest stand, at a location roughly, 4.5 km from the radiosonde launch site. The shear stress values obtained by means of the newly proposed stability function were in slightly better agreement with the eddy correlation values than those obtained by means of a Businger-Dyer type stability function. The general robustness of the profile method can be attributed in part to prior knowledge of the regional surface roughness (z ₀=1.2 m) and the momentum displacement height (d ₀=6.0 m), which were determined from neutral wind profile analysis. The 100 m drag coefficient for the unstable conditions above this broken forest surface was found to beu _* ² /V ₁₀₀ ² =0.0173.     

Topographic influence on surface winds

Using data collected during 1975–1976 from a meteorological network operating in the vicinity of the Columbia Generating Site approximately 8 km south of Portage, Wis., the influence of the Baraboo Hills on the surface wind field is determined. Half-hour means of wind speed and direction measured at 9 m at three sites were used to compute divergence and vorticity using Bellamy's method. The data were grouped into 18 sectors each 20 deg wide and averages computed for each quantity. Results indicate that for wind directions perpendicular to the eastern edge of the Baraboo Hills, the surface (9m) wind field is significantly perturbed up to 4 km from the bluffs. The largest convergence of 2.1 × 10^–4 s^–1 occurs with 160 deg wind direction and the largest divergence of 1.2 × 10^–4 s^–1 with 290 deg wind direction. The maximum anticyclonic vorticity was 1.6 × 10^–4 s^–1 at 210 deg and the maximum cyclonic vorticity was 1.6 × 10^–4 s^–1 at 330 deg.     

Development of an eye-wall like structure in a tropical cyclone model simulation

Structural changes during the intensification of a tropical storm into a hurricane in a numerical simulation are examined. A 10 layer primitive equation model that employs a horizontal grid spacing of 20 km over 4400 × 4400 km area is integrated. An elongated band in vertical motion over the storm area intensifies slowly during the first few hours. In the upper troposphere high pressures arise due to condensational heating. Between 8–12 h strong outflow winds develop in the upper troposphere due to the increased pressure gradients. Strong divergence occurs in the outflow wind region, and a large increase in the vertical motion, condensational heating and intensification rate of the storm ensues. Between 12–24 h the elongated band of the storm stage transforms into an eye-wall like structure, and the tropical storm intensifies into a hurricane. Regions with negative moist potential vorticity appear in the high troposphere. Widening of area of condensation and slanting of the convergence area occurs with height in the high level negative moist potential vorticity regions. Results suggest that the formation of anvil clouds in some cases may be due to the development of slantwise convection on the outer periphery of a hurricane's eye-wall.     

Gustiness and coherent structure under weak wind period in atmospheric boundary layer

Statistical analysis of turbulent and gusty characteristics in the atmospheric boundary layer under weak wind period has been carried out.The data used in the analysis were from the multilevel ultrasonic anemometer-thermometers at 47 m,120 m,and 280 m levels on Beijing 325 m meteorological tower.The time series of 3D atmospheric velocity were analyzed by using conventional Fourier spectral analysis and decompose into three parts:basic mean flow(period 10 min),gusty disturbances(1 min period 10 min)and turbulence fluctuations(period 1 min).The results show that under weak mean wind condition:1)the gusty disturbances are the most strong fluctuations,contribute about 60% kinetic energy of eddy kinetic energy and 80% downward flux of momentum,although both the eddy kinetic energy and momentum transport are small in comparison with those in strong mean wind condition;2)the gusty wind disturbances are anisotropic;3)the gusty wind disturbances have obviously coherent structure,and their horizontal and vertical component are negatively correlated and make downward transport of momentum more effectively;4)the friction velocities related to turbulence and gusty wind are approximately constant with height in the surface layer.     

Observations and Modelling of the Atmospheric Boundary Layer Over Sea-Ice in a Svalbard Fjord

Sonic anemometer and profile mast measurements made in Wahlenbergfjorden, Svalbard Arctic archipelago, in May 2006 and April 2007 were employed to study the atmospheric boundary layer over sea-ice. The turbulent surface fluxes of momentum and sensible heat were calculated using eddy correlation and gradient methods. The results showed that the literature-based universal functions underestimated turbulent mixing in strongly stable conditions. The validity of the Monin-Obukhov similarity theory was questionable for cross-fjord flow directions and in the presence of mesoscale variability or topographic effects. The aerodynamic roughness length showed a dependence on the wind direction. The mean roughness length for along-fjord wind directions was (2.4 ± 2.6) × 10⁻⁴ m, whereas that for cross-fjord directions was (5.4 ± 2.8) × 10⁻³ m. The thermal stratification and turbulent fluxes were affected by the synoptic situation with large differences between the 2 years. Channelling effects and drainage flows occurred especially during a weak large-scale flow. The study periods were simulated applying the Weather Research and Forecasting (WRF) model with 1-km horizontal resolution in the finest domain. The results for the 2-m air temperature and friction velocity were good, but the model failed to reproduce the spatial variability in wind direction between measurement sites 3 km apart. The model suggested that wind shear above the stable boundary layer provided a non-local source for the turbulence observed.     

Gravity Wave Momentum Fluxes and Surface Pressure Drags due to the Pyrénées: Variation with Numerical Model Horizontal Resolution

Summary. ?A hydrostatic numerical model is used to simulate the lee wave event IOP3 (0000 GMT to 1200 GMT 15^th October 1990) from the PYREX mountain experiment. Results from integrations at different horizontal resolutions are used to investigate the effect on surface pressure drag and the vertical flux of horizontal momentum due to orographically forced gravity waves. In particular, results showing the dependence on resolution of the partitioning between resolved and parametrized wave drag and fluxes are presented. With the model horizontal gridlength changing from 50 km to 10 km the majority of wave momentum flux changes from being parametrized to becoming resolved. More significantly, there is a change in the profile of flux with height. At 50 km resolution the largest inferred mean flow decelerations are at lower stratospheric level due to the parametrization scheme. At 10 km resolution this is shifted, with less deceleration high up and more wave deceleration lower down in the troposphere. Numerical weather prediction models are now beginning to take account of such low level drag with beneficial results. Received March 2, 1999/Revised July 15, 1999     

MOMENTUM BUDGET AND SHELTER MECHANISM OF BOUNDARY-LAYER FLOW NEAR A SHELTERBELT

We use a nonhydrostatic shelterbelt boundary-layer turbulence model with Mellor–Yamada second-order closure to evaluate quantitatively the dynamic processes of surface boundary-layer flow perturbed by shelterbelts of different densities and to understand the shelter mechanism. We first analyze the drag exerted on air by shelterbelts of different densities, a root cause of any shelter function, and the resulting wind reduction. The results show that the effectiveness of a shelter is determined not only by its total drag but also by the distribution of the drag-generated momentum deficit in the sheltered area, and that medium-dense shelterbelts have the maximum shelter effect. We also analyze the horizontal momentum budget and find that the shelter mechanism is the product of several processes. The results reveal that strong vertical mean transport and the pressure gradient also play important roles in shelter efficiency. The pressure perturbation caused by the shelter extends far downstream of the shelter, and combines with advective transport to provide the larger shelter efficiency of medium-dense shelterbelts. We finally analyze the changes of perturbed pressure, turbulence, and vertical velocity with shelterbelt density to further clarify the shelter mechanism.     

Equilibrium Atmospheric Boundary-Layer Flow: Computational Fluid Dynamics Simulation with Balanced Forces

Forcing relationships in steady, neutrally stratified atmospheric boundary-layer (ABL) flow are thoroughly analyzed. The ABL flow can be viewed as balanced between a forcing and a drag term. The drag term results from turbulent stress divergence, and above the ABL, both the drag and the forcing terms vanish. In computational wind engineering applications, the ABL flow is simulated not by directly specifying a forcing term in the ABL but by specifying boundary conditions for the simulation domain. Usually, these include the inflow boundary and the top boundary conditions. This ‘boundary-driven’ ABL flow is dynamically different from its real counterpart, and this is the major reason that the simulated boundary-driven ABL flow does not maintain horizontal homogeneity. Here, first a dynamical approach is proposed to develop a neutrally stratified equilibrium ABL flow. Computational fluid dynamics (CFD) software (Fluent 6.3) with the standard \(k\) – \(\varepsilon \) turbulence model is employed, and by applying a driving force profile, steady equilibrium ABL flows are simulated by the model. Profiles of wind speed and turbulent kinetic energy (TKE) derived using this approach are reasonable in comparison with the conventional logarithmic law and with observational data respectively. Secondly, the equilibrium ABL profiles apply as inflow conditions to simulate the boundary-driven ABL flow. Simulated properties between the inlet and the outlet sections across a fetch of 10 km are compared. Although profiles of wind speed, TKE, and its dissipation rate are consistently satisfactory under higher wind conditions, a deviation of TKE and its dissipation rate between the inlet and outlet are apparent (7–8 %) under lower wind-speed conditions (2 m s \(^{-1}\) at 10 m). Furthermore, the simulated surface stress systematically decreases in the downwind direction. A redistribution of the pressure field is also found in the simulation domain, which provides a different driving pattern from the realistic case in the ABL.     

Wind shear at tilted inversions

Vertical wind shear at a temperature inversion can be caused by baroclinicity associated with a tilt of the inversion. Four observational cases of tilted inversions are presented. The tilts on horizontal scales of 20–100 km range from 2–10 × 10^-3 and the vertical wind shear is between 1 and 25 m/s per 100 m. In general, there is remarkable agreement between observed and geostrophic wind shear.The observations show that the inversion tilt is particularly strong at the edges of mesoscale cloud fields. The Richardson number can reach subcritical values. Cloud fields may be surrounded by a cyclonically rotating wind field and cloud gaps by an anticyclonically rotating wind field.     

Cross‐shore variations of near‐surface wind velocity and atmospheric turbulence at the land‐sea boundary during CASP

Abstract

Airborne measurements of mean wind velocity and turbulence in the atmospheric boundary layer under wintertime conditions of cold offshore advection suggest that at a height of 50 m the mean wind speed increases with offshore distance by roughly 20% over a horizontal scale of order 10 km. Similarly, the vertical gust velocity and turbulent kinetic energy decay on scales of order 3.5 km by factors of 1.5 and 3.2, respectively. The scale of cross‐shore variations in the vertical fluxes of heat and downwind momentum is also 10 km, and the momentum flux is found to be roughly constant to 300 m, whereas the heat flux decreases with height. The stability parameter, z/L (where z = 50 m and L is the local Monin‐Obukhov length), is generally small over land but may reach order one over the warm ocean. The magnitude and horizontal length scales associated with the offshore variations in wind speed and turbulence are reasonably consistent with model results for a simple roughness change, but a more sophisticated model is required to interpret the combined effects of surface roughness and heat flux contrasts between land and sea.

Comparisons between aircraft and profile‐adjusted surface measurements of wind speed indicate that Doppler biases of 1–2 m s^?1 in the aircraft data caused by surface motions must be accounted for. In addition, the wind direction measurements of the Minimet anemometer buoy deployed in CASP are found to be in error by 25 ± 5°, possibly due to a misalignment of the anemometer vane. The vertical fluxes of heat and momentum show reasonably good agreement with surface estimates based on the Minimet data.     

Micrometeorological Observations of a Microburst in Southern Finland

On the afternoon of 3 July 2004 in Hyytiälä (Juupajoki, Finland), convective cells produced a strong downburst causing forest damage. The SMEAR II field station, situated near the damage site, enabled a unique micrometeorological analysis of a microburst with differences above and inside the canopy. At the time of the event, a squall line associated with a cold front was crossing Hyytiälä with a reflectivity maximum in the middle of the squall line. A bow echo, rear-inflow notch, and probable mesovortex were observed in radar data. The bow echo moved west-north-west, and its apex travelled just north of Hyytiälä. The turbulence data were analysed at two locations above the forest canopy and at one location at sub-canopy. At 1412 EET (Eastern European Time, UTC+2), the horizontal and vertical wind speed increased and the wind veered, reflecting the arrival of a gust front. At the same time, the carbon dioxide concentration increased due to turbulent mixing, the temperature decreased due to cold air flow from aloft and aerosol particle concentration decreased due to rain scavenging. An increase in the number concentration of ultra-fine particles (< 10 nm) was detected, supporting the new particle formation either from cloud outflow or due to rain. Five minutes after the gust front (1417 EET), strong horizontal and downward vertical wind speed gusts occurred with maxima of 22 and 15 m s^?1, respectively, reflecting the microburst. The turbulence spectra before, during and after the event were consistent with traditional turbulence spectral theory.     

A Simple Parametrization for the Concentration Variance Dissipation in a Lagrangian Single-Particle Model

A new quasi-analytical mixed-layer model is formulated describing the evolution of the convective atmospheric boundary layer (ABL) during cold-air outbreaks (CAO) over polar oceans downstream of the marginal sea-ice zones. The new model is superior to previous ones since it predicts not only temperature and mixed-layer height but also the height-averaged horizontal wind components. Results of the mixed-layer model are compared with dropsonde and aircraft observations carried out during several CAOs over the Fram Strait and also with results of a 3D non-hydrostatic (NH3D) model. It is shown that the mixed-layer model reproduces well the observed ABL height, temperature, low-level baroclinicity and its influence on the ABL wind speed. The mixed-layer model underestimates the observed ABL temperature only by about 10 %, most likely due to the neglect of condensation and subsidence. The comparison of the mixed-layer and NH3D model results shows good agreement with respect to wind speed including the formation of wind-speed maxima close to the ice edge. It is concluded that baroclinicity within the ABL governs the structure of the wind field while the baroclinicity above the ABL is important in reproducing the wind speed. It is shown that the baroclinicity in the ABL is strongest close to the ice edge and slowly decays further downwind. Analytical solutions demonstrate that the \(\mathrm{e}\)-folding distance of this decay is the same as for the decay of the difference between the surface temperature of open water and of the mixed-layer temperature. This distance characterizing cold-air mass transformation ranges from 450 to 850 km for high-latitude CAOs.