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.

Pre-Columbian landscape impact and agriculture in the Monumental Mound region of the Llanos de Moxos, lowland Bolivia

We present a multiproxy study of land use by a pre-Columbian earth mounds culture in the Bolivian Amazon. The Monumental Mounds Region (MMR) is an archaeological sub-region characterized by hundreds of pre-Columbian habitation mounds associated with a complex network of canals and causeways, and situated in the forest–savanna mosaic of the Llanos de Moxos. Pollen, phytolith, and charcoal analyses were performed on a sediment core from a large lake (14 km2), Laguna San José (14°56.97′S, 64°29.70′W).We found evidence of high levels of anthropogenic burning from AD 400 to AD 1280, corroborating dated occupation layers in two nearby excavated habitation mounds. The charcoal decline pre-dates the arrival of Europeans by at least 100 yr, and challenges the notion that the mounds culture declined because of European colonization. We show that the surrounding savanna soils were sufficiently fertile to support crops, and the presence of maize throughout the record shows that the area was continuously cultivated despite land-use change at the end of the earthmounds culture. We suggest that burning was largely confined to the savannas, rather than forests, and that pre-Columbian deforestation was localized to the vicinity of individual habitation mounds, whereas the inter-mound areas remained largely forested.

Estimation of sea surface temperature from the Spinning Enhanced Visible and Infra Red Imager, improved using numerical weather prediction

Most of the operational Sea Surface Temperature (SST) products derived from satellite infrared radiometry use multi-spectral algorithms. They show, in general, reasonable performances with root mean square (RMS) residuals around 0.5 K when validated against buoy measurements, but have limitations, particularly a component of the retrieval error that relates to such algorithms' limited ability to cope with the full variability of atmospheric absorption and emission. We propose to use forecast atmospheric profiles and a radiative transfer model to simulate the algorithmic errors of multi-spectral algorithms. In the practical case of SST derived from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard Meteosat Second Generation (MSG), we demonstrate that simulated algorithmic errors do explain a significant component of the actual errors observed for the non linear (NL) split window algorithm in operational use at the Centre de Météorologie Spatiale (CMS). The simulated errors, used as correction terms, reduce significantly the regional biases of the NL algorithm as well as the standard deviation of the differences with drifting buoy measurements. The availability of atmospheric profiles associated with observed satellite-buoy differences allows us to analyze the origins of the main algorithmic errors observed in the SEVIRI field of view: a negative bias in the inter-tropical zone, and a mid-latitude positive bias. We demonstrate how these errors are explained by the sensitivity of observed brightness temperatures to the vertical distribution of water vapour, propagated through the SST retrieval algorithm.

The reversibility of sea level rise

During the last century, global climate has been warming, and projections indicate that such a warming is likely to continue over coming decades. Most of the extra heat is stored in the ocean, resulting in thermal expansion of seawater and global mean sea level rise. Previous studies have shown that after CO2 emissions cease or CO2 concentration is stabilized, global mean surface air temperature stabilizes or decreases slowly, but sea level continues to rise. Using idealized CO2 scenario simulations with a hierarchy of models including an AOGCM and a step-response model, the authors show how the evolution of thermal expansion can be interpreted in terms of the climate energy balance and the vertical profile of ocean warming. Whereas surface temperature depends on cumulative CO2 emissions, sea level rise due to thermal expansion depends on the time profile of emissions. Sea level rise is smaller for later emissions, implying that targets to limit sea level rise would need to refer to the rate of emissions, not only to the time integral. Thermal expansion is in principle reversible, but to halt or reverse it quickly requires the radiative forcing to be reduced substantially, which is possible on centennial time scales only by geoengineering. If it could be done, the results indicate that heat would leave the ocean more readily than it entered, but even if thermal expansion were returned to zero, the geographical pattern of sea level would be altered. Therefore, despite any aggressive CO2 mitigation, regional sea level change is inevitable.

Vertical distribution and radiative effects of mineral dust and biomass burning aerosol over West Africa during DABEX

This paper presents measurements of the vertical distribution of aerosol extinction coefficient over West Africa during the Dust and Biomass-burning Aerosol Experiment (DABEX)/African Monsoon Multidisciplinary Analysis dry season Special Observing Period Zero (AMMA-SOP0). In situ aircraft measurements from the UK FAAM aircraft have been compared with two ground-based lidars (POLIS and ARM MPL) and an airborne lidar on an ultralight aircraft. In general, mineral dust was observed at low altitudes (up to 2 km), and a mixture of biomass burning aerosol and dust was observed at altitudes of 2–5 km. The study exposes difficulties associated with spatial and temporal variability when intercomparing aircraft and ground measurements. Averaging over many profiles provided a better means of assessing consistent errors and biases associated with in situ sampling instruments and retrievals of lidar ratios. Shortwave radiative transfer calculations and a 3-year simulation with the HadGEM2-A climate model show that the radiative effect of biomass burning aerosol was somewhat sensitive to the vertical distribution of aerosol. In particular, when the observed low-level dust layer was included in the model, the absorption of solar radiation by the biomass burning aerosols increased by 10%. We conclude that this absorption enhancement was caused by the dust reflecting solar radiation up into the biomass burning aerosol layer. This result illustrates that the radiative forcing of anthropogenic absorbing aerosol can be sensitive to the presence of natural aerosol species.

Black carbon vertical profiles strongly affect its radiative forcing uncertainty

The impact of black carbon (BC) aerosols on the global radiation balance is not well constrained. Here twelve global aerosol models are used to show that at least 20% of the present uncertainty in modeled BC direct radiative forcing (RF) is due to diversity in the simulated vertical profile of BC mass. Results are from phases 1 and 2 of the global aerosol model intercomparison project (AeroCom). Additionally, a significant fraction of the variability is shown to come from high altitudes, as, globally, more than 40% of the total BC RF is exerted above 5 km. BC emission regions and areas with transported BC are found to have differing characteristics. These insights into the importance of the vertical profile of BC lead us to suggest that observational studies are needed to better characterize the global distribution of BC, including in the upper troposphere.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

The influence of Antarctica on the momentum budget of the southern extratropics

The antarctic plateau acts as a strong heat sink for the global climate, cooling the atmosphere and radiating energy to space. A cold dense atmospheric boundary layer is formed. Strong surface winds are formed as the boundary layer drains off the plateau. These drainage winds and the eddy fluxes necessary to maintain them are analysed in a general circulation model (GCM). The drainage flow is well represented in the GCM. The associated mean meridional circulation is analysed in isentropic coordinates. The momentum budget over Antarctica reveals a balance between the Eliassen-Palm flux convergence and the Coriolis torque exerted by the mean meridional mass flux. Both vertical and horizontal components of the Eliassen-Palm flux contribute, the vertical component being the greater.

Modeling the summertime evolution of sea-ice melt ponds

1] We present a mathematical model describing the summer melting of sea ice. We simulate the evolution of melt ponds and determine area coverage and total surface ablation. The model predictions are tested for sensitivity to the melt rate of unponded ice, enhanced melt rate beneath the melt ponds, vertical seepage, and horizontal permeability. The model is initialized with surface topographies derived from laser altimetry corresponding to first-year sea ice and multiyear sea ice. We predict that there are large differences in the depth of melt ponds and the area of coverage between the two types of ice. We also find that the vertical seepage rate and the melt rate of unponded ice are important in determining the total surface ablation and area covered by melt ponds.

Reduced-order dynamics of the Martian atmospheric dynamics

In this paper we explore the possibility of deriving low-dimensional models of the dynamics of the Martian atmosphere. The analysis consists of a Proper Orthogonal Decomposition (POD) of the atmospheric streamfunction after first decomposing the vertical structure with a set of eigenmodes. The vertical modes were obtained from the quasi-geostrophic vertical structure equation. The empirical orthogonal functions (EOFs) were optimized to represent the atmospheric total energy. The total energy was used as the criterion to retain those modes with large energy content and discard the rest. The principal components (PCs) were analysed by means of Fourier analysis, so that the dominant frequencies could be identified. It was possible to observe the strong influence of the diurnal cycle and to identify the motion and vacillation of baroclinic waves.

Bollène-2002 experiment: radar quantitative precipitation estimation in the Cévennes–Vivarais region, France

The Bollène-2002 Experiment was aimed at developing the use of a radar volume-scanning strategy for conducting radar rainfall estimations in the mountainous regions of France. A developmental radar processing system, called Traitements Régionalisés et Adaptatifs de Données Radar pour l’Hydrologie (Regionalized and Adaptive Radar Data Processing for Hydrological Applications), has been built and several algorithms were specifically produced as part of this project. These algorithms include 1) a clutter identification technique based on the pulse-to-pulse variability of reflectivity Z for noncoherent radar, 2) a coupled procedure for determining a rain partition between convective and widespread rainfall R and the associated normalized vertical profiles of reflectivity, and 3) a method for calculating reflectivity at ground level from reflectivities measured aloft. Several radar processing strategies, including nonadaptive, time-adaptive, and space–time-adaptive variants, have been implemented to assess the performance of these new algorithms. Reference rainfall data were derived from a careful analysis of rain gauge datasets furnished by the Cévennes–Vivarais Mediterranean Hydrometeorological Observatory. The assessment criteria for five intense and long-lasting Mediterranean rain events have proven that good quantitative precipitation estimates can be obtained from radar data alone within 100-km range by using well-sited, well-maintained radar systems and sophisticated, physically based data-processing systems. The basic requirements entail performing accurate electronic calibration and stability verification, determining the radar detection domain, achieving efficient clutter elimination, and capturing the vertical structure(s) of reflectivity for the target event. Radar performance was shown to depend on type of rainfall, with better results obtained with deep convective rain systems (Nash coefficients of roughly 0.90 for point radar–rain gauge comparisons at the event time step), as opposed to shallow convective and frontal rain systems (Nash coefficients in the 0.6–0.8 range). In comparison with time-adaptive strategies, the space–time-adaptive strategy yields a very significant reduction in the radar–rain gauge bias while the level of scatter remains basically unchanged. Because the Z–R relationships have not been optimized in this study, results are attributed to an improved processing of spatial variations in the vertical profile of reflectivity. The two main recommendations for future work consist of adapting the rain separation method for radar network operations and documenting Z–R relationships conditional on rainfall type.

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.

Simulation of the global ENSO–Tropical cyclone teleconnection by a high-resolution coupled general circulation model

This study assesses the influence of the El Niño–Southern Oscillation (ENSO) on global tropical cyclone activity using a 150-yr-long integration with a high-resolution coupled atmosphere–ocean general circulation model [High-Resolution Global Environmental Model (HiGEM); with N144 resolution: ~90 km in the atmosphere and ~40 km in the ocean]. Tropical cyclone activity is compared to an atmosphere-only simulation using the atmospheric component of HiGEM (HiGAM). Observations of tropical cyclones in the International Best Track Archive for Climate Stewardship (IBTrACS) and tropical cyclones identified in the Interim ECMWF Re-Analysis (ERA-Interim) are used to validate the models. Composite anomalies of tropical cyclone activity in El Niño and La Niña years are used. HiGEM is able to capture the shift in tropical cyclone locations to ENSO in the Pacific and Indian Oceans. However, HiGEM does not capture the expected ENSO–tropical cyclone teleconnection in the North Atlantic. HiGAM shows more skill in simulating the global ENSO–tropical cyclone teleconnection; however, variability in the Pacific is overpronounced. HiGAM is able to capture the ENSO–tropical cyclone teleconnection in the North Atlantic more accurately than HiGEM. An investigation into the large-scale environmental conditions, known to influence tropical cyclone activity, is used to further understand the response of tropical cyclone activity to ENSO in the North Atlantic and western North Pacific. The vertical wind shear response over the Caribbean is not captured in HiGEM compared to HiGAM and ERA-Interim. Biases in the mean ascent at 500 hPa in HiGEM remain in HiGAM over the western North Pacific; however, a more realistic low-level vorticity in HiGAM results in a more accurate ENSO–tropical cyclone teleconnection.

Impact of atmospheric transport on the evolution of microphysical and optical properties of Saharan dust

Saharan dust affects the climate by altering the radiation balance and by depositing minerals to the Atlantic Ocean. Both are dependent on particle size. We present aircraft measurements comprising 42 profiles of size distribution (0.1–300 µm), representing freshly uplifted dust, regional aged dust, and dust in the Saharan Air Layer (SAL) over the Canary Islands. The mean effective diameter of dust in SAL profiles is 4.5 µm smaller than that in freshly uplifted dust, while the vertical structure changes from a low shallow layer (0–1.5 km) to a well-mixed deep Saharan dust layer (0–5 km). Size distributions show a loss of 60 to 90% of particles larger than 30 µm 12 h after uplift. The single scattering albedo (SSA) increases from 0.92 to 0.94 to 0.95 between fresh, aged, and SAL profiles: this is enough to alter heating rates by 26%. Some fresh dust close to the surface shows SSA as low as 0.85