黄土高原复杂地形上边界层低空急流对近地层湍流的影响

利用中尺度气象数值模式（Weather Research and Forecasting Model,WRF）模拟风场,结合兰州大学半干旱气候与环境观测站（Semi-Arid Climate and Environment Observatory of Lanzhou University,SACOL）湍流观测资料,分析了黄土高原复杂地形上稳定边界层低空急流对近地层湍流活动的影响.黄土高原复杂地形上稳定边界层低空急流的形成与地形作用引发的局地环流有关.低空急流对近地层湍流活动有强烈影响,剪切作用使小尺度湍涡活动加剧,湍动能增大,同时非平稳运动被压制.低空急流发生时,观测数据有87.3%是弱稳定情形（梯度理查森数小于0.25）;而无低空急流时,对应时段的观测表明65.4%属于强稳定层结（梯度理查森数大于0.3）,非平稳运动造成湍流功率谱在低频端迅速增大.与无低空急流和弱低空急流情形相比,强低空急流发生时,近地层湍动能增大1倍,湍动能在垂直方向上的传递增大1个量级,且方向向下,约为-3 × 10^-3 m³·s^-3,湍流在上层产生并向下传递.     

Stable-boundary-layer regimes from the perspective of the low-level jet

This paper reviews results from two field studies of the nocturnal stable atmospheric boundary layer (SBL) over the Great Plains of the United States. Data from a scanning remote-sensing system, a High-Resolution Doppler Lidar (HRDL), provided measurements of mean and turbulent wind components at high spatial and temporal resolution through the lowest 500–1000 m of the atmosphere. This data set has allowed the characteristics of the low-level jet (LLJ) maximum (speed, height, direction) to be documented through entire nights. LLJs form after sunset and produce strong shear in the layer below the LLJ maximum or nose, which is a source of turbulence and mixing in the SBL. Simultaneous HRDL measurements of turbulence quantities related to turbulence kinetic energy (TKE) has allowed the turbulence in the subjet layer to be related to LLJ properties. Turbulence structure was found to be a function of the bulk stability of the subjet layer. For the strong-LLJ (> 15 m s⁻¹), weakly stable cases the strength of the turbulence is proportional to the strength of the LLJ. For these cases with nearly continuous turbulence in the subjet layer, low-level jet scaling, in which lengths are scaled by the LLJ height and velocity variables are scaled by the LLJ speed, was found to be appropriate. For the weak-wind (< 5 m s⁻¹ in the lowest 200 m), very stable boundary layer (vSBL), the boundary layer was found to be very shallow (sometimes < 10 m deep), and turbulent fluxes between the earth’s surface and the atmosphere were found to be essentially shut down. For more intermediate wind speeds and stabilities, the SBL shows varying degrees of intermittency due to various mechanisms, including shearinstability and other gravity waves, density currents, and other mesoscale disturbances.     

Nocturnal basin low-level jets: an integrated study

Low-level jets (LLJs) are a very common feature in the nocturnal stably stratified boundary layer. Many factors can intervene in their generation, linked basically to effects of baroclinity. A special kind of low-level jets is composed by the nocturnal katabatic and basin flows, generated over terrain slopes. A study of observed LLJs in the Duero Basin is shown here, combining observational data and modelling experiments. Normalized in respect to the maximum wind height, the dynamic characteristics of the jets are similar: a two-layer system, with a stably stratified layer below the jet maximum and a near neutral layer above, with a very stable layer separating them at the level of the wind maximum. There is vertical mixing above and below the jet, and the connection between these layers takes place occasionally in a very turbulent manner.     

On the generation,propagation and radiation of magneto–acoustic–gravity waves with application to stars

Abstract

The study of the mechanisms controlling the stratification in closed fluid regions is an important branch of geophysical fluid dynamics. Part of this subject can be handled with a simple linear model, consisting of a buoyancy layer at the non-horizontal boundaries of a container and an advective-diffusive interior coupled by volume continuity. The model is valid under the following conditions: firstly, the buoyancy-frequency characterizing the solution must be sufficiently large to give rise to a flow pattern of boundary layer type and, secondly, the non-horizontal walls must not have too large thermal conductivity.

The main purpose of the present paper is to summarise previous work done by the authors in this field and to present some consequences of their theory not previously discussed.

Three important cases are discussed; certain stationary solutions, the decay of a given stratification and the build up of a stratification in a homogeneous fluid. The experimental results concerning the afore-mentioned cases are presented.     

Convection Regime of the Development of an Impulse–Buoyant Transient Jet Submerged in a Homogeneous Fluid

An intergral model of a transient vertical impulse–buoyant jet is suggested. The model contains a universal equation describing the propagation of the upper boundary of the convection front depending on the strength of the point source of buoyancy and momentum. The convective regime of jet propagation is considered, which includes a class of self-similar solutions corresponding to the buoyancy sources, whose strength varies with time following power and exponential laws. The obtained numerical solutions are compared with available experimental data on the profiles of vertical velocity and buoyancy on the jet axis.     

Structure of numerically simulated katabatic and anabatic flows along steep slopes

Direct numerical simulation (DNS) is applied to investigate properties of katabatic and anabatic flows along thermally perturbed (in terms of surface buoyancy flux) sloping surfaces in the absence of rotation. Numerical experiments are conducted for homogeneous surface forcings over infinite planar slopes. The simulated flows are the turbulent analogs of the Prandtl (1942) one-dimensional laminar slope flow. The simulated flows achieve quasi-steady periodic regimes at large times, with turbulent fluctuations being modified by persistent low-frequency oscillatory motions with frequency equal to the product of the ambient buoyancy frequency and the sine of the slope angle. These oscillatory wave-type motions result from interactions between turbulence and ambient stable stratification despite the temporal constancy of the surface buoyant forcing. The structure of the mean-flow fields and turbulence statistics in simulated slope flows is analyzed. An integral dynamic similarity constraint for steady slope/wall flows forced by surface buoyancy flux is derived and quantitatively verified against the DNS data.     

Richardson number profiles in laboratory experiments applied to shallow seas

Abstract

A unified analysis is given of the critical conditions for the onset of stratification due to either a vertical or a horizontal buoyancy flux, with tidal or wind stirring.

The critical conditions for the onset of stratification with a horizontal buoyancy flux are found to be of the form of ratios of the tidal slope, or wind setup, to the equivalent surface slope due to the lateral density gradient. These ratios, which are easily determined from sea data, indicate that the profiles of critical flux Richardson Number, averaged over the stirring cycle, are similar to those inferred from the laboratory experiments of Hopfinger and Linden (1982) in which there is zero mean shear turbulence with a stabilising buoyancy flux, and also that the efficiency for the conversion of kinetic energy to potential energy for tidal stirring is similar to that for wind stirring.

The observed much greater efficiency for wind stirring, compared with tidal stirring with a vertical buoyancy flux, is also consistent with the existence of flux Richardson Number profiles in the sea similar to those occurring in the corresponding laboratory experiments. Using the solution of the turbulent kinetic energy equation for the water column, the relative importance of the production of turbulent kinetic energy, and its diffusion by turbulence are assessed, and the critical conditions for the onset of stratification with a vertical buoyancy flux are shown to reduce the classical Simpson—Hunter form.     

A comprehensive physical pattern of land-air dynamic and thermal structure on the Qinghai-Xizang Plateau

According to the boundary layer observations of three stations (Garze, Damxung and Qamdu) and relevant earth satellite, radiosonde and surface observations during the intensive observational period (IOP) of the second Tibetan (Qinghai-Xizang) Plateau Experiment of atmospheric science (TIPEX), the land-air physical process and dynamic model on the Tibetan Plateau were comprehensively analyzed in this study. The dynamic characteristics of boundary layer and the rules of turbulent motion on the plateau were illustrated. The characteristics of distributions of wind speed and direction with mutiple-layer structure and deep convective mixed layer on the plateau, the strong buoyancy effect in turbulent motion on the plateau on which the air density is obviously smaller than on the plain, and the Ekman spiral and its dynamic pump effect of the plateau deep boundary layer have been found. The local static distribution of water vapor and the horizontal advection of water vapor in the plateau boundary layer were studied. The abnomal thermodynamic structure on the plateau surface and boundary layer, including the plateau strong radiation phenomenon and strong heating source characteristics of the middle plateau, was also analyzed. The authors synthesized the above dynamic and thermodynamic structures of both surface and boundary layers on the plateau and posed the comprehensive physical model of the turbulence and convective mixture mechanism on the plateau boundary layer. The characteristics of formation, development and movement for convective cloud cluster over the plateau influencing floods in the Yangtze River area of China were studied. The conceptual model of dynamic and thermodynamic structures of turbulent motion and convective plume related to the frequent occurrence of "pop-corn-like" cloud system is given as well.     

On the theory of thin rotating jets: A quasi‐geostrophic time dependent model

The time‐dependent meandering in a thin baroclinic jet over bottom topography is discussed in the quasi‐geostrophic approximation. The motion of the axis of the jet is taken to be vertically coherent and the axis itself is defined as inextensible. The motion is examined from a frame of reference moving with the axis but fixed at an arbitrary longitude in terms of an open ocean spatial initial value problem. The velocities of the axis and of the jet are quasi‐geostrophic, and vorticity conservation for the first non‐geostrophic components constrains the evolution of the axis and gives a path equation. The spatial linearized stability problem is studied and the jet is found to be baroclinically unstable to path disturbances of sufficiently high frequency which amplify downstream. An exact solution is obtained to the nonlinear path equation over a flat bottom with no ß‐effect. The evolution of the path of these unstable meanders is such that the path closes itself and forms rings (at which point the analysis breaks down). It is proposed that the baroclinic jet processes studied here are fundamental to the dynamics of Gulf Stream meandering and isolated eddy production.     

Turbulent entrainment at an inversion

Three separate conditions controlling the structure of the inversion are discussed. The rate of change of boundary layer temperature and the vertical gradient of potential temperature in the overlying air set the maximum limit to possible rise rates. In transient conditions, the rise of the inversion layer is sometimes dominated by turbulent transport and erosion, however this cannot continue indefinitely. The horizontal advection above the inversion can also change the relative temperature increase at the inversion, and so may partly determine inversion history.In this paper the theory of the plume-connective boundary layer has been extended to allow the plumes to erode material from the overlying stable air, and so permit study of how the boundary layer could deepen by eroding upward into the capping inversion. Only clear-air processes without condensation are included in the model described here.When the mechanism of convection is examined, it is demonstrated that, because of energetic considerations, the rate of erosion, by buoyancy generated turbulence of the overlying air down across the inversion, will not usually determine the rate of rise of the inversion, except briefly. Other effects often dominate in determining the average erosion rate. Furthermore, because of the limitations set by the buoyancy of the entrained air, sharp inversions, exceeding about 1°C, cannot be formed without the action of latent heat from evaporating droplets. An intermediate buffer layer forms in the inversion zone whenever the erosion rate is high enough to produce a temperature increase above the boundary layer, sufficient to inhibit direct entrainment.     

Numerical simulation of 2D buoyant jets in ice-covered and temperature-stratified water

A two-dimensional (2D) unsteady simulation model is applied to the problem of a submerged warm water discharge into a stratified lake or reservoir with an ice cover. Numerical simulations and analyses are conducted to gain insight into large-scale convective recirculation and flow processes in a cold waterbody induced by a buoyant jet. Jet behaviors under various discharge temperatures are captured by directly modeling flow and thermal fields. Flow structures and processes are described by the simulated spatial and temporal distributions of velocity and temperature in various regions: deflection, recirculation, attachment, and impingement. Some peculiar hydrothermal and dynamic features, e.g. reversal of buoyancy due to the dilution of a warm jet by entraining cold ambient water, are identified and examined. Simulation results show that buoyancy is the most important factor controlling jet behavior and mixing processes. The inflow boundary is treated as a liquid wall from which the jet is offset. Similarity and difference in effects of boundaries perpendicular and parallel to flow, and of buoyancy on jet attachment and impingement, are discussed. Symmetric flow configuration is used to de-emphasize the Coanda effect caused by offset.     

Comparison of methods for modelling the behaviour of bubbles produced by marine seismic airguns

Three models for the dynamics of seismic airgun‐generated bubbles and their associated far‐field signals are developed and compared with geophysical data. The first model of an airgun‐generated bubble uses a spherical approximation, the second is an approximate Lagrangian model which allows for small deformations from a spherical shape, whilst the final model is an axisymmetric boundary‐integral method which permits the bubble to evolve into highly non‐spherical geometries. The boundary‐integral method also allows both geometric interference and strong dynamic interactions in multi‐bubble studies. When comparing the spherical model to experimental data there are three apparent, significant differences: the magnitude of the primary pressure peak, which is greater in the model; the subsequent decay of the pressure peaks and motion – the experimental data demonstrating greater decay and a slower rise rate; and the frequency of oscillation, which is slower in the experimental data. It is believed that the first discrepancy is due to the initial stages of expansion where the compressed air is forced to sparge through the airgun ports. The other differences indicate that there is some other energy‐loss mechanism which is not accounted for in the spherical bubble model. Non‐spherical bubble behaviour is investigated through the use of two different deformable many‐bubble codes and their predictions are compared with the spherical model and experimental data. The Lagrangian model predicts the formation of a buoyancy‐driven liquid jet on the first collapse of a typical airgun bubble; however, the model breaks down when the bubble becomes significantly deformed, due to a low‐order spherical‐harmonic approximation for the potential. The axisymmetric boundary‐integral code models the jet shape accurately and it is found that these bubbles evolve to toroidal geometries when the jet impacts on the opposite surface of the bubble. This highly non‐spherical behaviour is readily observed on high‐speed films of airgun bubbles, and is one key source of energy loss; it damps the pulsations of the bubble and slows its rise speed. Inter‐bubble interactions are investigated using the two deformable bubble models, and the predictions are compared to field data. It was found that as the bubbles approach each other, their periods of oscillation increase in accordance with observations, and jets are formed in the direction of motion upon collapse.     

Some laboratory experiments with a rotating,differentially heated fluid annulus

Summary With a view to modelling more realistically certain large-scale meteorological and oceanographic flows, some experiments are described in which temperature and velocity fields are measured in a rotating, differentially heated fluid annulus, and their dependence upon the imposed boundary conditions is investigated. It is demonstrated that with suitable construction of the annulus walls the strength of the zonal baroclinic motion may be chosen independently of the basic density stratification. Most of the measurements described are for the symmetric flow regime. In the experiments it is found that certain aspects of the observed temperature and flow fields — in particular the basic stratification and the slope of the isotherms — agree reasonably well with theoretical estimates.     

Vortex pairs and power station cooling water

Cooling water discharged from power stations in the U.K. is frequently released from an outlet in an estuary or the sea. The warm water forms a thermal plume which is slightly buoyant and which spreads horizontally over the water surface while mixing vertically downwards with the cooler ambient water. In this paper, the possibility of vortex pair production at the cooling water outlet is considered as a mechanism contributing to this spreading of the warm water.The motion of a vortex pair contained between two rigid plane boundaries is an idealization of the flow between the water surface and the sea bed. The resulting motion is calculated from potential theory and viscous effects are neglected. The problem of deciding what strength to assign the vortices is discussed and specific consideration of shear and buoyancy at the outlet is detailed. It is observed that bifurcation of the vortex pair is determined by the initial position of the vortices and is unlikely to occur in conditions relevant to U.K. power station discharges.It is calculated that, in the absence of turbulence, the motion of such vortex pairs would result in horizontal spreading of the warm water which is greater than that observed at site surveys. It is concluded that turbulence in the ambient receiving water is sufficient to destroy vortices produced by the discharge during the early stages of the plume development.     

On the magnitude of upwelling fluxes in shelf-break canyons

A hydrodynamic model is employed to derive the magnitude of on-shelf fluxes through a shelf-break canyon for a wide range of canyon sizes and ambient oceanic conditions. Predicted canyon-upwelling fluxes are of the order of 0.05–0.1 Sv (1 Sv=1 million m³/s), being several orders of magnitude greater than upslope fluxes in the bottom Ekman layer on the ambient continental slope. On the basis of ∼150 simulations conducted, a bulk formula of upwelling flux in a submarine canyon is derived. For typical conditions, the upwelling flux varies quadratically with forcing strength (speed of incident flow), linearly with canyon depth, and is inversely proportional to the buoyancy frequency of the density stratification inside the canyon. Other parameters such as density stratification above shelf-break depth and bottom friction are found to have minor influences on the resultant canyon-upwelling flux.     

A dynamo model with double diffusive convection for Mercury's core

A recent dynamo model for Mercury assumes that the upper part of the planet's fluid core is thermally stably stratified because the temperature gradient at the core–mantle boundary is subadiabatic. Vigorous convection driven by a superadiabatic temperature gradient at the boundary of a growing solid inner core and by the associated release of light constituents takes place in a deep sub-layer and powers a dynamo. These models have been successful at explaining the observed weak global magnetic field at Mercury's surface. They have been based on the concept of codensity, which combines thermal and compositional sources of buoyancy into a single variable by assuming the same diffusivity for both components. Actual diffusivities in planetary cores differ by a large factor. To overcome the limitation of the codensity model, we solve two separate transport equations with different diffusivities in a double diffusive dynamo model for Mercury. When temperature and composition contribute comparable amounts to the buoyancy force, we find significant differences to the codensity model. In the double diffusive case convection penetrates the upper layer with a net stable density stratification in the form of finger convection. Compared to the codensity model, this enhances the poloidal magnetic field in the nominally stable layer and outside the core, where it becomes too strong compared to observation. Intense azimuthal flow in the stable layer generates a strong axisymmetric toroidal field. We find in double diffusive models a surface magnetic field of the observed strength when compositional buoyancy plays an inferior role for driving the dynamo, which is the case when the sulphur concentration in Mercury's core is only a fraction of a percent.     

A three-layered spectral model with application to wind-induced motion in the presence of stratification and a bottom slope

The main features of a new three-dimensional spectral model for the motion of a stratified sea are described. The model is applied to determine wind-induced elevations and currents over an idealised shelf-edge slope during the summer season of thermal stratification. A wind stress field spanning the slope is applied in the form of a rectangular pulse through time, the direction of the stress being parallel to the shelf edge. It is shown that strong vertical motion occurs near the top of the slope in the wind region; inertial currents are excited. Outside the wind region, long-period waves trapped over the slope and propagating along it may be identified.     

Numerical simulation of the influence of aerosol radiation effect on urban boundary layer

With the intensification of pollution and urbanization, the aerosol radiation effect continues to play an important role in the urban boundary layer. In this paper, a winter pollution process in Beijing has been taken as an example, and a new aerosol vertical profile in the radiative parameterization scheme within the Weather Forecast Research and Forecasting (WRF) model has been updated to study the effect of aerosols on radiation and the boundary layer. Furthermore, the interactions among aerosols, urbanization, and planetary boundary layer (PBL) meteorology were discussed through a series of numerical experiments. The results show the following: (1) The optimization improves the performance of the model in simulating the distribution features of air temperature, humidity, and wind in Beijing. (2) The aerosols reduce the surface temperature by reducing solar radiation and increasing the temperature in the upper layer by absorbing or backscattering solar radiation. The changes in the PBL temperature lead to more stable atmospheric stratification, reducing the energy transfer from the surface and the height of the boundary layer. (3) With the increase in the aerosol optical depth, the atmospheric stratification most likely becomes stable over rural areas, most likely becomes stable over suburb areas, and has great difficultly becoming stable over urban areas. Aerosol radiative forcing, underlying urban surfaces, and the interaction between them are the main factors that affect the changes in the meteorological elements in the PBL.     

A Model Study of the Strong and Weak Wind,Stably Stratified Nocturnal Boundary Layer: Influence of Gentle Slopes

With the exception of intermittency and waves, a brief review of the observed and modeled mean structure of the nocturnal boundary layer (NBL) is presented. The effect of gentle slopes on strong and weak wind NBL was investigated here using a one-dimensional model, with a simple correction term to account for the slope effects, identical to the one used by Brost and Wyngaard (1978). The study indicates that the wind profiles, temperature profiles and surface layer turbulence characteristics are extremely sensitive to the imposed geostrophic wind when small slopes are present especially for light winds. This is due to the complex interaction between the buoyancy driven slope flow and the imposed geostrophic wind that in turn influence the shear generation of turbulence. Finally, the current issues in the modeling of weak wind boundary layer are discussed.