干湿地表的湍流特征及其对深对流影响的大涡模拟

利用大涡模式模拟了对流边界层结构演变以及深对流触发过程.通过改变鲍恩比的敏感性试验研究不同大气初始状况下湿润和干旱下垫面湍流特征及其对深对流触发过程的影响.结果表明:干旱下垫面的混合层干而暖,厚度较大;湿润下垫面相反.由于地表感热通量对热力湍流形成的作用更大,干旱下垫面上湍流混合和夹卷作用更强,使得水汽和相当位温在边界...     

Observations and Large-Eddy Simulations of Entrainment in the Sheared Sahelian Boundary Layer

At the top of the planetary boundary layer, the entrainment of air, which incorporates dry and warm air from the free troposphere into the boundary layer, is a key process for exchanges with the free troposphere since it controls the growth of the boundary layer. Here, we focus on the semi-arid boundary layer where the entrainment process is analyzed using aircraft observations collected during the African Monsoon Multidisciplinary Analysis experiment and large-eddy simulations. The role of the entrainment is specifically enhanced in this region where very large gradients at the planetary boundary-layer top can be found due to the presence of the moist, cold monsoon flow on which the dry, warm Harmattan flow is superimposed. A first large-eddy simulation is designed based on aircraft observations of 5 June 2006 during the transition period between dry conditions and the active monsoon phase. The simulation reproduces the boundary-layer development and dynamics observed on this day. From this specific case, sensitivity tests are carried out to cover a range of conditions observed during seven other flights made in the same transition period in order to describe the entrainment processes in detail. The combination of large-eddy simulations and observations allows us to test the parametrization of entrainment in a mixed-layer model with zero-order and first-order approximations for the entrainment zone. The latter representation of the entrainment zone gives a better fit with the conditions encountered in the Sahelian boundary layer during the transition period because large entrainment thicknesses are observed. The sensitivity study also provides an opportunity to highlight the contribution of shear stress and scalar jumps at the top of the boundary layer in the entrainment process, and to test a relevant parametrization published in the recent literature for a mixed-layer model.     

An evaluation of neutral and convective planetary boundary-layer parameterizations relative to large eddy simulations

This paper compares a number of one-dimensional closure models for the planetary boundary layer (PBL) that are currently in use in large-scale atmospheric models. Using the results of a large-eddy simulation (LES) model as the standard of comparison, the PBL models are evaluated over a range of stratifications from free convective to neutral and a range of surface shear stresses. Capping inversion strengths for the convective cases range from weakly to strongly capped. Six prototypical PBL models are evaluated in this study, which focuses on the accuracy of the boundary-layer fluxes of momentum, heat, and two passive scalars. One scalar mimics humidity and the other is a top-down scalar entrained into the boundary layer from above. A set of measures based on the layer-averaged differences of these fluxes from the LES solutions is developed. In addition to the methodological framework and suite of LES solutions, the main result of the evaluation is the recognition that all of the examined PBL parameterizations have difficulty reproducing the entrainment at the top of the PBL, as given by the LES, in most parameter regimes. Some of the PBL models are relatively accurate in their entrainment flux in a subset of parameter regimes. The sensitivity of the PBL models to vertical resolution is explored, and substantive differences are observed in the performance of the PBL models, relative to LES, at low resolution typical of large scale atmospheric models.     

A mixed-layer model for regional evaporation

A simple mixed-layer model is developed to describe evaporation into a convective planetary boundary layer (PBL). The model comprises volume budget equations for temperature and humidity, equations to describe transport through the surface layer which is treated as part of the lower boundary, and equations to describe entrainment at the top of the PBL. The ground surface is modelled as a canopy resistance. The model was integrated with canopy resistance, surface-layer resistance and available energy, (R _n – G), input as given functions of time, and the simulated PBL was allowed to grow into an atmosphere with known temperature and humidity profiles.Two variants of the mixed-layer model were tested using data from the KNMI tower site at Cabauw in the Netherlands. These variants differed only in the formulation of entrainment: one used a formulation developed by Driedonks (1982) while the other was a simpler formulation. Simulated evaporation agreed very well with observations irrespective of which entrainment formulation was used, despite discrepancies between simulated and observed PBL height growth which were sometimes quite large for the simpler formulation. Sensitivity analysis of the model confirms that good PBL height-growth predictions are not always a prerequisite for good evaporation predictions.     

Evidence of a Mixed-Layer Dynamics Contribution to the Entrainment Process

The internal thermal boundary layer developing over the Mediterranean during a cold-air outbreak associated with a Tramontane event has been studied by means of airborne lidar, in situ sensors, and a modelling approach that consisted of nesting the University of Washington (UW) planetary boundary-layer (PBL) model in an advective zero-order jump model. This approach bypasses some of the deficiencies associated with each model: the absence of the dynamics in the mixed layer for the zero-order jump model and the lack of an inversion at the PBL top for the UW PBL model. Particular attention is given to the parameterization of the entrainment flux at the PBL top. Values of the entrainment closure parameter derived with the model when matching PBL structure observations are much lower than those derived with standard zero-order jump models. They also are in good agreement with values measured in different meteorological situations by other studies. This improvement is a result of the introduction of turbulent kinetic energy production in the mixed layer.     

Observations of Diurnal Cycles Over a West African Meridional Transect: Pre-Monsoon and Full-Monsoon Seasons

We document and characterize the climatology of the diurnal cycles encountered along a West African transect during the pre-monsoon and full-monsoon periods. The meridional gradient in low-level properties is fundamental for the monsoon dynamics and here, for the first time, it is studied based on a large set of observations from the African Monsoon Multidisciplinary Analysis (AMMA) campaign. A detailed analysis of surface energy budget, boundary-layer structures and cloud occurrence is carried out to investigate the diurnal cycles of the low levels. A relatively weak meridional gradient of net radiation is observed during the pre-monsoon period, and a large gradient in sensible heat flux is found over the transect with values increasing from south to north. This, as well as the boundary-layer structures, partly explains the large contrasts in the diurnal amplitude of potential temperature and specific humidity along the transect. During the monsoon period, the atmospheric regimes drastically change involving strong interactions between the surface, atmosphere and clouds. The maximum in net radiation is shifted northwards, towards the Sahel, which potentially has a significant impact on the monsoon circulation. The sensible heat flux is considerably reduced and the diurnal amplitude is strongly damped, while the daytime boundary-layer growth decreases significantly in the Sahel related to changes in the balance of boundary-layer processes. These results highlight the contrasted diurnal cycle regimes encountered over West Africa under dry, moist and wet conditions. They provide observationally-based diagnostics to investigate the ability of models to handle the representation of the diurnal cycle over land.     

Relation between equilbrium evaporation and the saturation pressure budget

The limiting surface Bowen ratios are calculated which maintain mixed-layer saturation pressure and relative humidity, both with and without boundary-layer entrainment. The equations are formally the same as those of Culf (1994) for equilibrium evaporation, but differ numerically because the coefficients are calculated at the mixed-layer saturation temperature (at the lifting condensation level) rather than at the surface temperature. The diurnal cycle over land is used to illustrate the value of these constraints.     

On the Effect of Surface Heat-Flux Heterogeneities on the Mixed-Layer-Top Entrainment

We used a set of large-eddy simulations to investigate the effect of one-dimensional stripe-like surface heat-flux heterogeneities on mixed-layer top entrainment. The profiles of sensible heat flux and the temporal evolution of the boundary-layer depth revealed decreased entrainment for small heat-flux amplitudes and increased entrainment for large heat-flux amplitudes, compared to the homogeneously-heated mixed layer. For large heat-flux amplitudes the largest entrainment was observed for patch sizes in the order of the boundary-layer depth, while for significantly smaller or larger patch sizes entrainment was similar as in the homogeneous case. In order to understand the underlying physics of this impact, a new approach was developed to infer local information on entrainment by means of the local flux divergence. We found an entrainment maximum over the centre of the stronger heated surface patch, where thermal energy is accumulated by the secondary circulation (SC) that was induced by the surface heterogeneity. Furthermore, we observed an entrainment maximum over the less heated patch as well, which we suppose is to be linked to the SC-induced horizontal flow convergence at the top of the convective boundary layer (CBL). For small heat-flux amplitudes a counteracting effect dominates that decreases entrainment, which we suppose is the horizontal advection of cold air in the lower, and warm air in the upper, CBL by the SC, stabilizing the CBL and thus weakening thermal convection. Moreover, we found that a mean wind can reduce the heterogeneity-induced impact on entrainment. If the flow is aligned perpendicular to the border between the differentially-heated patches, the SC and thus its impact on entrainment vanishes due to increased horizontal mixing, even for moderate wind speeds. However, if the flow is directed parallel to the border between the differentially-heated patches, the SC and thus its impact on entrainment persists.     

On the onset of entrainment instability over the ocean

The marine atmospheric boundary layer is characterized by cool temperatures and high humidity. Clouds are observed over most of the oceans. It is generally accepted that these overcast cloud decks break up into scattered fragments due to cloud-top entrainment instability. That is, if the air above the boundary layer is sufficiently cool and dry relative to cloud top, the buoyancy flux will be directed upwards and entrainment can occur freely.A boundary-layer model is used to test the sensitivity of the model atmosphere to the various processes which promote the onset of cloud-top entrainment instability. It is found that the transition from a solid cloud deck to scattered cumulus clouds depends on a rate process. The cloud cover is sensitive to mesoscale variations in sea surface temperature only if the cloud-top inversion is sufficiently weak.     

Development of a Single-Column Model in RegCM4 and Its Preliminary Application for Evaluating PBL Schemes in Simulating the Dry Convection Boundary Layer

A single-column model (SCM) is developed in the regional climate model RegCM4. The evolution of a dry convection boundary layer (DCBL) is used to evaluate this SCM. Moreover, four planetary boundary layer (PBL) schemes, namely the Holtslag-Boville scheme (HB), Yonsei University scheme (YSU), and two University of Washington schemes (UW01, Grenier-Bretherton-McCaa scheme and UW09, Bretherton-Park scheme), are compared by using the SCM approach. A large-eddy simulation (LES) of the DCBL is performed as a benchmark to examine how well a PBL parameterization scheme reproduces the LES results, and several diagnostic outputs are compared to evaluate the schemes. The results show that the SCM is proper constructed. In general, with the DCBL case, the YSU scheme performs best for reproducing the LES results, which include well-mixed features and vertical sensible heat fluxes; the simulated wind speed, turbulent kinetic energy, entrainment flux, and height of the entrainment zone are all underestimated in the UW09; the UW01 has all those biases of the UW09 but larger, and the simulated potential temperature is not well mixed; the HB is the least skillful scheme, by which the PBL height, entrainment flux, height of the entrainment zone, and the vertical gradients within the mixed layer are all overestimated, and a inversion layer near the top of the surface layer is wrongly simulated.Although more cases and further testing are required, these simulations show encouraging results towards the use of this SCM framework for evaluating the simulated physical processes by the RegCM4.     

Improving Non-local Parameterization of the Convective Boundary Layer

We propose improvements in the “non-local” parameterization scheme of the convective boundary layer. The countergradient terms for components of the momentum fluxes are introduced in a form analogous to those for other scalars. The scheme also includes explicit expressions for entrainment fluxes of momentum, temperature, and humidity. A simplified procedure for calculating the boundary-layer height is proposed, consisting of two steps: the evaluation of the convection level, followed by the assessment of the depth of the interfacial layer.     

Equatorial Indian Ocean evaporation estimates from operational meteorological satellites and some inferences in the context of monsoon onset and activity

The estimation of evaporation from the sea surface is not yet achieved adequately by remote sensing techniques, in general. However, for approximate averaged estimates over moderate space and time scales over a specific tropical region, e.g., weekly values over the Indian Ocean as needed in monsoon moisture diagnosis, it may be possible to extrapolate satellite wind and humidity data to the ocean surface and then use bulk aerodynamic parameterization for estimating evaporation. In the present investigation, GOES low-level cloud winds and TIROS-N moisture profiles over the Indian Ocean are extrapolated to the ocean surface. The planetary boundary layer (PBL) wind shear is obtained over different sub-regions and periods during the monsoon season, by reference to objectively analysed fields. These shear values are applied to GOES satellite winds to obtain sea-surface winds. The humidity extrapolation was based on (i) an exponential fit for water vapour density and (ii) a vertical distribution of relative humidity approximately proportional to atmospheric pressure. The exchange coefficient is varied slightly depending on wind speed and boundary-layer stability inferred approximately from TIROS-N sea surface temperatures and temperature profiles. The evaporation estimate as based on these satellite parameters is assessed by comparison with ships' surface observations. Sensible heat exchange is also estimated and assessed. Some inferences based on these estimates are also presented, in relation to monsoon onset and activity.     

Detection of Entrainment Influences on Surface-Layer Measurements and Extension of Monin–Obukhov Similarity Theory

We present a method to detect influences of boundary-layer processes on surface-layer measurements, using statistics and spectra of surface-layer variables only. We validated our detection method with boundary-layer measurements. Furthermore, we confirm that Monin–Obukhov similarity functions fit well to temperature-variance data obtained at two different homogeneous surfaces. However, we found that humidity variance measurements deviate from the universal functions above one of the two studied surfaces for days on which entrained air reached the surface layer. These results confirm that Monin–Obukhov similarity theory should be used with care in the analysis of surface-layer data. Finally, we propose the use of an extra term in flux-variance relations that depends on the entrainment ratio for humidity and on the boundary-layer height. If boundary-layer measurements are not available, we show how the entrainment ratio for humidity can be approximated from the skewness of the humidity distribution.     

Responses of global monsoon and seasonal cycle of precipitation to precession and obliquity forcing

In this paper, the response of global monsoon to changes in orbital forcing is investigated using a coupled atmosphere–ocean general circulation model with an emphasis on relative roles of precession and obliquity changes. When precession decreases, there are inter-hemispheric asymmetric responses in monsoonal precipitation, featuring a significant increase over most parts of the Northern Hemisphere (NH) monsoon regions and a decrease over the Southern Hemisphere (SH) monsoon regions. In contrast, when obliquity increases, global monsoon is enhanced except for the American monsoon. Dynamic effects (caused by changes in winds with humidity unchanged) dominate the monsoonal precipitation response to both precession and obliquity forcing, while thermodynamic effects (caused by changes in humidity with winds unchanged) is related to the northward extension of the North African summer monsoon. During minimum precession, the seasonal cycle of tropical precipitation is advanced with respect to the maximum precession. The rainfall increase in the transitional season (April-June in the NH and October-December in the SH) is dominated by the dynamic component. From an energetics perspective, the southward (northward) cross-equatorial energy transport during April-June (October-December) corresponds to a northward (southward) shift of tropical precipitation, which results in a seasonal advance in the migration of tropical precipitation. Nonetheless, there is no significant change in the seasonal cycle in response to obliquity forcing.

     

Atmospheric Boundary-Layer Dynamics with Constant Bowen Ratio

Motivated by the observation that the diurnal evolution of sensible and latent heat fluxes tends to maintain a constant Bowen ratio, we derive approximate solutions of the ordinary differential equations of a simplified atmospheric boundary-layer (ABL) model. Neglecting the early morning transition, the potential temperature and specific humidity of the mixed layer are found to be linearly related to the ABL height. Similar behaviour is followed by the inversion strengths of temperature and humidity at the top of the ABL. The potential temperature of the mixed layer depends on the entrainment parameter and the free-atmosphere temperature lapse rate, while the specific humidity also depends on the free-atmosphere humidity lapse rate and the Bowen ratio. The temporal dynamics appear only implicitly in the evolution of the height of the boundary layer, which in turn depends on the time-integrated surface sensible heat flux. Studying the limiting behaviour of the Bowen ratio for very low and very large values of net available energy, we also show how the tendency to maintain constant Bowen ratio during midday hours stems from its relative insensitivity to the atmospheric conditions for large values of net available energy. The analytical expression for the diurnal evolution of the ABL obtained with constant Bowen ratio is simple and provides a benchmark for the results of more complex models.     

Frequency of Boundary-Layer-Top Fluctuations in Convective and Stable Conditions Using Laser Remote Sensing

The planetary boundary-layer (PBL) height is determined with high temporal and altitude resolution from lidar backscatter profiles. Then, the frequencies of daytime thermal updrafts and downdrafts and of nighttime gravity waves are obtained applying a fast Fourier transform on the temporal fluctuation of the PBL height. The principal frequency components of each spectrum are related to the dominant processes occurring at the daytime and nighttime PBL top. Two groups of cases are selected for the study: one group combines daytime cases, measured in weak horizontal wind conditions and dominated by convection. The cases show higher updraft and downdraft frequencies for the shallow, convective boundary layer and lower frequencies for a deep PBL. For cases characterized by strong horizontal winds, the frequencies directly depend on the wind speed. The temporal variation of the PBL height is determined also in the likely presence of lee waves. For nighttime cases, the main frequency components in the spectra do not show a real correlation with the nocturnal PBL height. Altitude fluctuations of the top of the nocturnal boundary layer are observed even though the boundary layer is statically stable. These oscillations are associated with the wind shear effect and with buoyancy waves at the PBL top.     

A bulk mass flux convection scheme for climate model: description and moisture sensitivity

A convection scheme for climate model is developed based on Tiedtke’s (Mon Weather Rev 117:1779–1800, 1989) bulk mass flux framework and is evaluated with observational data and cloud resolving model simulation data. The main differences between the present parameterization and Tiedtke’s parameterization are the convection trigger, fractional entrainment and detrainment rate formulations, and closure method. Convection is triggered if the vertical velocity of a rising parcel is positive at the level at which the parcel is saturated. The fractional entrainment rate depends on the vertical velocity and buoyancy of the parcel as well as the environmental relative humidity. For the fractional detrainment rate, a linear decrease in the updraft mass flux above maximum buoyancy level is assumed. In the closure method, the cloud base mass flux is determined by considering both cloud layer instability and subcloud layer turbulent kinetic energy as controlling factors in the strength of the convection. The convection scheme is examined in a single column framework as well as using a general circulation model. The present bulk mass flux (BMF) scheme is compared with a simplified Relaxed Arakawa-Schubert (RAS) scheme. In contrast to the RAS, which specifies the cloud top, cloud top height in BMF depends on environmental properties, by considering the conditions of both the parcel and its environment in a fractional entrainment and detrainment rate formulations. As a result, BMF shows improved sensitivity in depth and strength of convection on environmental humidity compared to RAS, by strengthening coupling between cloud and environment. When the mid to lower troposphere is dry, the cloud resolving model and BMF produce cloud top around the dry layer and moisten the layer. In the framework of general circulation model, enhanced coupling between convection and environmental humidity in BMF results in improved representation of eastward propagating intraseasonal variability in the tropics—the Madden-Julian oscillation.     

Three-dimensional numerical study of turbulence in an entraining mixed layer

The mean structure calculated by a three-dimensional numerical model of a heated planetary boundary layer, in simulation of DAY 33 of the Australian Wangara data, has been previously described. The present study supplements it by describing properties of the calculated turbulence.A major finding is the importance of entrainment upon turbulence statistics relating to specific humidity, relative to those for potential temperature. The variances, skewness and spectra of velocity, temperature and humidity are presented, as are budget equations for kinetic energy, temperature and humidity variances and heat/moisture fluxes. These are interpreted with regard to the relative importance of the surface flux vs the flux due to entrainment at the top of the mixed layer, and in regard to the structure which would occur if the entrainment were to vanish.The Rotte-type closure assumption is tested for the correlation between the pressure fluctuation and the vertical gradient of vertical velocity, potential temperature, or specific humidity, and found to be qualitatively correct except near the top of the mixed layer.NCAR is sponsored by the National Science Foundation (U.S.A.).     

Resolved Versus Parametrized Boundary-Layer Plumes. Part I: A Parametrization-Oriented Conditional Sampling in Large-Eddy Simulations

A conditional sampling based on the combination of a passive tracer emitted at the surface and thermodynamic variables is proposed to characterise organized structures in large-eddy simulations of cloud-free and cloudy boundary layers. The sampling is evaluated against more traditional sampling of dry thermals or clouds. It enables the characterization of convective updrafts from the surface to the top of the boundary layer (or the top of cumulus clouds), describing in particular the transition from the sub-cloud to the cloud layer, and retrieves plume characteristics, entrainment and detrainment rates, variances and fluxes. This sampling is used to analyze the contribution of boundary-layer thermals to vertical fluxes and variances.     

Obtaining Potential Virtual Temperature Profiles,Entrainment Fluxes,and Spectra from Mini Unmanned Aerial Vehicle Data

We present a simple but effective small unmanned aerial vehicle design that is able to make high-resolution temperature and humidity measurements of the atmospheric boundary layer. The air model used is an adapted commercial design, and is able to carry all the instrumentation (barometer, temperature and humidity sensor, and datalogger) required for such measurements. It is fitted with an autopilot that controls the plane’s ascent and descent in a spiral to 1800 m above ground. We describe the results obtained on three different days when the plane, called Aerolemma-3, flew continuously throughout the day. Surface measurements of the sensible virtual heat flux made simultaneously allowed the calculation of all standard convective turbulence scales for the boundary layer, as well as a rigorous test of existing models for the entrainment flux at the top of the boundary layer, and for its growth. A novel approach to calculate the entrainment flux from the top-down, bottom-up model of Wynagaard and Brost is used. We also calculated temperature fluctuations by means of a spectral high-pass filter, and calculated their spectra. Although the time series are small, tapering proved ineffective in this case. The spectra from the untapered series displayed a consistent −5/3 behaviour, and from them it was possible to calculate a dimensionless dissipation function, which exhibited the expected similarity behaviour against boundary-layer bulk stability. The simplicity, ease of use and economy of such small aircraft make us optimistic about their usefulness in boundary-layer research.