The GCM–Reality Intercomparison Project for SPARC (GRIPS): scientific issues and initial results

To investigate the effects of the middle atmosphere on climate, the World Climate Research Programme is supporting the project "Stratospheric Processes and their Role in Climate" (SPARC). A central theme of SPARC, to examine model simulations of the coupled troposphere—middle atmosphere system, is being performed through the initiative called GRIPS (GCM—Reality Intercomparison Project for SPARC). In this paper, an overview of the objectives of GRIPS is given. Initial activities include an assessment of the performance of middle atmosphere climate models, and preliminary results from this evaluation are presented here. It is shown that although all 13 models evaluated represent most major features of the mean atmospheric state, there are deficiencies in the magnitude and location of the features, which cannot easily be traced to the formulation (resolution or the parameterizations included) of the models. Most models show a cold bias in all locations, apart from the tropical tropopause region where they can be either too warm or too cold. The strengths and locations of the major jets are often misrepresented in the models. Looking at three—dimensional fields reveals, for some models, more severe deficiencies in the magnitude and positioning of the dominant structures (such as the Aleutian high in the stratosphere), although undersampling might explain some of these differences from observations. All the models have shortcomings in their simulations of the present—day climate, which might limit the accuracy of predictions of the climate response to ozone change and other anomalous forcing.

On the 'downward control' of extratropical diabatic circulations by eddy-induced mean zonal forces

The situation considered is that of a zonally symmetric model of the middle atmosphere subject to a given quasi-steady zonal force F̄, conceived to be the result of irreversible angular momentum transfer due to the upward propagation and breaking of Rossby and gravity waves together with any other dissipative eddy effects that may be relevant. The model's diabatic heating is assumed to have the qualitative character of a relaxation toward some radiatively determined temperature field. To the extent that the force F̄ may be regarded as given, and the extratropical angular momentum distribution is realistic, the extratropical diabatic mass flow across a given isentropic surface may be regarded as controlled exclusively by the F̄ distribution above that surface (implying control by the eddy dissipation above that surface and not, for instance, by the frequency of tropopause folding below). This “downward control” principle expresses a critical part of the dynamical chain of cause and effect governing the average rate at which photochemical products like ozone become available for folding into, or otherwise descending into, the extratropical troposphere. The dynamical facts expressed by the principle are also relevant, for instance, to understanding the seasonal-mean rate of upwelling of water vapor to the summer mesopause, and the interhemispheric differences in stratospheric tracer transport. The robustness of the principle is examined when F̄ is time-dependent. For a global-scale, zonally symmetric diabatic circulation with a Brewer-Dobson-like horizontal structure given by the second zonally symmetric Hough mode, with Rossby height HR = 13 km in an isothermal atmosphere with density scale height H = 7 km, the vertical partitioning of the unsteady part of the mass circulation caused by fluctuations in F̄ confined to a shallow layer LF̄ is always at least 84% downward. It is 90% downward when the force fluctuates sinusoidally on twice the radiative relaxation timescale and 95% if five times slower. The time-dependent adjustment when F̄ is changed suddenly is elucidated, extending the work of Dickinson (1968), when the atmosphere is unbounded above and below. Above the forcing, the adjustment is characterized by decay of the meridional mass circulation cell at a rate proportional to the radiative relaxation rate τr−1 divided by {1 + (4H2/HR2)}. This decay is related to the boundedness of the angular momentum that can be taken up by the finite mass of air above LF̄ without causing an ever-increasing departure from thermal wind balance. Below the forcing, the meridional mass circulation cell penetrates downward at a speed τr−1 HR2/H. For the second Hough mode, the time for downward penetration through one density scale height is about 6 days if the radiative relaxation time is 20 days, the latter being representative of the lower stratosphere. At any given altitude, a steady state is approached. The effect of a rigid lower boundary on the time-dependent adjustment is also considered. If a frictional planetary boundary layer is present then a steady state is ultimately approached everywhere, with the mass circulation extending downward from LF̄ and closing via the boundary layer. Satellite observations of temperature and ozone are used in conjunction with a radiative transfer scheme to estimate the altitudes from which the lower stratospheric diabatic vertical velocity is controlled by the effective F̄ in the real atmosphere. The data appear to indicate that about 80% of the effective control is usually exerted from below 40 km but with significant exceptions up to 70 km (in the high latitude southern hemispheric winter). The implications for numerical modelling of chemical transport are noted.

The importance of surface pressure changes in the response of the atmosphere to zonally-symmetric thermal and mechanical forcing

he classical problem of the response of a balanced, axisymmetric vortex to thermal and mechanical forcing is re-examined, paying special attention to the lower boundary condition. The correct condition is DΦ/Dt = 0, where Φ is the geopotential and D/Dt the material derivative, which explicitly accounts for a mass redistribution as part of the mean-flow response. This redistribution is neglected when using the boundary condition Dp/Dt = 0, which has conventionally been applied in this problem. It is shown that applying the incorrect boundary condition, and thereby ignoring the surface pressure change, leads to a zonal wind acceleration δū/δt that is too strong, especially near the surface. The effect is significant for planetary-scale forcing even when applied at tropopause level. A comparison is made between the mean-flow evolution in a baroclinic life-cycle, as simulated in a fully nonlinear, primitive-equation model, and that predicted by using the simulated eddy fluxes in the zonally-symmetric response problem. Use of the correct lower boundary condition is shown to lead to improved agreement.

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.

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.

Tracing troposphere-to-stratosphere transport above a mid-latitude deep convective system

Within the project SPURT (trace gas measurements in the tropopause region) a variety of trace gases have been measured in situ in order to investigate the role of dynamical and chemical processes in the extra-tropical tropopause region. In this paper we report on a flight on 10 November 2001 leading from Hohn, Germany (52�N) to Faro, Portugal (37�N) through a strongly developed deep stratospheric intrusion. This streamer was associated with a large convective system over the western Mediterranean with potentially significant troposphere-to-stratosphere transport. Along major parts of the flight we measured unexpectedly high NOy mixing ratios. Also H2O mixing ratios were significantly higher than stratospheric background levels confirming the extraordinary chemical signature of the probed air masses in the interior of the streamer. Backward trajectories encompassing the streamer enable to analyze the origin and physical characteristics of the air masses and to trace troposphere-to-stratosphere transport. Near the western flank of the streamer features caused by long range transport, such as tropospheric filaments characterized by sudden drops in the O3 and NOy mixing ratios and enhanced CO and H2O can be reconstructed in great detail using the reverse domain filling technique. These filaments indicate a high potential for subsequent mixing with the stratospheric air. At the south-western edge of the streamer a strong gradient in the NOy and the O3 mixing ratios coincides very well with a sharp gradient in potential vorticity in the ECMWF fields. In contrast, in the interior of the streamer the observed highly elevated NOy and H2O mixing ratios up to a potential temperature level of 365K and potential vorticity values of maximum 10 PVU cannot be explained in terms of resolved troposphere-to-stratosphere transport along the backward trajectories. Also mesoscale simulations with a High Resolution Model reveal no direct evidence for convective H2O injection up to this level. Elevated H2O mixing ratios in the ECMWF and HRM are seen only up to about tropopause height at 340 hPa and 270 hPa, respectively, well below flight altitude of about 200 hPa. However, forward tracing of the convective influence as identified by satellite brightness temperature measurements and counts of lightning strokes shows that during this part of the flight the aircraft was closely following the border of an air mass which was heavily impacted by convective activity over Spain and Algeria. This is evidence that deep convection at mid-latitudes may have a large impact on the tracer distribution of the lowermost stratosphere reaching well above the thunderstorms anvils as claimed by recent studies using cloud-resolving models.

Jet characterization in the upper troposphere/lower stratosphere (UTLS): applications to climatology and transport studies

A method of classifying the upper tropospheric/lower stratospheric (UTLS) jets has been developed that allows satellite and aircraft trace gas data and meteorological fields to be efficiently mapped in a jet coordinate view. A detailed characterization of multiple tropopauses accompanies the jet characterization. Jet climatologies show the well-known high altitude subtropical and lower altitude polar jets in the upper troposphere, as well as a pattern of concentric polar and subtropical jets in the Southern Hemisphere, and shifts of the primary jet to high latitudes associated with blocking ridges in Northern Hemisphere winter. The jet-coordinate view segregates air masses differently than the commonly-used equivalent latitude (EqL) coordinate throughout the lowermost stratosphere and in the upper troposphere. Mapping O3 data from the Aura Microwave Limb Sounder (MLS) satellite and the Winter Storms aircraft datasets in jet coordinates thus emphasizes different aspects of the circulation compared to an EqL-coordinate framework: the jet coordinate reorders the data geometrically, thus highlighting the strong PV, tropopause height and trace gas gradients across the subtropical jet, whereas EqL is a dynamical coordinate that may blur these spatial relationships but provides information on irreversible transport. The jet coordinate view identifies the concentration of stratospheric ozone well below the tropopause in the region poleward of and below the jet core, as well as other transport features associated with the upper tropospheric jets. Using the jet information in EqL coordinates allows us to study trace gas distributions in regions of weak versus strong jets, and demonstrates weaker transport barriers in regions with less jet influence. MLS and Atmospheric Chemistry Experiment-Fourier Transform Spectrometer trace gas fields for spring 2008 in jet coordinates show very strong, closely correlated, PV, tropopause height and trace gas gradients across the jet, and evidence of intrusions of stratospheric air below the tropopause below and poleward of the subtropical jet; these features are consistent between instruments and among multiple trace gases. Our characterization of the jets is facilitating studies that will improve our understanding of upper tropospheric trace gas evolution.

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.

Satellite observations and modeling of transport in the upper troposphere through the lower mesosphere during the 2006 major stratospheric sudden warming

An unusually strong and prolonged stratospheric sudden warming (SSW) in January 2006 was the first major SSW for which globally distributed long-lived trace gas data are available covering the upper troposphere through the lower mesosphere. We use Aura Microwave Limb Sounder (MLS), Atmospheric Chemistry Experiment-Fourier Transform Spectrometer (ACE-FTS) data, the SLIMCAT Chemistry Transport Model (CTM), and assimilated meteorological analyses to provide a comprehensive picture of transport during this event. The upper tropospheric ridge that triggered the SSW was associated with an elevated tropopause and layering in trace gas profiles in conjunction with stratospheric and tropospheric intrusions. Anomalous poleward transport (with corresponding quasi-isentropic troposphere-to-stratosphere exchange at the lowest levels studied) in the region over the ridge extended well into the lower stratosphere. In the middle and upper stratosphere, the breakdown of the polar vortex transport barrier was seen in a signature of rapid, widespread mixing in trace gases, including CO, H2O, CH4 and N2O. The vortex broke down slightly later and more slowly in the lower than in the middle stratosphere. In the middle and lower stratosphere, small remnants with trace gas values characteristic of the pre-SSW vortex lingered through the weak and slow recovery of the vortex. The upper stratospheric vortex quickly reformed, and, as enhanced diabatic descent set in, CO descended into this strong vortex, echoing the fall vortex development. Trace gas evolution in the SLIMCAT CTM agrees well with that in the satellite trace gas data from the upper troposphere through the middle stratosphere. In the upper stratosphere and lower mesosphere, the SLIMCAT simulation does not capture the strong descent of mesospheric CO and H2O values into the reformed vortex; this poor CTM performance in the upper stratosphere and lower mesosphere results primarily from biases in the diabatic descent in assimilated analyses.

Measurements of NO, NOy, N2O, and O3 during SPURT: implications for transport and chemistry in the lowermost stratosphere

We present measurements of NO, NOy, O3, and N2O within the lowermost stratosphere (LMS) over Europe obtained during the SPURT project. The measurements cover all seasons between November 2001 and July 2003. They span a broad band of latitudes from 30° N to 75° N and a potential temperature range from 290 to 380 K. The measurements represent a comprehensive data set of these tracers and reveal atmospheric transport processes that influence tracer distributions in the LMS. Median mixing ratios of stratospheric tracers in equivalent latitude-potential temperature coordinates show a clear seasonal cycle related to the Brewer-Dobson circulation, with highest values in spring and lowest values in autumn. Vertical tracer profiles show strong gradients at the extratropical tropopause, suggesting that vertical (cross-isentropic) mixing is reduced above the tropopause. Pronounced meridional gradients in the tracer mixing ratios are found on potential temperature surfaces in the LMS. This suggests strongly reduced mixing along isentropes. Concurrent large gradients in static stability in the vertical direction, and of PV in the meridional direction, suggest the presence of a mixing barrier. Seasonal cycles were found in the correlation slopes ΔO3/ΔN2O and ΔNOy/ΔN2O well above the tropopause. Absolute slope values are smallest in spring indicating chemically aged stratospheric air originating from high altitudes and latitudes. Larger values were measured in summer and autumn suggesting that a substantial fraction of air takes a "short-cut" from the tropical tropopause region into the extratropical LMS. The seasonal change in the composition of the LMS has direct implications for the ozone chemistry in this region. Comparisons of measured NO with the critical NO value at which net ozone production changes from negative to positive, imply ozone production up to 20 K above the local tropopause in spring, up to 30 K in summer, and up to 40 K in autumn. Above these heights, and in winter, net ozone production is negative.

Seasonal cycles and variability of O3 and H2O in the UT/LMS during SPURT

Airborne high resolution in situ measurements of a large set of trace gases including ozone (O3) and total water (H2O) in the upper troposphere and the lowermost stratosphere (UT/LMS) have been performed above Europe within the SPURT project. SPURT provides an extensive data coverage of the UT/LMS in each season within the time period between November 2001 and July 2003. In the LMS a distinct spring maximum and autumn minimum is observed in O3, whereas its annual cycle in the UT is shifted by 2–3 months later towards the end of the year. The more variable H2O measurements reveal a maximum during summer and a minimum during autumn/winter with no phase shift between the two atmospheric compartments. For a comprehensive insight into trace gas composition and variability in the UT/LMS several statistical methods are applied using chemical, thermal and dynamical vertical coordinates. In particular, 2-dimensional probability distribution functions serve as a tool to transform localised aircraft data to a more comprehensive view of the probed atmospheric region. It appears that both trace gases, O3 and H2O, reveal the most compact arrangement and are best correlated in the view of potential vorticity (PV) and distance to the local tropopause, indicating an advanced mixing state on these surfaces. Thus, strong gradients of PV seem to act as a transport barrier both in the vertical and the horizontal direction. The alignment of trace gas isopleths reflects the existence of a year-round extra-tropical tropopause transition layer. The SPURT measurements reveal that this layer is mainly affected by stratospheric air during winter/spring and by tropospheric air during autumn/summer. Normalised mixing entropy values for O3 and H2O in the LMS appear to be maximal during spring and summer, respectively, indicating highest variability of these trace gases during the respective seasons.

Highly resolved observations of trace gases in the lowermost stratosphere and upper troposphere from the SPURT project: an overview

During SPURT (Spurenstofftransport in der Tropopausenregion, trace gas transport in the tropopause region) we performed measurements of a wide range of trace gases with different lifetimes and sink/source characteristics in the northern hemispheric upper troposphere (UT) and lowermost stratosphere (LMS). A large number of in-situ instruments were deployed on board a Learjet 35A, flying at altitudes up to 13.7 km, at times reaching to nearly 380 K potential temperature. Eight measurement campaigns (consisting of a total of 36 flights), distributed over all seasons and typically covering latitudes between 35° N and 75° N in the European longitude sector (10° W–20° E), were performed. Here we present an overview of the project, describing the instrumentation, the encountered meteorological situations during the campaigns and the data set available from SPURT. Measurements were obtained for N2O, CH4, CO, CO2, CFC12, H2, SF6, NO, NOy, O3 and H2O. We illustrate the strength of this new data set by showing mean distributions of the mixing ratios of selected trace gases, using a potential temperature-equivalent latitude coordinate system. The observations reveal that the LMS is most stratospheric in character during spring, with the highest mixing ratios of O3 and NOy and the lowest mixing ratios of N2O and SF6. The lowest mixing ratios of NOy and O3 are observed during autumn, together with the highest mixing ratios of N2O and SF6 indicating a strong tropospheric influence. For H2O, however, the maximum concentrations in the LMS are found during summer, suggesting unique (temperature- and convection-controlled) conditions for this molecule during transport across the tropopause. The SPURT data set is presently the most accurate and complete data set for many trace species in the LMS, and its main value is the simultaneous measurement of a suite of trace gases having different lifetimes and physical-chemical histories. It is thus very well suited for studies of atmospheric transport, for model validation, and for investigations of seasonal changes in the UT/LMS, as demonstrated in accompanying and elsewhere published studies.

Model simulations and aircraft measurements of vertical, seasonal and latitudinal O3 and CO distributions over Europe

During a series of 8 measurement campaigns within the SPURT project (2001-2003), vertical profiles of CO and O3 have been obtained at subtropical, middle and high latitudes over western Europe, covering the troposphere and lowermost stratosphere up to ~14 km altitude during all seasons. The seasonal and latitudinal variation of the measured trace gas profiles are compared to simulations with the chemical transport model MATCH. In the troposphere reasonable agreement between observations and model predictions is achieved for CO and O3, in particular at subtropical and mid-latitudes, while the model overestimates (underestimates) CO (O3 in the lowermost stratosphere particularly at high latitudes, indicating too strong simulated bi-directional exchange across the tropopause. By the use of tagged tracers in the model, long-range transport of Asian air masses is identified as the dominant source of CO pollution over Europe in the free troposphere.

Determination of eddy-diffusivity in the lowermost stratosphere

We present a 2D-advection-diffusion model that simulates the main transport pathways influencing tracer distributions in the lowermost stratosphere (LMS). The model describes slow diabatic descent of aged stratospheric air, vertical (cross-isentropic) and horizontal (along isentropes) diffusion within the LMS and across the tropopause using equivalent latitude and potential temperature coordinates. Eddy diffusion coefficients parameterize the integral effect of dynamical processes leading to small scale turbulence and mixing. They were specified by matching model simulations to observed CO distributions. Interestingly, the model suggests mixing across isentropes to be more important than horizontal mixing across surfaces of constant equivalent latitude, shining new light on the interplay between various transport mechanisms in the LMS. The model achieves a good description of the small scale tracer features at the tropopause with squared correlation coefficients R2 = 0.72…0.94.

Seasonality and extent of extratropical TST derived from in-situ CO measurements during SPURT

We present airborne in-situ trace gas measurements which were performed on eight campaigns between November 2001 and July 2003 during the SPURT-project (SPURenstofftransport in der Tropopausenregion, trace gas transport in the tropopause region). The measurements on a quasi regular basis allowed an overview of the seasonal variations of the trace gas distribution in the tropopause region over Europe from 35°-75°N to investigate the influence of transport and mixing across the extratropical tropopause on the lowermost stratosphere. From the correlation of CO and O3 irreversible mixing of tropospheric air into the lowermost stratosphere is identified. The CO distribution indicates that transport and subsequent mixing of tropospheric air across the extratropical tropopause predominantly affects a layer, which closely follows the shape of the local tropopause. In addition, the seasonal cycle of CO2 illustrates the strong coupling of that layer to the extratropical troposphere. Both, horizontal gradients of CO on isentropes as well as the CO-O3-distribution in the lowermost stratosphere reveal that the influence of quasi-horizontal transport and subsequent mixing weakens with distance from the local tropopause. The mixing layer extends to about 25 K in potential temperature above the local tropopause exhibiting only a weak seasonality. However, at large distances from the tropopause a significant influence of tropospheric air is still evident. The relation between N2O and CO2 indicates that a significant contribution of air originating from the tropical tropopause contributes to the background air in the extratropical lowermost stratosphere.