A numerical case study of a cold air outbreak over East Asia during late spring

Summary A case study of a late spring cold air outbreak associated with heavy rainfall and snow storms over East Asia is carried out with a numerical model. This was the most severe case over northern China in 50 years. In the lower troposphere, the evolution of the edge anticyclone and the northerly low-level orographic jet east of Qinghai-Tibetan Plateau are well simulated in the ECMWF limited area model but the associated cold surge is under predicted. In the experiments with envelope and valley filled orography, results are improved. Neither cold surge nor strong frontogenesis is simulated in a no latent heating experiment, indicating that the feedback of latent heat release is quite essential to the burst of cold air in this case.With 9 Figures     

Structure of the marine atmospheric boundary layer during two cold air outbreaks of varying intensities: GALE 86

Aircraft, surface, upper air and satellite measurements have been used to observe the evolution and growth of the convective Marine Atmospheric Boundary Layer (MABL) offshore of North Carolina in close proximity to the Gulf Stream, during the intense cold air outbreak of 28 January 1986 and the moderate event of 12 February 1986, as part of the Genesis of Atlantic Lows Experiment (GALE). Air mass modification processes, driven primarily by the ocean-atmosphere exchanges of surface turbulent sensible and latent heat fluxes, caused the overlying air mass to warm and moisten as it advected over the warmer waters of the eastern United States continental shelf. Maximum observed near-surface total heat fluxes were 1045 and 811 W·m^–2 over the core of the Gulf Stream, for 28 January and 12 February 1986, respectively. The observed changes in the overlying air mass occurred almost instantaneously as the ambient flow traversed different underlying SST conditions.The turbulent structure showed a buoyancy-dominated MABL below approximately 0.8z/h. However, shear was also observed to be an important production term above 0.8z/h and below 0.1z/h for the 28 January 1986 event. Dissipation of turbulent kinetic energy was the dominant destruction term in the budgets, but vertical transport of energy was a strong contributor below 0.5z/h, above which this term became a source of turbulent energy. Additionally, the normalized standard deviations of the horizontal velocity components showed a near-equal contribution to the turbulence, while the vertical velocity components displayed the characteristic mid-layer maximum profile observed for a convective, well-mixed boundary layer.     

A numerical study of the convective boundary layer

Computations of the buoyantly unstable Ekman layer are performed at low Reynolds number. The turbulent fields are obtained directly by solving the three-dimensional time-dependent Navier-Stokes equations (using the Boussinesq approximation to account for buoyancy effects), and no turbulence model is needed. Two levels of heating are considered, one quite vigorous, the other more moderate. Statistics for the vigorously heated case are found to agree reasonably well with laboratory, field, and large-eddy simulation results, when Deardorff's mixed-layer scaling is used. No indication of large-scale longitudinal roll cells is found in this convection-dominated flow, for which the inversion height to Obukhov length scale ratio –z _i /L _*=26. However, when heating is more moderate (so that –z _i /L _*=2), evidence of coherent rolls is present. About 10% of the total turbulent kinetic energy and turbulent heat flux, and 20% of the Reynolds shear stress, are estimated to be a direct consequence of the observed cells.     

DMS oxidation and turbulent transport in the marine boundary layer: A numerical study

DMS oxidation in the marine boundary layer has been simulated with a mesoscale meteorological model including detailed physical parameterizations. The impact of vertical turbulent transport on the DMS and SO₂ diurnal cycles with and without in-cloud SO₂ oxidation has been studied in a one-dimensional version of the model and compared to results obtained with a zero-dimensional box model. Initialisation has been done using balanced values issued from the imposed sea-air fluxes, dry deposition fluxes and chemical source/sink terms. Particular emphasis has been put on the important role played by evolving vertical mixing in the marine boundary layer.     

Observations of atmospheric boundary layer evolution above the Gulf Stream frontal zone

Ship borne measurements of atmospheric boundary layer (ABL) parameters, sea-surface temperature and radar signals are analyzed to reveal the effects of the ABL transformation above the Gulf Stream temperature frontal zone. It was found that local changes in vertical gradients of wind speed and air temperature are well correlated with sub-mesoscale (~ 10 km) sea surface temperature variations. These effects are accompanied by appropriate variations in surface wind stresses that were identified from microwave backscatter.For steady atmospheric conditions the same effects were observed on spatial scales of 100 km, demonstrating positive radar signal contrast of the Gulf Stream warm waters with respect to surrounding Sargasso sea and shelf water areas. A simplified model of the ABL, accounting for an effect of spatial inhomogeneity by introducing an internal boundary layer, is used to analyze field observations. The model is able to reproduce both sub-mesoscale and mesoscale ABL evolution.     

A numerical study of the marine stratocumulus cloud layer

A one-dimensional grid-level model including longwave radiative transfer and a level-4 second-order turbulent transfer closure which contains prognostic equations for turbulent quantities, is used to study the physics and dynamics of inversion-capped marine stratocumulus clouds.A set of numerical experiments had been performed to examined the role of sea surface temperature, large-scale vertical velocity, wind speed, and vertical wind shear in the formation and the structure of low-level clouds. For a given sea surface and geostrophic wind speed, stratocumulus clouds can grow higher with smaller large-scale subsidence as less dry air entrains into the cloud. Clouds grow higher with higher sea surface temperature for a given geostrophic wind speed and large-scale subsidence as a result of enhanced moist convection. In high wind speeds, the entire cloud deck is lifted up because of larger surface energy flux. In the budget studies of the turbulent kinetic energy (TKE), the buoyancy term is a major source term when the wind speed and the vertical shear are small across the inversion top. When the wind speed and the vertical wind shear across the inversion top become large, the mixed layer is decoupled into a cloud and a subcloud layer. In the TKE budget studies, the shear generation term becomes an important term in the budgets of the TKE and the variance of vertical velocity.     

Numerical simulation of the development of the convective boundary layer during a cold air outbreak

During cold air outbreaks, cold and stably stratified air masses are advected from land or ice surfaces over a warmer sea surface. Due to the heating from below, a convective boundary layer develops. For small fetches, convection is organized in the form of horizontal roll vortices, which at greater distances join in a zone with open or closed cells. The formation of the convective boundary layer, and the associated roll vortices, is simulated with a numerical model and results are compared with observations obtained during the MASEX experiment off the east coast of the United States. To validate the model, a comparison with a one-dimensional mixed-layer model is also made, with special attention given to the exact representation of the observed initial and boundary conditions. Comparisons between model results and observations show good qualitative and quantitative correspondence in mean temperature and heat flux profiles respectively at different distances from the coast. Maximum values of vertical velocity are well reproduced. Turbulent kinetic energy is found to be concentrated in the small updraft regions of the rolls, which is in accordance with observations from the MASEX-experiment.     

Numerical simulation of January 28 cold air outbreak during GALE part II: The mesoscale circulation and marine boundary layer

A two-dimensional (2-D) mesoscale numerical model is applied to simulate the January 28 cold-air outbreak over the Gulf Stream region during the Intensive Observation Period-2 (IOP-2) of the 1986 Genesis of Atlantic Lows Experiment (GALE). The model utilizes a turbulence closure which involves the turbulent kinetic energy (TKE) and dissipation () equations and combines the level 2.5 formulations of Mellor and Yamada (1982) for better determination of the eddy Prandtl number.The modeled marine boundary layer (MBL) is in good agreement with the observations (Wayland and Raman, 1989) showing a low-level jet west of the Gulf Stream warm core and a constrained boundary layer due to the middle-level (2–4.5 km) stable layer. The MBL-induced single cloud and rain band first appears east of the Gulf Stream boundary, and then moves offshore at the speed of the circulation front. The front, however, moves slightly slower than the ambient flow. Removal of the tropopause does not influence the low-level circulation and the movement of the front. The speed of the front is slightly larger in the baroclinic downshear flow than in the barotropic flow. The results also indicate that the observed high cloud streets propagating downwind of the Gulf Stream may be related to upper-level baroclinic lee waves triggered by an elevated density mountain. The density mountain waves, however, become evanescent as the baroclinity (which gives a larger Scorer parameter) is removed.The modeled 2-D circulation systems are found to be sensitive to differing eddy Prandtl numbers, in contrast to the 1-D model results presented in Part I. Sensitivities become increasingly important as the clouds begin to interact with the MBL. A constant eddy Prandtl number of unity produces a more slantwise convection compared to that by the level 2.5 case. Cloud development is stronger in slantwise convection than in upright convection. The fastest development of clouds can be explained in terms of the conditional symmetric instability (CSI), which begins as the MBL baroclinity becomes sufficiently large.     

A case study of warm air advection over a melting snow surface

When relatively warm, moist air moves over a snow surface, sensible heat and moisture are extracted from its lower layers and used to melt the snow. The depth of the cooled layer depends on horizontal wind speeds and the presence of high vertical wind shear. The mechanism for air mass modification appears to be turbulent mixing.     

夜间城市边界层发展的数值研究

本文使用非定常非线性二维数值模式,研究夜间气流流经城市热岛上空引起的风场、温度场和垂直涡旋扩散参数的调整以及城市热岛环流的发展演变.     

A numerical study of cloud streets in the planetary boundary layer

A two-dimensional numerical model is used to study the influence of small non-precipitating clouds on horizontal roll vortices in the planetary boundary layer. The model explicitly represents the large-scale two-dimensional motions whilst small-scale eddies are parameterized by a buoyancy dependent mixing-length hypothesis. It is applied to conditions corresponding to an observed case of cloud street formation.     

A numerical study of daily transitions in the convective boundary layer

We examine daily (morning–afternoon) transitions in the atmospheric boundary layer based on large-eddy simulations. Under consideration are the effects of the stratification at the top of the mixed layer and of the wind shear. The results describe the transitory behaviour of temperature and wind velocity, their second moments, the boundary-layer height Z _m (defined by the maximum of the potential temperature gradient) and its standard deviation σ _m , the mixed-layer height z _i (defined by the minimum of the potential temperature flux), entrainment velocity W _e, and the entrainment flux H _i . The entrainment flux and the entrainment velocity are found to lag slightly in time with respect to the surface temperature flux. The simulations imply that the atmospheric values of velocity variances, measured at various instants during the daytime, and normalized in terms of the actual convective scale w_*, are not expected to collapse to a single curve, but to produce a significant scatter of observational points. The measured values of the temperature variance, normalized in terms of the actual convective scale Θ_*, are expected to form a single curve in the mixed layer, and to exhibit a considerable scatter in the interfacial layer.     

A numerical study on the meso-scale pollutant dispersion over a sloped surface in the stable boundary layer

A numerical model is used to study the properties of pollutant dispersion over a large uniformly-sloped surface in the stable atmospheric boundary layer. By simulating the structure of boundary layer flow to improve the advective wind field and vertical eddy diffusivities included in the advection-diffusion equa-tion, this numerical model permits an estimation of the distribution of pollutant concentration for more real-istic atmospheric diffusion conditions.     

南宁平流雾并伴有毛毛雨天气的边界层特征

通过对2005-2011年间南宁出现的7次平流雾并同时伴有毛毛雨天气过程的温度、湿度和风场等边界层特征进行统计分析,得出在广西地面处于高压后部,锋面北退,500hPa为宽广的西南气流控制和850hPa强锋区的形势下,当逆温层下降到1200m以下,厚度400-800m,逆温层上下温差不大和水汽处在准饱和状态时有利于南宁出现雨雾并存的天气.     

Measurements of peroxyacetylnitrate in the marine boundary layer over the Atlantic

The atmospheric concentration of peroxyacetylnitrate (PAN) was measured during a cruise of the R.S. Polarstern from Bremerhaven (Germany) to Rio Grande do Sul (Brazil) in September/ October 1988. The measurements were made in-situ by a combination of electron capture gaschromatography with a cryogenic preconcentration step. The theoretical lower limit of detection (3) was 0.4 ppt. The mixing ratios of PAN varied by more than three orders of magnitude from 2000 ppt in the English Channel to less than 0.4 ppt south of the Azores (38° N). South of 35° N, PAN levels were below the detection limit, except at 30–31° S off the eastern coast of South America. Here, PAN mixing ratios of 10 to 100 ppt were detected in continentally influenced air masses. Detectable levels of PAN were mostly observed in air masses of continental or high northern origin. Changes in the wind directions were usually associated with substantial changes in the PAN mixing ratios.     

Structure of the marine boundary layer over north western Indian Ocean during 1983 summer monsoon

The spatial variability of the structure of the lower troposphere over the northwestern Indian Ocean for the period 12th July to 2nd September, 1983 has been studied using upper air data collected during the first scientific cruise of ORV Sagar Kanya.An analysis of thermodynamic structure and kinematics of the marine boundary layer for different zonal and meridional sections revealed the following features: (a) Temperature and humidity inversions were generally absent over the study area except over a few locations in the western region; (b) Large-scale subsidence was found over the central equatorial Indian Ocean; (c) The convective activity over the western Indian Ocean was found to be moderately suppressed as compared to the eastern region; (d) The zonal and meridional components of winds along the equator and 10° N zonal section exhibited a mirror-image-like distribution.     

A numerical study of the influence of an air temperature-inversion layer and a seawater density-jump layer on the structure of interacting boundary layers

The effects of an air-temperature inversion in the atmosphere and a seawater density jump in the ocean on the structure of the atmospheric and oceanic boundary layers are studied by use of a coupled model. The numerical model consists of a closed system of equations for velocities, turbulent kinetic energy, turbulent exchange coefficient, local turbulent length scale, and stratification expressions for both air and sea boundary layers. The effects of the temperature inversion and the density jump are incorporated into the equations of turbulent kinetic energy of the atmosphere and ocean by a parameterization. A series of numerical experiments was conducted to determine the effects of various strengths of the inversion layer and surface heat fluxes in the atmosphere and of the density-jump layer in the ocean on the structure of the interacting boundary layers.The numerical results show that the temperature inversion in the atmosphere and density jump in the ocean have strong influences on turbulent structure [especially on the turbulent exchange coefficient (TEC) and turbulent kinetic energy (TKE)] and on air-sea interaction characteristics. Maxima of TKE and TEC strongly decrease with increasing strength of the inversion layer, and they disappear for strong inversions in the atmosphere. Certain strengths (density differences between the upper and the lower layers) of the density-jump layer in the ocean (₂ 0.1 g/cm³) produce double maxima in TEC-profiles and TKE-profiles in the ocean. The magnitudes of air-sea interaction characteristics such as geostrophic drag coefficient, and surface drift current increase with increasing strength of the density-jump layer in the ocean. The density-jump layer plays the role of a barrier that limits vertical mixing in the ocean. The numerical results agree well with available observed data and accepted quantitive understanding of the influences of a temperature inversion layer and a density-jump layer on the interacting atmospheric and oceanic boundary layers.     

坡地斜压稳定大气边界层结构的数值研究

本文用数值模式较为全面地探讨了斜坡地形斜压稳定大气边界层的结构,并与平坦地形的结果进行了比较,给出了两者之间很有意思的差别.     

Factors influencing the marine boundary layer during a cold-air outbreak

The model for the cloud-topped marine boundary layer during a cold air outbreak developed by Stage and Businger (1981a) is used in conjunction with a test profile based on a fall outbreak episode over Lake Ontario to study factors influencing marine boundary-layer evolution. Sensitivity tests are done which show changes in layer evolution resulting from variation of wind speed, radiative sky temperature, water surface temperature, humidity of the shoreline sounding and divergence. The behavior of the layer in the presence of a region of cold-water upwelling near the shore is also investigated. It is found that the main effect of the upwelling region is to delay modification of the boundary-layer air.     

海面粗糙度可变时定常非中性大气边界层数值模式

本文结出一个海面粗糙度可变时定常非中性大气边界层数值模式,其中粗糙度是摩擦速度的函数,湍流交换系数是湍流热通量和风切变的函数,而热通量在海气温差已知时由近地层通量梯度关系求出。解得的风和若干边界层特征参数符合观测及理论考虑。因此,只要知道地转风、纬度及海气温差即可用本模式计算海面大气边界层中不同高度的风。