Quantification of chemical and physical processes influencing ozone during long-range transport using a trajectory ensemble

During long-range transport, many distinct processes – including photochemistry, deposition, emissions and mixing – contribute to the transformation of air mass composition. Partitioning the effects of different processes can be useful when considering the sensitivity of chemical transformation to, for example, a changing environment or anthropogenic influence. However, transformation is not observed directly, since mixing ratios are measured, and models must be used to relate changes to processes. Here, four cases from the ITCT-Lagrangian 2004 experiment are studied. In each case, aircraft intercepted a distinct air mass several times during transport over the North Atlantic, providing a unique dataset and quantifying the net changes in composition from all processes. A new framework is presented to deconstruct the change in O3 mixing ratio (Δ O3) into its component processes, which were not measured directly, taking into account the uncertainty in measurements, initial air mass variability and its time evolution. The results show that the net chemical processing (Δ O3chem) over the whole simulation is greater than net physical processing (Δ O3phys) in all cases. This is in part explained by cancellation effects associated with mixing. In contrast, each case is in a regime of either net photochemical destruction (lower tropospheric transport) or production (an upper tropospheric biomass burning case). However, physical processes influence O3 indirectly through addition or removal of precursor gases, so that changes to physical parameters in a model can have a larger effect on Δ O3chem than Δ O3phys. Despite its smaller magnitude, the physical processing distinguishes the lower tropospheric export cases, since the net photochemical O3 change is −5 ppbv per day in all three cases. Processing is quantified using a Lagrangian photochemical model with a novel method for simulating mixing through an ensemble of trajectories and a background profile that evolves with them. The model is able to simulate the magnitude and variability of the observations (of O3, CO, NOy and some hydrocarbons) and is consistent with the time-average OH following air-masses inferred from hydrocarbon measurements alone (by Arnold et al., 2007). Therefore, it is a useful new method to simulate air mass evolution and variability, and its sensitivity to process parameters.

Mechanisms for tropical upwelling in the stratosphere

The dynamics of the tropical upwelling branch of the stratospheric Brewer–Dobson circulation are examined, with a particular focus on the role of the middle-atmosphere Hadley circulation. Upwelling is examined in terms of both the diabatic circulation and Lagrangian trajectories using a zonally symmetric balance model. The behavior of the wave-driven circulation in the presence of angular momentum redistribution by the Hadley circulation is also considered. The results of the zonally symmetric model are compared with fields from a middle-atmosphere GCM. It is found that the Hadley circulation makes a significant contribution to annual mean tropical upwelling at the upwelling maximum in the vicinity of the stratopause, and can account for most of the annual mean upwelling seen in the GCM. In the mid- to lower stratosphere, the role of the Hadley circulation is much weaker and wave drag appears to be required to explain the observed upwelling, although the Hadley circulation makes a nonnegligible contribution to the annual cycle of the upwelling. Subtropical wave drag can produce annual mean upwelling through a nonlinear mechanism; viscosity is not required. However, the magnitude of the observed upwelling suggests that wave drag must penetrate quite close to the equator.

On the nature of large-scale mixing in the stratosphere and mesosphere

Studies of tracer transport in the stratosphere have shown that adiabatic quasi-horizontal tracer evolution is controlled primarily by the large-scale low-frequency component of the flow. This behavior is consistent with the concept of chaotic advection, wherein the Eulerian velocity field is spatially coherent and temporally quasi-regular on timescales over which the Lagrangian evolution is chaotic. In this study, winds from a middle atmosphere general circulation model (the Canadian Middle Atmosphere Model) are used to compare and contrast the nature of tracer evolution in the stratosphere and mesosphere. It is found that the concept of chaotic advection is relevant in the stratosphere but not in the mesosphere. The explanation for this behavior is the increased strength of gravity wave activity in the mesosphere as compared with the stratosphere, which leads to shallower kinetic energy spectra on synoptic scales and a much shorter Eulerian correlation time. The shallower kinetic energy spectra imply that tracer evolution in the mesosphere is spectrally local, in contrast with the spectrally nonlocal regime that prevails in the stratosphere. This means that tracer advection calculations in the mesosphere are controlled primarily by the gravity wave spectrum and are intrinsically resolution dependent.

A closer look at chaotic advection in the stratosphere: part I: geometric structure

The relevance of chaotic advection to stratospheric mixing and transport is addressed in the context of (i) a numerical model of forced shallow-water flow on the sphere, and (ii) a middle-atmosphere general circulation model. It is argued that chaotic advection applies to both these models if there is suitable large-scale spatial structure in the velocity field and if the velocity field is temporally quasi-regular. This spatial structure is manifested in the form of “cat’s eyes” in the surf zone, such as are commonly seen in numerical simulations of Rossby wave critical layers; by analogy with the heteroclinic structure of a temporally aperiodic chaotic system the cat’s eyes may be thought of as an “organizing structure” for mixing and transport in the surf zone. When this organizing structure exists, Eulerian and Lagrangian autocorrelations of the velocity derivatives indicate that velocity derivatives decorrelate more rapidly along particle trajectories than at fixed spatial locations (i.e., the velocity field is temporally quasi-regular). This phenomenon is referred to as Lagrangian random strain.

A general method for finding extremal states of Hamiltonian dynamical systems, with applications to perfect fluids

In addition to the Hamiltonian functional itself, non-canonical Hamiltonian dynamical systems generally possess integral invariants known as ‘Casimir functionals’. In the case of the Euler equations for a perfect fluid, the Casimir functionals correspond to the vortex topology, whose invariance derives from the particle-relabelling symmetry of the underlying Lagrangian equations of motion. In a recent paper, Vallis, Carnevale & Young (1989) have presented algorithms for finding steady states of the Euler equations that represent extrema of energy subject to given vortex topology, and are therefore stable. The purpose of this note is to point out a very general method for modifying any Hamiltonian dynamical system into an algorithm that is analogous to those of Vallis etal. in that it will systematically increase or decrease the energy of the system while preserving all of the Casimir invariants. By incorporating momentum into the extremization procedure, the algorithm is able to find steadily-translating as well as steady stable states. The method is applied to a variety of perfect-fluid systems, including Euler flow as well as compressible and incompressible stratified flow.

An exact local conservation theorem for finite-amplitude disturbances to non-parallel shear flows, with remarks on Hamiltonian structure and on Arnol'd's stability theorems

Disturbances of arbitrary amplitude are superposed on a basic flow which is assumed to be steady and either (a) two-dimensional, homogeneous, and incompressible (rotating or non-rotating) or (b) stably stratified and quasi-geostrophic. Flow over shallow topography is allowed in either case. The basic flow, as well as the disturbance, is assumed to be subject neither to external forcing nor to dissipative processes like viscosity. An exact, local ‘wave-activity conservation theorem’ is derived in which the density A and flux F are second-order ‘wave properties’ or ‘disturbance properties’, meaning that they are O(a2) in magnitude as disturbance amplitude a [rightward arrow] 0, and that they are evaluable correct to O(a2) from linear theory, to O(a3) from second-order theory, and so on to higher orders in a. For a disturbance in the form of a single, slowly varying, non-stationary Rossby wavetrain, $\overline{F}/\overline{A}$ reduces approximately to the Rossby-wave group velocity, where (${}^{-}$) is an appropriate averaging operator. F and A have the formal appearance of Eulerian quantities, but generally involve a multivalued function the correct branch of which requires a certain amount of Lagrangian information for its determination. It is shown that, in a certain sense, the construction of conservable, quasi-Eulerian wave properties like A is unique and that the multivaluedness is inescapable in general. The connection with the concepts of pseudoenergy (quasi-energy), pseudomomentum (quasi-momentum), and ‘Eliassen-Palm wave activity’ is noted. The relationship of this and similar conservation theorems to dynamical fundamentals and to Arnol'd's nonlinear stability theorems is discussed in the light of recent advances in Hamiltonian dynamics. These show where such conservation theorems come from and how to construct them in other cases. An elementary proof of the Hamiltonian structure of two-dimensional Eulerian vortex dynamics is put on record, with explicit attention to the boundary conditions. The connection between Arnol'd's second stability theorem and the suppression of shear and self-tuning resonant instabilities by boundary constraints is discussed, and a finite-amplitude counterpart to Rayleigh's inflection-point theorem noted

Mean motions induced by baroclinic instability in a jet

A study is made of the zonal-mean motions induced by a growing baroclinic wave in several contexts, under the framework of three different analysis schemes: the conventional Eulerian mean (EM), the transformed Eulerian mean (TEM), and the generalized Lagrangian mean (GLM). The effect of meridional shear in the initial jet on these induced mean motions is considered by treating the instability problem in the context of the two-layer model. The conceptual simplicity of the TEM formulation is shown to be useful in diagnosing the dynamics of instability, much as it has been found helpful in many problems of wave, mean-flow interaction. In addition, it is found that the TEM vertical velocity is a very good indicator of the GLM vertical velocity. However, the GLM meridional velocity is always convergent towards the centre of instability activity, and is not at all well represented by the nondivergent TEM meridional velocity. In comparing the results with Uryu's (1979) calculation of the GLM circulation induced by a growing Eady wave, it is found that the inclusion of meridional jet shear in the present work leads to some strikingly different effects in the GLM zonal wind acceleration. In the case of pure baroclinic instability treated by Uryu, the Eulerian and Stokes accelerations nearly cancel each other in the centre of the channel, leaving a weak Lagrangian acceleration opposed to the Eulerian one. In the more general case of mixed baroclinic-barotropic instability, however, the Eulerian and Stokes accelerations can reinforce one another, leading to a very strong Lagrangian zonal wind

Evidence for changes in stratospheric transport and mixing over the past three decades based on multiple data sets and tropical leaky pipe analysis

Variability in the strength of the stratospheric Lagrangian mean meridional or Brewer-Dobson circulation and horizontal mixing into the tropics over the past three decades are examined using observations of stratospheric mean age of air and ozone. We use a simple representation of the stratosphere, the tropical leaky pipe (TLP) model, guided by mean meridional circulation and horizontal mixing changes in several reanalyses data sets and chemistry climate model (CCM) simulations, to help elucidate reasons for the observed changes in stratospheric mean age and ozone. We find that the TLP model is able to accurately simulate multiyear variability in ozone following recent major volcanic eruptions and the early 2000s sea surface temperature changes, as well as the lasting impact on mean age of relatively short-term circulation perturbations. We also find that the best quantitative agreement with the observed mean age and ozone trends over the past three decades is found assuming a small strengthening of the mean circulation in the lower stratosphere, a moderate weakening of the mean circulation in the middle and upper stratosphere, and a moderate increase in the horizontal mixing into the tropics. The mean age trends are strongly sensitive to trends in the horizontal mixing into the tropics, and the uncertainty in the mixing trends causes uncertainty in the mean circulation trends. Comparisons of the mean circulation and mixing changes suggested by the measurements with those from a recent suite of CCM runs reveal significant differences that may have important implications on the accurate simulation of future stratospheric climate.

Transport timescales and tracer properties in the extratropical UTLS

A comprehensive evaluation of seasonal backward trajectories initialized in the northern hemisphere lowermost stratosphere (LMS) has been performed to investigate the factors that determine the temporal and spatial structure of troposphere-to-stratosphere-transport (TST) and it's impact on the LMS. In particular we explain the fundamental role of the transit time since last TST (tTST) for the chemical composition of the LMS. According to our results the structure of the LMS can be characterized by a layer with tTST<40 days forming a narrow band around the local tropopause. This layer extends about 30 K above the local dynamical tropopause, corresponding to the extratropical tropopause transition layer (ExTL) as identified by CO. The LMS beyond this layer shows a relatively well defined separation as marked by an aprupt transition to longer tTST indicating less frequent mixing and a smaller fraction of tropospheric air. Thus the LMS constitutes a region of two well defined regimes of tropospheric influence. These can be characterized mainly by different transport times from the troposphere and different fractions of tropospheric air. Carbon monoxide (CO) mirrors this structure of tTST due to it's finite lifetime on the order of three months. Water vapour isopleths, on the other hand, do not uniquely indicate TST and are independent of tTST, but are determined by the Lagrangian Cold Point (LCP) of air parcels. Most of the backward trajectories from the LMS experienced their LCP in the tropics and sub-tropics, and TST often occurs 20 days after trajectories have encountered their LCP. Therefore, ExTL properties deduced from CO and H2O provide totally different informations on transport and particular TST for the LMS.

Explosive development of winter storm Xynthia over the subtropical North Atlantic Ocean

In winter of 2009–2010 south-western Europe was hit by several destructive windstorms. The most important was Xynthia (26–28 February 2010), which caused 64 reported casualties and was classified as the 2nd most expensive natural hazard event for 2010 in terms of economic losses. In this work we assess the synoptic evolution, dynamical characteristics and the main impacts of storm Xynthia, whose genesis, development and path were very uncommon. Wind speed gusts observed at more than 500 stations across Europe are evaluated as well as the wind gust field obtained with a regional climate model simulation for the entire North Atlantic and European area. Storm Xynthia was first identified on 25 February around 30° N, 50° W over the subtropical North Atlantic Ocean. Its genesis occurred on a region characterized by warm and moist air under the influence of a strong upper level wave embedded in the westerlies. Xynthia followed an unusual SW–NE path towards Iberia, France and central Europe. The role of moist air masses on the explosive development of Xynthia is analysed by considering the evaporative sources. A lagrangian model is used to identify the moisture sources, sinks and moisture transport associated with the cyclone during its development phase. The main supply of moisture is located over an elongated region of the subtropical North Atlantic Ocean with anomalously high SST, confirming that the explosive development of storm Xynthia had a significant contribution from the subtropics.

Solutions of the fully compressible semi-geostrophic system

The fully compressible semi-geostrophic system is widely used in the modelling of large-scale atmospheric flows. In this paper, we prove rigorously the existence of weak Lagrangian solutions of this system, formulated in the original physical coordinates. In addition, we provide an alternative proof of the earlier result on the existence of weak solutions of this system expressed in the so-called geostrophic, or dual, coordinates. The proofs are based on the optimal transport formulation of the problem and on recent general results concerning transport problems posed in the Wasserstein space of probability measures.

Using coordinated observations in polarized white light and Faraday rotation to probe the spatial position and magnetic field of an interplanetary sheath

Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons and therefore depend on both viewing geometry and electron density. The Faraday rotation (FR) of radio waves from extragalactic pulsars and quasars, which arises due to the presence of such solar wind features, depends on the line-of-sight magnetic field component B ∥ and the electron density. To understand coordinated WL and FR observations of CMEs, we perform forward magnetohydrodynamic modeling of an Earth-directed shock and synthesize the signatures that would be remotely sensed at a number of widely distributed vantage points in the inner heliosphere. Removal of the background solar wind contribution reveals the shock-associated enhancements in WL and FR. While the efficiency of Thomson scattering depends on scattering angle, WL radiance I decreases with heliocentric distance r roughly according to the expression Ir –3. The sheath region downstream of the Earth-directed shock is well viewed from the L4 and L5 Lagrangian points, demonstrating the benefits of these points in terms of space weather forecasting. The spatial position of the main scattering site r sheath and the mass of plasma at that position M sheath can be inferred from the polarization of the shock-associated enhancement in WL radiance. From the FR measurements, the local B ∥sheath at r sheath can then be estimated. Simultaneous observations in polarized WL and FR can not only be used to detect CMEs, but also to diagnose their plasma and magnetic field properties.

Aircraft routing with minimal climate impact: the REACT4C climate cost function modelling approach (V1.0)

In addition to CO2, the climate impact of aviation is strongly influenced by non-CO2 emissions, such as nitrogen oxides, influencing ozone and methane, and water vapour, which can lead to the formation of persistent contrails in ice-supersaturated regions. Because these non-CO2 emission effects are characterised by a short lifetime, their climate impact largely depends on emission location and time; that is to say, emissions in certain locations (or times) can lead to a greater climate impact (even on the global average) than the same emission in other locations (or times). Avoiding these climate-sensitive regions might thus be beneficial to climate. Here, we describe a modelling chain for investigating this climate impact mitigation option. This modelling chain forms a multi-step modelling approach, starting with the simulation of the fate of emissions released at a certain location and time (time-region grid points). This is performed with the chemistry–climate model EMAC, extended via the two submodels AIRTRAC (V1.0) and CONTRAIL (V1.0), which describe the contribution of emissions to the composition of the atmosphere and to contrail formation, respectively. The impact of emissions from the large number of time-region grid points is efficiently calculated by applying a Lagrangian scheme. EMAC also includes the calculation of radiative impacts, which are, in a second step, the input to climate metric formulas describing the global climate impact of the emission at each time-region grid point. The result of the modelling chain comprises a four-dimensional data set in space and time, which we call climate cost functions and which describes the global climate impact of an emission at each grid point and each point in time. In a third step, these climate cost functions are used in an air traffic simulator (SAAM) coupled to an emission tool (AEM) to optimise aircraft trajectories for the North Atlantic region. Here, we describe the details of this new modelling approach and show some example results. A number of sensitivity analyses are performed to motivate the settings of individual parameters. A stepwise sanity check of the results of the modelling chain is undertaken to demonstrate the plausibility of the climate cost functions.

Application of complex extreme learning machine to multiclass classification problems with high dimensionality: a THz spectra classification problem

We extend extreme learning machine (ELM) classifiers to complex Reproducing Kernel Hilbert Spaces (RKHS) where the input/output variables as well as the optimization variables are complex-valued. A new family of classifiers, called complex-valued ELM (CELM) suitable for complex-valued multiple-input–multiple-output processing is introduced. In the proposed method, the associated Lagrangian is computed using induced RKHS kernels, adopting a Wirtinger calculus approach formulated as a constrained optimization problem similarly to the conventional ELM classifier formulation. When training the CELM, the Karush–Khun–Tuker (KKT) theorem is used to solve the dual optimization problem that consists of satisfying simultaneously smallest training error as well as smallest norm of output weights criteria. The proposed formulation also addresses aspects of quaternary classification within a Clifford algebra context. For 2D complex-valued inputs, user-defined complex-coupled hyper-planes divide the classifier input space into four partitions. For 3D complex-valued inputs, the formulation generates three pairs of complex-coupled hyper-planes through orthogonal projections. The six hyper-planes then divide the 3D space into eight partitions. It is shown that the CELM problem formulation is equivalent to solving six real-valued ELM tasks, which are induced by projecting the chosen complex kernel across the different user-defined coordinate planes. A classification example of powdered samples on the basis of their terahertz spectral signatures is used to demonstrate the advantages of the CELM classifiers compared to their SVM counterparts. The proposed classifiers retain the advantages of their ELM counterparts, in that they can perform multiclass classification with lower computational complexity than SVM classifiers. Furthermore, because of their ability to perform classification tasks fast, the proposed formulations are of interest to real-time applications.

The slowly evolving background state of the atmosphere

The theory of wave–mean flow interaction requires a partition of the atmospheric flow into a notional background state and perturbations to it. Here, a background state, known as the Modified Lagrangian Mean (MLM), is defined as the zonally symmetric state obtained by requiring that every potential vorticity (PV) contour lying within an isentropic layer encloses the same mass and circulation as in the full flow. For adiabatic and frictionless flow, these two integral properties are time-invariant and the MLM state is a steady solution of the primitive equations. The time dependence in the adiabatic flow is put into the perturbations, which can be described by a wave-activity conservation law that is exact even at large amplitude. Furthermore, the effects of non-conservative processes on wave activity can be calculated from the conservation law. A new method to calculate the MLM state is introduced, where the position of the lower boundary is obtained as part of the solution. The results are illustrated using Northern Hemisphere ERA-Interim data. The MLM state evolves slowly, implying that the net non-conservative effects are weak. Although ‘adiabatic eddy fluxes’ cannot affect the MLM state, the effects of Rossby-wave breaking, PV filamentation and subsequent dissipation result in sharpening of the polar vortex edge and meridional shifts in the MLM zonal flow, both at tropopause level and on the winter stratospheric vortex. The rate of downward migration of wave activity during stratospheric sudden warmings is shown to be given by the vertical scale associated with polar vortex tilt divided by the time-scale for wave dissipation estimated from the wave-activity conservation law. Aspects of troposphere–stratosphere interaction are discussed. The new framework is suitable to examine the climate and its interactions with disturbances, such as midlatitude storm tracks, and makes a clean partition between adiabatic and non-conservative processes.