Doppler lidar measurements of oriented planar ice crystals falling from supercooled and glaciated layer clouds

The properties of planar ice crystals settling horizontally have been investigated using a vertically pointing Doppler lidar. Strong specular reflections were observed from their oriented basal facets, identified by comparison with a second lidar pointing 4° from zenith. Analysis of 17 months of continuous high-resolution observations reveals that these pristine crystals are frequently observed in ice falling from mid-level mixed-phase layer clouds (85% of the time for layers at −15 °C). Detailed analysis of a case study indicates that the crystals are nucleated and grow rapidly within the supercooled layer, then fall out, forming well-defined layers of specular reflection. From the lidar alone the fraction of oriented crystals cannot be quantified, but polarimetric radar measurements confirmed that a substantial fraction of the crystal population was well oriented. As the crystals fall into subsaturated air, specular reflection is observed to switch off as the crystal faces become rounded and lose their faceted structure. Specular reflection in ice falling from supercooled layers colder than −22 °C was also observed, but this was much less pronounced than at warmer temperatures: we suggest that in cold clouds it is the small droplets in the distribution that freeze into plates and produce specular reflection, whilst larger droplets freeze into complex polycrystals. The lidar Doppler measurements show that typical fall speeds for the oriented crystals are ≈ 0.3 m s−1, with a weak temperature correlation; the corresponding Reynolds number is Re ∼ 10, in agreement with light-pillar measurements. Coincident Doppler radar observations show no correlation between the specular enhancement and the eddy dissipation rate, indicating that turbulence does not control crystal orientation in these clouds. Copyright © 2010 Royal Meteorological Society

Parameter estimation using a particle method: inferring mixing coefficients from sea level observations

This paper presents a first attempt to estimate mixing parameters from sea level observations using a particle method based on importance sampling. The method is applied to an ensemble of 128 members of model simulations with a global ocean general circulation model of high complexity. Idealized twin experiments demonstrate that the method is able to accurately reconstruct mixing parameters from an observed mean sea level field when mixing is assumed to be spatially homogeneous. An experiment with inhomogeneous eddy coefficients fails because of the limited ensemble size. This is overcome by the introduction of local weighting, which is able to capture spatial variations in mixing qualitatively. As the sensitivity of sea level for variations in mixing is higher for low values of mixing coefficients, the method works relatively well in regions of low eddy activity.

Flow structure and variability in the subtropical Indian Ocean: Instability of the South Indian Ocean Countercurrent

The origin of the eddy variability around the 25°S band in the Indian Ocean is investigated. We have found that the surface circulation east of Madagascar shows an anticyclonic subgyre bounded to the south by eastward flow from southwest Madagascar, and to the north by the westward flowing South Equatorial Current (SEC) between 15° and 20°S. The shallow, eastward flowing South Indian Ocean Countercurrent (SICC) extends above the deep reaching, westward flowing SEC to 95°E around the latitude of the high variability band. Applying a two-layer model reveals that regions of large vertical shear along the SICC-SEC system are baroclinically unstable. Estimates of the frequencies (3.5–6 times/year) and wavelengths (290–470 km) of the unstable modes are close to observations of the mesoscale variability derived from altimetry data. It is likely then that Rossby wave variability locally generated in the subtropical South Indian Ocean by baroclinic instability is the origin of the eddy variability around 25°S as seen, for example, in satellite altimetry.

Quantifying mesoscale variability in ocean transient tracer fields

A powerful way to test the realism of ocean general circulation models is to systematically compare observations of passive tracer concentration with model predictions. The general circulation models used in this way cannot resolve a full range of vigorous mesoscale activity (on length scales between 10–100 km). In the real ocean, however, this activity causes important variability in tracer fields. Thus, in order to rationally compare tracer observations with model predictions these unresolved fluctuations (the model variability error) must be estimated. We have analyzed this variability using an eddy‐resolving reduced‐gravity model in a simple midlatitude double‐gyre configuration. We find that the wave number spectrum of tracer variance is only weakly sensitive to the distribution of (large scale slowly varying) tracer sources and sinks. This suggests that a universal passive tracer spectrum may exist in the ocean. We estimate the spectral shape using high‐resolution measurements of potential temperature on an isopycnal in the upper northeast Atlantic Ocean, finding a slope near k −1.7 between 10 and 500 km. The typical magnitude of the variance is estimated by comparing tracer simulations using different resolutions. For CFC‐ and tritium‐type transient tracers the peak magnitude of the model variability saturation error may reach 0.20 for scales shorter than 100 km. This is of the same order as the time mean saturation itself and well over an order of magnitude greater than the instrumental uncertainty.

Water use of a bioenergy plantation increases in a future high CO2 world

Fast-growing poplar trees may in future be used as a source of renewable energy for heat, electricity and biofuels such as bioethanol. Water use in Populus x euramericana (clone I214), following long-term exposure to elevated CO2 in the POPFACE (poplar free-air carbon dioxide enrichment) experiment, is quantified here. Stomatal conductance was measured and, during two measurement campaigns made before and after coppicing, whole-tree water use was determined using heat-balance sap-flow gauges, first validated using eddy covariance measurements of latent heat flux. Water use was determined by the balance between leaf-level reductions in stomatal conductance and tree-level stimulations in transpiration. Reductions in stomatal conductance were found that varied between 16 and 39% relative to ambient air. Whole-tree sap flow was increased in plants growing under elevated CO2, on average, by 12 and 23%, respectively, in the first and in the second measurement campaigns. These results suggest that future CO2 concentrations may result in an increase in seasonal water use in fast-growing, short-rotation Populus plantations.

On the structure of Langmuir turbulence

The Stokes drift induced by surface waves distorts turbulence in the wind-driven mixed layer of the ocean, leading to the development of streamwise vortices, or Langmuir circulations, on a wide range of scales. We investigate the structure of the resulting Langmuir turbulence, and contrast it with the structure of shear turbulence, using rapid distortion theory (RDT) and kinematic simulation of turbulence. Firstly, these linear models show clearly why elongated streamwise vortices are produced in Langmuir turbulence, when Stokes drift tilts and stretches vertical vorticity into horizontal vorticity, whereas elongated streaky structures in streamwise velocity fluctuations (u) are produced in shear turbulence, because there is a cancellation in the streamwise vorticity equation and instead it is vertical vorticity that is amplified. Secondly, we develop scaling arguments, illustrated by analysing data from LES, that indicate that Langmuir turbulence is generated when the deformation of the turbulence by mean shear is much weaker than the deformation by the Stokes drift. These scalings motivate a quantitative RDT model of Langmuir turbulence that accounts for deformation of turbulence by Stokes drift and blocking by the air–sea interface that is shown to yield profiles of the velocity variances in good agreement with LES. The physical picture that emerges, at least in the LES, is as follows. Early in the life cycle of a Langmuir eddy initial turbulent disturbances of vertical vorticity are amplified algebraically by the Stokes drift into elongated streamwise vortices, the Langmuir eddies. The turbulence is thus in a near two-component state, with suppressed and . Near the surface, over a depth of order the integral length scale of the turbulence, the vertical velocity (w) is brought to zero by blocking of the air–sea interface. Since the turbulence is nearly two-component, this vertical energy is transferred into the spanwise fluctuations, considerably enhancing at the interface. After a time of order half the eddy decorrelation time the nonlinear processes, such as distortion by the strain field of the surrounding eddies, arrest the deformation and the Langmuir eddy decays. Presumably, Langmuir turbulence then consists of a statistically steady state of such Langmuir eddies. The analysis then provides a dynamical connection between the flow structures in LES of Langmuir turbulence and the dominant balance between Stokes production and dissipation in the turbulent kinetic energy budget, found by previous authors.

Baroclinic stationary waves in aquaplanet models

An aquaplanet model is used to study the nature of the highly persistent low-frequency waves that have been observed in models forced by zonally symmetric boundary conditions. Using the Hayashi spectral analysis of the extratropical waves, the authors find that a quasi-stationary wave 5 belongs to a wave packet obeying a well-defined dispersion relation with eastward group velocity. The components of the dispersion relation with k ≥ 5 baroclinically convert eddy available potential energy into eddy kinetic energy, whereas those with k < 5 are baroclinically neutral. In agreement with Green’s model of baroclinic instability, wave 5 is weakly unstable, and the inverse energy cascade, which had been previously proposed as a main forcing for this type of wave, only acts as a positive feedback on its predominantly baroclinic energetics. The quasi-stationary wave is reinforced by a phase lock to an analogous pattern in the tropical convection, which provides further amplification to the wave. It is also found that the Pedlosky bounds on the phase speed of unstable waves provide guidance in explaining the latitudinal structure of the energy conversion, which is shown to be more enhanced where the zonal westerly surface wind is weaker. The wave’s energy is then trapped in the waveguide created by the upper tropospheric jet stream. In agreement with Green’s theory, as the equator-to-pole SST difference is reduced, the stationary marginally stable component shifts toward higher wavenumbers, while wave 5 becomes neutral and westward propagating. Some properties of the aquaplanet quasi-stationary waves are found to be in interesting agreement with a low frequency wave observed by Salby during December–February in the Southern Hemisphere so that this perspective on low frequency variability, apart from its value in terms of basic geophysical fluid dynamics, might be of specific interest for studying the earth’s atmosphere.

Controls of carbon dioxide concentrations and fluxes above central London

Eddy-covariance measurements of carbon dioxide fluxes were taken semi-continuously between October 2006 and May 2008 at 190 m height in central London (UK) to quantify emissions and study their controls. Inner London, with a population of 8.2 million (~5000 inhabitants per km2) is heavily built up with 8% vegetation cover within the central boroughs. CO2 emissions were found to be mainly controlled by fossil fuel combustion (e.g. traffic, commercial and domestic heating). The measurement period allowed investigation of both diurnal patterns and seasonal trends. Diurnal averages of CO2 fluxes were found to be highly correlated to traffic. However changes in heating-related natural gas consumption and, to a lesser extent, photosynthetic activity that controlled the seasonal variability. Despite measurements being taken at ca. 22 times the mean building height, coupling with street level was adequate, especially during daytime. Night-time saw a higher occurrence of stable or neutral stratification, especially in autumn and winter, which resulted in data loss in post-processing. No significant difference was found between the annual estimate of net exchange of CO2 for the expected measurement footprint and the values derived from the National Atmospheric Emissions Inventory (NAEI), with daytime fluxes differing by only 3%. This agreement with NAEI data also supported the use of the simple flux footprint model which was applied to the London site; this also suggests that individual roughness elements did not significantly affect the measurements due to the large ratio of measurement height to mean building height.

On the response of the Antarctic Circumpolar Current transport to climate change in coupled climate models

The CMIP3 (IPCC AR4) models show a consistent intensification and poleward shift of the westerly winds over the Southern Ocean during the 21st century. However, the responses of the Antarctic Circumpolar Currents (ACC) show great diversity in these models, with many even showing reductions in transport. To obtain some understanding of diverse responses in the ACC transport, we investigate both external atmospheric and internal oceanic processes that control the ACC transport responses in these models. While the strengthened westerlies act to increase the tilt of isopycnal surfaces and hence the ACC transport through Ekman pumping effects, the associated changes in buoyancy forcing generally tend to reduce the surface meridional density gradient. The steepening of isopycnal surfaces induced by increased wind forcing leads to enhanced (parameterized) eddy-induced transports that act to reduce the isopycnal slopes. There is also considerable narrowing of the ACC that tends to reduce the ACC transport, caused mainly by the poleward shifts of the subtropical gyres and to a lesser extent by the equatorward expansions of the subpolar gyres in some models. If the combined effect of these retarding processes is larger than that of enhanced Ekman pumping, the ACC transport will be reduced. In addition, the effect of Ekman pumping on the ACC is reduced in weakly stratified models. These findings give insight into the reliability of IPCC-class model predictions of the Southern Ocean circulation, and into the observed decadal-scale steady ACC transport.

Spin-up and adjustment of the Antarctic Circumpolar Current and global pycnocline

A theory is presented for the adjustment of the Antarctic Circumpolar Current (ACC) and global pycnocline to a sudden and sustained change in wind forcing. The adjustment timescale is controlled by the mesoscale eddy diffusivity across the ACC, the mean width of the ACC, the surface area of the ocean basins to the north, and deep water formation in the North Atlantic. In particular, northern sinking may have the potential to shorten the timescale and reduce its sensitivity to Southern Ocean eddies, but the relative importance of northern sinking and Southern Ocean eddies cannot be determined precisely, largely due to limitations in the parameterization of northern sinking. Although it is clear that the main processes that control the adjustment timescale are those which counteract the deepening of the global pycnocline, the theory also suggests that the timescale can be subtly modified by wind forcing over the ACC and global diapycnal mixing. Results from calculations with a reduced-gravity model compare well with the theory. The multidecadal-centennial adjustment timescale implies that long observational time series will be required to detect dynamic change in the ACC due to anthropogenic forcing. The potential role of Southern Ocean mesoscale eddy activity in determining both the equilibrium state of the ACC and the timescale over which it adjusts suggests that the response to anthropogenic forcing may be different in coupled ocean-atmosphere climate models that parameterize and resolve mesoscale eddies.

Further investigation of the impact of idealised continents and SST distributions on the Northern Hemisphere storm tracks

Building on studies by Brayshaw et al. (2009, 2011) of the basic ingredients of the North Atlantic storm track (land-sea contrast, orography and SST), this article investigates the impact of Eurasian topography and Pacific SST anomalies on the North Pacific and Atlantic storm tracks through a hierarchy of atmospheric GCM simulations using idealised boundary conditions in the HadGAM1 model. The Himalaya-Tibet mountain complex is found to play a crucial role in shaping the North Pacific storm track. The northward deflection of the westerly flow around northern Tibet generates an extensive pool of very cold air in the north-eastern tip of the Asian continent, which strengthens the meridional temperature gradient and favours baroclinic growth in the western Pacific. The Kuroshio SST front is also instrumental in strengthening the Pacific storm track through its impact on near-surface baroclinicity, while the warm waters around Indonesia tend to weaken it through the impact on baroclinicity of stationary Rossby waves propagating poleward from the convective heating regions. Three mechanisms by which the Atlantic storm track may be affected by changes in the boundary conditions upstream of the Rockies are discussed. In the model configuration used here, stationary Rossby waves emanating from Tibet appear to weaken the North Atlantic storm track substantially, whereas those generated over the cold waters off Peru appear to strengthen it. Changes in eddy-driven surface winds over the Pacific generally appear to modify the flow over the Rocky Mountains, leading to consistent modifications in the Atlantic storm track. The evidence for each of these mechanisms is, however, ultimately equivocal in these simulations.

Tropical Cyclones in a hieararchy of climate models of increasing resolution

Tropical Cyclone (TC) is normally not studied at the individual level with Global Climate Models (GCMs), because the coarse grid spacing is often deemed insufficient for a realistic representation of the basic underlying processes. GCMs are indeed routinely deployed at low resolution, in order to enable sufficiently long integrations, which means that only large-scale TC proxies are diagnosed. A new class of GCMs is emerging, however, which is capable of simulating TC-type vortexes by retaining a horizontal resolution similar to that of operational NWP GCMs; their integration on the latest supercomputers enables the completion of long-term integrations. The UK-Japan Climate Collaboration and the UK-HiGEM projects have developed climate GCMs which can be run routinely for decades (with grid spacing of 60 km) or centuries (with grid spacing of 90 km); when coupled to the ocean GCM, a mesh of 1/3 degrees provides eddy-permitting resolution. The 90 km resolution model has been developed entirely by the UK-HiGEM consortium (together with its 1/3 degree ocean component); the 60 km atmospheric GCM has been developed by UJCC, in collaboration with the Met Office Hadley Centre.

Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain

Scintillometry is an established technique for determining large areal average sensible heat fluxes. The scintillometer measurement is related to sensible heat flux via Monin–Obukhov similarity theory, which was developed for ideal homogeneous land surfaces. In this study it is shown that judicious application of scintillometry over heterogeneous mixed agriculture on undulating topography yields valid results when compared to eddy covariance (EC). A large aperture scintillometer (LAS) over a 2.4 km path was compared with four EC stations measuring sensible (H) and latent (LvE) heat fluxes over different vegetation (cereals and grass) which when aggregated were representative of the LAS source area. The partitioning of available energy into H and LvE varied strongly for different vegetation types, with H varying by a factor of three between senesced winter wheat and grass pasture. The LAS derived H agrees (one-to-one within the experimental uncertainty) with H aggregated from EC with a high coefficient of determination of 0.94. Chronological analysis shows individual fields may have a varying contribution to the areal average sensible heat flux on short (weekly) time scales due to phenological development and changing soil moisture conditions. Using spatially aggregated measurements of net radiation and soil heat flux with H from the LAS, the areal averaged latent heat flux (LvELAS) was calculated as the residual of the surface energy balance. The regression of LvELAS against aggregated LvE from the EC stations has a slope of 0.94, close to ideal, and demonstrates that this is an accurate method for the landscape-scale estimation of evaporation over heterogeneous complex topography.

Wind-direction effects on urban-type flows

Practically all extant work on flows over obstacle arrays, whether laboratory experiments or numerical modelling, is for cases where the oncoming wind is normal to salient faces of the obstacles. In the field, however, this is rarely the case. Here, simulations of flows at various directions over arrays of cubes representing typical urban canopy regions are presented and discussed. The computations are of both direct numerical simulation and large-eddy simulation type. Attention is concentrated on the differences in the mean flow within the canopy region arising from the different wind directions and the consequent effects on global properties such as the total surface drag, which can change very significantly—by up to a factor of three in some circumstances. It is shown that for a given Reynolds number the typical viscous forces are generally a rather larger fraction of the pressure forces (principally the drag) for non-normal than for normal wind directions and that, dependent on the surface morphology, the average flow direction deep within the canopy can be largely independent of the oncoming wind direction. Even for regular arrays of regular obstacles, a wind direction not normal to the obstacle faces can in general generate a lateral lift force (in the direction normal to the oncoming flow). The results demonstrate this and it is shown how computations in a finite domain with the oncoming flow generated by an appropriate forcing term (e.g. a pressure gradient) then lead inevitably to an oncoming wind direction aloft that is not aligned with the forcing term vector.

The wind in the willows: flows in forest canopies in complex terrain

Forest canopies are important components of the terrestrial carbon budget, which has motivated a worldwide effort, FLUXNET, to measure CO2 exchange between forests and the atmosphere. These measurements are difficult to interpret and to scale up to estimate exchange across a landscape. Here we review the effects of complex terrain on the mean flow, turbulence, and scalar exchange in canopy flows, as exemplified by adjustment to forest edges and hills, including the effects of stable stratification. We focus on the fundamental fluid mechanics, in which developments in theory, measurements, and modeling, particularly through large-eddy simulation, are identifying important processes and providing scaling arguments. These developments set the stage for the development of predictive models that can be used in combination with measurements to estimate exchange at the landscape scale.