A Lagrangian model of air-mass photochemistry and mixing using a trajectory ensemble: the Cambridge Tropospheric Trajectory model of Chemistry And Transport (CiTTyCAT) version 4.2

A Lagrangian model of photochemistry and mixing is described (CiTTyCAT, stemming from the Cambridge Tropospheric Trajectory model of Chemistry And Transport), which is suitable for transport and chemistry studies throughout the troposphere. Over the last five years, the model has been developed in parallel at several different institutions and here those developments have been incorporated into one "community" model and documented for the first time. The key photochemical developments include a new scheme for biogenic volatile organic compounds and updated emissions schemes. The key physical development is to evolve composition following an ensemble of trajectories within neighbouring air-masses, including a simple scheme for mixing between them via an evolving "background profile", both within the boundary layer and free troposphere. The model runs along trajectories pre-calculated using winds and temperature from meteorological analyses. In addition, boundary layer height and precipitation rates, output from the analysis model, are interpolated to trajectory points and used as inputs to the mixing and wet deposition schemes. The model is most suitable in regimes when the effects of small-scale turbulent mixing are slow relative to advection by the resolved winds so that coherent air-masses form with distinct composition and strong gradients between them. Such air-masses can persist for many days while stretching, folding and thinning. Lagrangian models offer a useful framework for picking apart the processes of air-mass evolution over inter-continental distances, without being hindered by the numerical diffusion inherent to global Eulerian models. The model, including different box and trajectory modes, is described and some output for each of the modes is presented for evaluation. The model is available for download from a Subversion-controlled repository by contacting the corresponding authors.

Examination of wind storms over Central Europe with respect to circulation weather types and NAO phases

The occurrence of wind storms in Central Europe is investigated with respect to large-scale atmospheric flow and local wind speeds in the investigation area. Two different methods of storm identification are applied for Central Europe as the target region: one based on characteristics of large-scale flow (circulation weather types, CWT) and the other on the occurrence of extreme wind speeds. The identified events are examined with respect to the NAO phases and CWTs under which they occur. Pressure patterns, wind speeds and cyclone tracks are investigated for storms assigned to different CWTs. Investigations are based on ERA40 reanalysis data. It is shown that about 80% of the storm days in Central Europe are connected with westerly flow and that Central European storm events primarily occur during a moderately positive NAO phase, while strongly positive NAO phases (6.4% of all days) account for more than 20% of the storms. A storm occurs over Central Europe during about 10% of the days with a strong positive NAO index. The most frequent pathway of cyclone systems associated with storms over Central Europe leads from the North Atlantic over the British Isles, North Sea and southern Scandinavia into the Baltic Sea. The mean intensity of the systems typically reaches its maximum near the British Isles. Differences between the characteristics for storms identified from the CWT identification procedure (gale days, based on MSLP fields) and those from extreme winds at Central European grid points are small, even though only 70% of the storm days agree. While most storms occur during westerly flow situations, specific characteristics of storms during the other CWTs are also considered. Copyright © 2009 Royal Meteorological Society

Changes in storm track and cyclone activity in three SRES ensemble experiments with the ECHAM5/MPI-OM1 GCM

Synoptic activity over the Northern Hemisphere is evaluated in ensembles of ECHAM5/MPI-OM1 simulations for recent climate conditions (20C) and for three climate scenarios (following SRES A1B, A2, B1). A close agreement is found between the simulations for present day climate and the respective results from reanalysis. Significant changes in the winter mid-tropospheric storm tracks are detected in all three scenario simulations. Ensemble mean climate signals are rather similar, with particularly large activity increases downstream of the Atlantic storm track over Western Europe. The magnitude of this signal is largely dependent on the imposed change in forcing. However, differences between individual ensemble members may be large. With respect to the surface cyclones, the scenario runs produce a reduction in cyclonic track density over the mid-latitudes, even in the areas with increasing mid-tropospheric activity. The largest decrease in track densities occurs at subtropical latitudes, e.g., over the Mediterranean Basin. An increase of cyclone intensities is detected for limited areas (e.g., near Great Britain and Aleutian Isles) for the A1B and A2 experiments. The changes in synoptic activity are associated with alterations of the Northern Hemisphere circulation and background conditions (blocking frequencies, jet stream). The North Atlantic Oscillation index also shows increased values with enhanced forcing. With respect to the effects of changing synoptic activity, the regional change in cyclone intensities is accompanied by alterations of the extreme surface winds, with increasing values over Great Britain, North and Baltic Seas, as well as the areas with vanishing sea ice, and decreases over much of the subtropics.

Analysis of frequency and intensity of European winter storm events from a multi-model perspective, at synoptic and regional scales

This study focuses on the analysis of winter (October-November-December-January-February-March; ONDJFM) storm events and their changes due to increased anthropogenic greenhouse gas concentrations over Europe. In order to assess uncertainties that are due to model formulation, 4 regional climate models (RCMs) with 5 high resolution experiments, and 4 global general circulation models (GCMs) are considered. Firstly, cyclone systems as synoptic scale processes in winter are investigated, as they are a principal cause of the occurrence of extreme, damage-causing wind speeds. This is achieved by use of an objective cyclone identification and tracking algorithm applied to GCMs. Secondly, changes in extreme near-surface wind speeds are analysed. Based on percentile thresholds, the studied extreme wind speed indices allow a consistent analysis over Europe that takes systematic deviations of the models into account. Relative changes in both intensity and frequency of extreme winds and their related uncertainties are assessed and related to changing patterns of extreme cyclones. A common feature of all investigated GCMs is a reduced track density over central Europe under climate change conditions, if all systems are considered. If only extreme (i.e. the strongest 5%) cyclones are taken into account, an increasing cyclone activity for western parts of central Europe is apparent; however, the climate change signal reveals a reduced spatial coherency when compared to all systems, which exposes partially contrary results. With respect to extreme wind speeds, significant positive changes in intensity and frequency are obtained over at least 3 and 20% of the European domain under study (35–72°N and 15°W–43°E), respectively. Location and extension of the affected areas (up to 60 and 50% of the domain for intensity and frequency, respectively), as well as levels of changes (up to +15 and +200% for intensity and frequency, respectively) are shown to be highly dependent on the driving GCM, whereas differences between RCMs when driven by the same GCM are relatively small.

Dependence of winter precipitation over Portugal on NAO and baroclinic wave activity

The relationship between winter (DJF) rainfall over Portugal and the variable large scale circulation is addressed. It is shown that the poles of the sea level pressure (SLP) field variability associated with rainfall variability are shifted about 15° northward with respect to those used in standard definitions of the North Atlantic Oscillation (NAO). It is suggested that the influence of NAO on rainfall dominantly arises from the associated advection of humidity from the Atlantic Ocean. Rainfall is also related to different aspects of baroclinic wave activity, the variability of the latter quantity in turn being largely dependent on the NAO. A negative NAO index (leading to increased westerly surface geostrophic winds into Portugal) is associated with an increased number of deep (ps<980 hPa) surface lows over the central North Atlantic and of intermediate (980

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.

Radiative-dynamical climatology of the first-generation Canadian middle atmosphere model

The Canadian Middle Atmosphere Modelling (MAM) project is a collaboration between thé Atmospheric Environment Service (AES) of Environment Canada and several Canadian universities. Its goal is thé development of a comprehensive General Circulation Model of the troposphere-stratosphere-mesosphere System, starting from the AES/CCCma third-generation atmospheric General Circulation Model. This paper describes the basic features of the first-generation Canadian MAM and some aspects of its radiative-dynamical climatology. Standard first-order mean diagnostics are presented for monthly means and for the annual cycle of zonal-mean winds and temperatures. The mean meridional circulation is examined, and comparison is made between thé steady diabatic, downward controlled, and residual stream functions. It is found that downward control holds quite well in the monthly mean through most of the middle atmosphere, even during equinoctal periods. The relative roles of different drag processes in determining the mean downwelling over the wintertime polar middle stratosphere is examined, and the vertical structure of the drag is quantified.

Lower-stratospheric radiative damping and polar-night jet oscillation events

The effect of stratospheric radiative damping time scales on stratospheric variability and on stratosphere–troposphere coupling is investigated in a simplified global circulation model by modifying the vertical profile of radiative damping in the stratosphere while holding it fixed in the troposphere. Perpetual-January conditions are imposed, with sinusoidal topography of zonal wavenumber 1 or 2. The depth and duration of the simulated sudden stratospheric warmings closely track the lower-stratospheric radiative time scales. Simulations with the most realistic profiles of radiative damping exhibit extended time-scale recoveries analogous to polar-night jet oscillation (PJO) events, which are observed to follow sufficiently deep stratospheric warmings. These events are characterized by weak lower-stratospheric winds and enhanced stability near the tropopause, which persist for up to 3 months following the initial warming. They are obtained with both wave-1 and wave-2 topography. Planetary-scale Eliassen–Palm (EP) fluxes entering the vortex are also suppressed, which is in agreement with observed PJO events. Consistent with previous studies, the tropospheric jets shift equatorward in response to the warmings. The duration of the shift is closely correlated with the period of enhanced stability. The magnitude of the shift in these runs, however, is sensitive only to the zonal wavenumber of the topography. Although the shift is sustained primarily by synoptic-scale eddies, the net effect of the topographic form drag and the planetary-scale fluxes is not negligible; they damp the surface wind response but enhance the vertical shear. The tropospheric response may also reduce the generation of planetary waves, further extending the stratospheric dynamical time scales.

Dynamical control of the mesosphere by orographic and non-orographic gravity wave drag during the extended northern winters of 2006 and 2009

A version of the Canadian Middle Atmosphere Model (CMAM) that is nudged toward reanalysis data up to 1 hPa is used to examine the impacts of parameterized orographic and non-orographic gravity wave drag (OGWD and NGWD) on the zonal-mean circulation of the mesosphere during the extended northern winters of 2006 and 2009 when there were two large stratospheric sudden warmings. The simulations are compared to Aura Microwave Limb Sounder (MLS) observations of mesospheric temperature, carbon monoxide (CO) and derived zonal winds. The control simulation, which uses both OGWD and NGWD, is shown to be in good agreement with MLS. The impacts of OGWD and NGWD are assessed using simulations in which those sources of wave drag are removed. In the absence of OGWD the mesospheric zonal winds in the months preceding the warmings are too strong, causing increased mesospheric NGWD, which drives excessive downwelling, resulting in overly large lower mesospheric values of CO prior to the warming. NGWD is found to be most important following the warmings when the underlying westerlies are too weak to allow much vertical propagation of the orographic gravity waves to the mesosphere. NGWD is primarily responsible for driving the circulation that results in the descent of CO from the thermosphere following the warmings. Zonal mean mesospheric winds and temperatures in all simulations are shown to be strongly constrained by (i.e. slaved to) the stratosphere. Finally, it is demonstrated that the responses to OGWD and NGWD are non-additive due to their dependence and influence on the background winds and temperatures.

What processes drive the ocean heat transport?

The ocean contributes to regulating the Earth’s climate through its ability to transport heat from the equator to the poles. In this study we use long simulations of an ocean model to investigate whether the heat transport is carried primarily by wind-driven gyres or whether it is dominated by deep circulations associated with abyssal mixing and high latitude convection. The heat transport is computed as a function of temperature classes. In the Pacific and Indian ocean, the bulk of the heat transport is associated with wind-driven gyres confined to the thermocline. In the Atlantic, the thermocline gyres account for only 40% of the total heat transport. The remaining 60% is associated with a circulation reaching down to cold waters below the thermocline. Using a series of sensitivity experiments, we show that this deep heat transport is primarily set by the strength and patterns of surface winds and only secondarily by diabatic processes at high latitudes in the North Atlantic. Abyssal mixing below 2000 m has hardly any impact on ocean heat transport. A major implication is that the role of the ocean in regulating Earth’s climate strongly depends on how surface winds change across different climates in both hemispheres at low and high latitudes.

The dependence of southern ocean meridional overturning on wind stress

An eddy-resolving numerical model of a zonal flow, meant to resemble the Antarctic Circumpolar Current, is described and analyzed using the framework of J. Marshall and T. Radko. In addition to wind and buoyancy forcing at the surface, the model contains a sponge layer at the northern boundary that permits a residual meridional overturning circulation (MOC) to exist at depth. The strength of the residual MOC is diagnosed for different strengths of surface wind stress. It is found that the eddy circulation largely compensates for the changes in Ekman circulation. The extent of the compensation and thus the sensitivity of the MOC to the winds depend on the surface boundary condition. A fixed-heat-flux surface boundary severely limits the ability of the MOC to change. An interactive heat flux leads to greater sensitivity. To explain the MOC sensitivity to the wind strength under the interactive heat flux, transformed Eulerian-mean theory is applied, in which the eddy diffusivity plays a central role in determining the eddy response. A scaling theory for the eddy diffusivity, based on the mechanical energy balance, is developed and tested; the average magnitude of the diffusivity is found to be proportional to the square root of the wind stress. The MOC sensitivity to the winds based on this scaling is compared with the true sensitivity diagnosed from the experiments.

Optimising the FAMOUS climate model: inclusion of global carbon cycling

FAMOUS fills an important role in the hierarchy of climate models, both explicitly resolving atmospheric and oceanic dynamics yet being sufficiently computationally efficient that either very long simulations or large ensembles are possible. An improved set of carbon cycle parameters for this model has been found using a perturbed physics ensemble technique. This is an important step towards building the "Earth System" modelling capability of FAMOUS, which is a reduced resolution, and hence faster running, version of the Hadley Centre Climate model, HadCM3. Two separate 100 member perturbed parameter ensembles were performed; one for the land surface and one for the ocean. The land surface scheme was tested against present-day and past representations of vegetation and the ocean ensemble was tested against observations of nitrate. An advantage of using a relatively fast climate model is that a large number of simulations can be run and hence the model parameter space (a large source of climate model uncertainty) can be more thoroughly sampled. This has the associated benefit of being able to assess the sensitivity of model results to changes in each parameter. The climatologies of surface and tropospheric air temperature and precipitation are improved relative to previous versions of FAMOUS. The improved representation of upper atmosphere temperatures is driven by improved ozone concentrations near the tropopause and better upper level winds.

Idealised simulations of sting-jet cyclones

An idealised modelling study of sting-jet cyclones is presented. Sting jets are descending mesoscale jets that occur in some extratropical cyclones and produce localised regions of strong low-level winds in the frontal fracture region. Moist baroclinic lifecycle (LC1) simulations are performed with modifications to produce cyclones resembling observed sting-jet cyclones. A sting jet exists in the idealised control cyclone with similar characteristics to the sting jet in a simulation of windstorm Gudrun (a confirmed sting-jet case). Unlike in windstorm Gudrun, a low-level layer of strong moist static stability prohibits the descent of the strong winds from above the boundary layer to the surface in the idealised case. Conditional symmetric instability (CSI) exists in the cloud head and dissipates as the sting jet leaves the cloud head and descends. The descending, initially moist, sting-jet trajectories consistently have negative or near-zero saturated moist potential vorticity but moist static stability and inertial stability, consistent with CSI release; the moist static stability becomes negative during the period of most rapid descent, by which time the air is relatively dry implying conditional instability release is unlikely. Sensitivity experiments show that the existence of the sting jet is robust to changes in the initial state, and that the initial tropospheric static stability significantly impacts the descent rate of the sting jet. Inertial and conditional instability are probably being released in the experiment with the weakest initial static stability. This suggests that sting jets can arise through the release of all three instabilities associated with negative saturated moist potential vorticity. While evaporative cooling occurs along the sting-jet trajectories, a sensitivity experiment with evaporation effects turned off shows no significant change to the wind strength or descent rate of the sting jet implying that instability release is the dominant sting-jet driving mechanism.

Wave activity in the tropical tropopause layer in seven reanalysis and four chemistry climate model data sets

Sub-seasonal variability including equatorial waves significantly influence the dehydration and transport processes in the tropical tropopause layer (TTL). This study investigates the wave activity in the TTL in 7 reanalysis data sets (RAs; NCEP1, NCEP2, ERA40, ERA-Interim, JRA25, MERRA, and CFSR) and 4 chemistry climate models (CCMs; CCSRNIES, CMAM, MRI, and WACCM) using the zonal wave number-frequency spectral analysis method with equatorially symmetric-antisymmetric decomposition. Analyses are made for temperature and horizontal winds at 100 hPa in the RAs and CCMs and for outgoing longwave radiation (OLR), which is a proxy for convective activity that generates tropopause-level disturbances, in satellite data and the CCMs. Particular focus is placed on equatorial Kelvin waves, mixed Rossby-gravity (MRG) waves, and the Madden-Julian Oscillation (MJO). The wave activity is defined as the variance, i.e., the power spectral density integrated in a particular zonal wave number-frequency region. It is found that the TTL wave activities show significant difference among the RAs, ranging from ∼0.7 (for NCEP1 and NCEP2) to ∼1.4 (for ERA-Interim, MERRA, and CFSR) with respect to the averages from the RAs. The TTL activities in the CCMs lie generally within the range of those in the RAs, with a few exceptions. However, the spectral features in OLR for all the CCMs are very different from those in the observations, and the OLR wave activities are too low for CCSRNIES, CMAM, and MRI. It is concluded that the broad range of wave activity found in the different RAs decreases our confidence in their validity and in particular their value for validation of CCM performance in the TTL, thereby limiting our quantitative understanding of the dehydration and transport processes in the TTL.

Helga's dune: 40 years of dune dynamics in the Namib Desert

The dynamics of a small linear dune on the northern margin of the Namib Sand Sea have been monitored using erosion pins placed at the dune tip since 1969. GPS measurements of these pins enabled estimation of the rates of advance and lateral migration of the dune. The average rate of advance of the dune tip over the period 1969–2012 was 1.99 m yr–1 towards 015°. Rates of advance and lateral movement varied over the period of monitoring, with a decrease in rates of advance by a factor of 50%, but an increase in the rate of lateral movement. Changes in dune behavior appear to be related to changes in wind regime and the vegetation cover of the interdune area, as a result of increased rainfall in recent years. This study demonstrates the dynamic nature of the tip of this dune and its sensitivity to changes in winds and sand supply.