Potential Temperature Profiles Research Articles

Parameterizations of the daytime friction velocity, temperature scale, and upslope flow over gently inclined terrain in calm synoptic conditions

A set of new parameterizations for the friction velocity and temperature scale over gently sloped terrain and in calm synoptic conditions are theoretically derived. The friction velocity is found to be proportional to the product of the square root of the total accumulated heating in the boundary layer and the sinusoidal function of the slope angle, while the temperature scale is proportional to the product of the boundary layer depth, the sinusoidal function of the slope angle and the potential temperature gradient in the free atmosphere. Using the new friction velocity parameterization, together with a parameterization of eddy diffusivity and an initial potential temperature profile around sunrise, an improved parameterization for the thermally induced upslope flow profile is derived by solving the Prandtl equations. The upslope flow profile is found to be simply proportional to the friction velocity.

Doppler Lidar Estimation of Mixing Height Using Turbulence, Shear, and Aerosol Profiles

Abstract The concept of boundary layer mixing height for meteorology and air quality applications using lidar data is reviewed, and new algorithms for estimation of mixing heights from various types of lower-tropospheric coherent Doppler lidar measurements are presented. Velocity variance profiles derived from Doppler lidar data demonstrate direct application to mixing height estimation, while other types of lidar profiles demonstrate relationships to the variance profiles and thus may also be used in the mixing height estimate. The algorithms are applied to ship-based, high-resolution Doppler lidar (HRDL) velocity and backscattered-signal measurements acquired on the R/V Ronald H. Brown during Texas Air Quality Study (TexAQS) 2006 to demonstrate the method and to produce mixing height estimates for that experiment. These combinations of Doppler lidar–derived velocity measurements have not previously been applied to analysis of boundary layer mixing height—over the water or elsewhere. A comparison of the results to those derived from ship-launched, balloon-radiosonde potential temperature and relative humidity profiles is presented.

Parameterization of PBL Processes in an Atmospheric General Circulation Model: Description and Preliminary Assessment

Abstract This paper presents the basic features of a newly developed planetary boundary layer (PBL) parameterization, and the performance assessment of a version of the University of California, Los Angeles (UCLA), Atmospheric General Circulation Model (AGCM) to which the parameterization is incorporated. The UCLA AGCM traditionally uses a framework in which a sigma-type vertical coordinate for the PBL shares a coordinate surface with the free atmosphere at the PBL top. This framework facilitates an explicit representation of processes concentrated near the PBL top, which is crucially important especially for predicting PBL clouds. In the new framework, multiple layers are introduced between the PBL top and earth’s surface, allowing for predictions of the vertical profiles of potential temperature, total water mixing ratio, and horizontal winds within the PBL. The vertically integrated “bulk” turbulent kinetic energy (TKE) is also predicted for the PBL. The PBL-top mass entrainment is determined through an equation including the effects of TKE and the radiative and evaporative cooling processes concentrated near the PBL top. The surface fluxes are determined from an aerodynamic formula in which the velocity scale depends both on the square root of TKE and the grid-scale PBL velocity at the lowermost model layer. The turbulent fluxes within the PBL are determined through an approach that includes the effects of both large convective and small diffusive eddies. AGCM simulations with the new formulation of PBL are analyzed with a focus on the seasonal and diurnal variations. The simulated seasonal cycle of stratocumulus over the eastern oceans is realistic, as are the diurnal cycles of the PBL depth and precipitation over land. The simulated fluxes of latent heat, momentum, and shortwave radiation at the ocean surface and baroclinic activity in the middle latitudes show significant improvements over the previous versions of the AGCM based on the single-layer PBL.

Surface Heterogeneity Effects on Regional-Scale Fluxes in Stable Boundary Layers: Surface Temperature Transitions

Abstract Large-eddy simulation, with recently developed dynamic subgrid-scale models, is used to study the effect of heterogeneous surface temperature distributions on regional-scale turbulent fluxes in the stable boundary layer (SBL). Simulations are performed of a continuously turbulent SBL with surface heterogeneity added in the form of streamwise transitions in surface temperature. Temperature differences between patches of 6 and 3 K are explored with patch length scales ranging from one-half to twice the equivalent homogeneous boundary layer height. The surface temperature heterogeneity has important effects on the mean wind speed and potential temperature profiles as well as on the surface heat flux distribution. Increasing the difference between the patch temperatures results in decreased magnitude of the average surface heat flux, with a corresponding increase in the mean potential temperature in the boundary layer. The simulation results are also used to test existing models for average surface fluxes over heterogeneous terrain. The tested models fail to fully represent the average turbulent heat flux, with models that break the domain into homogeneous subareas grossly underestimating the heat flux magnitude over patches with relatively colder surface temperatures. Motivated by these results, a new parameterization based on local similarity theory is proposed. The new formulation is found to correct the bias over the cold patches, resulting in improved average surface heat flux calculations.

Rainfall regimes associated with banded convection in the Cévennes-Vivarais area

It is well known that relief plays an important part in the triggering and enhancement of rainfall. Our study is focussed on the southeastern part of France where the topography and the specific meteorological conditions lead to an important hydrometeorological potential risk. From a climatologic point of view, we consider that these rainfall events arise from two kinds of convection: a deep convection essentially governed by synoptic conditions and where the relief has little direct impact and a shallow convection which is strongly controlled by the air circulation within the relief that leads to banded rain patterns. This study aims at understanding and analysing the atmospheric variables that control banded convection in this region. The methodology is first based on the exploration of the meteorological and rainfall data bases in order to identify the meteorological characteristics associated with this convection. Our results show that banded convection events can be characterized by specific vertical profiles of the dynamics variables (wind velocity, shear) and thermodynamic variables (humidity and wet bulb potential temperature profiles, stratification). We thus propose a generic sounding that is used as an input to the MesoNH meteorological model. Preliminary simulation results show a banded organisation of rainfall, which confirms our selection and allow for sensibility studies and further investigations of meteorological characteristics associated with this particular precipitation pattern.

Effects of water table dynamics on regional climate: A case study over east Asian monsoon area

Groundwater is an important component of the hydrologic cycle, and its anomaly will result in variations of soil moisture, water, and energy balances between the land surface and atmosphere, which ultimately influence climate. In this study, we implement a groundwater model into the regional climate model RegCM3, which is called RegCM3_GW, and investigate the effects of water table dynamics on regional climate. Numerical experiments by RegCM3_GW and RegCM3 over the east Asian monsoon area show that incorporating the water table dynamics into the regional climate model reduces the systematic biases of the simulated precipitation by 38.5% and 39.8% over semiarid and humid regions, respectively, and increases the bias slightly by 5.6% over semihumid regions. To seek the reasons for the differences of simulated precipitation, we analyze the atmospheric water vapor budget and the local water cycle among the water table, soil moisture, evapotranspiration (ET), and convective precipitation. It is found that the top and root zone soil layers become wetter and enhance the bare soil evaporation but do not always increase the transpiration. Because of the variations of each ET's component, the obvious enhancements of ET occur in semiarid regions and contribute to more instable profiles of pseudoequivalent potential temperature. The atmospheric moisture budget analysis indicates that the recycling rate and precipitation efficiency increase greatly over semiarid regions, which presents a local aquifer‐atmosphere feedback, while the variations of atmospheric water vapor transport control the development of precipitation over semihumid and humid regions. Therefore, the effects of water table dynamics on regional climate consist of the local aquifer‐atmosphere interaction and the changes of circulation originated from ambient aquifer‐atmosphere interaction, and the latter factor plays an important role in the monsoon area. Sensitivity of the results to a change in convection parameterization is also explored.

Thermodynamic Properties of Mesoscale Convective Systems Observed during BAMEX

Abstract Dropsonde observations from the Bow Echo and Mesoscale Convective Vortex Experiment (BAMEX) are used to document the spatiotemporal variability of temperature, moisture, and wind within mesoscale convective systems (MCSs). Onion-type sounding structures are found throughout the stratiform region of MCSs, but the temperature and moisture variability is large. Composite soundings were constructed and statistics of thermodynamic variability were generated within each subregion of the MCS. The calculated air vertical velocity helped identify subsaturated downdrafts. It was found that lapse rates within the cold pool varied markedly throughout the MCS. Layered wet-bulb potential temperature profiles seem to indicate that air within the lowest several kilometers comes from a variety of source regions. It was also found that lapse-rate transitions across the 0°C level were more common than isothermal, melting layers. The authors discuss the implications these findings have and how they can be used to validate future high-resolution numerical simulations of MCSs.

Seasonal variation of cloud activity and atmospheric profiles over the eastern part of the Tibetan Plateau

Cumulus activity over the Tibetan Plateau has a clear pattern of seasonal progression. Activity occurs with significant frequency in the eastern part of the plateau from April to mid‐May. From mid‐May to mid‐June the frequency of cumulus activity drastically decreases and then increases after mid‐June. Causes of such definite seasonal progression are revealed by in situ radiosonde observation data. In early spring the vertical profile of potential temperature is almost uniform from the surface to 7000 m above sea level around noon, and atmospheric stratification is vertically neutral for dry convection. Under such a condition, dry convection easily occurs and develops cumulus cloud that is frequently observed. From mid‐May to mid‐June the atmosphere has conditionally unstable stratification, but the atmosphere is relatively dry and hardly saturated and cumulus activity occurs less frequently. Though the atmosphere is also conditionally unstable after mid‐June, the increase in total precipitable water content allows easy atmospheric saturation and cumulus activity recommences. At the same time, satellite data and reanalysis data show, in the eastern part of the Tibetan Plateau, a consistent correspondence of the seasonal variation of cumulus activity, total precipitable water content, and vertical instability of the atmosphere. Results of synoptic‐scale analysis indicate the enhancement of cumulus convection in the northern part of India and the Tibetan High cause a moisture increase over the Tibetan Plateau and leads to the recommencement of frequent cumulus activity over the plateau.

Estimates of the height of the boundary layer using SODAR and rawinsoundings in Amazonia

During the LBA campaign in Amazonia 2002, simultaneous measurements were made of the boundary layer using different instruments (rawinsoundings and SODAR). The profiles of potential temperature and humidity were used to estimates the height of the boundary layer using 3 different techniques. The SODAR's measurements did not capture the shallow morning boundary layer observed at the profiles.

Comparison of nocturnal inversion characteristics obtained by sodar and microwave temperature profiler

Sodar echogram is a usual means to estimate the height of temperature inversion. Our data of long-term synchronous measurements with sodar and microwave profiler in urban and rural areas showed an intricate relation between intensity of backscattering and potential temperature profiles. Qualitative relationships between echogram pattern, wind speed and temperature profiles were obtained. Distributions of temperature gradients in rural and urban measurements for summer and winter time were calculated.

Simulations of the Urban Planetary Boundary Layer in an Arid Metropolitan Area

AbstractA modified version of the fifth-generation Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model (MM5) was applied to the arid Phoenix, Arizona, metropolitan region. The ability of the model to simulate characteristics of the summertime urban planetary boundary layer (PBL) was tested by comparing model results with observations from two field campaigns conducted in May/June 1998 and June 2001. The modified MM5 included a refined land use/cover classification and updated land use data for Phoenix and bulk approaches of characteristics of the urban surface energy balance. PBL processes were simulated by a version of MM5’s Medium-Range Forecast Model (MRF) scheme that was enhanced by new surface flux and nonlocal mixing approaches. Simulated potential temperature profiles were tested against radiosonde data, indicating that the modified MRF scheme was able to simulate vertical mixing and the evolution and height of the PBL with good accuracy and better than the original MRF scheme except in the late afternoon. During both simulation periods, it is demonstrated that the modified MM5 simulated near-surface air temperatures and wind speeds in the urban area consistently and considerably better than the standard MM5 and that wind direction simulations were improved slightly.

The YAK-AEROSIB transcontinental aircraft campaigns: new insights on the transport of CO2, CO and O3 across Siberia

Two airborne campaigns were carried out to measure the tropospheric concentrations and variability of CO2, CO and O3 over Siberia. In order to quantify the influence of remote and regional natural and anthropogenic sources, we analysed a total of 52 vertical profiles of these species collected in April and September 2006, every ∼200 km and up to 7 km altitude. CO2 and CO concentrations were high in April 2006 (respectively 385–390 ppm CO2 and 160–200 ppb CO) compared to background values. CO concentrations up to 220 ppb were recorded above 3.5 km over eastern Siberia, with enhancements in 500–1000 m thick layers. The presence of CO enriched air masses resulted from a quick frontal uplift of a polluted air mass exposed to northern China anthropogenic emissions and to fire emissions in northern Mongolia. A dominant Asian origin for CO above 4 km (71.0%) contrasted with a dominant European origin below this altitude (70.9%) was deduced both from a transport model analysis, and from the contrasted ΔCO/ΔCO2 ratio vertical distribution. In September 2006, a significant O3 depletion (∼ –30 ppb) was repeatedly observed in the boundary layer, as diagnosed from virtual potential temperature profiles and CO2 gradients, compared to the free troposphere aloft, suggestive of a strong O3 deposition over Siberian forests.

Measurements and Mesoscale Modeling of Autumnal Vertical Ozone Profiles in Southern Taiwan

Vertical measurements of ozone were made using a tethered balloon at the Linyuan site in Kaohsiung County, southern Taiwan. Ozone was monitored at altitudes of 0, 100, 300, 500, and 1000 m from November 23 to 25 in 2005. The potential temperature profiles revealed a stable atmosphere during the study period, largely because of the dominance of the high-pressure system and nocturnal radiation cooling close to the surface. The mixing height was low (50-300 m), particularly in the late night and early morning. The surface ozone concentrations that were predicted using TAPM (The Air Pollution Model) were high (33.7-119 ppbv) in the daytime (10:00-16:00) and were low (10-40 ppbv) at other times; the predictions of which were consistent with the observations. The simulated surface ozone concentrations reveal that costal lands typically had higher ozone concentrations than those inland, because most industrial parks are located in or close to the boundaries of Kaohsiung City. Both measurements and simulations indicate that daytime ozone concentrations decreased quickly with increasing height at altitudes below 300 m; while nighttime ozone concentrations were lower at low altitudes (50 to 300 m) than at higher altitudes, partly because of dry deposition and titration of surface ozone by the near-surface nitrogen oxides (NO(subscript x)) and partly because of the existence of the residual layer above the stable nocturnal boundary layer. The simulations show a good correlation between the maximum daytime surface ozone concentration and average nighttime ozone concentration above the nocturnal boundary layer.

On the Scale-dependence of the Gradient Richardson Number in the Residual Layer

We present results of a technique for examining the scale-dependence of the gradient Richardson number, Ri, in the nighttime residual layer. The technique makes use of a series of high-resolution, in situ, vertical profiles of wind speed and potential temperature obtained during CASES-99 in south-eastern Kansas, U.S.A. in October 1999. These profiles extended from the surface, through the nighttime stable boundary layer, and well into the residual layer. Analyses of the vertical gradients of both wind speed, potential temperature and turbulence profiles over a wide range of vertical scale sizes are used to estimate profiles of the local Ri and turbulence structure as a function of scale size. The utility of the technique lies both with the extensive height range of the residual layer as well as with the fact that the sub-metre resolution of the raw profiles enables a metre-by-metre ‘sliding’ average of the scale-dependent Richardson number values over hundreds of metres vertically. The results presented here show that small-scale turbulence is a ubiquitous and omnipresent feature of the residual layer, and that the region is dynamic and highly variable, exhibiting persistent turbulent structure on vertical scales of a few tens of metres or less. Furthermore, these scales are comparable to the scales over which the Ri is less than or equal to the critical value of Ric of 0.25, although turbulence is also shown to exist in regions with significantly larger Ri values, an observation at least consistent with the concept of hysteresis in turbulence generation and maintenance. Insofar as the important scale sizes are comparable to or smaller than the resolution of current models, it follows that, in order to resolve the observed details of small Ri values and the concomitant turbulence generation, future models need to be capable of significantly higher resolutions.

Atmospheric CO<sub>2</sub> modeling at the regional scale: an intercomparison of 5 meso-scale atmospheric models

Abstract. Atmospheric CO2 modeling in interaction with the surface fluxes, at the regional scale is developed within the frame of the European project CarboEurope-IP and its Regional Experiment component. In this context, five meso-scale meteorological models at 2 km resolution participate in an intercomparison exercise. Using a common experimental protocol that imposes a large number of rules, two days of the CarboEurope Regional Experiment Strategy (CERES) campaign are simulated. A systematic evaluation of the models is done in confrontation with the observations, using statistical tools and direct comparisons. Thus, temperature and relative humidity at 2 m, wind direction, surface energy and CO2 fluxes, vertical profiles of potential temperature as well as in-situ CO2 concentrations comparisons between observations and simulations are examined. These comparisons reveal a cold bias in the simulated temperature at 2 m, the latent heat flux is often underestimated. Nevertheless, the CO2 concentrations heterogeneities are well captured by most of the models. This intercomparison exercise shows also the models ability to represent the meteorology and carbon cycling at the synoptic and regional scale in the boundary layer, but also points out some of the major shortcomings of the models.

Observations of the Effects of Atmospheric Stability on Turbulence Statistics Deep within an Urban Street Canyon

Abstract Data obtained in downtown Oklahoma City, Oklahoma, during the Joint Urban 2003 atmospheric dispersion study have been analyzed to investigate the effects of upstream atmospheric stability on turbulence statistics in an urban core. The data presented include turbulent heat and momentum fluxes at various vertical and horizontal locations in the lower 30% of the street canyon. These data have been segregated into three broad stability classification regimes: stable (z/L > 0.2), neutral (−0.2 < z/L < 0.2), and unstable (z/L < −0.2) based on upstream measurements of the Monin–Obukhov length scale L. Most of the momentum-related turbulence statistics were insensitive to upstream atmospheric stability, while the energy-related statistics (potential temperatures and kinematic heat fluxes) were more sensitive. In particular, the local turbulence intensity inside the street canyon varied little with atmospheric stability but always had large magnitudes. Measurements of turbulent momentum fluxes indicate the existence of regions of upward transport of high horizontal momentum fluid near the ground that is associated with low-level jet structures for all stabilities. The turbulent kinetic energy normalized by a local shear stress velocity collapses the data well and shows a clear repeatable pattern that appears to be stability invariant. The magnitude of the normalized turbulent kinetic energy increases rapidly as the ground is approached. This behavior is a result of a much more rapid drop in the correlation between the horizontal and vertical velocities than in the velocity variances. This lack of correlation in the turbulent momentum fluxes is consistent with previous work in the literature. It was also observed that the mean potential temperatures almost always decrease with increasing height in the street canyon and that the vertical heat fluxes are always positive regardless of upstream atmospheric stability. In addition, mean potential temperature profiles are slightly more unstable during the unstable periods than during the neutral or stable periods. The magnitudes of all three components of the heat flux and the variability of the heat fluxes decrease with increasing atmospheric stability. In addition, the cross-canyon and along-canyon heat fluxes are as large as the vertical component of the heat fluxes in the lower portion of the canyon.

Role of land-surface temperature feedback on model performance for the stable boundary layer

At present a variety of boundary-layer schemes is in use in numerical models and often a large variation of model results is found. This is clear from model intercomparisons, such as organized within the GEWEX Atmospheric Boundary Layer Study (GABLS). In this paper we analyze how the specification of the land-surface temperature affects the results of a boundary-layer scheme, in particular for stable conditions. As such we use a well established column model of the boundary layer and we vary relevant parameters in the turbulence scheme for stable conditions. By doing so, we can reproduce the outcome for a variety of boundary-layer models. This is illustrated with the original set-up of the second GABLS intercomparison study using prescribed geostrophic winds and land-surface temperatures as inspired by (but not identical to) observations of CASES-99 for a period of more than two diurnal cycles. The model runs are repeated using a surface temperature that is calculated with a simple land-surface scheme. In the latter case, it is found that the range of model results in stable conditions is reduced for the sensible heat fluxes, and the profiles of potential temperature and wind speed. However, in the latter case the modelled surface temperatures are rather different than with the original set-up, which also impacts on near-surface air temperature and wind speed. As such it appears that the model results in stable conditions are strongly influenced by non-linear feedbacks in which the magnitude of the geostrophic wind speed and the related land-surface temperature play an important role.

Analysis of large‐scale conditions associated with convection over the Indian monsoon region

AbstractThis paper examines the thermodynamics of the atmosphere in relation to occurrence of convective rainfall over the Indian region. Various thermodynamical and kinematic parameters such as convective available potential energy (CAPE), convective inhibition energy (CINE), equivalent potential temperature and relative vorticity field are computed based on model analysis field of the limited area model of India Meteorological Department. The data period of this study is from 15 March 2001 to 28 February 2002. The study shows that the spatial distribution of CAPE in pre‐monsoon and monsoon season exhibits a belt of high CAPE along the east coast extending north up to Gangetic West Bengal and another zone of high CAPE along the southwest coast. During post‐monsoon and winter seasons, the belt of maximum CAPE shifts over the extreme southern Peninsula. The annual variation of CAPE over all stations is bi‐modal in nature with twin peaks each at the beginning and at the end of the monsoon season. Magnitude of CINE on the other hand is nearly zero during the monsoon season. The result shows that presence of strong thermodynamic environment is not sufficient for the occurrence of deep convection. Factors like proper dynamic conditions play a very important role in controlling the occurrence of deep convection. The surface warming during the pre‐monsoon season destabilizes the atmosphere, but large‐scale moist convection, generally, does not take place until CINE becomes close to zero (in the monsoon season). Preceding the monsoon season, the land sea contrast heating strengthens circulation in the lower atmosphere, which in turn leads to the moistening and destabilization of lapse rate. This process starts about a month prior to the actual start of the monsoon rainfall, allowing moisture to build up to lead to the widespread monsoon rainfall. During winter months, rainfall occurs over northwest India in association with high negative value of CINE and over southern Peninsula in association with high value of CAPE. Substantial changes are noticed in the equivalent potential temperature (θe) and vorticity profiles between days with and without rainfall over northwest India (Delhi). The results suggest that the increase in moisture in association with the occurrence of rainfall over Delhi is of advective origin. Where as, other regions (Kolkata, Mumbai and Thiruvananthapuram) display relatively smaller changes in the vertical profiles of θe, indicating an in situ presence of moisture. Copyright © 2007 Royal Meteorological Society

Similarity theory and calculation of turbulent fluxes at the surface for the stably stratified atmospheric boundary layer

In this paper we revise the similarity theory for the stably stratified atmospheric boundary layer (ABL), formulate analytical approximations for the wind velocity and potential temperature profiles over the entire ABL, validate them against large-eddy simulation and observational data, and develop an improved surface flux calculation technique for use in operational models.

Comparison between Backscatter Lidar and Radiosonde Measurements of the Diurnal and Nocturnal Stratification in the Lower Troposphere

AbstractA collection of boundary layer heights has been derived from measurements performed by a ground-based backscatter lidar in Neuchâtel, Switzerland (47.000°N, 6.967°E, 485 m ASL). A dataset of 98 cases have been collected during 2 yr. From these data, 61 are noon and 37 are midnight cases. The following two different schemes were used to retrieve the mixed layer depth and the height of the residual layer from the measurements: the gradient and variance methods. The obtained values were compared with those derived from the potential temperature profiles as computed from radiosonde data. For nocturnal cases, the height of the first aerosol layer above the residual layer was also compared to the corresponding potential temperature value. Correlation coefficients between lidar and radiosonde in both convective and stable conditions are between 0.88 and 0.97.