Refine
Year of publication
- 2016 (61) (remove)
Document Type
- Article (61) (remove)
Language
- English (61)
Has Fulltext
- yes (61)
Is part of the Bibliography
- no (61)
Keywords
- Atmospheric chemistry (3)
- Andes (1)
- BSRN (1)
- Bayesian network (1)
- Biogeochemistry (1)
- COSMO-CLM (1)
- Clausius–Clapeyron scaling (1)
- Climate-change impacts (1)
- Cloud properties (1)
- Convection-permitting climate model (1)
Institute
- Geowissenschaften (61) (remove)
This study reports and discusses new radiometric ages, petrographical and volcanological observations and whole rock geochemical data of the rocks of the Rudnik Mts. volcano-intrusive complex. The complex hosts a Pb-Zn-Ag deposit and belongs to the Serbo-Macedonian metallogenetic belt. Two distinct igneous events are distignuished. The first occurred >30 Ma and was mainly characterized by extrusive and shallow intrusive dacites and andesites and was unrelated to mineralization. The second igneous event occurred <23 Ma and was highly heterogeneous in terms of volcanic products and petrographic varieties, but with predominance of quartzlatites. The dacite-andesites (first event) and the quartzlatites (second event) are geochemically similar and display a calc-alkaline affinity and highly incompatible element enriched patterns on spider diagrams, but the younger quartzlatites are richer in K2O, Rb and Ba and poorer in Sr. This is taken as evidence that mixing between an ultrapotassic lamprophyre/lamproite magma and an acid calc-alkaline (dacite-like) magma was essential petrogenetic processes during the second event. The proposed simplified volcanological model suggests that this mixing was responsible for triggering strongly explosive volcanic activity as well as for providing conditions for active hydrothermal and mineralization processes. The observed link between a specific magmatic phase and ore deposit formation can be a general phennomenon in the Balkans, and must be addressed by further and more advanced studies.
Atmospheric aerosols and their effect on clouds are thought to be important for anthropogenic radiative forcing of the climate, yet remain poorly understood1. Globally, around half of cloud condensation nuclei originate from nucleation of atmospheric vapours2. It is thought that sulfuric acid is essential to initiate most particle formation in the atmosphere3,4, and that ions have a relatively minor role5. Some laboratory studies, however, have reported organic particle formation without the intentional addition of sulfuric acid, although contamination could not be excluded6,7. Here we present evidence for the formation of aerosol particles from highly oxidized biogenic vapours in the absence of sulfuric acid in a large chamber under atmospheric conditions. The highly oxygenated molecules (HOMs) are produced by ozonolysis of α-pinene. We find that ions from Galactic cosmic rays increase the nucleation rate by one to two orders of magnitude compared with neutral nucleation. Our experimental findings are supported by quantum chemical calculations of the cluster binding energies of representative HOMs. Ion-induced nucleation of pure organic particles constitutes a potentially widespread source of aerosol particles in terrestrial environments with low sulfuric acid pollution.
About half of present-day cloud condensation nuclei originate from atmospheric nucleation, frequently appearing as a burst of new particles near midday1. Atmospheric observations show that the growth rate of new particles often accelerates when the diameter of the particles is between one and ten nanometres2,3. In this critical size range, new particles are most likely to be lost by coagulation with pre-existing particles4, thereby failing to form new cloud condensation nuclei that are typically 50 to 100 nanometres across. Sulfuric acid vapour is often involved in nucleation but is too scarce to explain most subsequent growth5,6, leaving organic vapours as the most plausible alternative, at least in the planetary boundary layer7,8,9,10. Although recent studies11,12,13 predict that low-volatility organic vapours contribute during initial growth, direct evidence has been lacking. The accelerating growth may result from increased photolytic production of condensable organic species in the afternoon2, and the presence of a possible Kelvin (curvature) effect, which inhibits organic vapour condensation on the smallest particles (the nano-Köhler theory)2,14, has so far remained ambiguous. Here we present experiments performed in a large chamber under atmospheric conditions that investigate the role of organic vapours in the initial growth of nucleated organic particles in the absence of inorganic acids and bases such as sulfuric acid or ammonia and amines, respectively. Using data from the same set of experiments, it has been shown15 that organic vapours alone can drive nucleation. We focus on the growth of nucleated particles and find that the organic vapours that drive initial growth have extremely low volatilities (saturation concentration less than 10−4.5 micrograms per cubic metre). As the particles increase in size and the Kelvin barrier falls, subsequent growth is primarily due to more abundant organic vapours of slightly higher volatility (saturation concentrations of 10−4.5 to 10−0.5 micrograms per cubic metre). We present a particle growth model that quantitatively reproduces our measurements. Furthermore, we implement a parameterization of the first steps of growth in a global aerosol model and find that concentrations of atmospheric cloud concentration nuclei can change substantially in response, that is, by up to 50 per cent in comparison with previously assumed growth rate parameterizations.
The growth of freshly formed aerosol particles can be the bottleneck in their survival to cloud condensation nuclei. It is therefore crucial to understand how particles grow in the atmosphere. Insufficient experimental data has impeded a profound understanding of nano-particle growth under atmospheric conditions. Here we study nano-particle growth in the CLOUD (Cosmics Leaving OUtdoors Droplets) chamber, starting from the formation of molecular clusters. We present measured growth rates at sub-3 nm sizes with different atmospherically relevant concentrations of sulphuric acid, water, ammonia and dimethylamine. We find that atmospheric ions and small acid-base clusters, which are not generally accounted for in the measurement of sulphuric acid vapour, can participate in the growth process, leading to enhanced growth rates. The availability of compounds capable of stabilizing sulphuric acid clusters governs the magnitude of these effects and thus the exact growth mechanism. We bring these observations into a coherent framework and discuss their significance in the atmosphere.
The abiotic and biotic drivers of rapid diversification in Andean bellflowers (Campanulaceae)
(2016)
The tropical Andes of South America, the world's richest biodiversity hotspot, are home to many rapid radiations. While geological, climatic, and ecological processes collectively explain such radiations, their relative contributions are seldom examined within a single clade. We explore the contribution of these factors by applying a series of diversification models that incorporate mountain building, climate change, and trait evolution to the first dated phylogeny of Andean bellflowers (Campanulaceae: Lobelioideae). Our framework is novel for its direct incorporation of geological data on Andean uplift into a macroevolutionary model. We show that speciation and extinction are differentially influenced by abiotic factors: speciation rates rose concurrently with Andean elevation, while extinction rates decreased during global cooling. Pollination syndrome and fruit type, both biotic traits known to facilitate mutualisms, played an additional role in driving diversification. These abiotic and biotic factors resulted in one of the fastest radiations reported to date: the centropogonids, whose 550 species arose in the last 5 million yr. Our study represents a significant advance in our understanding of plant evolution in Andean cloud forests. It further highlights the power of combining phylogenetic and Earth science models to explore the interplay of geology, climate, and ecology in generating the world's biodiversity.
Evaluation of radiation components in a global freshwater model with station-based observations
(2016)
In many hydrological models, the amount of evapotranspired water is calculated using the potential evapotranspiration (PET) approach. The main driver of several PET approaches is net radiation, whose downward components are usually obtained from meteorological input data, whereas the upward components are calculated by the model itself. Thus, uncertainties can be large due to both the input data and model assumptions. In this study, we compare the radiation components of the WaterGAP Global Hydrology Model, driven by two meteorological input datasets and two radiation setups from ERA-Interim reanalysis. We assess the performance with respect to monthly observations provided by the Baseline Surface Radiation Network (BSRN) and the Global Energy Balance Archive (GEBA). The assessment is done for the global land area and specifically for energy/water limited regions. The results indicate that there is no optimal radiation input throughout the model variants, but standard meteorological input datasets perform better than those directly obtained by ERA-Interim reanalysis for the key variable net radiation. The low number of observations for some radiation components, as well as the scale mismatch between station observations and 0.5° × 0.5° grid cell size, limits the assessment.
Binary nucleation of sulphuric acid-water particles is expected to be an important process in the free troposphere at low temperatures. SAWNUC (Sulphuric Acid Water Nucleation) is a model of binary nucleation that is based on laboratory measurements of the binding energies of sulphuric acid and water in charged and neutral clusters. Predictions of SAWNUC are compared for the first time comprehensively with experimental binary nucleation data from the CLOUD chamber at European Organization for Nuclear Research. The experimental measurements span a temperature range of 208–292 K, sulphuric acid concentrations from 1·106 to 1·109 cm−3, and distinguish between ion-induced and neutral nucleation. Good agreement, within a factor of 5, is found between the experimental and modeled formation rates for ion-induced nucleation at 278 K and below and for neutral nucleation at 208 and 223 K. Differences at warm temperatures are attributed to ammonia contamination which was indicated by the presence of ammonia-sulphuric acid clusters, detected by an Atmospheric Pressure Interface Time of Flight (APi-TOF) mass spectrometer. APi-TOF measurements of the sulphuric acid ion cluster distributions (math formula with i = 0, 1, ..., 10) show qualitative agreement with the SAWNUC ion cluster distributions. Remaining differences between the measured and modeled distributions are most likely due to fragmentation in the APi-TOF. The CLOUD results are in good agreement with previously measured cluster binding energies and show the SAWNUC model to be a good representation of ion-induced and neutral binary nucleation of sulphuric acid-water clusters in the middle and upper troposphere.
Chemical reduction-oxidation mechanisms within mantle rocks link to the terrestrial carbon cycle by influencing the depth at which magmas can form, their composition, and ultimately the chemistry of gases released into the atmosphere. The oxidation state of the uppermost mantle has been widely accepted to be unchanged over the past 3800 m.y., based on the abundance of redox-sensitive elements in greenstone belt–associated samples of different ages. However, the redox signal in those rocks may have been obscured by their complex origins and emplacement on continental margins. In contrast, the source and processes occurring during decompression melting at spreading ridges are relatively well constrained. We retrieve primary redox conditions from metamorphosed mid-oceanic ridge basalts (MORBs) and picrites of various ages (ca. 3000–550 Ma), using V/Sc as a broad redox proxy. Average V/Sc values for Proterozoic suites (7.0 ± 1.4, 2σ, n = 6) are similar to those of modern MORB (6.8 ± 1.6), whereas Archean suites have lower V/Sc (5.2 ± 0.4, n = 5). The lower Archean V/Sc is interpreted to reflect both deeper melt extraction from the uppermost mantle, which becomes more reduced with depth, and an intrinsically lower redox state. The pressure-corrected oxygen fugacity (expressed relative to the fayalite-magnetite-quartz buffer, ΔFMQ, at 1 GPa) of Archean sample suites (ΔFMQ –1.19 ± 0.33, 2σ) is significantly lower than that of post-Archean sample suites, including MORB (ΔFMQ –0.26 ± 0.44). Our results imply that the reducing Archean atmosphere was in equilibrium with Earth’s mantle, and further suggest that magmatic gases crossed the threshold that allowed a build-up in atmospheric O2 levels ca. 3000 Ma, accompanied by the first “whiffs” of oxygen in sediments of that age.
When assessing global water resources with hydrological models, it is essential to know the methodological uncertainties in the water resources estimates. The study presented here quantifies effects of the uncertainty in the spatial and temporal patterns of meteorological variables on water balance components at the global, continental and grid cell scale by forcing the global hydrological model WaterGAP 2.2 (ISI-MIP 2.1) with five state-of-the-art climate forcing input data-sets. While global precipitation over land during 1971–2000 varies between 103 500 and 111 000 km3 yr−1, global river discharge varies between 39 200 and 42 200 km3 yr−1. Temporal trends of global wa- ter balance components are strongly affected by the uncertainty in the climate forcing (except human water abstractions), and there is a need for temporal homogenization of climate forcings (in particular WFD/WFDEI). On about 10–20 % of the global land area, change of river discharge between two consecutive 30 year periods was driven more strongly by changes of human water use including dam construction than by changes in precipitation. This number increases towards the end of the 20th century due to intensified human water use and dam construction. The calibration approach of WaterGAP against observed long-term average river discharge reduces the impact of climate forcing uncertainty on estimated river discharge significantly. Different homgeneous climate forcings lead to a variation of Q of only 1.6 % for the 54 % of global land area that are constrained by discharge observations, while estimated renewable water resources in the remaining uncalibrated regions vary by 18.5 %. Uncertainties are especially high in Southeast Asia where Global Runoff Data Centre (GRDC) data availability is very sparse. By sharing already available discharge data, or installing new streamflow gauging stations in such regions, water balance uncertainties could be reduced which would lead to an improved assessment of the world’s water resources.
Ice nucleating particles over the Eastern Mediterranean measured by unmanned aircraft systems
(2016)
During an intensive field campaign on aerosol, clouds and ice nucleation in the Eastern Mediterranean in April 2016, we have measured the abundance of ice nucleating particles (INP) in the lower troposphere from unmanned aircraft systems (UAS). Aerosol samples were collected by miniaturized electrostatic precipitators onboard the UAS at altitudes up to 2.5 km. The number of INP in these samples, which are active in the deposition and condensation modes at temperatures from −20 to −30 ◦C, were analyzed immediately after collection on site using the ice nucleus counter FRIDGE. During the one month campaign we encountered a series of Saharan dust plumes that traveled at several kilometers altitude. Here we present INP data from 42 individual flights, together with aerosol number concentrations, observations of lidar backscattering, dust concentrations derived by the dust transport model DREAM (Dust Regional Atmospheric Model), and results from scanning electron microscopy. The effect of the dust plumes is reflected by the coincidence of INP with the particulate mass (PM), the lidar signal and with the predicted dust mass of the model. This suggests that mineral dust or a constituent related to dust was a major contributor to the ice nucleating properties of the aerosol. Peak concentrations of above 100 INP std.l -1 were measured at −30 ◦C. The INP concentration in elevated plumes was on average a factor of 10 higher than at ground level. Since desert dust is transported for long distances over wide areas of the globe predominantly at several km altitude we conclude that INP measurements at ground level may be of limited significance for the situation at the level of cloud formation.
Recent analysis of long-term balloon-borne measurements of Antarctic stratospheric condensation nuclei (CN) and temperature combined with global model calculations showed the wide extent of a mid stratospheric layer of new particles. Here the nucleation model SAWNUC is used to derive Antarctic stratospheric gaseous sulfuric acid profiles and to investigate the nucleation process of this CN layer. The sulfuric acid profiles were derived for an altitude range of 18-32 km between July and October by simulating air parcel trajectories that descend inside the polar vortex and calculating the sulfuric acid amount that reproduces the observations. The derived sulfuric acid concentrations (volume mixing ratios) are of the order of magnitude of 104 cm-3 (10-14) in July. In the following months the concentrations increase to about 107 cm-3 (10 exp -11) in October. They depend strongly on the temperature because a given temperature leaves only a small sulfuric acid range to reproduce the observed magnitude of CN. Ion-induced nucleation occurs. However, while it dominates nucleation at higher temperatures it has no significant influence on the nucleation rates at lower temperatures. Preexisting particles significantly reduce nucleation at sulfuric acid mixing ratios below 1 ppt. First estimates of sulfuric acid production rates range from 0.5 to about 500 molecules cm-3 s-1. A production mechanism for gaseous sulfuric acid during the Antarctic winter seems to be necessary to fully explain the observations. The derived sulfuric acid profiles compare well with midlatitude and Arctic sulfuric acid concentrations.
Small-scale thermal upwellings under the northern East African Rift from S travel time tomography
(2016)
There is a long-standing debate over how many and what types of plumes underlie the East African Rift and whether they do or do not drive its extension and consequent magmatism and seismicity. Here we present a new tomographic study of relative teleseismic S and SKS residuals that expands the resolution from previous regional studies below the northern East African Rift to image structure from the surface to the base of the transition zone. The images reveal two low-velocity clusters, below Afar and west of the Main Ethiopian Rift, that extend throughout the upper mantle and comprise several smaller-scale (about 100 km diameter), low-velocity features. These structures support those of our recent P tomographic study below the region. The relative magnitude of S to P residuals is around 3.5, which is consistent with a predominantly thermal nature of the anomalies. The S and P velocity anomalies in the low-velocity clusters can be explained by similar excess temperatures in the range of 100–200°C, consistent with temperatures inferred from other seismic, geochemical, and petrological studies. Somewhat stronger VS anomalies below Afar than west of the Main Ethiopian Rift may include an expression of volatiles and/or melt in this region. These results, together with a comparison with previous larger-scale tomographic models, indicate that these structures are likely small-scale upwellings with mild excess temperatures, rising from a regional thermal boundary layer at the base of the upper mantle.
The Fire Modeling Intercomparison Project (FireMIP), phase 1: experimental and analytical protocols
(2016)
The important role of fire in regulating vegetation community composition and contributions to emissions of greenhouse gases and aerosols make it a critical component of dynamic global vegetation models and Earth system models. Over two decades of development, a wide variety of model structures and mechanisms have been designed and incorporated into global fire models, which have been linked to different vegetation models. However, there has not yet been a systematic examination of how these different strategies contribute to model performance. Here we describe the structure of the first phase of the Fire Model Intercomparison Project (FireMIP), which for the first time seeks to systematically compare a number of models. By combining a standardized set of input data and model experiments with a rigorous comparison of model outputs to each other and to observations, we will improve the understanding of what drives vegetation fire, how it can best be simulated, and what new or improved observational data could allow better constraints on model behavior. Here we introduce the fire models used in the first phase of FireMIP, the simulation protocols applied, and the benchmarking system used to evaluate the models. The works published in this journal are distributed under the Creative Commons Attribution 3.0 License. This license does not affect the Crown copy-right work, which is re-usable under the Open Government Licence (OGL). The Creative Commons Attribution 3.0 License and the OGL are interoperable and do not conflict with, reduce, or limit each other.
In this study, we aim to reconstruct a relevant and new database of monthly zonal mean distribution of carbon dioxide (CO2) at global scale extending from the upper-troposphere (UT) to stratosphere (S). This product can be used for model and satellite validation in the UT/S, as a prior for inversion modelling and mainly to analyse a plausible feature of the stratospherictropospheric exchange as well as the stratospheric circulation and its variability. To do so, we investigate the ability of a Lagrangian trajectory model guided by ERA-Interim reanalysis to construct the CO2 abundance in the UT/S. From 10 year backward trajectories and tropospheric observations of CO2, we reconstruct upper-tropospheric and stratospheric CO2 over the period 2000–2010. The inter-comparisons of the reconstructed CO2 with mid-latitude vertical profiles measured by balloon samples as well as quasi-horizontal air samples from ER-2 aircraft during SOLVE and CONTRAIL campaigns exhibit a remarkable agreement. That demonstrates the potential of Lagrangian model to reconstruct CO2 in the UT/S. The zonal mean distribution exhibits relatively large CO2 in the tropical stratosphere due to the seasonal variation of the tropical upwelling of Brewer-Dobson circulation. During winter and spring, the tropical pipe is relatively isolated but is less confined during summer and autumn so that high CO2 values are more readily transported out of the tropics to the mid- and high latitude stratosphere. The shape of the vertical profiles suggests that relatively high CO2 above 20 km altitude mainly enter the stratosphere through tropical upwelling. CO2 mixing ratio is relatively low in the polar and tropical regions above 25 km. On average the CO2 mixing ratio decreases with altitude by 6-8 ppmv from the UT to stratosphere (e.g. up to 35 km) and is nearly constant with altitude.
The fractional release factor (FRF) gives information on the amount of a halocarbon that is released at some point in the stratosphere from its source form to the inorganic form, which can harm the ozone layer through catalytic reactions. The quantity is of major importance because it directly affects the calculation of the Ozone Depletion Potential (ODP). To apply FRF in this context, steady-state values are needed, thus representing a molecular property for a given atmospheric situation. In particular, these values should be independent of the tropospheric trends of the respective halogenated trace gases.
We analyzed the temporal evolution of FRF from ECHAM/MESSy Atmospheric Chemistry (EMAC) model simulations for several halocarbons and nitrous oxide between 1965–2011 on different mean age levels and found that the current formulation of FRF yields highly time-dependent values. We show that this is caused by the way that the tropospheric trend is handled in the current calculation method of FRF.
Taking into account chemical loss in the calculation of stratospheric mixing ratios reduces the time-dependence in correlations of different tracers. Therefore we implemented a loss term in the formulation of FRF and applied the parameterization of a "mean arrival time" to our data set.
We find that the time-dependence in FRF can almost be compensated by applying a new trend correction in the calculation of FRF. We suggest that this new method should be used to calculate time-independent FRF, which can then be used e.g. for the calculation of ODP
Desert dust is one of the most abundant ice nucleating particle types in the atmosphere. Traditionally, clay minerals were assumed to determine the ice nucleation ability of desert dust and constituted the focus of ice nucleation studies over several decades. Recently some feldspar species were identified to be ice active at much higher temperatures than clay minerals, redirecting studies to investigate the contribution of feldspar to ice nucleation on desert dust. However, so far no study has shown the atmospheric relevance of this mineral phase.
For this study four dust samples were collected after airborne transport in the troposphere from the Sahara to different locations (Crete, the Peloponnese, Canary Islands, and the Sinai Peninsula). Additionally, 11 dust samples were collected from the surface from nine of the biggest deserts worldwide. The samples were used to study the ice nucleation behavior specific to different desert dusts. Furthermore, we investigated how representative surface-collected dust is for the atmosphere by comparing to the ice nucleation activity of the airborne samples. We used the IMCA-ZINC setup to form droplets on single aerosol particles which were subsequently exposed to temperatures between 233 and 250 K. Dust particles were collected in parallel on filters for offline cold-stage ice nucleation experiments at 253–263 K. To help the interpretation of the ice nucleation experiments the mineralogical composition of the dusts was investigated. We find that a higher ice nucleation activity in a given sample at 253 K can be attributed to the K-feldspar content present in this sample, whereas at temperatures between 238 and 245 K it is attributed to the sum of feldspar and quartz content present. A high clay content, in contrast, is associated with lower ice nucleation activity. This confirms the importance of feldspar above 250 K and the role of quartz and feldspars determining the ice nucleation activities at lower temperatures as found by earlier studies for monomineral dusts. The airborne samples show on average a lower ice nucleation activity than the surface-collected ones. Furthermore, we find that under certain conditions milling can lead to a decrease in the ice nucleation ability of polymineral samples due to the different hardness and cleavage of individual mineral phases causing an increase of minerals with low ice nucleation ability in the atmospherically relevant size fraction. Comparison of our data set to an existing desert dust parameterization confirms its applicability for climate models. Our results suggest that for an improved prediction of the ice nucleation ability of desert dust in the atmosphere, the modeling of emission and atmospheric transport of the feldspar and quartz mineral phases would be key, while other minerals are only of minor importance.
Desert dust is one of the most abundant ice nucleating particle types in the atmosphere. Tra ditionally, clay minerals were assumed to determine the ice nucleation ability of desert dust and constituted the focus of ice nucleation studies. Only recently some feldspar species were identified to be ice-active at much higher temperatures than clay minerals, redirecting studies to investigate the contribution of feldspar to ice nucleation on desert dust. However, so far no study has shown the atmospheric relevance of this mineral phase.
For this study four dust samples were collected after airborne transport in the troposphere from the Sahara to different locations (Crete, the Peloponnese, Canary Islands and the Sinai Peninsula). Additionally, eleven dust samples were collected from the surface from nine of the biggest deserts worldwide. The samples were used to study the ice nucleation behavior specific to different desert dusts. Furthermore we investigated how representative ice nucleation on surface-collected dust is for that in the atmosphere by comparing to the ice nucleation activity of the airborne samples. We used the IMCA-ZINC set-up to form droplets on single aerosol particles which were subsequently exposed to temperatures between 233 - 250 K. Dust particles were collected in parallel on filters for offline cold stage ice nucleation experiments at 253 - 263 K. To help the interpretation of the results from the ice nucleation experiments the mineralogical composition of the dusts was investigated.We found that a higher ice nucleation activity in a given sample can be attributed at 253 K to the K-feldspar content present in this sample whereas at temperatures between 238 - 245 K it is attributed to the sum of feldspar and quartz content present. A high clay content on the other hand is associated with a lower ice nucleation activity of a sample. This confirms the importance of feldspar at T > 250 K and the role of quartz and feldspars determining the ice nucleation activities at lower T as found by earlier studies for monomineral dust surrogates. Furthermore, we find that milling may lead to a decrease in the ice nucleation ability of polymineral samples due to a change in mineralogical composition in the atmospherically relevant size fraction arising from the different hardness and cleavage of individual mineral phases. Comparison of our comprehensive data set to an existing desert dust parameterization confirms its applicability for climate models. Our results suggest that for an improved prediction of the ice nucleation ability of desert dust in the atmosphere, the modelling of emission and atmospheric transport of the feldspar and quartz mineral phases would be key while other minerals are only of minor importance.
In-situ single particle analysis of ice particle residuals (IPR) and out-of-cloud aerosol particles was conducted by means of laser ablation mass spectrometry during the intensive INUIT-JFJ/CLACE campaign at the high alpine research station Jungfraujoch (3580 m a.s.l.) in January/February 2013. During the four week campaign more than 70000 out-of-cloud aerosol particles and 595 IPR were analyzed covering a particle size diameter range from 100 nm to 3 μm. The IPR were sampled during 273 hours while the station was covered by mixed-phase clouds at ambient temperatures between -27 °C and -6 °C. The identification of particle types is based on laboratory studies of different types of biological, mineral and anthropogenic aerosol particles. As outcome instrument specific marker peaks for the different investigated particle types were obtained and applied to the field data. The results show that the sampled IPR contain a larger relative amount of natural, primary aerosol, like soil dust (13 %) and minerals (11 %), in comparison to out-of-cloud aerosol particles (2 % and <1 %, respectively). Additionally, anthropogenic aerosol particles, like particles from industrial emissions and lead-containing particles, were found to be more abundant in the IPR than in the out-of-cloud aerosol. The out of-cloud aerosol contained a large fraction of aged particles (30 %, including organic material and secondary inorganics), whereas this particle type was much less abundant (3 %) in the IPR. In a selected subset of the data where a direct comparison between out-of-cloud aerosol particles and IPR in air masses with similar origin was possible, a pronounced enhancement of biological particles was found in the IPR.
Size-resolved long-term measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations and hygroscopicity were conducted at the remote Amazon Tall Tower Observatory (ATTO) in the central Amazon Basin over a 1-year period and full seasonal cycle (March 2014–February 2015). The measurements provide a climatology of CCN properties characteristic of a remote central Amazonian rain forest site.
The CCN measurements were continuously cycled through 10 levels of supersaturation (S = 0.11 to 1.10 %) and span the aerosol particle size range from 20 to 245 nm. The mean critical diameters of CCN activation range from 43 nm at S = 1.10 % to 172 nm at S = 0.11 %. The particle hygroscopicity exhibits a pronounced size dependence with lower values for the Aitken mode (κAit = 0.14 ± 0.03), higher values for the accumulation mode (κAcc = 0.22 ± 0.05), and an overall mean value of κmean = 0.17 ± 0.06, consistent with high fractions of organic aerosol.
The hygroscopicity parameter, κ, exhibits remarkably little temporal variability: no pronounced diurnal cycles, only weak seasonal trends, and few short-term variations during long-range transport events. In contrast, the CCN number concentrations exhibit a pronounced seasonal cycle, tracking the pollution-related seasonality in total aerosol concentration. We find that the variability in the CCN concentrations in the central Amazon is mostly driven by aerosol particle number concentration and size distribution, while variations in aerosol hygroscopicity and chemical composition matter only during a few episodes.
For modeling purposes, we compare different approaches of predicting CCN number concentration and present a novel parametrization, which allows accurate CCN predictions based on a small set of input data.
Size-resolved long-term measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations as well as hygroscopicity were conducted at the remote Amazon Tall Tower Observatory (ATTO) in the central Amazon Basin over a one-year period and full seasonal cycle (March 2014 - February 2015). The presented measurements provide a climatology of CCN properties for a characteristic central Amazonian rain forest site.
The CCN measurements were continuously cycled through 10 levels of supersaturation (S = 0.11 to 1.10 %) and span the aerosol particle size range from 20 to 245 nm. The observed mean critical diameters of CCN activation range from 43 nm at S = 1.10 % to 172 nm at S = 0.11 %. The particle hygroscopicity exhibits a pronounced size dependence with lower values for the Aitken mode (κAit = 0.14 ± 0.03), elevated values for the accumulation mode (κAcc = 0.22 ± 0.05), and an overall mean value of κmean = 0.17 ± 0.06, consistent with high fractions of organic aerosol.
The hygroscopicity parameter κ exhibits remarkably little temporal variability: no pronounced diurnal cycles, weak seasonal trends, and few short-term variations during long-range transport events. In contrast, the CCN number concentrations exhibit a pronounced seasonal cycle, tracking the pollution-related seasonality in total aerosol concentration. We find that the variability in the CCN concentrations in the central Amazon is mostly driven by aerosol particle number concentration and size distribution, while variations in aerosol hygroscopicity and chemical composition matter only during a few episodes.
For modelling purposes, we compare different approaches of predicting CCN number concentration and present a novel parameterization, which allows accurate CCN predictions based on a small set of input data.
Profiles of CFC-11 (CCl3F) and CFC-12 (CCl2F2) of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) aboard the European satellite Envisat have been retrieved from versions MIPAS/4.61 to MIPAS/4.62 and MIPAS/5.02 to MIPAS/5.06 level-1b data using the scientific level-2 processor run by Karlsruhe Institute of Technology (KIT), Institute of Meteorology and Climate Research (IMK) and Consejo Superior de Investigaciones Científicas (CSIC), Instituto de Astrofísica de Andalucía (IAA). These profiles have been compared to measurements taken by the balloon-borne cryosampler, Mark IV (MkIV) and MIPAS-Balloon (MIPAS-B), the airborne MIPAS-STRatospheric aircraft (MIPAS-STR), the satellite-borne Atmospheric Chemistry Experiment Fourier transform spectrometer (ACE-FTS) and the High Resolution Dynamic Limb Sounder (HIRDLS), as well as the ground-based Halocarbon and other Atmospheric Trace Species (HATS) network for the reduced spectral resolution period (RR: January 2005–April 2012) of MIPAS. ACE-FTS, MkIV and HATS also provide measurements during the high spectral resolution period (full resolution, FR: July 2002–March 2004) and were used to validate MIPAS CFC-11 and CFC-12 products during that time, as well as profiles from the Improved Limb Atmospheric Spectrometer, ILAS-II. In general, we find that MIPAS shows slightly higher values for CFC-11 at the lower end of the profiles (below ∼ 15 km) and in a comparison of HATS ground-based data and MIPAS measurements at 3 km below the tropopause. Differences range from approximately 10 to 50 pptv ( ∼ 5–20 %) during the RR period. In general, differences are slightly smaller for the FR period. An indication of a slight high bias at the lower end of the profile exists for CFC-12 as well, but this bias is far less pronounced than for CFC-11 and is not as obvious in the relative differences between MIPAS and any of the comparison instruments. Differences at the lower end of the profile (below ∼ 15 km) and in the comparison of HATS and MIPAS measurements taken at 3 km below the tropopause mainly stay within 10–50 pptv (corresponding to ∼ 2–10 % for CFC-12) for the RR and the FR period. Between ∼ 15 and 30 km, most comparisons agree within 10–20 pptv (10–20 %), apart from ILAS-II, which shows large differences above ∼ 17 km. Overall, relative differences are usually smaller for CFC-12 than for CFC-11. For both species – CFC-11 and CFC-12 – we find that differences at the lower end of the profile tend to be larger at higher latitudes than in tropical and subtropical regions. In addition, MIPAS profiles have a maximum in their mixing ratio around the tropopause, which is most obvious in tropical mean profiles. Comparisons of the standard deviation in a quiescent atmosphere (polar summer) show that only the CFC-12 FR error budget can fully explain the observed variability, while for the other products (CFC-11 FR and RR and CFC-12 RR) only two-thirds to three-quarters can be explained. Investigations regarding the temporal stability show very small negative drifts in MIPAS CFC-11 measurements. These instrument drifts vary between ∼ 1 and 3 % decade−1. For CFC-12, the drifts are also negative and close to zero up to ∼ 30 km. Above that altitude, larger drifts of up to ∼ 50 % decade−1 appear which are negative up to ∼ 35 km and positive, but of a similar magnitude, above.
The organic rich Livello Bonarelli formed as a result of oxygen deficiency and carbonate dissolution in the oceans during the Cenomanian/Turonian (C/T) transition. During this Ocean Anoxic Event 2 (OAE2), a combination of factors caused increased productivity, incomplete decomposition of organic matter and widespread deposition of black shales. Although these sediments are extensively studied, the exact extent, cause, timing and duration of oceanic anoxia are debated (Sinton and Duncan, 1997; Mitchell et al., 2008). Contrasting causal mechanisms have been suggested, including stratification of the water column (Lanci et al., 2010) versus intensification of the hydrological cycle driving a dynamic ocean circulation (Trabucho-Alexandre et al., 2010). Studies on trace-elemental and (radiogenic) isotope compositions of Cenomanian marine successions have suggested a volcanic origin of OAE2, by delivering nutrients to the semi-enclosed proto-North Atlantic (Zheng et al., 2013, and references therein; Du Vivier et al., 2014). Deciphering the importance of volcanic and oceanographic processes requires tight constraints on their relative timing. Regularly occurring black cherts and shales below the Livello Bonarelli demonstrate that oceanic conditions in the Umbria-Marche Basin were punctuated by episodes of regional anoxia from the mid-Cenomanian onwards. Their hierarchical stacking pattern suggests an orbital control on the deposition of organic rich horizons (Mitchell et al., 2008; Lanci et al., 2010). Stable carbon isotope data reveal that long-term 15 variations in eccentricity paced the carbon cycle (Sprovieri et al., 2013) and sea level changes (Voigt et al., 2006) of the Late Cretaceous. Here we investigate the role of orbital forcing on climate and the carbon cycle, and, specifically, on organic-rich sedimentation prior, during, and after OAE2.
We also explore the potential for establishing an anchored astrochronology for the C/T interval in Europe. Recent improvements in the astronomical solution (La2011; Laskar et al., 2011b) and in the intercalibration of radiometric and astronomical dating techniques (Kuiper et al., 2008; Renne et al., 2013) allow the extension of the astronomical time scale into the Cretaceous. The C/T boundary in the Western Interior (USA) has been dated at 93.90 ± 0.15 Ma by intercalibration of radio-isotopic and astrochronologic time scales (Meyers et al., 2012b). Also, reinterpretation of proxy records spanning the C/T interval seems to resolve discrepancies in reported durations of the OAE2 (Sageman et al., 2006; Meyers et al., 2012a). The well-documented Italian rhythmic successions, reference sections for climatic processes in the Tethyan realm, need to be tied in with the absolute time scale. Biostratigraphic correlation to radioisotopically-dated ash beds in the Western Interior is complicated by the provinciality of faunas and floras. However, δ13C stratigraphy provides a reliable correlation tool (Gale et al., 2005) and we present a new 40Ar/39Ar age for the Thatcher bentonite from the Western Interior occurring within the mid-Cenomanian δ13C event (MCE). This study integrates the well-developed cyclostratigraphy from the Umbria-Marche Basin with radioisotopic ages from the Western Interior and derives a numerical timescale for this critical interval in Earth’s history.
The oceans at the time of the Cenomanian–Turonian transition were abruptly perturbed by a period of bottom-water anoxia. This led to the brief but widespread deposition of black organic-rich shales, such as the Livello Bonarelli in the Umbria–Marche Basin (Italy). Despite intensive studies, the origin and exact timing of this event are still debated. In this study, we assess leading hypotheses about the inception of oceanic anoxia in the Late Cretaceous greenhouse world by providing a 6 Myr long astronomically tuned timescale across the Cenomanian–Turonian boundary. We procure insights into the relationship between orbital forcing and the Late Cretaceous carbon cycle by deciphering the imprint of astronomical cycles on lithologic, physical properties, and stable isotope records, obtained from the Bottaccione, Contessa and Furlo sections in the Umbria–Marche Basin. The deposition of black shales and cherts, as well as the onset of oceanic anoxia, is related to maxima in the 405 kyr cycle of eccentricity-modulated precession. Correlation to radioisotopic ages from the Western Interior (USA) provides unprecedented age control for the studied Italian successions. The most likely tuned age for the base of the Livello Bonarelli is 94.17 ± 0.15 Ma (tuning 1); however, a 405 kyr older age cannot be excluded (tuning 2) due to uncertainties in stratigraphic correlation, radioisotopic dating, and orbital configuration. Our cyclostratigraphic framework suggests that the exact timing of major carbon cycle perturbations during the Cretaceous may be linked to increased variability in seasonality (i.e. a 405 kyr eccentricity maximum) after the prolonged avoidance of seasonal extremes (i.e. a 2.4 Myr eccentricity minimum). Volcanism is probably the ultimate driver of oceanic anoxia, but orbital periodicities determine the exact timing of carbon cycle perturbations in the Late Cretaceous. This unites two leading hypotheses about the inception of oceanic anoxia in the Late Cretaceous greenhouse world.
We present the first 3-D model of seismic P and S velocities in the crust and uppermost mantle beneath the Gulf of Aqaba and surrounding areas based on the results of passive travel time tomography. The tomographic inversion was performed based on travel time data from ∼ 9000 regional earthquakes provided by the Egyptian National Seismological Network (ENSN), and this was complemented with data from the International Seismological Centre (ISC). The resulting P and S velocity patterns were generally consistent with each other at all depths. Beneath the northern part of the Red Sea, we observed a strong high-velocity anomaly with abrupt limits that coincide with the coastal lines. This finding may indicate the oceanic nature of the crust in the Red Sea, and it does not support the concept of gradual stretching of the continental crust. According to our results, in the middle and lower crust, the seismic anomalies beneath the Gulf of Aqaba seem to delineate a sinistral shift (∼ 100 km) in the opposite flanks of the fault zone, which is consistent with other estimates of the left-lateral displacement in the southern part of the Dead Sea Transform fault. However, no displacement structures were visible in the uppermost lithospheric mantle.
Biomass burning impacts vegetation dynamics, biogeochemical cycling, atmospheric chemistry, and climate, with sometimes deleterious socio-economic impacts. Under future climate projections it is often expected that the risk of wildfires will increase. Our ability to predict the magnitude and geographic pattern of future fire impacts rests on our ability to model fire regimes, using either well-founded empirical relationships or process-based models with good predictive skill. While a large variety of models exist today, it is still unclear which type of model or degree of complexity is required to model fire adequately at regional to global scales. This is the central question underpinning the creation of the Fire Model Intercomparison Project (FireMIP), an international initiative to compare and evaluate existing global fire models against benchmark data sets for present-day and historical conditions. In this paper we review how fires have been represented in fire-enabled dynamic global vegetation models (DGVMs) and give an overview of the current state of the art in fire-regime modelling. We indicate which challenges still remain in global fire modelling and stress the need for a comprehensive model evaluation and outline what lessons may be learned from FireMIP.
Reconstructions of biomass burning from sediment-charcoal records to improve data–model comparisons
(2016)
The location, timing, spatial extent, and frequency of wildfires are changing rapidly in many parts of the world, producing substantial impacts on ecosystems, people, and potentially climate. Paleofire records based on charcoal accumulation in sediments enable modern changes in biomass burning to be considered in their long-term context. Paleofire records also provide insights into the causes and impacts of past wildfires and emissions when analyzed in conjunction with other paleoenvironmental data and with fire models. Here we present new 1000-year and 22 000-year trends and gridded biomass burning reconstructions based on the Global Charcoal Database version 3 (GCDv3), which includes 736 charcoal records (57 more than in version 2). The new gridded reconstructions reveal the spatial patterns underlying the temporal trends in the data, allowing insights into likely controls on biomass burning at regional to global scales. In the most recent few decades, biomass burning has sharply increased in both hemispheres but especially in the north, where charcoal fluxes are now higher than at any other time during the past 22 000 years. We also discuss methodological issues relevant to data–model comparisons and identify areas for future research. Spatially gridded versions of the global data set from GCDv3 are provided to facilitate comparison with and validation of global fire simulations.
MIPAS-Envisat is a satellite-borne sensor which measured vertical profiles of a wide range of trace gases from 2002 to 2012 using IR emission spectroscopy. We present geophysical validation of the MIPAS-Envisat operational retrieval (version 6.0) of N2O, CH4, CFC-12, and CFC-11 by the European Space Agency (ESA). The geophysical validation data are derived from measurements of samples collected by a cryogenic whole air sampler flown to altitudes of up to 34 km by means of large scientific balloons. In order to increase the number of coincidences between the satellite and the balloon observations, we applied a trajectory matching technique. The results are presented for different time periods due to a change in the spectroscopic resolution of MIPAS in early 2005. Retrieval results for N2O, CH4, and CFC-12 show partly good agreement for some altitude regions, which differs for the periods with different spectroscopic resolution. The more recent low spectroscopic resolution data above 20 km altitude show agreement with the combined uncertainties, while there is a tendency of the earlier high spectral resolution data set to underestimate these species above 25 km. The earlier high spectral resolution data show a significant overestimation of the mixing ratios for N2O, CH4, and CFC-12 below 20 km. These differences need to be considered when using these data. The CFC-11 results from the operation retrieval version 6.0 cannot be recommended for scientific studies due to a systematic overestimation of the CFC-11 mixing ratios at all altitudes.
Influence of sea surface roughness length parameterization on Mistral and Tramontane simulations
(2016)
The Mistral and Tramontane are mesoscale winds in southern France and above the Western Mediterranean Sea. They are phenomena well suited for studying channeling effects as well as atmosphere–land/ocean processes. This sensitivity study deals with the influence of the sea surface roughness length parameterizations on simulated Mistral and Tramontane wind speed and wind direction. Several simulations with the regional climate model COSMO-CLM were performed for the year 2005 with varying values for the Charnock parameter α. Above the western Mediterranean area, the simulated wind speed and wind direction pattern on Mistral days changes depending on the parameterization used. Higher values of α lead to lower simulated wind speeds. In areas, where the simulated wind speed does not change much, a counterclockwise rotation of the simulated wind direction is observed.
We present the characterization and application of a new gas chromatography time-of-flight mass spectrometry instrument (GC-TOFMS) for the quantitative analysis of halocarbons in air samples. The setup comprises three fundamental enhancements compared to our earlier work (Hoker et al., 2015): (1) full automation, (2) a mass resolving power R = m/Δm of the TOFMS (Tofwerk AG, Switzerland) increased up to 4000 and (3) a fully accessible data format of the mass spectrometric data. Automation in combination with the accessible data allowed an in-depth characterization of the instrument. Mass accuracy was found to be approximately 5 ppm in mean after automatic recalibration of the mass axis in each measurement. A TOFMS configuration giving R = 3500 was chosen to provide an R-to-sensitivity ratio suitable for our purpose. Calculated detection limits are as low as a few femtograms by means of the accurate mass information. The precision for substance quantification was 0.15 % at the best for an individual measurement and in general mainly determined by the signal-to-noise ratio of the chromatographic peak. Detector non-linearity was found to be insignificant up to a mixing ratio of roughly 150 ppt at 0.5 L sampled volume. At higher concentrations, non-linearities of a few percent were observed (precision level: 0.2 %) but could be attributed to a potential source within the detection system. A straightforward correction for those non-linearities was applied in data processing, again by exploiting the accurate mass information. Based on the overall characterization results, the GC-TOFMS instrument was found to be very well suited for the task of quantitative halocarbon trace gas observation and a big step forward compared to scanning, quadrupole MS with low mass resolving power and a TOFMS technique reported to be non-linear and restricted by a small dynamical range.
Amines are potentially important for atmospheric new particle formation, but their concentrations are usually low with typical mixing ratios in the pptv range or even smaller. Therefore, the demand for highly sensitive gas-phase amine measurements has emerged in the last several years. Nitrate chemical ionization mass spectrometry (CIMS) is routinely used for the measurement of gasphase sulfuric acid in the sub-pptv range. Furthermore, extremely low volatile organic compounds (ELVOCs) can be detected with a nitrate CIMS. In this study we demonstrate that a nitrate CIMS can also be used for the sensitive measurement of dimethylamine (DMA, (CH3/2NH) using the NO3−•(HNO3)1 − 2• (DMA) cluster ion signal. Calibration measurements were made at the CLOUD chamber during two different measurement campaigns. Good linearity between 0 and ~120 pptv of DMA as well as a sub-pptv detection limit of 0.7 pptv for a 10 min integration time are demonstrated at 278K and 38% RH.
We present a compact and versatile cryofocusing– thermodesorption unit, which we developed for quantitative analysis of halogenated trace gases in ambient air. Possible applications include aircraft-based in situ measurements, in situ monitoring and laboratory operation for the analysis of flask samples. Analytes are trapped on adsorptive material cooled by a Stirling cooler to low temperatures (e.g. -80°C) and subsequently desorbed by rapid heating of the adsorptive material (e.g. 200°C). The set-up involves neither the exchange of adsorption tubes nor any further condensation or refocusing steps. No moving parts are used that would require vacuum insulation. This allows for a simple and robust design. Reliable operation is ensured by the Stirling cooler, which neither contains a liquid refrigerant nor requires refilling a cryogen. At the same time, it allows for significantly lower adsorption temperatures compared to commonly used Peltier elements. We use gas chromatography – mass spectrometry (GC–MS) for separation and detection of the preconcentrated analytes after splitless injection. A substance boiling point range of approximately -80 to +150°C and a substance mixing ratio range of less than 1 ppt (pmol mol−1)to more than 500 ppt in preconcentrated sample volumes of 0.1 to 10 L of ambient air is covered, depending on the application and its analytical demands. We present the instrumental design of the preconcentration unit and demonstrate capabilities and performance through the examination of analyte breakthrough during adsorption, repeatability of desorption and analyte residues in blank tests. Examples of application are taken from the analysis of flask samples collected at Mace Head Atmospheric Research Station in Ireland using our laboratory GC–MS instruments and by data obtained during a research flight with our in situ aircraft instrument GhOSTMS (Gas chromatograph for the Observation of Tracers – coupled with a Mass Spectrometer).
We present a compact and versatile cryofocusing–thermodesorption unit, which we developed for quantitative analysis of halogenated trace gases in ambient air. Possible applications include aircraft-based in situ measurements, in situ monitoring and laboratory operation for the analysis of flask samples. Analytes are trapped on adsorptive material cooled by a Stirling cooler to low temperatures (e.g. −80 °C) and subsequently desorbed by rapid heating of the adsorptive material (e.g. +200 °C). The set-up involves neither the exchange of adsorption tubes nor any further condensation or refocusing steps. No moving parts are used that would require vacuum insulation. This allows for a simple and robust design. Reliable operation is ensured by the Stirling cooler, which neither contains a liquid refrigerant nor requires refilling a cryogen. At the same time, it allows for significantly lower adsorption temperatures compared to commonly used Peltier elements. We use gas chromatography – mass spectrometry (GC–MS) for separation and detection of the preconcentrated analytes after splitless injection. A substance boiling point range of approximately −80 to +150 °C and a substance mixing ratio range of less than 1 ppt (pmol mol−1) to more than 500 ppt in preconcentrated sample volumes of 0.1 to 10 L of ambient air is covered, depending on the application and its analytical demands. We present the instrumental design of the preconcentration unit and demonstrate capabilities and performance through the examination of analyte breakthrough during adsorption, repeatability of desorption and analyte residues in blank tests. Examples of application are taken from the analysis of flask samples collected at Mace Head Atmospheric Research Station in Ireland using our laboratory GC–MS instruments and by data obtained during a research flight with our in situ aircraft instrument GhOST-MS (Gas chromatograph for the Observation of Tracers – coupled with a Mass Spectrometer).
Knowledge about mass discrimination effects in a chemical ionization mass spectrometer (CIMS) is crucial for quantifying, e.g., the recently discovered extremely low volatile organic compounds (ELVOCs) and other compounds for which no calibration standard exists so far. Here, we present a simple way of estimating mass discrimination effects of a nitrate-based chemical ionization atmospheric pressure interface time-of-flight (CI-APi-TOF) mass spectrometer. Characterization of the mass discrimination is achieved by adding different perfluorinated acids to the mass spectrometer in amounts sufficient to deplete the primary ions significantly. The relative transmission efficiency can then be determined by comparing the decrease of signals from the primary ions and the increase of signals from the perfluorinated acids at higher masses. This method is in use already for PTR-MS; however, its application to a CI-APi-TOF brings additional difficulties, namely clustering and fragmentation of the measured compounds, which can be treated with statistical analysis of the measured data, leading to selfconsistent results. We also compare this method to a transmission estimation obtained with a setup using an electrospray ion source, a high-resolution differential mobility analyzer and an electrometer, which estimates the transmission of the instrument without the CI source. Both methods give different transmission curves, indicating non-negligible mass discrimination effects of the CI source. The absolute transmission of the instrument without the CI source was estimated with the HR-DMA method to plateau between the m=z range of 127 and 568 Th at around 1.5 %; however, for the CI source included, the depletion method showed a steady increase in relative transmission efficiency from the m=z range of the primary ion (mainly at 62 Th) to around 550 Th by a factor of around 5. The main advantages of the depletion method are that the instrument is used in the same operation mode as during standard measurements and no knowledge of the absolute amount of the measured substance is necessary, which results in a simple setup.
AirCore-HR: a high resolution column sampling to enhance the vertical description of CH₄ and CO₂
(2016)
An original and innovative sampling system called AirCore was presented by NOAA in 2010 (Karion et al., 2010). It consists of a long (> 100 m) and narrow (< 1 cm) stainless steel tube that can retain a profile of atmospheric air. The captured air sample has then to be analyzed with a gas analyzer for trace mole fraction. In this study, we introduce a new AirCore aiming at improved resolution along the vertical with the objectives to: (i) better capture the vertical distribution of CO2 and CH4, (ii) provide a tool to compare AirCores and validate the estimated vertical resolution achieved by AirCores. This AirCore-HR (high resolution) consists of a 300 m tube, combining 200 m of 1/8 in. (3.175 mm) tube and a 100 m of 1/4 in. (6.35 mm) tube. This new configuration allows to achieve a vertical resolution of 300 m up to 15 km and better than 500 m up to 22 km (if analysis of the retained sample is performed within 3 hours). The AirCore-HR was flown for the first time during the annual StratoScience campaign from CNES in August 2014 from Timmins (Ontario, Canada). High-resolution vertical profiles of CO2 and CH4 up to 25 km were successfully retrieved. These profiles revealed well defined transport structures in the troposphere (also seen in CAMS-ECMWF high resolution forecasts of CO2 and CH4 profiles) and captured the decrease of CO2 and CH4 in the stratosphere. The multi-instruments gondola from the flight carried two other low-resolution AirCore-GUF that allowed to perform direct comparisons and study the underlying processing method used to convert the sample of air to greenhouse gases vertical profiles. In particular, degrading the AirCore-HR derived profiles to the low resolution of AirCore-GUF yields an excellent match between both sets of CH4 profiles, and shows a good consistency between vertical structures of CO2 and CH4. These results fully validate the theoretical vertical resolution achievable by AirCores. Finally, the uncertainties associated with the measurements are assessed, yielding an average uncertainty below 3 ppb for CH4 and 0.25 ppm for CO2 with the major source of uncertainty coming from the potential loss of air sample on the ground and the choice of the starting and ending point of the collected air sample inside the tube. In an ideal case where the sample would be fully retained, it would be possible to know precisely the pressure at which air was sampled last and thus to improve the overall uncertainty to about 0.1 ppm for CO2 and 2 ppb for CH4.
Recently significant advances have been made in the collection, detection and characterization of ice nucleating particles (INPs). Ice nuclei are particles that facilitate the heterogeneous formation of ice within the atmospheric aerosol by lowering the free energy barrier to spontaneous nucleation and growth of ice from atmospheric water and/or vapor. The Frankfurt isostatic diffusion chamber (FRankfurt Ice nucleation Deposition freezinG Experiment: FRIDGE) is an INP collection and offline detection system that has become widely deployed and shows additional potential for ambient measurements. Since its initial development FRIDGE has gone through several iterations and improvements. Here we describe improvements that have been made in the collection and analysis techniques. We detail the uncertainties inherent in the measurement method and suggest a systematic method of error analysis for FRIDGE measurements. Thus what is presented herein should serve as a foundation for the dissemination of all current and future measurements using FRIDGE instrumentation.
Biomass burning impacts vegetation dynamics, biogeochemical cycling, atmospheric chemistry, and climate, with sometimes deleterious socio-economic impacts. Under future climate projections it is often expected that the risk of wildfires will increase. Our ability to predict the magnitude and geographic pattern of future fire impacts rests on our ability to model fire regimes, either using well-founded empirical relationships or process-based models with good predictive skill. A large variety of models exist today and it is still unclear which type of model or degree of complexity is required to model fire adequately at regional to global scales. This is the central question underpinning the creation of the Fire Model Intercomparison Project - FireMIP, an international project to compare and evaluate existing global fire models against benchmark data sets for present-day and historical conditions. In this paper we summarise the current state-of-the-art in fire regime modelling and model evaluation, and outline what lessons may be learned from FireMIP.
We discuss applications of a recently developed method for model reduction based on linear response theory of weakly coupled dynamical systems. We apply the weak coupling method to simple stochastic differential equations with slow and fast degrees of freedom. The weak coupling model reduction method results in general in a non-Markovian system; we therefore discuss the Markovianization of the system to allow for straightforward numerical integration. We compare the applied method to the equations obtained through homogenization in the limit of large timescale separation between slow and fast degrees of freedom. We numerically compare the ensemble spread from a fixed initial condition, correlation functions and exit times from a domain. The weak coupling method gives more accurate results in all test cases, albeit with a higher numerical cost.
When assessing global water resources with hydrological models, it is essential to know about methodological uncertainties. The values of simulated water balance components may vary due to different spatial and temporal aggregations, reference periods, and applied climate forcings, as well as due to the consideration of human water use, or the lack thereof. We analyzed these variations over the period 1901–2010 by forcing the global hydrological model WaterGAP 2.2 (ISIMIP2a) with five state-of-the-art climate data sets, including a homogenized version of the concatenated WFD/WFDEI data set. Absolute values and temporal variations of global water balance components are strongly affected by the uncertainty in the climate forcing, and no temporal trends of the global water balance components are detected for the four homogeneous climate forcings considered (except for human water abstractions). The calibration of WaterGAP against observed long-term average river discharge Q significantly reduces the impact of climate forcing uncertainty on estimated Q and renewable water resources. For the homogeneous forcings, Q of the calibrated and non-calibrated regions of the globe varies by 1.6 and 18.5 %, respectively, for 1971–2000. On the continental scale, most differences for long-term average precipitation P and Q estimates occur in Africa and, due to snow undercatch of rain gauges, also in the data-rich continents Europe and North America. Variations of Q at the grid-cell scale are large, except in a few grid cells upstream and downstream of calibration stations, with an average variation of 37 and 74 % among the four homogeneous forcings in calibrated and non-calibrated regions, respectively. Considering only the forcings GSWP3 and WFDEI_hom, i.e., excluding the forcing without undercatch correction (PGFv2.1) and the one with a much lower shortwave downward radiation SWD than the others (WFD), Q variations are reduced to 16 and 31 % in calibrated and non-calibrated regions, respectively. These simulation results support the need for extended Q measurements and data sharing for better constraining global water balance assessments. Over the 20th century, the human footprint on natural water resources has become larger. For 11–18% of the global land area, the change of Q between 1941–1970 and 1971–2000 was driven more strongly by change of human water use including dam construction than by change in precipitation, while this was true for only 9–13 % of the land area from 1911–1940 to 1941–1970.
Coccolithophores are an abundant phytoplankton group that exhibit remarkable diversity in their biology, ecology, and calcitic exoskeletons (coccospheres). Their extensive fossil record is testament to their important biogeochemical role and is a valuable archive of biotic responses to environmental change stretching back over 200 million years. However, to realise the potential of this archive requires an understanding of the physiological processes that underpin coccosphere architecture. Using culturing experiments on four modern coccolithophore species (Calcidiscus leptoporus, Calcidiscus quadriperforatus, Helicosphaera carteri and Coccolithus braarudii) from three long-lived families, we investigate how coccosphere architecture responds to shifts from exponential (rapid cell division) to stationary (slowed cell division) growth phases as cell physiology reacts to nutrient depletion. These experiments reveal statistical differences in cell size and the number of coccoliths per cell between these two growth phases, specifically that cells in exponential-phase growth are typically smaller with fewer coccoliths, whereas cells experiencing growth-limiting nutrient depletion have larger coccosphere sizes and greater numbers of coccoliths per cell. Although the exact numbers are species-specific, these growthphase
shifts in coccosphere geometry demonstrate that the core physiological responses of cells to nutrient depletion results in increased cell sizes and coccoliths per cell across four different coccolithophore families (Calcidiscaceae, Coccolithaceae, Isochrysidaceae, Helicosphaeraceae), a representative diversity of this phytoplankton group. Building on this, direct comparison of coccosphere geometries in modern and fossil coccolithophores enables a proxy for growth phase to be developed that allows growth responses to environmental change to be investigated throughout their evolutionary history. Our data also shows that changes in growth rate and coccoliths per cell associated with growth-phase shifts can substantially alter cellular calcite production. Coccosphere geometry is therefore a valuable tool for accessing growth information in the fossil record that provides unprecedented insights into the biotic responses of species to environmental change and its potential biogeochemical consequences.
We present the prototype of a regional climate system model based on the COSMO-CLM regional climate model coupled with several model components, analyze the performance of the couplings and present a strategy to find an optimum configuration with respect to computational costs and time to solution.
The OASIS3-MCT coupler is used to couple COSMO-CLM with two land surface models (CLM and VEG3D), a regional ocean model for the Mediterranean Sea (NEMO-MED12), two ocean models for the North and Baltic Sea (NEMO-NORDIC and TRIMNP+CICE) and the atmospheric component of an earth system model (MPI-ESM). We present a unified OASIS3-MCT interface which handles all couplings in a similar way, minimizes the model source code modifications and describes the physics and numerics of the couplings. Furthermore, we discuss solutions for specific regional coupling problems like handling of different domains, multiple usage of MCT interpolation library and efficient exchange of 3D fields.
A series of real-case simulations over Europe has been conducted and the computational performance of the couplings has been analyzed. The usage of the LUCIA tool of the OASIS3-MCT coupler enabled separation of the direct costs of: coupling, load imbalance and additional computations. The resulting limits for time to solution and costs are shown and the potential of further improvement of the computational efficiency is summarized for each coupling.
It was found that the OASIS3-MCT coupler keeps the direct coupling costs of communication and horizontal interpolation small in comparison with the costs of the additional computations and load imbalance for all investigated couplings. For the first time this could be demonstrated for an exchange of approximately 450 2D fields per time step necessary for the atmosphere-atmosphere coupling between COSMO-CLM and MPI-ESM.
A procedure for finding an optimum configuration for each of the couplings was developed considering the time to solution and costs of the simulations. The optimum configurations are presented for sequential and concurrent coupling layouts. The procedure applied can be regarded as independent on the specific coupling layout and coupling details.
Convection-permitting climate model are promising tools for improved representation of extremes, but the number of regions for which these models have been evaluated are still rather limited to make robust conclusions. In addition, an integrated interpretation of near-surface characteristics (typically temperature and precipitation) together with cloud properties is limited. The objective of this paper is to comprehensively evaluate the performance of a ‘state-of-the-art’ regional convection-permitting climate model for a mid-latitude coastal region with little orographic forcing. For this purpose, an 11-year integration with the COSMO-CLM model at Convection-Permitting Scale (CPS) using a grid spacing of 2.8 km was compared with in-situ and satellite-based observations of precipitation, temperature, cloud properties and radiation (both at the surface and the top of the atmosphere). CPS clearly improves the representation of precipitation, in especially the diurnal cycle, intensity and spatial distribution of hourly precipitation. Improvements in the representation of temperature are less obvious. In fact the CPS integration overestimates both low and high temperature extremes. The underlying cause for the overestimation of high temperature extremes was attributed to deficiencies in the cloud properties: The modelled cloud fraction is only 46 % whereas a cloud fraction of 65 % was observed. Surprisingly, the effect of this deficiency was less pronounced at the radiation balance at the top of the atmosphere due to a compensating error, in particular an overestimation of the reflectivity of clouds when they are present. Overall, a better representation of convective precipitation and a very good representation of the daily cycle in different cloud types were demonstrated. However, to overcome remaining deficiencies, additional efforts are necessary to improve cloud characteristics in CPS. This will be a challenging task due to compensating deficiencies that currently exist in ‘state-of-the-art’ models, yielding a good representation of average climate conditions. In the light of using the CPS models to study climate change it is necessary that these deficiencies are addressed in future research.
This study aims to assess the skill of regional climate models (RCMs) at reproducing the climatology of Mediterranean cyclones. Seven RCMs are considered, five of which were also coupled with an oceanic model. All simulations were forced at the lateral boundaries by the ERA-Interim reanalysis for a common 20-year period (1989–2008). Six different cyclone tracking methods have been applied to all twelve RCM simulations and to the ERA-Interim reanalysis in order to assess the RCMs from the perspective of different cyclone definitions. All RCMs reproduce the main areas of high cyclone occurrence in the region south of the Alps, in the Adriatic, Ionian and Aegean Seas, as well as in the areas close to Cyprus and to Atlas mountains. The RCMs tend to underestimate intense cyclone occurrences over the Mediterranean Sea and reproduce 24–40 % of these systems, as identified in the reanalysis. The use of grid nudging in one of the RCMs is shown to be beneficial, reproducing about 60 % of the intense cyclones and keeping a better track of the seasonal cycle of intense cyclogenesis. Finally, the most intense cyclones tend to be similarly reproduced in coupled and uncoupled model simulations, suggesting that modeling atmosphere–ocean coupled processes has only a weak impact on the climatology and intensity of Mediterranean cyclones.
This paper is a contribution to the special issue on Med-CORDEX, an international coordinated initiative dedicated to the multi-component regional climate modelling (atmosphere, ocean, land surface, river) of the Mediterranean under the umbrella of HyMeX, CORDEX, and Med-CLIVAR and coordinated by Samuel Somot, Paolo Ruti, Erika Coppola, Gianmaria Sannino, Bodo Ahrens, and Gabriel Jordà.
In this study we investigate the scaling of precipitation extremes with temperature in the Mediterranean region by assessing against observations the present day and future regional climate simulations performed in the frame of the HyMeX and MED-CORDEX programs. Over the 1979–2008 period, despite differences in quantitative precipitation simulation across the various models, the change in precipitation extremes with respect to temperature is robust and consistent. The spatial variability of the temperature–precipitation extremes relationship displays a hook shape across the Mediterranean, with negative slope at high temperatures and a slope following Clausius–Clapeyron (CC)-scaling at low temperatures. The temperature at which the slope of the temperature–precipitation extreme relation sharply changes (or temperature break), ranges from about 20 °C in the western Mediterranean to <10 °C in Greece. In addition, this slope is always negative in the arid regions of the Mediterranean. The scaling of the simulated precipitation extremes is insensitive to ocean–atmosphere coupling, while it depends very weakly on the resolution at high temperatures for short precipitation accumulation times. In future climate scenario simulations covering the 2070–2100 period, the temperature break shifts to higher temperatures by a value which is on average the mean regional temperature change due to global warming. The slope of the simulated future temperature–precipitation extremes relationship is close to CC-scaling at temperatures below the temperature break, while at high temperatures, the negative slope is close, but somewhat flatter or steeper, than in the current climate depending on the model. Overall, models predict more intense precipitation extremes in the future. Adjusting the temperature–precipitation extremes relationship in the present climate using the CC law and the temperature shift in the future allows the recovery of the temperature–precipitation extremes relationship in the future climate. This implies negligible regional changes of relative humidity in the future despite the large warming and drying over the Mediterranean. This suggests that the Mediterranean Sea is the primary source of moisture which counteracts the drying and warming impacts on relative humidity in parts of the Mediterranean region.
Mistral and tramontane wind speed and wind direction patterns in regional climate simulations
(2016)
The Mistral and Tramontane are important wind phenomena that occur over southern France and the northwestern Mediterranean Sea. Both winds travel through constricting valleys before flowing out towards the Mediterranean Sea. The Mistral and Tramontane are thus interesting phenomena, and represent an opportunity to study channeling effects, as well as the interactions between the atmosphere and land/ocean surfaces. This study investigates Mistral and Tramontane simulations using five regional climate models with grid spacing of about 50 km and smaller. All simulations are driven by ERA-Interim reanalysis data. Spatial patterns of surface wind, as well as wind development and error propagation along the wind tracks from inland France to offshore during Mistral and Tramontane events, are presented and discussed. To disentangle the results from large-scale error sources in Mistral and Tramontane simulations, only days with well simulated large-scale sea level pressure field patterns are evaluated. Comparisons with the observations show that the large-scale pressure patterns are well simulated by the considered models, but the orographic modifications to the wind systems are not well simulated by the coarse-grid simulations (with a grid spacing of about 50 km), and are reproduced slightly better by the higher resolution simulations. On days with Mistral and/or Tramontane events, most simulations underestimate (by 13 % on average) the wind speed over the Mediterranean Sea. This effect is strongest at the lateral borders of the main flow—the flow width is underestimated. All simulations of this study show a clockwise wind direction bias over the sea during Mistral and Tramontane events. Simulations with smaller grid spacing show smaller biases than their coarse-grid counterparts.
We report on HCFC-22 data acquired by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) in the reduced spectral resolution nominal observation mode. The data cover the period from January 2005 to April 2012 and the altitude range from the upper troposphere (above cloud top altitude) to about 50 km. The profile retrieval was performed by constrained nonlinear least squares fitting of modelled spectra to the measured limb spectral radiances. The spectral ν4-band at 816.5 ± 13 cm−1 was used for the retrieval. A Tikhonov-type smoothing constraint was applied to stabilise the retrieval. In the lower stratosphere, we find a global volume mixing ratio of HCFC-22 of about 185 pptv in January 2005. The rate of linear growth in the lower latitudes lower stratosphere was about 6 to 7 pptv year−1 in the period 2005–2012. The profiles obtained were compared with ACE-FTS satellite data v3.5, as well as with MkIV balloon profiles and cryosampler balloon measurements. Between 13 and 22 km, average agreement within −3 to +5 pptv (MIPAS – ACE) with ACE-FTS v3.5 profiles is demonstrated. Agreement with MkIV solar occultation balloon-borne measurements is within 10–20 pptv below 30 km and worse above, while in situ cryosampler balloon measurements are systematically lower over their full altitude range by 15–50 pptv below 24 km and less than 10 pptv above 28 km. MIPAS HCFC-22 time series below 10 km altitude are shown to agree mostly well to corresponding time series of near-surface abundances from the NOAA/ESRL and AGAGE networks, although a more pronounced seasonal cycle is obvious in the satellite data. This is attributed to tropopause altitude fluctuations and subsidence of polar winter stratospheric air into the troposphere. A parametric model consisting of constant, linear, quasi-biennial oscillation (QBO) and several sine and cosine terms with different periods has been fitted to the temporal variation of stratospheric HCFC-22 for all 10°-latitude/1-to-2-km-altitude bins. The relative linear variation was always positive, with relative increases of 40–70 % decade−1 in the tropics and global lower stratosphere, and up to 120 % decade−1 in the upper stratosphere of the northern polar region and the southern extratropical hemisphere. Asian HCFC-22 emissions have become the major source of global upper tropospheric HCFC-22. In the upper troposphere, monsoon air, rich in HCFC-22, is instantaneously mixed into the tropics. In the middle stratosphere, between 20 and 30 km, the observed trend is inconsistent with the trend at the surface (corrected for the age of stratospheric air), hinting at circulation changes. There exists a stronger positive trend in HCFC-22 in the Southern Hemisphere and a more muted positive trend in the Northern Hemisphere, implying a potential change in the stratospheric circulation over the observation period.
Under certain conditions, secondary organic aerosol (SOA) particles can exist in the atmosphere in an amorphous solid or semi-solid state. To determine their relevance to processes such as ice nucleation or chemistry occurring within particles requires knowledge of the temperature and relative humidity (RH) range for SOA to exist in these states. In the Cosmics Leaving Outdoor Droplets (CLOUD) experiment at The European Organisation for Nuclear Research (CERN), we deployed a new in situ optical method to detect the viscous state of α-pinene SOA particles and measured their transition from the amorphous highly viscous state to states of lower viscosity. The method is based on the depolarising properties of laboratory-produced non-spherical SOA particles and their transformation to non-depolarising spherical particles at relative humidities near the deliquescence point. We found that particles formed and grown in the chamber developed an asymmetric shape through coagulation. A transition to a spherical shape was observed as the RH was increased to between 35 % at −10 °C and 80 % at −38 °C, confirming previous calculations of the viscosity-transition conditions. Consequently, α-pinene SOA particles exist in a viscous state over a wide range of ambient conditions, including the cirrus region of the free troposphere. This has implications for the physical, chemical, and ice-nucleation properties of SOA and SOA-coated particles in the atmosphere.
There are strong indications that particles containing secondary organic aerosol (SOA) exhibit amorphous solid or semi-solid phase states in the atmosphere. This may facilitate heterogeneous ice nucleation and thus influence cloud properties. However, experimental ice nucleation studies of biogenic SOA are scarce. Here, we investigated the ice nucleation ability of viscous SOA particles.
The SOA particles were produced from the ozone initiated oxidation of α-pinene in an aerosol chamber at temperatures in the range from −38 to −10 °C at 5–15 % relative humidity with respect to water to ensure their formation in a highly viscous phase state, i.e. semi-solid or glassy. The ice nucleation ability of SOA particles with different sizes was investigated with a new continuous flow diffusion chamber. For the first time, we observed heterogeneous ice nucleation of viscous α-pinene SOA for ice saturation ratios between 1.3 and 1.4 significantly below the homogeneous freezing limit. The maximum frozen fractions found at temperatures between −39.0 and −37.2 °C ranged from 6 to 20 % and did not depend on the particle surface area. Global modelling of monoterpene SOA particles suggests that viscous biogenic SOA particles are indeed present in regions where cirrus cloud formation takes place. Hence, they could make up an important contribution to the global ice nucleating particle budget.
Cloud microphysical processes involving the ice phase in tropospheric clouds are among the major uncertainties in cloud formation, weather, and general circulation models. The detection of aerosol particles, liquid droplets, and ice crystals, especially in the small cloud particle-size range below 50 μm, remains challenging in mixed phase, often unstable environments. The Cloud Aerosol Spectrometer with Polarization (CASPOL) is an airborne instrument that has the ability to detect such small cloud particles and measure the variability in polarization state of their backscattered light. Here we operate the versatile Cosmics Leaving OUtdoor Droplets (CLOUD) chamber facility at the European Organization for Nuclear Research (CERN) to produce controlled mixed phase and other clouds by adiabatic expansions in an ultraclean environment, and use the CASPOL to discriminate between different aerosols, water, and ice particles. In this paper, optical property measurements of mixed-phase clouds and viscous secondary organic aerosol (SOA) are presented. We report observations of significant liquid–viscous SOA particle polarization transitions under dry conditions using CASPOL. Cluster analysis techniques were subsequently used to classify different types of particles according to their polarization ratios during phase transition. A classification map is presented for water droplets, organic aerosol (e.g., SOA and oxalic acid), crystalline substances such as ammonium sulfate, and volcanic ash. Finally, we discuss the benefits and limitations of this classification approach for atmospherically relevant concentrations and mixtures with respect to the CLOUD 8–9 campaigns and its potential contribution to tropical troposphere layer analysis.
The CARIBIC (Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container) passenger aircraft observatory performed in situ measurements at 10–12 km altitude in the South Asian summer monsoon anticyclone between June and September 2008. These measurements enable us to investigate this atmospheric region (which so far has mostly been observed from satellites) using the broad suite of trace gases and aerosol particles measured by CARIBIC. Elevated levels of a variety of atmospheric pollutants (e.g. carbon monoxide, total reactive nitrogen oxides, aerosol particles, and several volatile organic compounds) were recorded. The measurements provide detailed information about the chemical composition of air in different parts of the monsoon anticyclone, particularly of ozone precursors. While covering a range of 3500 km inside the monsoon anticyclone, CARIBIC observations show remarkable consistency, i.e. with distinct latitudinal patterns of trace gases during the entire monsoon period.
Using the CARIBIC trace gas and aerosol particle measurements in combination with the Lagrangian particle dispersion model FLEXPART, we investigated the characteristics of monsoon outflow and the chemical evolution of air masses during transport. The trajectory calculations indicate that these air masses originated mainly from South Asia and mainland Southeast Asia. Estimated photochemical ages of the air were found to agree well with transport times from a source region east of 90–95° E. The photochemical ages of the air in the southern part of the monsoon anticyclone were systematically younger (less than 7 days) and the air masses were mostly in an ozone-forming chemical mode. In its northern part the air masses were older (up to 13 days) and had unclear ozone formation or destruction potential. Based on analysis of forward trajectories, several receptor regions were identified. In addition to predominantly westward transport, we found evidence for efficient transport (within 10 days) to the Pacific and North America, particularly during June and September, and also of cross-tropopause exchange, which was strongest during June and July. Westward transport to Africa and further to the Mediterranean was the main pathway during July.
The first concerted multi-model intercomparison of halogenated very short-lived substances (VSLS) has been performed, within the framework of the ongoing Atmospheric Tracer Transport Model Intercomparison Project (TransCom). Eleven global models or model variants participated, simulating the major natural bromine VSLS, bromoform (CHBr3) and dibromomethane (CH2Br2), over a 20-year period (1993-2012). The overarching goal of TransCom-VSLS was to provide a reconciled model estimate of the stratospheric source gas injection (SGI) of bromine from these gases, to constrain the current measurement-derived range, and to investigate inter-model differences
due to emissions and transport processes. Models ran with standardised idealised chemistry, to isolate differences due to transport, and we investigated the sensitivity of results to a range of VSLS emission inventories. Models were tested in their ability to reproduce the observed seasonal and spatial distribution of VSLS at the surface, using measurements from NOAA’s long-term global monitoring network, and in the tropical troposphere, using recent aircraft measurements - including high altitude observations from the NASA Global Hawk platform.
The models generally capture the seasonal cycle of surface CHBr3 and CH2Br2 well, with a strong model measurement correlation (r ≥0.7) and a low sensitivity to the choice of emission inventory, at most sites. In a given model, the absolute model-measurement agreement is highly sensitive to the choice of emissions and inter-model differences are also apparent, even when using the same inventory, highlighting the challenges faced in evaluating such inventories at the global scale. Across the ensemble, most consistency is found within the tropics where most of the models (8 out of 11) achieve optimal agreement to surface CHBr3 observations using the lowest of the three CHBr3 emission inventories tested (similarly, 8 out of 11 models for CH2Br2). In general, the models are able to reproduce well observations of CHBr3 and CH2Br2 obtained in the tropical tropopause layer (TTL) at various locations throughout the Pacific. Zonal variability in VSLS loading in the TTL is generally consistent among models, with CHBr3 (and to a lesser extent CH2Br2) most elevated over the tropical West Pacific during boreal winter. The models also indicate the Asian Monsoon during boreal summer to be an important pathway for VSLS reaching the stratosphere, though the strength of this signal varies considerably among models.
We derive an ensemble climatological mean estimate of the stratospheric bromine SGI from CHBr3 and CH2Br2 of 2.0 (1.2-2.5) ppt, ∼57% larger than the best estimate from the most re- cent World Meteorological Organization (WMO) Ozone Assessment Report. We find no evidence for a long-term, transport-driven trend in the stratospheric SGI of bromine over the simulation period. However, transport-driven inter-annual variability in the annual mean bromine SGI is of the order of a ±5%, with SGI exhibiting a strong positive correlation with ENSO in the East Pacific