Natural and anthropogenic methane fluxes in Eurasia: a mesoscale quantification by generalized atmospheric inversion

Eight surface observation sites providing quasi-continuous measurements of atmospheric methane mixingratios have been operated since the mid-2000’s in Siberia. For the first time in a single work, we assimilate 1 year of these in situ observations in an atmospheric inversion. Our objective is to quantify methane surface fluxes from anthropogenic and wetland sources at the mesoscale in the Siberian lowlands for the year 2010. To do so, we first inquire about the way the inversion uses the observations and the way the fluxes are constrained by the observation sites. As atmospheric inver- sions at the mesoscale suffer from mis-quantified sources of uncertainties, we follow recent innovations in inversion techniques and use a new inversion approach which quantifies the uncertainties more objectively than the previous inversion systems. We find that, due to errors in the representation of the atmospheric transport and redundant pieces of information, only one observation every few days is found valuable by the inversion. The remaining high-resolution quasi-continuous signal is representative of very local emission patterns difficult to analyse with a mesoscale system. An analysis of the use of information by the inversion also reveals that the observation sites constrain methane emissions within a radius of 500 km. More observation sites than the ones currently in operation are then necessary to constrain the whole Siberian lowlands. Still, the fluxes within the constrained areas are quantified with objectified uncertainties. Finally, the tolerance intervals for posterior methane fluxes are of roughly 20 % (resp. 50 %) of the fluxes for anthropogenic (resp. wetland) sources. About 50–70 % of Siberian lowlands emissions are constrained by the inversion on average on an annual basis. Extrapolating the figures on the constrained areas to the whole Siberian lowlands, we find a regional methane budget of 5–28 TgCH4 for the year 2010, i.e. 1–5 % of the global methane emissions. As very few in situ observations are available in the region of interest, observations of methane total columns from the Greenhouse Gas Observing SATellite (GOSAT) are tentatively used for the evaluation of the inversion results, but they exhibit only a marginal signal from the fluxes within the region of interest.

Evaluating atmospheric methane inversion model results for Pallas, northern Finland

A state-of-the-art inverse model, CarbonTracker Data Assimilation Shell (CTDAS), was used to optimize estimates of methane (CH4) surface fluxes using atmospheric observations of CH4 as a constraint. The model consists of the latest version of the TM5 atmospheric chemistry-transport model and an ensemble Kalman filter based data assimilation system. The model was constrained by atmospheric methane surface concentrations, obtained from the World Data Centre for Greenhouse Gases (WDCGG). Prior methane emissions were specified for five sources: biosphere, anthropogenic, fire, termites and ocean, of which bio-sphere and anthropogenic emissions were optimized. Atmospheric CH 4 mole fractions for 2007 from northern Finland calculated from prior and optimized emissions were compared with observations. It was found that the root mean squared errors of the posterior esti - mates were more than halved. Furthermore, inclusion of NOAA observations of CH 4 from weekly discrete air samples collected at Pallas improved agreement between posterior CH 4 mole fraction estimates and continuous observations, and resulted in reducing optimized biosphere emissions and their uncertainties in northern Finland.

Accumulation rates and predominant atmospheric sources of natural and anthropogenic Hg and Pb on the Faroe Islands

A monolith representing 5420 14C yr of peat accumulation was collected from a blanket bog at Myrarnar, Faroe Islands. The maximum Hg concentration (498 ng/g at a depth of 4.5 cm) coincides with the maximum concentration of anthropogenic Pb (111 μg/g). Age dating of recent peat accumulation using 210Pb (CRS model) shows that the maxima in Hg and Pb concentrations occur at AD 1954 ± 2. These results, combined with the isotopic composition of Pb in that sample (206Pb/207Pb = 1.1720 ± 0.0017), suggest that coal burning was the dominant source of both elements. From the onset of peat accumulation (ca. 4286 BC) until AD 1385, the ratios Hg/Br and Hg/Se were constant (2.2 ± 0.5 × 10-4 and 8.5 ± 1.8 × 10-3, respectively). Since then, Hg/Br and Hg/Se values have increased, also reaching their maxima in AD 1954. The age date of the maximum concentrations of anthropogenic Hg and Pb in the Faroe Islands is consistent with a previous study of peat cores from Greenland and Denmark (dated using the atmospheric bomb pulse curve of 14C), which showed maximum concentrations in AD 1953. The average rate of atmospheric Hg accumulation from 1520 BC to AD 1385 was 1.27 ± 0.38 μg/m2/yr. The Br and Se concentrations and the background Hg/Br and Hg/Se ratios were used to calculate the average rate of natural Hg accumulation for the same period, 1.32 ± 0.36 μg/m2/yr and 1.34 ± 0.29 μg/m2/yr, respectively. These fluxes are similar to the preanthropogenic rates obtained using peat cores from Switzerland, southern Greenland, southern Ontario, Canada, and the northeastern United States. Episodic volcanic emissions and the continual supply of marine aerosols to the Faroe Islands, therefore, have not contributed significantly to the Hg inventory or the Hg accumulation rates, relative to these other areas. The maximum rate of Hg accumulation was 34 μg/m2/yr. The greatest fluxes of anthropogenic Hg accumulation calculated using Br and Se, respectively, were 26 and 31 μg/m2/yr. The rate of atmospheric Hg accumulation in 1998 (16 μg/m2/yr) is comparable to the values recently obtained by atmospheric transport modeling for Denmark, the Faroe Islands, and Greenland.

Impact of oceanic circulation changes on atmospheric δ¹³ CO₂

δ¹³ CO₂ measured in Antarctic ice cores provides constraints on oceanic and terrestrial carbon cycle processes linked with millennial-scale changes in atmospheric CO₂. However, the interpretation of δ¹³ CO₂ is not straight-forward. Using carbon isotope-enabled versions of the LOVECLIM and Bern3D models, we perform a set of sensitivity experiments in which the formation rates of North Atlantic Deep Water (NADW), North Pacific Deep Water (NPDW), Antarctic Bottom Water (AABW), and Antarctic Intermediate Water (AAIW) are varied. We study the impact of these circulation changes on atmospheric δ¹³ CO₂ as well as on the oceanic δ¹³ CO₂ distribution. In general, we find that the formation rates of AABW, NADW, NPDW, and AAIW are negatively correlated with changes in δ¹³ CO₂: namely, strong oceanic ventilation decreases atmospheric δ¹³ CO₂. However, since large-scale oceanic circulation reorganizations also impact nutrient utilization and the Earth’s climate, the relationship between atmospheric δ¹³ CO₂ levels and ocean ventilation rate is not unequivocal. In both models atmospheric δ¹³ CO₂ is very sensitive to changes in AABW formation rates: increased AABW formation enhances the transport of low δ¹³ CO₂ waters to the surface and decreases atmospheric δ¹³ CO₂. By contrast, the impact of NADW changes on atmospheric δ¹³ CO₂ is less robust and might be model dependent. This results from complex interplay between global climate, carbon cycle, and the formation rate of NADW, a water body characterized by relatively high δ¹³ CO₂.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

Zeeman effect in atmospheric O₂ measured by ground-based microwave radiometry

In this work we study the Zeeman effect on stratospheric O₂ using ground-based microwave radiometer measurements. The interaction of the Earth magnetic field with the oxygen dipole leads to a splitting of O₂ energy states, which polarizes the emission spectra. A special campaign was carried out in order to measure this effect in the oxygen emission line centered at 53.07 GHz. Both a fixed and a rotating mirror were incorporated into the TEMPERA (TEMPERature RAdiometer) in order to be able to measure under different observational angles. This new configuration allowed us to change the angle between the observational path and the Earth magnetic field direction. Moreover, a high-resolution spectrometer (1 kHz) was used in order to measure for the first time the polarization state of the radiation due to the Zeeman effect in the main isotopologue of oxygen from ground-based microwave measurements. The measured spectra showed a clear polarized signature when the observational angles were changed, evidencing the Zeeman effect in the oxygen molecule. In addition, simulations carried out with the Atmospheric Radiative Transfer Simulator (ARTS) allowed us to verify the microwave measurements showing a very good agreement between model and measurements. The results suggest some interesting new aspects for research of the upper atmosphere.

Hygroscopic growth of atmospheric aerosol particles based on active remote sensing and radiosounding measurements: selected cases in southeastern Spain

A new methodology based on combining active and passive remote sensing and simultaneous and collocated radiosounding data to study the aerosol hygroscopic growth effects on the particle optical and microphysical properties is presented. The identification of hygroscopic growth situations combines the analysis of multispectral aerosol particle backscatter coefficient and particle linear depolarization ratio with thermodynamic profiling of the atmospheric column. We analyzed the hygroscopic growth effects on aerosol properties, namely the aerosol particle backscatter coefficient and the volume concentration profiles, using data gathered at Granada EARLINET station. Two study cases, corresponding to different aerosol loads and different aerosol types, are used for illustrating the potential of this methodology. Values of the aerosol particle backscatter coefficient enhancement factors range from 2.1 ± 0.8 to 3.9 ± 1.5, in the ranges of relative humidity 60–90 and 40–83%, being similar to those previously reported in the literature. Differences in the enhancement factor are directly linked to the composition of the atmospheric aerosol. The largest value of the aerosol particle backscatter coefficient enhancement factor corresponds to the presence of sulphate and marine particles that are more affected by hygroscopic growth. On the contrary, the lowest value of the enhancement factor corresponds to an aerosol mixture containing sulphates and slight traces of mineral dust. The Hänel parameterization is applied to these case studies, obtaining results within the range of values reported in previous studies, with values of the γ exponent of 0.56 ± 0.01 (for anthropogenic particles slightly influenced by mineral dust) and 1.07 ± 0.01 (for the situation dominated by anthropogenic particles), showing the convenience of this remote sensing approach for the study of hygroscopic effects of the atmospheric aerosol under ambient unperturbed conditions. For the first time, the retrieval of the volume concentration profiles for these cases using the Lidar Radiometer Inversion Code (LIRIC) allows us to analyze the aerosol hygroscopic growth effects on aerosol volume concentration, observing a stronger increase of the fine mode volume concentration with increasing relative humidity.

Comparison of co-located independent ground-based middle atmospheric wind and temperature measurements with numerical weather prediction models

High-resolution, ground-based and independent observations including co-located wind radiometer, lidar stations, and infrasound instruments are used to evaluate the accuracy of general circulation models and data-constrained assimilation systems in the middle atmosphere at northern hemisphere midlatitudes. Systematic comparisons between observations, the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analyses including the recent Integrated Forecast System cycles 38r1 and 38r2, the NASA’s Modern-Era Retrospective Analysis for Research and Applications (MERRA) reanalyses, and the free-running climate Max Planck Institute–Earth System Model–Low Resolution (MPI-ESM-LR) are carried out in both temporal and spectral dom ains. We find that ECMWF and MERRA are broadly consistent with lidar and wind radiometer measurements up to ~40 km. For both temperature and horizontal wind components, deviations increase with altitude as the assimilated observations become sparser. Between 40 and 60 km altitude, the standard deviation of the mean difference exceeds 5 K for the temperature and 20 m/s for the zonal wind. The largest deviations are observed in winter when the variability from large-scale planetary waves dominates. Between lidar data and MPI-ESM-LR, there is an overall agreement in spectral amplitude down to 15–20 days. At shorter time scales, the variability is lacking in the model by ~10 dB. Infrasound observations indicate a general good agreement with ECWMF wind and temperature products. As such, this study demonstrates the potential of the infrastructure of the Atmospheric Dynamics Research Infrastructure in Europe project that integrates various measurements and provides a quantitative understanding of stratosphere-troposphere dynamical coupling for numerical weather prediction applications.

High-resolution AMS (super 14) C dating of post-bomb peat archives of atmospheric pollutants.

Peat deposits in Greenland and Denmark were investigated to show that high-resolution dating of these archives of atmospheric deposition can be provided for the last 50 years by radiocarbon dating using the atmospheric bomb pulse. (super 14) C was determined in macrofossils from sequential one cm slices using accelerator mass spectrometry (AMS). Values were calibrated with a general-purpose curve derived from annually averaged atmospheric (super 14) CO (sub 2) values in the northernmost northern hemisphere (NNH, 30 degrees -90 degrees N). We present a through review of (super 14) C bomb-pulse data from the NNH including our own measurements made in tree rings and seeds from Arizona as well as other previously published data. We show that our general-purpose calibration curve is valid for the whole NNH producing accurate dates within 1-2 years. In consequence, (super 14) C AMS can precisely date individual points in recent peat deposits within the range of the bomb-pulse (from the mid-1950s on). Comparing the (super 14) C AMS results with the customary dating method for recent peat profiles by (super 210) Pb, we show that the use of (super 137) Cs to validate and correct (super 210) Pb dates proves to be more problematic than previously supposed. As a unique example of our technique, we show how this chronometer can be applied to identify temporal changes in Hg concentrations from Danish and Greenland peat cores.

Covariation of deep Southern Ocean oxygenation and atmospheric CO2 through the last ice age

No single mechanism can account for the full amplitude of past atmospheric carbon dioxide (CO2) concentration variability over glacial–interglacial cycles. A build-up of carbon in the deep ocean has been shown to have occurred during the Last Glacial Maximum. However, the mechanisms responsible for the release of the deeply sequestered carbon to the atmosphere at deglaciation, and the relative importance of deep ocean sequestration in regulating millennial-timescale variations in atmospheric CO2 concentration before the Last Glacial Maximum, have remained unclear. Here we present sedimentary redox-sensitive trace metal records from the Antarctic Zone of the Southern Ocean that provide a reconstruction of transient changes in deep ocean oxygenation and, by inference, respired carbon storage throughout the last glacial cycle. Our data suggest that respired carbon was removed from the abyssal Southern Ocean during the Northern Hemisphere cold phases of the deglaciation, when atmospheric CO2 concentration increased rapidly, reflecting—at least in part— a combination of dwindling iron fertilization by dust and enhanced deep ocean ventilation. Furthermore, our records show that the observed covariation between atmospheric CO2 concentration and abyssal Southern Ocean oxygenation was maintained throughout most of the past 80,000 years. This suggests that on millennial timescales deep ocean circulation and iron fertilization in the Southern Ocean played a consistent role in modifying atmospheric CO2 concentration.

Biological and physical controls in the Southern Ocean on past millennial-scale atmospheric CO2 changes

Millennial-scale climate changes during the last glacial period and deglaciation were accompanied by rapid changes in atmospheric CO2 that remain unexplained. While the role of the Southern Ocean as a ’control valve’ on ocean–atmosphere CO2 exchange has been emphasized, the exact nature of this role, in particular the relative contributions of physical (for example, ocean dynamics and air–sea gas exchange) versus biological processes (for example, export productivity), remains poorly constrained. Here we combine reconstructions of bottom-water [O2], export production and 14C ventilation ages in the sub-Antarctic Atlantic, and show that atmospheric CO2 pulses during the last glacial- and deglacial periods were consistently accompanied by decreases in the biological export of carbon and increases in deep-ocean ventilation via southern-sourced water masses. These findings demonstrate how the Southern Ocean’s ’organic carbon pump’ has exerted a tight control on atmospheric CO2, and thus global climate, specifically via a synergy of both physical and biological processes.