A review of urban roughness sublayer turbulence

It is becoming increasingly important that we can understand and model flow processes in urban areas. Applications such as weather forecasting, air quality and sustainable urban development rely on accurate modelling of the interface between an urban surface and the atmosphere above. This review gives an overview of current understanding of turbulence generated by an urban surface up to a few building heights, the layer called the roughness sublayer (RSL). High quality datasets are also identified which can be used in the development of suitable parameterisations of the urban RSL. Datasets derived from physical and numerical modelling, and full-scale observations in urban areas now exist across a range of urban-type morphologies (e.g. street canyons, cubes, idealised and realistic building layouts). Results show that the urban RSL depth falls within 2 – 5 times mean building height and is not easily related to morphology. Systematic perturbations away from uniform layouts (e.g. varying building heights) have a significant impact on RSL structure and depth. Considerable fetch is required to develop an overlying inertial sublayer, where turbulence is more homogeneous, and some authors have suggested that the “patchiness” of urban areas may prevent inertial sublayers from developing at all. Turbulence statistics suggest similarities between vegetation and urban canopies but key differences are emerging. There is no consensus as to suitable scaling variables, e.g. friction velocity above canopy vs. square root of maximum Reynolds stress, mean vs. maximum building height. The review includes a summary of existing modelling practices and highlights research priorities.

Note: A balloon-borne accelerometer technique for measuring atmospheric turbulence

A weather balloon and its suspended instrument package behave like a pendulum with a moving pivot. This dynamical system is exploited here for the detection of atmospheric turbulence. By adding an accelerometer to the instrument package, the size of the swings induced by atmospheric turbulence can be measured. In test flights, strong turbulence has induced accelerations greater than 5g, where g = 9.81 m s−2. Calibration of the accelerometer data with a vertically orientated lidar has allowed eddy dissipation rate values of between 10−3 and 10−2 m2 s−3 to be derived from the accelerometer data. The novel use of a whole weather balloon and its adapted instrument package can be used as a new instrument to make standardized in situ measurements of turbulence.

Rapid changes in the redox conditions of the western Tethys Ocean during the early Aptian oceanic anoxic event

The early Aptian (125 to 121 Ma) records an episode of severe environmental change including a major perturbation of the carbon cycle, an oceanic anoxic event (OAE 1a, 122.5 Ma), a platform drowning episode and a biocalcification crisis. We propose to trace changes in the oxygenation state of the ocean during the early Aptian anoxic event using the redox-sensitive trace-element (RSTE) distribution, phosphorus accumulation rates (PARs) and organic-matter characterization in three different basins of the western Tethys. The following sections have been investigated: Gorgo a Cerbara (central Italy) in the Umbria Marche basin, Glaise (SE France) in the Vocontian basin and Cassis/La Bédoule (SE France) located in the Provencal basin. In the Gorgo a Cerbara section, RSTE distributions show a low background level along the main part of the section, contrasted by different maxima in concentrations within the Selli level. In the Glaise section, the Goguel level displays a weak increase in RSTE contents coeval with moderate TOC values. At Cassis/La Bédoule, no significant RSTE enrichments have been observed in sediments equivalent to the Selli level. These differences in the records of the geochemical proxies of the Selli level or its equivalent indicate the deposition under different redox conditions, probably related to the paleogeography. Our data indicate the development of anoxic–euxinic conditions in the deeper part of the Tethys during OAE 1a, whereas in the shallower environments, conditions were less reducing. Moreover, at Gorgo a Cerbara, the Selli level is characterized by rapid changes in the intensity of reducing conditions in the water column. Ocean eutrophication seems to be a major factor in the development and the persistence of anoxia as suggested by the PAR evolution. Higher PAR values at the onset of OAE 1a suggest an increase in nutrient input, whereas the return to lower values through the first part of the OAE 1a interval may be related to the weakened capacity to retain P in the sedimentary reservoir due to bottom-water oxygen depletion. This general pattern is contrasted by the data of Gorgo a Cerbara, where the sediments deposited during the OAE 1a interval show P-enrichments (mainly authigenic P). This is associated with maxima in TOC values and Corg:Ptot ratios, suggesting that a part of the remobilized P was trapped in the sediments and as such prevented from returning to the water column.

The expression of the Cenomanian-Turonian oceanic anoxic event in Tibet

The Gongzha section of Tibet, China is located at the northern margin of the Indian Plate (SE Tethys) and is characterized by hemipelagic grey marls and marly limestones, light grey limestones and silty limestones, but no organic-rich sediments. High-resolution biostratigraphy reveals an expanded Cenomanian–Turonian (CT) boundary interval and the δ13C record includes the main features of the classical positive carbon-isotope excursion that characterizes the CT oceanic anoxic event. The biotic response inferred from the foraminifera suggests that oxic to dysoxic conditions prevailed, except for a short interval marked by peak abundance of Heterohelix that indicates a significantly dysoxic environment during the δ13C “b” peak excursion. The overall decreasing trend in redox-sensitive trace elements (RSTE) during the maximum δ13C excursion confirms the absence of significant longer-lasting anoxia in the Gongzha section. Enrichments in RSTE are linked to phases of increased detrital input. Chemical weathering indices suggest that the upper Cenomanian sediments accumulated under an increasingly hot and humid climate that culminated near the CT boundary. In the early Turonian lower weathering indices suggest a warm, drier climatic regime with reduced continental runoff. Phosphorus mass-accumulation rates show a significant peak at the onset of the positive δ13C excursion, followed by a decrease up to the basal Turonian. This pattern is positively correlated with the long-term decrease in detrital index as also observed in numerous other CT boundary sections (e.g., Eastbourne, Pueblo, and Whadi El Ghaib, Sinaï). Long-term phosphorus accumulation in the Gongzha section is therefore associated with changes in detrital input. The overall decreased detrital input can be explained by the increasingly remote continental sources due to the major transgression at the end of Cenomanian, coupled with changes in continental weathering intensity linked to increasingly more arid climate conditions.

Observational constraints on atmospheric and oceanic cross-equatorial heat transports: revisiting the precipitation asymmetry problem in climate models

Satellite based top-of-atmosphere (TOA) and surface radiation budget observations are combined with mass corrected vertically integrated atmospheric energy divergence and tendency from reanalysis to infer the regional distribution of the TOA, atmospheric and surface energy budget terms over the globe. Hemispheric contrasts in the energy budget terms are used to determine the radiative and combined sensible and latent heat contributions to the cross-equatorial heat transports in the atmosphere (AHT_EQ) and ocean (OHT_EQ). The contrast in net atmospheric radiation implies an AHT_EQ from the northern hemisphere (NH) to the southern hemisphere (SH) (0.75 PW), while the hemispheric difference in sensible and latent heat implies an AHT_EQ in the opposite direction (0.51 PW), resulting in a net NH to SH AHT_EQ (0.24 PW). At the surface, the hemispheric contrast in the radiative component (0.95 PW) dominates, implying a 0.44 PW SH to NH OHT_EQ. Coupled model intercomparison project phase 5 (CMIP5) models with excessive net downward surface radiation and surface-to-atmosphere sensible and latent heat transport in the SH relative to the NH exhibit anomalous northward AHT_EQ and overestimate SH tropical precipitation. The hemispheric bias in net surface radiative flux is due to too much longwave surface radiative cooling in the NH tropics in both clear and all-sky conditions and excessive shortwave surface radiation in the SH subtropics and extratropics due to an underestimation in reflection by clouds.

Atmospheric and oceanic contributions to irreducible forecast uncertainty of Arctic surface climate

Uncertainty of Arctic seasonal to interannual predictions arising from model errors and initial state uncertainty has been widely discussed in the literature, whereas the irreducible forecast uncertainty (IFU) arising from the chaoticity of the climate system has received less attention. However, IFU provides important insights into the mechanisms through which predictability is lost, and hence can inform prioritization of model development and observations deployment. Here, we characterize how internal oceanic and surface atmospheric heat fluxes contribute to IFU of Arctic sea ice and upper ocean heat content in an Earth system model by analyzing a set of idealized ensemble prediction experiments. We find that atmospheric and oceanic heat flux are often equally important for driving unpredictable Arctic-wide changes in sea ice and surface water temperatures, and hence contribute equally to IFU. Atmospheric surface heat flux tends to dominate Arctic-wide changes for lead times of up to a year, whereas oceanic heat flux tends to dominate regionally and on interannual time scales. There is in general a strong negative covariance between surface heat flux and ocean vertical heat flux at depth, and anomalies of lateral ocean heat transport are wind-driven, which suggests that the unpredictable oceanic heat flux variability is mainly forced by the atmosphere. These results are qualitatively robust across different initial states, but substantial variations in the amplitude of IFU exist. We conclude that both atmospheric variability and the initial state of the upper ocean are key ingredients for predictions of Arctic surface climate on seasonal to interannual time scales.

The impact of oceanic heat transport on the atmospheric circulation

A general circulation model of intermediate complexity with an idealized Earth-like aquaplanet setup is used to study the impact of changes in the oceanic heat transport on the global atmospheric circulation. Focus is on the atmospheric mean meridional circulation and global thermodynamic properties. The atmosphere counterbalances to a large extent the imposed changes in the oceanic heat transport, but, nonetheless, significant modifications to the atmospheric general circulation are found. Increasing the strength of the oceanic heat transport up to 2.5 PW leads to an increase in the global mean near-surface temperature and to a decrease in its equator-to-pole gradient. For stronger transports, the gradient is reduced further, but the global mean remains approximately constant. This is linked to a cooling and a reversal of the temperature gradient in the tropics. Additionally, a stronger oceanic heat transport leads to a decline in the intensity and a poleward shift of the maxima of both the Hadley and Ferrel cells. Changes in zonal mean diabatic heating and friction impact the properties of the Hadley cell, while the behavior of the Ferrel cell is mostly controlled by friction. The efficiency of the climate machine, the intensity of the Lorenz energy cycle and the material entropy production of the system decline with increased oceanic heat transport. This suggests that the climate system becomes less efficient and turns into a state of reduced entropy production as the enhanced oceanic transport performs a stronger large-scale mixing between geophysical fluids with different temperatures, thus reducing the available energy in the climate system and bringing it closer to a state of thermal equilibrium.

First-order turbulence closure for modelling complex canopy flows

Simple first-order closure remains an attractive way of formulating equations for complex canopy flows when the aim is to find analytic or simple numerical solutions to illustrate fundamental physical processes. Nevertheless, the limitations of such closures must be understood if the resulting models are to illuminate rather than mislead. We propose five conditions that first-order closures must satisfy then test two widely used closures against them. The first is the eddy diffusivity based on a mixing length. We discuss the origins of this approach, its use in simple canopy flows and extensions to more complex flows. We find that it satisfies most of the conditions and, because the reasons for its failures are well understood, it is a reliable methodology. The second is the velocity-squared closure that relates shear stress to the square of mean velocity. Again we discuss the origins of this closure and show that it is based on incorrect physical principles and fails to satisfy any of the five conditions in complex canopy flows; consequently its use can lead to actively misleading conclusions.

Langmuir turbulence and surface heating in the ocean surface boundary layer

This study uses large-eddy simulation to investigate the structure of the ocean surface boundary layer (OSBL) in the presence of Langmuir turbulence and stabilizing surface heat fluxes. The OSBL consists of a weakly stratified layer, despite a surface heat flux, above a stratified thermocline. The weakly stratified (mixed) layer is maintained by a combination of a turbulent heat flux produced by the wave-driven Stokes drift and downgradient turbulent diffusion. The scaling of turbulence statistics, such as dissipation and vertical velocity variance, is only affected by the surface heat flux through changes in the mixed layer depth. Diagnostic models are proposed for the equilibrium boundary layer and mixed layer depths in the presence of surface heating. The models are a function of the initial mixed layer depth before heating is imposed and the Langmuir stability length. In the presence of radiative heating, the models are extended to account for the depth profile of the heating.

Impact of the North Atlantic Oscillation on transatlantic flight routes and clear-air turbulence

The variation of wind-optimal transatlantic flight routes and their turbulence potential is investigated to understand how upper-level winds and large-scale flow patterns can affect the efficiency and safety of long-haul flights. In this study, the wind-optimal routes (WORs) that minimize the total flight time by considering wind variations are modeled for flights between John F. Kennedy International Airport (JFK) in New York, New York, and Heathrow Airport (LHR) in London, United Kingdom, during two distinct winter periods of abnormally high and low phases of North Atlantic Oscillation (NAO) teleconnection patterns. Eastbound WORs approximate the JFK–LHR great circle (GC) route following northerly shifted jets in the +NAO period. Those WORs deviate southward following southerly shifted jets during the −NAO period, because eastbound WORs fly closely to the prevailing westerly jets to maximize tailwinds. Westbound WORs, however, spread meridionally to avoid the jets near the GC in the +NAO period to minimize headwinds. In the −NAO period, westbound WORs are north of the GC because of the southerly shifted jets. Consequently, eastbound WORs are faster but have higher probabilities of encountering clear-air turbulence than westbound ones, because eastbound WORs are close to the jet streams, especially near the cyclonic shear side of the jets in the northern (southern) part of the GC in the +NAO (−NAO) period. This study suggests how predicted teleconnection weather patterns can be used for long-haul strategic flight planning, ultimately contributing to minimizing aviation’s impact on the environment

VELOCITY FIELD OF COMPRESSIBLE MAGNETOHYDRODYNAMIC TURBULENCE: WAVELET DECOMPOSITION AND MODE SCALINGS

We study compressible magnetohydrodynamic turbulence, which holds the key to many astrophysical processes, including star formation and cosmic-ray propagation. To account for the variations of the magnetic field in the strongly turbulent fluid, we use wavelet decomposition of the turbulent velocity field into Alfven, slow, and fast modes, which presents an extension of the Cho & Lazarian decomposition approach based on Fourier transforms. The wavelets allow us to follow the variations of the local direction of the magnetic field and therefore improve the quality of the decomposition compared to the Fourier transforms, which are done in the mean field reference frame. For each resulting component, we calculate the spectra and two-point statistics such as longitudinal and transverse structure functions as well as higher order intermittency statistics. In addition, we perform a Helmholtz-Hodge decomposition of the velocity field into incompressible and compressible parts and analyze these components. We find that the turbulence intermittency is different for different components, and we show that the intermittency statistics depend on whether the phenomenon was studied in the global reference frame related to the mean magnetic field or in the frame defined by the local magnetic field. The dependencies of the measures we obtained are different for different components of the velocity; for instance, we show that while the Alfven mode intermittency changes marginally with the Mach number, the intermittency of the fast mode is substantially affected by the change.

Blowout bifurcation and spatial mode excitation in the bubbling transition to turbulence

The transition to turbulence (spatio-temporal chaos) in a wide class of spatially extended dynamical system is due to the loss of transversal stability of a chaotic attractor lying on a homogeneous manifold (in the Fourier phase space of the system) causing spatial mode excitation Since the latter manifests as intermittent spikes this has been called a bubbling transition We present numerical evidences that this transition occurs due to the so called blowout bifurcation whereby the attractor as a whole loses transversal stability and becomes a chaotic saddle We used a nonlinear three-wave interacting model with spatial diffusion as an example of this transition (C) 2010 Elsevier B V All rights reserved

Low-dimensional chaos and wave turbulence in plasmas

We investigated drift-wave turbulence in the plasma edge of a small tokamak by considering solutions of the Hasegawa-Mima equation involving three interacting modes in Fourier space. The resulting low-dimensional dynamics presented periodic as well as chaotic evolution of the Fourier-mode amplitudes, and we performed the control of chaotic behaviour through the application of a fourth resonant wave of small amplitude.

Generation of Vortices and Observation of Quantum Turbulence in an Oscillating Bose-Einstein Condensate

We report on the experimental observation of vortex formation and production of tangled vortex distribution in an atomic BEC of (87)Rb atoms submitted to an external oscillatory perturbation. The oscillatory perturbations start by exciting quadrupolar and scissors modes of the condensate. Then regular vortices are observed finally evolving to a vortex tangle configuration. The vortex tangle is a signature of the presence of a turbulent regime in the cloud. We also show that this turbulent cloud has suppression of the aspect ratio inversion typically observed in quantum degenerate bosonic gases during free expansion.

Transition to quantum turbulence in finite-size superfluids

A novel concept of quantum turbulence in finite size superfluids, such as trapped bosonic atoms, is discussed. We have used an atomic (87)Rb Bose-Einstein condensate (BEC) to study the emergence of this phenomenon. In our experiment, the transition to the quantum turbulent regime is characterized by a tangled vortex lines formation, controlled by the amplitude and time duration of the excitation produced by an external oscillating field. A simple model is suggested to account for the experimental observations. The transition from the non-turbulent to the turbulent regime is a rather gradual crossover. But it takes place in a sharp enough way, allowing for the definition of an effective critical line separating the regimes. Quantum turbulence emerging in a finite-size superfluid may be a new idea helpful for revealing important features associated to turbulence, a more general and broad phenomenon. [GRAPHICS] Amplitude versus elapsed time diagram of magnetically excited BEC superfluid, presenting the evolution from the non-turbulent regime, with well separated vortices, to the turbulent regimes, with tangled vortices (C) 2011 by Astro Ltd. Published exclusively by WILEY-VCH Verlag GmbH & Co. KGaA