Historic records of organic compounds from a high Alpine glacier: influences of biomass burning, anthropogenic emissions, and dust transport

Historic records of α-dicarbonyls (glyoxal, methylglyoxal), carboxylic acids (C6–C12 dicarboxylic acids, pinic acid, p-hydroxybenzoic acid, phthalic acid, 4-methylphthalic acid), and ions (oxalate, formate, calcium) were determined with annual resolution in an ice core from Grenzgletscher in the southern Swiss Alps, covering the time period from 1942 to 1993. Chemical analysis of the organic compounds was conducted using ultra-high-performance liquid chromatography (UHPLC) coupled to electrospray ionization high-resolution mass spectrometry (ESI-HRMS) for dicarbonyls and long-chain carboxylic acids and ion chromatography for short-chain carboxylates. Long-term records of the carboxylic acids and dicarbonyls, as well as their source apportionment, are reported for western Europe. This is the first study comprising long-term trends of dicarbonyls and long-chain dicarboxylic acids (C6–C12) in Alpine precipitation. Source assignment of the organic species present in the ice core was performed using principal component analysis. Our results suggest biomass burning, anthropogenic emissions, and transport of mineral dust to be the main parameters influencing the concentration of organic compounds. Ice core records of several highly correlated compounds (e.g., p-hydroxybenzoic acid, pinic acid, pimelic, and suberic acids) can be related to the forest fire history in southern Switzerland. P-hydroxybenzoic acid was found to be the best organic fire tracer in the study area, revealing the highest correlation with the burned area from fires. Historical records of methylglyoxal, phthalic acid, and dicarboxylic acids adipic acid, sebacic acid, and dodecanedioic acid are comparable with that of anthropogenic emissions of volatile organic compounds (VOCs). The small organic acids, oxalic acid and formic acid, are both highly correlated with calcium, suggesting their records to be affected by changing mineral dust transport to the drilling site.

Ice-Core Based Assessment of Historical Anthropogenic Heavy Metal (Cd, Cu, Sb, Zn) Emissions in the Soviet Union

The development of strategies and policies aiming at the reduction of environmental exposure to air pollution requires the assessment of historical emissions. Although anthropogenic emissions from the extended territory of the Soviet Union (SU) considerably influenced concentrations of heavy metals in the Northern Hemisphere, Pb is the only metal with long-term historical emission estimates for this region available, whereas for selected other metals only single values exist. Here we present the first study assessing long-term Cd, Cu, Sb, and Zn emissions in the SU during the period 1935–1991 based on ice-core concentration records from Belukha glacier in the Siberian Altai and emission data from 12 regions in the SU for the year 1980. We show that Zn primarily emitted from the Zn production in Ust-Kamenogorsk (East Kazakhstan) dominated the SU heavy metal emission. Cd, Sb, Zn (Cu) emissions increased between 1935 and the 1970s (1980s) due to expanded non-ferrous metal production. Emissions of the four metals in the beginning of the 1990s were as low as in the 1950s, which we attribute to the economic downturn in industry, changes in technology for an increasing metal recovery from ores, the replacement of coal and oil by gas, and air pollution control.

Allowable carbon emissions lowered by multiple climate targets

Climate targets are designed to inform policies that would limit the magnitude and impacts of climate change caused by anthropogenic emissions of greenhouse gases and other substances. The target that is currently recognized by most world governments1 places a limit of two degrees Celsius on the global mean warming since preindustrial times. This would require large sustained reductions in carbon dioxide emissions during the twenty-first century and beyond2, 3, 4. Such a global temperature target, however, is not sufficient to control many other quantities, such as transient sea level rise5, ocean acidification6, 7 and net primary production on land8, 9. Here, using an Earth system model of intermediate complexity (EMIC) in an observation-informed Bayesian approach, we show that allowable carbon emissions are substantially reduced when multiple climate targets are set. We take into account uncertainties in physical and carbon cycle model parameters, radiative efficiencies10, climate sensitivity11 and carbon cycle feedbacks12, 13 along with a large set of observational constraints. Within this framework, we explore a broad range of economically feasible greenhouse gas scenarios from the integrated assessment community14, 15, 16, 17 to determine the likelihood of meeting a combination of specific global and regional targets under various assumptions. For any given likelihood of meeting a set of such targets, the allowable cumulative emissions are greatly reduced from those inferred from the temperature target alone. Therefore, temperature targets alone are unable to comprehensively limit the risks from anthropogenic emissions.

Identification of adequate neurally adjusted ventilatory assist (NAVA) during systematic increases in the NAVA level

Neurally adjusted ventilatory assist (NAVA) delivers airway pressure (P(aw)) in proportion to the electrical activity of the diaphragm (EAdi) using an adjustable proportionality constant (NAVA level, cm·H(2)O/μV). During systematic increases in the NAVA level, feedback-controlled down-regulation of the EAdi results in a characteristic two-phased response in P(aw) and tidal volume (Vt). The transition from the 1st to the 2nd response phase allows identification of adequate unloading of the respiratory muscles with NAVA (NAVA(AL)). We aimed to develop and validate a mathematical algorithm to identify NAVA(AL). P(aw), Vt, and EAdi were recorded while systematically increasing the NAVA level in 19 adult patients. In a multistep approach, inspiratory P(aw) peaks were first identified by dividing the EAdi into inspiratory portions using Gaussian mixture modeling. Two polynomials were then fitted onto the curves of both P(aw) peaks and Vt. The beginning of the P(aw) and Vt plateaus, and thus NAVA(AL), was identified at the minimum of squared polynomial derivative and polynomial fitting errors. A graphical user interface was developed in the Matlab computing environment. Median NAVA(AL) visually estimated by 18 independent physicians was 2.7 (range 0.4 to 5.8) cm·H(2)O/μV and identified by our model was 2.6 (range 0.6 to 5.0) cm·H(2)O/μV. NAVA(AL) identified by our model was below the range of visually estimated NAVA(AL) in two instances and was above in one instance. We conclude that our model identifies NAVA(AL) in most instances with acceptable accuracy for application in clinical routine and research.

Microgram-Level Radiocarbon Determination of Carbonaceous Particles in Firn and Ice Samples: Pretreatment and Oc/Ec Separation

Carbonaceous particles that comprise organic carbon (OC) and elemental carbon (EC) are of increasing interest in climate research because of their influence on the radiation balance of the Earth. The radiocarbon determination of particulate OC and EC extracted from ice cores provides a powerful tool to reconstruct the long-term natural and anthropogenic emissions of carbonaceous particles. However, this C-14-based source apportionment method has not been applied for the firn section, which is the uppermost part of Alpine glaciers with a typical thickness of up to 50 m. In contrast to glacier ice, firn samples are more easily contaminated through drilling and handling operations. In this study, an alternative decontamination method for firn samples consisting of chiselling off the outer parts instead of rinsing them was developed and verified. The obtained procedural blank of 2.8 +/- 0.8 mu g C for OC is a factor of 2 higher compared to the rinsing method used for ice, but still relatively low compared to the typical OC concentration in firn samples from Alpine glaciers. The EC blank of 0.3 +/- 0.1 mu g C is similar for both methods. For separation of OC and EC for subsequent C-14 analysis, a thermal-optical method instead of the purely thermal method was applied for the first time to firn and ice samples, resulting in a reduced uncertainty of both the mass and C-14 determination. OC and EC concentrations as well as their corresponding fraction of modern for firn and ice samples from Fiescherhorn and Jungfraujoch agree well with published results, validating the new method.

Three decades of global methane sources and sinks

Methane is an important greenhouse gas, responsible for about 20 of the warming induced by long-lived greenhouse gases since pre-industrial times. By reacting with hydroxyl radicals, methane reduces the oxidizing capacity of the atmosphere and generates ozone in the troposphere. Although most sources and sinks of methane have been identified, their relative contributions to atmospheric methane levels are highly uncertain. As such, the factors responsible for the observed stabilization of atmospheric methane levels in the early 2000s, and the renewed rise after 2006, remain unclear. Here, we construct decadal budgets for methane sources and sinks between 1980 and 2010, using a combination of atmospheric measurements and results from chemical transport models, ecosystem models, climate chemistry models and inventories of anthropogenic emissions. The resultant budgets suggest that data-driven approaches and ecosystem models overestimate total natural emissions. We build three contrasting emission scenarios � which differ in fossil fuel and microbial emissions � to explain the decadal variability in atmospheric methane levels detected, here and in previous studies, since 1985. Although uncertainties in emission trends do not allow definitive conclusions to be drawn, we show that the observed stabilization of methane levels between 1999 and 2006 can potentially be explained by decreasing-to-stable fossil fuel emissions, combined with stable-to-increasing microbial emissions. We show that a rise in natural wetland emissions and fossil fuel emissions probably accounts for the renewed increase in global methane levels after 2006, although the relative contribution of these two sources remains uncertain.

Comment on “The phase relation between atmospheric carbon dioxide and global temperature” Humlum et al. [Glob. Planet. Change 100: 51–69.]: Isotopes ignored

A recent study relying purely on statistical analysis of relatively short time series suggested substantial re-thinking of the traditional view about causality explaining the detected rising trend of atmospheric CO2 (atmCO2) concentrations. If these results are well-justified then they should surely compel a fundamental scientific shift in paradigms regarding both atmospheric greenhouse warming mechanism and global carbon cycle. However, the presented work suffers from serious logical deficiencies such as, 1) what could be the sink for fossil fuel CO2 emissions, if neither the atmosphere nor the ocean – as suggested by the authors – plays a role? 2) What is the alternative explanation for ocean acidification if the ocean is a net source of CO2 to the atmosphere? Probably the most provocative point of the commented study is that anthropogenic emissions have little influence on atmCO2 concentrations. The authors have obviously ignored the reconstructed and directly measured carbon isotopic trends of atmCO2 (both δ13C, and radiocarbon dilution) and the declining O2/N2 ratio, although these parameters provide solid evidence that fossil fuel combustion is the major source of atmCO2 increase throughout the Industrial Era.

Placing unprecedented recent fir growth in a European-wide and Holocene-long context

Forest decline played a pivotal role in motivating Europe's political focus on sustainability around 35 years ago. Silver fir (Abies alba) exhibited a particularly severe dieback in the mid-1970s, but disentangling biotic from abiotic drivers remained challenging because both spatial and temporal data were lacking. Here, we analyze 14 136 samples from living trees and historical timbers, together with 356 pollen records, to evaluate recent fir growth from a continent-wide and Holocene-long perspective. Land use and climate change influenced forest growth over the past millennium, whereas anthropogenic emissions of acidic sulfates and nitrates became important after about 1850. Pollution control since the 1980s, together with a warmer but not drier climate, has facilitated an unprecedented surge in productivity across Central European fir stands. Restricted fir distribution prior to the Mesolithic and again in the Modern Era, separated by a peak in abundance during the Bronze Age, is indicative of the long-term interplay of changing temperatures, shifts in the hydrological cycle, and human impacts that have shaped forest structure and productivity.

Fast parallel kriging-based stepwise uncertainty reduction with application to the identification of an excursion set

Stepwise uncertainty reduction (SUR) strategies aim at constructing a sequence of points for evaluating a function  f in such a way that the residual uncertainty about a quantity of interest progressively decreases to zero. Using such strategies in the framework of Gaussian process modeling has been shown to be efficient for estimating the volume of excursion of f above a fixed threshold. However, SUR strategies remain cumbersome to use in practice because of their high computational complexity, and the fact that they deliver a single point at each iteration. In this article we introduce several multipoint sampling criteria, allowing the selection of batches of points at which f can be evaluated in parallel. Such criteria are of particular interest when f is costly to evaluate and several CPUs are simultaneously available. We also manage to drastically reduce the computational cost of these strategies through the use of closed form formulas. We illustrate their performances in various numerical experiments, including a nuclear safety test case. Basic notions about kriging, auxiliary problems, complexity calculations, R code, and data are available online as supplementary materials.

An updated synthesis of the observed and projected impacts of climate change on the chemical, physical and biological processes in the oceans

The 5th Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC) states with very high certainty that anthropogenic emissions have caused measurable changes in the physical ocean environment. These changes are summarized with special focus on those that are predicted to have the strongest, most direct effects on ocean biological processes; namely, ocean warming and associated phenomena (including stratification and sea level rise) as well as deoxygenation and ocean acidification. The biological effects of these changes are then discussed for microbes (including phytoplankton), plants, animals, warm and cold-water corals, and ecosystems. The IPCC AR5 highlighted several areas related to both the physical and biological processes that required further research. As a rapidly developing field, there have been many pertinent studies published since the cut off dates for the AR5, which have increased our understanding of the processes at work. This study undertook an extensive review of recently published literature to update the findings of the AR5 and provide a synthesized review on the main issues facing future oceans. The level of detail provided in the AR5 and subsequent work provided a basis for constructing projections of the state of ocean ecosystems in 2100 under two the Representative Concentration Pathways RCP4.5 and 8.5. Finally the review highlights notable additions, clarifications and points of departure from AR5 provided by subsequent studies.

Water-use efficiency and transpiration across European forests during the Anthropocene

The Earth’s carbon and hydrologic cycles are intimately coupled by gas exchange through plant stomata1, 2, 3. However, uncertainties in the magnitude4, 5, 6 and consequences7, 8 of the physiological responses9, 10 of plants to elevated CO2 in natural environments hinders modelling of terrestrial water cycling and carbon storage11. Here we use annually resolved long-term δ13C tree-ring measurements across a European forest network to reconstruct the physiologically driven response of intercellular CO2 (Ci) caused by atmospheric CO2 (Ca) trends. When removing meteorological signals from the δ13C measurements, we find that trees across Europe regulated gas exchange so that for one ppmv atmospheric CO2 increase, Ci increased by ~0.76 ppmv, most consistent with moderate control towards a constant Ci/Ca ratio. This response corresponds to twentieth-century intrinsic water-use efficiency (iWUE) increases of 14 ± 10 and 22 ± 6% at broadleaf and coniferous sites, respectively. An ensemble of process-based global vegetation models shows similar CO2 effects on iWUE trends. Yet, when operating these models with climate drivers reintroduced, despite decreased stomatal opening, 5% increases in European forest transpiration are calculated over the twentieth century. This counterintuitive result arises from lengthened growing seasons, enhanced evaporative demand in a warming climate, and increased leaf area, which together oppose effects of CO2-induced stomatal closure. Our study questions changes to the hydrological cycle, such as reductions in transpiration and air humidity, hypothesized to result from plant responses to anthropogenic emissions.

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.

Continuous CO2/CH4/CO measurements (2012–2014) at Beromünster tall tower station in Switzerland

The understanding of the continental carbon budget is essential to predict future climate change. In order to quantify CO₂ and CH₄ fluxes at the regional scale, a measurement system was installed at the former radio tower in Beromünster as part of the Swiss greenhouse gas monitoring network (CarboCount CH). We have been measuring the mixing ratios of CO₂, CH₄ and CO on this tower with sample inlets at 12.5, 44.6, 71.5, 131.6 and 212.5 m above ground level using a cavity ring down spectroscopy (CRDS) analyzer. The first 2-year (December 2012–December 2014) continuous atmospheric record was analyzed for seasonal and diurnal variations and interspecies correlations. In addition, storage fluxes were calculated from the hourly profiles along the tower. The atmospheric growth rates from 2013 to 2014 determined from this 2-year data set were 1.78 ppm yr⁻¹, 9.66 ppb yr⁻¹ and and -1.27 ppb yr⁻¹ for CO₂, CH₄ and CO, respectively. After detrending, clear seasonal cycles were detected for CO₂ and CO, whereas CH₄ showed a stable baseline suggesting a net balance between sources and sinks over the course of the year. CO and CO₂ were strongly correlated (r² > 0.75) in winter (DJF), but almost uncorrelated in summer. In winter, anthropogenic emissions dominate the biospheric CO₂ fluxes and the variations in mixing ratios are large due to reduced vertical mixing. The diurnal variations of all species showed distinct cycles in spring and summer, with the lowest sampling level showing the most pronounced diurnal amplitudes. The storage flux estimates exhibited reasonable diurnal shapes for CO₂, but underestimated the strength of the surface sinks during daytime. This seems plausible, keeping in mind that we were only able to calculate the storage fluxes along the profile of the tower but not the flux into or out of this profile, since no Eddy covariance flux measurements were taken at the top of the tower.

Modeling the formation and aging of secondary organic aerosols during CalNex 2010

Four different literature parameterizations for the formation and evolution of urban secondary organic aerosol (SOA) frequently used in 3-D models are evaluated using a 0-D box model representing the Los Angeles metropolitan region during the California Research at the Nexus of Air Quality and Climate Change (CalNex) 2010 campaign. We constrain the model predictions with measurements from several platforms and compare predictions with particle- and gas-phase observations from the CalNex Pasadena ground site. That site provides a unique opportunity to study aerosol formation close to anthropogenic emission sources with limited recirculation. The model SOA that formed only from the oxidation of VOCs (V-SOA) is insufficient to explain the observed SOA concentrations, even when using SOA parameterizations with multi-generation oxidation that produce much higher yields than have been observed in chamber experiments, or when increasing yields to their upper limit estimates accounting for recently reported losses of vapors to chamber walls. The Community Multiscale Air Quality (WRF-CMAQ) model (version 5.0.1) provides excellent predictions of secondary inorganic particle species but underestimates the observed SOA mass by a factor of 25 when an older VOC-only parameterization is used, which is consistent with many previous model–measurement comparisons for pre-2007 anthropogenic SOA modules in urban areas. Including SOA from primary semi-volatile and intermediate-volatility organic compounds (P-S/IVOCs) following the parameterizations of Robinson et al. (2007), Grieshop et al. (2009), or Pye and Seinfeld (2010) improves model–measurement agreement for mass concentration. The results from the three parameterizations show large differences (e.g., a factor of 3 in SOA mass) and are not well constrained, underscoring the current uncertainties in this area. Our results strongly suggest that other precursors besides VOCs, such as P-S/IVOCs, are needed to explain the observed SOA concentrations in Pasadena. All the recent parameterizations overpredict urban SOA formation at long photochemical ages (3 days) compared to observations from multiple sites, which can lead to problems in regional and especially global modeling. However, reducing IVOC emissions by one-half in the model to better match recent IVOC measurements improves SOA predictions at these long photochemical ages. Among the explicitly modeled VOCs, the precursor compounds that contribute the greatest SOA mass are methylbenzenes. Measured polycyclic aromatic hydrocarbons (naphthalenes) contribute 0.7% of the modeled SOA mass. The amounts of SOA mass from diesel vehicles, gasoline vehicles, and cooking emissions are estimated to be 16–27, 35–61, and 19–35 %, respectively, depending on the parameterization used, which is consistent with the observed fossil fraction of urban SOA, 71(+-3) %. The relative contribution of each source is uncertain by almost a factor of 2 depending on the parameterization used. In-basin biogenic VOCs are predicted to contribute only a few percent to SOA. A regional SOA background of approximately 2.1 μgm-3 is also present due to the long-distance transport of highly aged OA, likely with a substantial contribution from regional biogenic SOA. The percentage of SOA from diesel vehicle emissions is the same, within the estimated uncertainty, as reported in previous work that analyzed the weekly cycles in OA concentrations (Bahreini et al., 2012; Hayes et al., 2013). However, the modeling work presented here suggests a strong anthropogenic source of modern carbon in SOA, due to cooking emissions, which was not accounted for in those previous studies and which is higher on weekends. Lastly, this work adapts a simple two-parameter model to predict SOA concentration and O/C from urban emissions. This model successfully predicts SOA concentration, and the optimal parameter combination is very similar to that found for Mexico City. This approach provides a computationally inexpensive method for predicting urban SOA in global and climate models. We estimate pollution SOA to account for 26 Tg yr-1 of SOA globally, or 17% of global SOA, one third of which is likely to be non-fossil.