Characterising the background errors for the boundary-layer capping inversion

The one-dimensional variational assimilation of vertical temperature information in the presence of a boundary-layer capping inversion is studied. For an optimal analysis of the vertical temperature profile, an accurate representation of the background error covariances is essential. The background error covariances are highly flow-dependent due to the variability in the presence, structure and height of the boundary-layer capping inversion. Flow-dependent estimates of the background error covariances are shown by studying the spread in an ensemble of forecasts. A forecast of the temperature profile (used as a background state) may have a significant error in the position of the capping inversion with respect to observations. It is shown that the assimilation of observations may weaken the inversion structure in the analysis if only magnitude errors are accounted for as is the case for traditional data assimilation methods used for operational weather prediction. The positional error is treated explicitly here in a new data assimilation scheme to reduce positional error, in addition to the traditional framework to reduce magnitude error. The distribution of the positional error of the background inversion is estimated for use with the new scheme.

Wave-driven wind jets in the marine atmospheric boundary layer

The interaction between ocean surface waves and the overlying wind leads to a transfer of momentum across the air–sea interface. Atmospheric and oceanic models typically allow for momentum transfer to be directed only downward, from the atmosphere to the ocean. Recent observations have suggested that momentum can also be transferred upward when long wavelength waves, characteristic of remotely generated swell, propagate faster than the wind speed. The effect of upward momentum transfer on the marine atmospheric boundary layer is investigated here using idealized models that solve the momentum budget above the ocean surface. A variant of the classical Ekman model that accounts for the wave-induced stress demonstrates that, although the momentum flux due to the waves penetrates only a small fraction of the depth of the boundary layer, the wind profile is profoundly changed through its whole depth. When the upward momentum transfer from surface waves sufficiently exceeds the downward turbulent momentum flux, then the near-surface wind accelerates, resulting in a low-level wave-driven wind jet. This increases the Coriolis force in the boundary layer, and so the wind turns in the opposite direction to the classical Ekman layer. Calculations of the wave-induced stress due to a wave spectrum representative of fast-moving swell demonstrate upward momentum transfer that is dominated by contributions from waves in the vicinity of the peak in the swell spectrum. This is in contrast to wind-driven waves whose wave-induced stress is dominated by very short wavelength waves. Hence the role of swell can be characterized by the inverse wave age based on the wave phase speed corresponding to the peak in the spectrum. For a spectrum of waves, the total momentum flux is found to reverse sign and become upward, from waves to wind, when the inverse wave age drops below the range 0.15–0.2, which agrees reasonably well with previously published oceanic observations.

Synoptic controls on boundary-layer characteristics

We report the characteristics of the three-dimensional, time evolving, atmospheric boundary layer that develops beneath an idealised, dry, baroclinic weather system. The boundary-layer structure is forced by thermal advection associated with the weather system. Large positive heat fluxes behind the cold front drive a vigorous convective boundary layer, whereas moderate negative heat fluxes in the warm sector between the cold and warm fronts generate shallow, stably stratified or neutral boundary layers. The forcing of the boundary-layer structure is quantified by forming an Eulerian mass budget integrated over the depth of the boundary layer. The mass budget indicates that tropospheric air is entrained into the boundary layer both in the vicinity of the high-pressure centre, and behind the cold front. It is then transported horizontally within the boundary layer and converges towards the cyclone's warm sector, whence it is ventilated out into the troposphere. This cycling of air is likely to be important for the ventilation of pollution out of the boundary layer, and for the transformation of the properties of large-scale air masses.

Controls on boundary layer ventilation: Boundary layer processes and large-scale dynamics

Midlatitude cyclones are important contributors to boundary layer ventilation. However, it is uncertain how efficient such systems are at transporting pollutants out of the boundary layer, and variations between cyclones are unexplained. In this study 15 idealized baroclinic life cycles, with a passive tracer included, are simulated to identify the relative importance of two transport processes: horizontal divergence and convergence within the boundary layer and large-scale advection by the warm conveyor belt. Results show that the amount of ventilation is insensitive to surface drag over a realistic range of values. This indicates that although boundary layer processes are necessary for ventilation they do not control the magnitude of ventilation. A diagnostic for the mass flux out of the boundary layer has been developed to identify the synoptic-scale variables controlling the strength of ascent in the warm conveyor belt. A very high level of correlation (R-2 values exceeding 0.98) is found between the diagnostic and the actual mass flux computed from the simulations. This demonstrates that the large-scale dynamics control the amount of ventilation, and the efficiency of midlatitude cyclones to ventilate the boundary layer can be estimated using the new mass flux diagnostic. We conclude that meteorological analyses, such as ERA-40, are sufficient to quantify boundary layer ventilation by the large-scale dynamics.

Evaluating forecasts of the evolution of the cloudy boundary layer using diurnal composites of radar and lidar observations

Observations of boundary-layer cloud have been made using radar and lidar at Chilbolton, Hampshire, UK. These have been compared with output from 7 different global and regional models. Fifty-five cloudy days have been composited to reveal the mean diurnal variation of cloud top and base heights, cloud thickness and liquid water path of the clouds. To enable like-for-like comparison between model and observations, the observations have been averaged on to the grid of each model. The composites show a distinct diurnal cycle in observed cloud; the cloud height exhibits a sinusoidal variation throughout the day with a maximum at around 1600 and a minimum at around 0700 UTC. This diurnal cycle is captured by six of the seven models analysed, although the models generally under-predict both cloud top and cloud base heights throughout the day. The two worst performing models in terms of cloud boundaries also have biases of around a factor of two in liquid water path; these were the only two models that did not include an explicit formulation for cloud-top entrainment.

Sea-breeze modification of the growth of a marine internal boundary layer

A numerical mesoscale model is used to make a high-resolution simulation of the marine boundary layer in the Persian Gulf, during conditions of offshore flow from Saudi Arabia. A marine internal boundary layer (MIBL) and a sea-breeze circulation (SBC) are found to co-exist. The sea breeze develops in the mid-afternoon, at which time its front is displaced several tens of kilometres offshore. Between the coast and the sea-breeze system, the MIBL that occurs is consistent with a picture described in the existing literature. However, the MIBL is perturbed by the SBC, the boundary layer deepening significantly seaward of the sea-breeze front. Our analysis suggests that this strong, localized deepening is not a direct consequence of frontal uplift, but rather that the immediate cause is the retardation of the prevailing, low-level offshore wind by the SBC. The simulated boundary-layer development can be accounted for by using a simple 1D Lagrangian model of growth driven by the surface heat flux. This model is obtained as a straightforward modification of an established MIBL analytic growth model.

Numerical modeling of the propagation environment in the atmospheric boundary layer over the Persian Gulf

Strong vertical gradients at the top of the atmospheric boundary layer affect the propagation of electromagnetic waves and can produce radar ducts. A three-dimensional, time-dependent, nonhydrostatic numerical model was used to simulate the propagation environment in the atmosphere over the Persian Gulf when aircraft observations of ducting had been made. A division of the observations into high- and low-wind cases was used as a framework for the simulations. Three sets of simulations were conducted with initial conditions of varying degrees of idealization and were compared with the observations taken in the Ship Antisubmarine Warfare Readiness/Effectiveness Measuring (SHAREM-115) program. The best results occurred with the initialization based on a sounding taken over the coast modified by the inclusion of data on low-level atmospheric conditions over the Gulf waters. The development of moist, cool, stable marine internal boundary layers (MIBL) in air flowing from land over the waters of the Gulf was simulated. The MIBLs were capped by temperature inversions and associated lapses of humidity and refractivity. The low-wind MIBL was shallower and the gradients at its top were sharper than in the high-wind case, in agreement with the observations. Because it is also forced by land–sea contrasts, a sea-breeze circulation frequently occurs in association with the MIBL. The size, location, and internal structure of the sea-breeze circulation were realistically simulated. The gradients of temperature and humidity that bound the MIBL cause perturbations in the refractivity distribution that, in turn, lead to trapping layers and ducts. The existence, location, and surface character of the ducts were well captured. Horizontal variations in duct characteristics due to the sea-breeze circulation were also evident. The simulations successfully distinguished between high- and low-wind occasions, a notable feature of the SHAREM-115 observations. The modeled magnitudes of duct depth and strength, although leaving scope for improvement, were most encouraging.

Idealized model for changes in equilibrium temperature, mixed layer depth, and boundary layer cloud over land in a doubled CO2 climate

An idealized equilibrium model for the undisturbed partly cloudy boundary layer (BL) is used as a framework to explore the coupling of the energy, water, and carbon cycles over land in midlatitudes and show the sensitivity to the clear‐sky shortwave flux, the midtropospheric temperature, moisture, CO2, and subsidence. The changes in the surface fluxes, the BL equilibrium, and cloud cover are shown for a warmer, doubled CO2 climate. Reduced stomatal conductance in a simple vegetation model amplifies the background 2 K ocean temperature rise to an (unrealistically large) 6 K increase in near‐surface temperature over land, with a corresponding drop of near‐surface relative humidity of about 19%, and a rise of cloud base of about 70 hPa. Cloud changes depend strongly on changes of mean subsidence; but evaporative fraction (EF) decreases. EF is almost uniquely related to mixed layer (ML) depth, independent of background forcing climate. This suggests that it might be possible to infer EF for heterogeneous landscapes from ML depth. The asymmetry of increased evaporation over the oceans and reduced transpiration over land increases in a warmer doubled CO2 climate.

Generation of inertia-gravity waves in the rotating thermal annulus by a localised boundary layer instability

Waves with periods shorter than the inertial period exist in the atmosphere (as inertia-gravity waves) and in the oceans (as Poincaré and internal gravity waves). Such waves owe their origin to various mechanisms, but of particular interest are those arising either from local secondary instabilities or spontaneous emission due to loss of balance. These phenomena have been studied in the laboratory, both in the mechanically-forced and the thermally-forced rotating annulus. Their generation mechanisms, especially in the latter system, have not yet been fully understood, however. Here we examine short period waves in a numerical model of the rotating thermal annulus, and show how the results are consistent with those from earlier laboratory experiments. We then show how these waves are consistent with being inertia-gravity waves generated by a localised instability within the thermal boundary layer, the location of which is determined by regions of strong shear and downwelling at certain points within a large-scale baroclinic wave flow. The resulting instability launches small-scale inertia-gravity waves into the geostrophic interior of the flow. Their behaviour is captured in fully nonlinear numerical simulations in a finite-difference, 3D Boussinesq Navier-Stokes model. Such a mechanism has many similarities with those responsible for launching small- and meso-scale inertia-gravity waves in the atmosphere from fronts and local convection.

Boundary layer dynamics over London, UK, as observed using Doppler lidar during REPARTEE-II

Urban boundary layers (UBLs) can be highly complex due to the heterogeneous roughness and heating of the surface, particularly at night. Due to a general lack of observations, it is not clear whether canonical models of boundary layer mixing are appropriate in modelling air quality in urban areas. This paper reports Doppler lidar observations of turbulence profiles in the centre of London, UK, as part of the second REPARTEE campaign in autumn 2007. Lidar-measured standard deviation of vertical velocity averaged over 30 min intervals generally compared well with in situ sonic anemometer measurements at 190 m on the BT telecommunications Tower. During calm, nocturnal periods, the lidar underestimated turbulent mixing due mainly to limited sampling rate. Mixing height derived from the turbulence, and aerosol layer height from the backscatter profiles, showed similar diurnal cycles ranging from c. 300 to 800 m, increasing to c. 200 to 850 m under clear skies. The aerosol layer height was sometimes significantly different to the mixing height, particularly at night under clear skies. For convective and neutral cases, the scaled turbulence profiles resembled canonical results; this was less clear for the stable case. Lidar observations clearly showed enhanced mixing beneath stratocumulus clouds reaching down on occasion to approximately half daytime boundary layer depth. On one occasion the nocturnal turbulent structure was consistent with a nocturnal jet, suggesting a stable layer. Given the general agreement between observations and canonical turbulence profiles, mixing timescales were calculated for passive scalars released at street level to reach the BT Tower using existing models of turbulent mixing. It was estimated to take c. 10 min to diffuse up to 190 m, rising to between 20 and 50 min at night, depending on stability. Determination of mixing timescales is important when comparing to physico-chemical processes acting on pollutant species measured simultaneously at both the ground and at the BT Tower during the campaign. From the 3 week autumnal data-set there is evidence for occasional stable layers in central London, effectively decoupling surface emissions from air aloft.

Determining the contribution of volcanic ash and boundary layer aerosol in backscatter lidar returns: a three‐component atmosphere approach

A solution of the lidar equation is discussed, that permits combining backscatter and depolarization measurements to quantitatively distinguish two different aerosol types with different depolarization properties. The method has been successfully applied to simultaneous observations of volcanic ash and boundary layer aerosol obtained in Exeter, United Kingdom, on 16 and 18 April 2010, permitting the contribution of the two aerosols to be quantified separately. First a subset of the atmospheric profiles is used where the two aerosol types belong to clearly distinguished layers, for the purpose of characterizing the ash in terms of lidar ratio and depolarization. These quantities are then used in a three‐component atmosphere solution scheme of the lidar equation applied to the full data set, in order to compute the optical properties of both aerosol types separately. On 16 April a thin ash layer, 100–400 m deep, is observed (average and maximum estimated ash optical depth: 0.11 and 0.2); it descends from ∼2800 to ∼1400 m altitude over a 6‐hour period. On 18 April a double ash layer, ∼400 m deep, is observed just above the morning boundary layer (average and maximum estimated ash optical depth: 0.19 and 0.27). In the afternoon the ash is entrained into the boundary layer, and the latter reaches a depth of ∼1800 m (average and maximum estimated ash optical depth: 0.1 and 0.15). An additional ash layer, with a very small optical depth, was observed on 18 April at an altitude of 3500–4000 m. By converting the lidar optical measurements using estimates of volcanic ash specific extinction, derived from other works, the observations seem to suggest approximate peak ash concentrations of ∼1500 and ∼1000 mg/m3,respectively, on the two observations dates.

RACORO extended-term, aircraft observations of boundary-layer clouds

A first-of-a-kind, extended-term cloud aircraft campaign was conducted to obtain an in-situ statistical characterization of continental boundary-layer clouds needed to investigate cloud processes and refine retrieval algorithms. Coordinated by the Atmospheric Radiation Measurement (ARM) Aerial Facility (AAF), the Routine AAF Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations (RACORO) field campaign operated over the ARM Southern Great Plains (SGP) site from 22 January to 30 June 2009, collecting 260 h of data during 59 research flights. A comprehensive payload aboard the Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft measured cloud microphysics, solar and thermal radiation, physical aerosol properties, and atmospheric state parameters. Proximity to the SGP's extensive complement of surface measurements provides ancillary data that supports modeling studies and facilitates evaluation of a variety of surface retrieval algorithms. The five-month duration enabled sampling a range of conditions associated with the seasonal transition from winter to summer. Although about two-thirds of the cloud flights occurred in May and June, boundary-layer cloud fields were sampled under a variety of environmental and aerosol conditions, with about 77% of the flights occurring in cumulus and stratocumulus. Preliminary analyses illustrate use of these data to analyze cloud-aerosol relationships, characterize the horizontal variability of cloud radiative impacts, and evaluate surface-based retrievals. We discuss how an extended-term campaign requires a simplified operating paradigm that is different from that used for typical, short-term, intensive aircraft field programs.