The interactive responses of water quality and hydrology to changes in multiple stressors, and implications for the long-term effective management of phosphorus

Soluble reactive phosphorus (SRP) plays a key role in eutrophication, a global problem decreasing habitat quality and in-stream biodiversity. Mitigation strategies are required to prevent SRP fluxes from exceeding critical levels, and must be robust in the face of potential changes in climate, land use and a myriad of other influences. To establish the longevity of these strategies it is therefore crucial to consider the sensitivity of catchments to multiple future stressors. This study evaluates how the water quality and hydrology of a major river system in the UK (the River Thames) respond to alterations in climate, land use and water resource allocations, and investigates how these changes impact the relative performance of management strategies over an 80-year period. In the River Thames, the relative contributions of SRP from diffuse and point sources vary seasonally. Diffuse sources of SRP from agriculture dominate during periods of high runoff, and point sources during low flow periods. SRP concentrations rose under any future scenario which either increased a) surface runoff or b) the area of cultivated land. Under these conditions, SRP was sourced from agriculture, and the most effective single mitigation measures were those which addressed diffuse SRP sources. Conversely, where future scenarios reduced flow e.g. during winters of reservoir construction, the significance of point source inputs increased, and mitigation measures addressing these issues became more effective. In catchments with multiple point and diffuse sources of SRP, an all-encompassing effective mitigation approach is difficult to achieve with a single strategy. In order to attain maximum efficiency, multiple strategies might therefore be employed at different times and locations, to target the variable nature of dominant SRP sources and pathways.

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.

Comments on some recent measurements of anomalously steep N2O and O3tracer spectra in the stratospheric surf zone

Recent aircraft measurements, primarily in the extratropics, of the horizontal variance of nitrous oxide (N2O) and ozone (O3) in the middle stratosphere indicate that horizontal spectra of the tracer variance scale nearly as k−2, where k is the spatial wavenumber along the aircraft flight track [Strahan and Mahlman, 1994; Bacmeister et al., 1996]. This spectral scaling has been regarded as inconsistent with the accepted picture of stratospheric tracer motion; large-scale quasi-two-dimensional tracer advection typically yields a k−1 scaling (i.e., the classical Batchelor spectrum). In this paper it is argued that the nearly k−2 scaling seen in the measurements is a natural outcome of quasi-two-dimensional filamentation of the polar vortex edge. The accepted picture of stratospheric tracer motion can thus be retained: no additional physical processes are needed to account for deviations from the Batchelor spectrum. Our argument is based on the finite lifetime of tracer filaments and on the “singularity spectrum” associated with a one-dimensional field composed of randomly spaced jumps in concentration.

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.

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.

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.

In-situ comparison of the NOy instruments flown in MOZAIC and SPURT

Two aircraft instruments for the measurement of total odd nitrogen (NOy) were compared side by side aboard a Learjet A35 in April 2003 during a campaign of the AFO2000 project SPURT (Spurengastransport in der Tropopausenregion). The instruments albeit employing the same measurement principle (gold converter and chemiluminescence) had different inlet configurations. The ECO-Physics instrument operated by ETH-Zürich in SPURT had the gold converter mounted outside the aircraft, whereas the instrument operated by FZ-Jülich in the European project MOZAIC III (Measurements of ozone, water vapour, carbon monoxide and nitrogen oxides aboard Airbus A340 in-service aircraft) employed a Rosemount probe with 80 cm of FEP-tubing connecting the inlet to the gold converter. The NOy concentrations during the flight ranged between 0.3 and 3 ppb. The two data sets were compared in a blind fashion and each team followed its normal operating procedures. On average, the measurements agreed within 7%, i.e. within the combined uncertainty of the two instruments. This puts an upper limit on potential losses of HNO3 in the Rosemount inlet of the MOZAIC instrument. Larger transient deviations were observed during periods after calibrations and when the aircraft entered the stratosphere. The time lag of the MOZAIC instrument observed in these instances is in accordance with the time constant of the MOZAIC inlet line determined in the laboratory for HNO3.

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.

Impact of vegetation and preferential source areas on global dust aerosol: results from a model study

[1] We present a model of the dust cycle that successfully predicts dust emissions as determined by land surface properties, monthly vegetation and snow cover, and 6-hourly surface wind speeds for the years 1982–1993. The model takes account of the role of dry lake beds as preferential source areas for dust emission. The occurrence of these preferential sources is determined by a water routing and storage model. The dust source scheme also explicitly takes into account the role of vegetation type as well as monthly vegetation cover. Dust transport is computed using assimilated winds for the years 1987–1990. Deposition of dust occurs through dry and wet deposition, where subcloud scavenging is calculated using assimilated precipitation fields. Comparison of simulated patterns of atmospheric dust loading with the Total Ozone Mapping Spectrometer satellite absorbing aerosol index shows that the model produces realistic results from daily to interannual timescales. The magnitude of dust deposition agrees well with sediment flux data from marine sites. Emission of submicron dust from preferential source areas are required for the computation of a realistic dust optical thickness. Sensitivity studies show that Asian dust source strengths are particularly sensitive to the seasonality of vegetation cover.

Elevated atmospheric CO2 and humidity delay leaf fall in Betula pendula, but not in Alnus glutinosa or Populus tremula × tremuloides

Context: Anthropogenic activity has increased the level of atmospheric CO2, which is driving an increase of global temperatures and associated changes in precipitation patterns. At Northern latitudes, one of the likely consequences of global warming is increased precipitation and air humidity. Aims: In this work, the effects of both elevated atmospheric CO2 and increased air humidity on trees commonly growing in northern European forests were assessed. Methods: The work was carried out under field conditions by using Free Air Carbon dioxide Enrichment (FACE) and Free Air Humidity Manipulation (FAHM) systems. Leaf litter fall was measured over 4 years (FACE) or 5 years (FAHM) to determine the effects of FACE and FAHM on leaf phenology. Results: Increasing air humidity delayed leaf litter fall in Betula pendula, but not in Populus tremula × tremuloides. Similarly, under elevated atmospheric CO2, leaf litter fall was delayed in Betula pendula, but not in Alnus glutinosa. Increased CO2 appeared to interact with periods of low precipitation in summer and high ozone levels during these periods to effect leaf fall. Conclusions: This work shows that increased CO2 and humidity delay leaf fall, but this effect is species specific.

SPARC Data Initiative: A comparison of ozone climatologies from international satellite limb sounders

A comprehensive quality assessment of the ozone products from 18 limb-viewing satellite instruments is provided by means of a detailed intercomparison. The ozone climatologies in form of monthly zonal mean time series covering the upper troposphere to lower mesosphere are obtained from LIMS, SAGE I/II/III, UARS-MLS, HALOE, POAM II/III, SMR, OSIRIS, MIPAS, GOMOS, SCIAMACHY, ACE-FTS, ACE-MAESTRO, Aura-MLS, HIRDLS, and SMILES within 1978–2010. The intercomparisons focus on mean biases of annual zonal mean fields, interannual variability, and seasonal cycles. Additionally, the physical consistency of the data is tested through diagnostics of the quasi-biennial oscillation and Antarctic ozone hole. The comprehensive evaluations reveal that the uncertainty in our knowledge of the atmospheric ozone mean state is smallest in the tropical and midlatitude middle stratosphere with a 1σ multi-instrument spread of less than ±5%. While the overall agreement among the climatological data sets is very good for large parts of the stratosphere, individual discrepancies have been identified, including unrealistic month-to-month fluctuations, large biases in particular atmospheric regions, or inconsistencies in the seasonal cycle. Notable differences between the data sets exist in the tropical lower stratosphere (with a spread of ±30%) and at high latitudes (±15%). In particular, large relative differences are identified in the Antarctic during the time of the ozone hole, with a spread between the monthly zonal mean fields of ±50%. The evaluations provide guidance on what data sets are the most reliable for applications such as studies of ozone variability, model-measurement comparisons, detection of long-term trends, and data-merging activities.

Reconciliation of halogen-induced ozone loss with the total-column ozone record

The observed depletion of the ozone layer from the 1980s onwards is attributed to halogen source gases emitted by human activities. However, the precision of this attribution is complicated by year-to-year variations in meteorology, that is, dynamical variability, and by changes in tropospheric ozone concentrations. As such, key aspects of the total-column ozone record, which combines changes in both tropospheric and stratospheric ozone, remain unexplained, such as the apparent absence of a decline in total-column ozone levels before 1980, and of any long-term decline in total-column ozone levels in the tropics. Here we use a chemistry–climate model to estimate changes in halogen-induced ozone loss between 1960 and 2010; the model is constrained by observed meteorology to remove the effects of dynamical variability, and driven by emissions of tropospheric ozone precursors to separate out changes in tropospheric ozone. We show that halogen-induced ozone loss closely followed stratospheric halogen loading over the studied period. Pronounced enhancements in ozone loss were apparent in both hemispheres following the volcanic eruptions of El Chichon and, in particular, Mount Pinatubo, which significantly enhanced stratospheric aerosol loads. We further show that approximately 40% of the long-term non-volcanic ozone loss occurred before 1980, and that long-term ozone loss also occurred in the tropical stratosphere. Finally, we show that halogen-induced ozone loss has declined by over 10% since stratospheric halogen loading peaked in the late 1990s, indicating that the recovery of the ozone layer is well underway.

The climate impact of past changes in halocarbons and CO2in the tropical UTLS region

A chemistry-climate model coupled to an ocean model is used to compare the climate impact of past (1960-2010) changes in concentrations of halocarbons with those of CO2 in the tropical upper troposphere and lower stratosphere. The halocarbon contribution to both upper troposphere warming and the associated increase in lower stratospheric upwelling is about 40% as large as that due to CO2. Trends in cold-point temperature and lower stratosphere water vapor are positive for both halocarbons and CO2, and are of about the same magnitude. Trends in lower stratosphere ozone are negative, due to the increased upwelling. These increases in water vapor and decreases in lower stratosphere ozone feed back on lower stratosphere temperature through radiative cooling. The radiative cooling from ozone is about a factor of two larger than that from water vapor in the vicinity of the cold-point tropopause, while water vapor dominates at heights above 50 hPa. For halocarbons this indirect radiative cooling more than offsets the direct radiative warming, and together with the adiabatic cooling accounts for the lack of a halocarbon-induced warming of the lower stratosphere. For CO2 the indirect cooling from increased water vapor and decreased ozone is of comparable magnitude to the direct warming from CO2 in the vicinity of the cold-point tropopause, and (together with the increased upwelling) lowers the height at which CO2 increases induce stratospheric cooling, thus explaining the relatively weak increase in cold-point temperature due to the CO2 increases.

Meteorology, air quality, and health in London: the ClearfLo project

The ClearfLo project provides integrated measurements of the meteorology, composition and particulate loading of London's urban atmosphere to improve predictive capability for air quality. Air quality and heat are strong health drivers and their accurate assessment and forecast are important in densely populated urban areas. However, the sources and processes leading to high concentrations of main pollutants such as ozone, nitrogen dioxide, and fine and coarse particulate matter in complex urban areas are not fully understood, limiting our ability to forecast air quality accurately. This paper introduces the ClearfLo project's interdisciplinary approach to investigate the processes leading to poor air quality and elevated temperatures. Within ClearfLo (www.clearflo.ac.uk), a large multi-institutional project funded by the UK Natural Environment Research Council (NERC), integrated measurements of meteorology, gaseous and particulate composition/loading within London's atmosphere were undertaken to understand the processes underlying poor air quality. Long-term measurement infrastructure installed at multiple levels (street and elevated), and at urban background, kerbside and rural locations were complemented with high-resolution numerical atmospheric simulations . Combining these (measurement/modeling) enhances understanding of seasonal variations in meteorology and composition together with the controlling processes. Two intensive observation periods (winter 2012 and summer Olympics 2012) focus upon the vertical structure and evolution of the urban boundary layer, chemical controls on nitrogen dioxide and ozone production, in particular the role of volatile organic compounds, and processes controlling the evolution, size, distribution and composition of particulate matter. The paper shows that mixing heights are deeper over London than in the rural surroundings and the seasonality of the urban boundary layer evolution controls when concentrations peak. The composition also reflects the seasonality of sources such as domestic burning and biogenic emissions.

Advances in understanding mineral dust and boundary layer processes over the Sahara from Fennec aircraft observations

The Fennec climate program aims to improve understanding of the Saharan climate system through a synergy of observations and modelling. We present a description of the Fennec airborne observations during 2011 and 2012 over the remote Sahara (Mauritania and Mali) and the advances in the understanding of mineral dust and boundary layer processes they have provided. Aircraft instrumentation aboard the UK FAAM BAe146 and French SAFIRE Falcon 20 is described, with specific focus on instrumentation specially developed and relevant to Saharan meteorology and dust. Flight locations, aims and associated meteorology are described. Examples and applications of aircraft measurements from the Fennec flights are presented, highlighting new scientific results delivered using a synergy of different instruments and aircraft. These include: (1) the first airborne measurement of dust particles sized up to 300 microns and associated dust fluxes in the Saharan atmospheric boundary layer (SABL), (2) dust uplift from the breakdown of the nocturnal low-level jet before becoming visible in SEVIRI satellite imagery, (3) vertical profiles of the unique vertical structure of turbulent fluxes in the SABL, (4) in-situ observations of processes in SABL clouds showing dust acting as CCN and IN at −15 °C, (5) dual-aircraft observations of the SABL dynamics, thermodynamics and composition in the Saharan heat low region (SHL), (6) airborne observations of a dust storm associated with a cold-pool (haboob) issued from deep convection over the Atlas, (7) the first airborne chemical composition measurements of dust in the SHL region with differing composition, sources (determined using Lagrangian backward trajectory calculations) and absorption properties between 2011 and 2012, (8) coincident ozone and dust surface area measurements suggest coarser particles provide a route for ozone depletion, (9) discrepancies between airborne coarse mode size distributions and AERONET sunphotometer retrievals under light dust loadings. These results provide insights into boundary layer and dust processes in the SHL region – a region of substantial global climatic importance.