550 Geowissenschaften
Refine
Year of publication
- 2016 (88) (remove)
Document Type
- Article (67)
- Doctoral Thesis (15)
- Contribution to a Periodical (3)
- Part of a Book (1)
- Master's Thesis (1)
- Working Paper (1)
Has Fulltext
- yes (88)
Is part of the Bibliography
- no (88)
Keywords
- Atmospheric chemistry (3)
- Anthropogene Klimaänderung (2)
- Modellierung (2)
- BSRN (1)
- Bayesian network (1)
- Biogeochemistry (1)
- Bochum (1)
- CCLM (1)
- CORDEX-EAS (1)
- COSMO-CLM (1)
Ecological networks are more sensitive to plant than to animal extinction under climate change
(2016)
Impacts of climate change on individual species are increasingly well documented, but we lack understanding of how these effects propagate through ecological communities. Here we combine species distribution models with ecological network analyses to test potential impacts of climate change on >700 plant and animal species in pollination and seed-dispersal networks from central Europe. We discover that animal species that interact with a low diversity of plant species have narrow climatic niches and are most vulnerable to climate change. In contrast, biotic specialization of plants is not related to climatic niche breadth and vulnerability. A simulation model incorporating different scenarios of species coextinction and capacities for partner switches shows that projected plant extinctions under climate change are more likely to trigger animal coextinctions than vice versa. This result demonstrates that impacts of climate change on biodiversity can be amplified via extinction cascades from plants to animals in ecological networks.
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.
Das Ziel dieser Arbeit ist die Untersuchung der stratosphärischen Meridionalzirkulation mit Hilfe von chemisch aktiven Spurengasen. Diese motiviert sich durch die Tatsache, dass der Klimawandel neben den viel erforschten Auswirkungen auf die Troposphäre, auch Reaktionen in der Stratosphäre zur Folge hat, welche bisher weit weniger tiefgehend untersucht wurden. Das macht die Stratosphäre zu einem aktuellen und frequentierten Forschungsgebiet der experimentellen und theoretischen Meteorologie. Neben vereinzelten hochaufgelösten in-situ Messungen und globalen Satellitendaten sind es hier vor allem globale numerische Klima-Chemiemodelle, die für Analysen genutzt werden. Für diese Arbeit wurden Daten des EMAC-Modells (engl.: ECHAM/MESSy Atmospheric Chemistry) ausgwertet, welche im Rahmen der ESCiMo (engl.: Earth System Chemistry integrated Modelling) Initiative vom MESSy-Konsortium (engl.: Modular Earth Submodel System) erstellt wurden. Die Zielsetzung dieser Arbeit war, ob sich etwaige Änderungen des stratosphärischen Transports anhand von modellierten, chemisch aktiven, idealisierten Spurengasen feststellen lassen. Idealisiert bedeutet hierbei, dass diese Gase ein konstantes Mischungsverhältnis am Erdboden aufweisen und den identischen chemischen Prozessen unterliegen wie die realistischen Tracer. Dies hat zur Folge, dass diese Spurengase somit nicht in das Strahlungsbudget des Modells rückkoppeln und ihre Verteilung nicht durch zeitliche troposphärische Trends beeinflusst wird. Zur Analyse des stratosphärischen Transports wurden die Differenzen der monatlich gemittelten Mischungsverhältnisse zweier Zeitpunkte der verschiedenen Substanzen im Vertikalprofil betrachtet und ausgewertet, wobei hier die photolytische Lebenszeit und die Zeitskala des Transports zu berücksichtigen war. Um die Saisonalität von Transport und Chemie zu berücksichtigen, wurden dazu die Monate März, Juni, September und Dezember analysiert.
Die Resultate zeigten, dass chemisch aktive Substanzen in der Tat geeignet sind Änderungen in der Dynamik festzustellen. So stellte sich heraus, dass mit einer allgemeinen Intensivierung der stratosphärischen Meridionalzirkulation im kommenden Jahrhundert gerechnet werden kann, wobei hiervon besonders die untere Stratosphäre betroffen ist. Eine Differenzierung welche Komponente der Zirkulation (Residualtransport oder bidirektionale quasi-horizontale Mischung) hierbei von übergeordneter Bedeutung ist, konnte nicht spezifiziert werden. Um abzuschätzen, ob sich die Änderung der Zirkulation durch Änderungen in den Mischungsverhältnissen von chemisch aktiven Substanzen mit Hilfe von direkten Messungen nachweisen lässt, wurde die atmosphärische Variabilität des Modells bestimmt und mit den Veränderungen dieser Mischungsverhältnisse verglichen. Es zeigte sich, dass diese modellierte atmosphärische Variabilität zum Teil deutlich größer war, als die Differenzen der Mischungsverhältnisse und so ohne eine Vielzahl von in-situ Messungen keine eindeutige Aussage zulassen. Um eine statistisch valide Aussage treffen zu können, müssen daher mehrere Messreihen innerhalb eines Monats durchgeführt werden. Zudem stellte sich heraus, dass der Monat Juni der bestmögliche Messzeitraum ist, da hier die natürliche Variabilität am geringsten ist. Zuletzt wurden die Spurengase mit vergleichsweise kleiner chemischer Lebenszeit auf normierten N2O-Isoplethen untersucht und die Verschiebung dieser Kurve zwischen den zwei Zeitpunkten analysiert. Die so gewonnenen Resultate ließen den Schluss zu, dass sich auf diese Weise die atmosphärische Variabilität reduzieren lässt und bei Nutzung mit experimentellen Daten eine zu den Tracer-Differenzen konsistente Aussage zulässt. So bestärkte diese Methode die These, dass sich der stratosphärische Transport innerhalb des 21. Jahrhunderts wahrscheinlich verstärken wird.
Background: Replicate population pairs that diverge in response to similar selective regimes allow for an investigation of (a) whether phenotypic traits diverge in a similar and predictable fashion, (b) whether there is gradual variation in phenotypic divergence reflecting variation in the strength of natural selection among populations, (c) whether the extent of this divergence is correlated between multiple character suites (i.e., concerted evolution), and (d) whether gradual variation in phenotypic divergence predicts the degree of reproductive isolation, pointing towards a role for adaptation as a driver of (ecological) speciation. Here, we use poeciliid fishes of the genera Gambusia and Poecilia that have repeatedly evolved extremophile lineages able to tolerate high and sustained levels of toxic hydrogen sulfide (H2S) to answer these questions.
Results: We investigated evolutionary divergence in response to H2S in Gambusia spp. (and to a lesser extent Poecilia spp.) using a multivariate approach considering the interplay of life history, body shape, and population genetics (nuclear miscrosatellites to infer population genetic differentiation as a proxy for reproductive isolation). We uncovered both shared and unique patterns of evolution: most extremophile Gambusia predictably evolved larger heads and offspring size, matching a priori predictions for adaptation to sulfidic waters, while variation in adult life histories was idiosyncratic. When investigating patterns for both genera (Gambusia and Poecilia), we found that divergence in offspring-related life histories and body shape were positively correlated across populations, but evidence for individual-level associations between the two character suites was limited, suggesting that genetic linkage, developmental interdependencies, or pleiotropic effects do not explain patterns of concerted evolution. We further found that phenotypic divergence was positively correlated with both environmental H2S-concentration and neutral genetic differentiation (a proxy for gene flow).
Conclusions: Our results suggest that higher toxicity exerts stronger selection, and that divergent selection appears to constrain gene flow, supporting a scenario of ecological speciation. Nonetheless, progress toward ecological speciation was variable, partially reflecting variation in the strength of divergent selection, highlighting the complexity of selective regimes even in natural systems that are seemingly governed by a single, strong selective agent.
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.
In light of the global sea-level rise and climate change of the 21th century, it is important to look back into the recent past in order to understand what the future might hold. A multi-proxy data set was compiled to evaluate the influence of geomorphological and environmental factors, such as antecedent topography, subsidence, sea level and climate, on reef, sand apron and lagoon development in modern carbonate platforms through the Holocene. Therefore, a combination of remote sensing and morphological data from 122 modern carbonate platforms and atolls in the Atlantic, Indian and Pacific Oceans were conducted, along with a case study from the oceanic (Darwinian) barrier-reef system of Bora Bora, French Polynesia, South Pacific.
The influence of antecedent topography and platform size as factors controlling Holocene sand apron development and extension in modern atolls and carbonate platforms is hypothesized. Antecedent topography describes the elevation and relief of the underlying Pleistocene topography (karst) and determines the distance from the sea floor to the rising postglacial sea level. Maximum lagoon depth and marginal reef thickness, when available in literature, were used as proxies for antecedent topography. Sand apron proportions of 122 atolls and carbonate platforms from the Atlantic, Indian and Pacific Oceans were quantified and correlated to maximum lagoon depth, total platform area and marginal reef thickness. This study shows that sand apron proportions increase with decreasing lagoon depths. Sand apron proportions also increase with decreasing platform area. The interaction of antecedent topography and Holocene sea-level rise is responsible for variations in accommodation space and at least determines the extension of the lateral expansion of sand aprons. In general, sand apron formation started when marginal reefs approached relative sea level. Spatial and regional variations in sea-level history let sand apron formation start earlier in the Indo-Pacific region (transgressive-regressive) than in the Western Atlantic Ocean (transgressive).
The influence of sea level, antecedent topography and subsidence of a volcanic island on late Quaternary reef development was evaluated based on six rotary core transects on the barrier and fringing reefs of Bora Bora. This study was designed to revalue the Darwinian model, the subsidence theory of reef development, which genetically connects fringing reef, barrier reef and atoll development by continuous subsidence of the volcanic basement. Postglacial sea-level rise, and to a minor degree subsidence, were identified as major factors controlling Holocene reef development in that they have created accommodation space and controlled reef architecture. Antecedent topography was also an important factor because the Holocene barrier reef is located on a Pleistocene barrier reef forming a topographic high. Pleistocene soil and basalt formed the pedestal of the fringing reef. Uranium-Thorium dating shows that barrier and fringing reefs developed contemporaneously during the Holocene.
In the barrier–reef lagoon of Bora Bora, the influence of environmental factors, such as sea level and climate, tsunamis and tropical cyclones controlling Holocene sediment dynamics was evaluated based on sedimentological, paleontological, geochronological and geochemical data. The lagoonal succession comprises mixed carbonate-siliciclastic sediments overlying peat and Pleistocene soil. The multi-proxy data set shows variations in grain-size, total organic carbon (proxy for primary productivity), Ca and Cl element intensities (proxies for carbonate availability and lagoonal salinity) during the mid-late Holocene. These patterns could result from event sedimentation during storms and correlate to event deposits found in nearby Tahaa, probably induced by elevated cyclone activity. Accordingly, elevated erosion and runoff from the volcanic island and lower lagoonal salinity would be a result of rainfall during repeated cyclone landfall. However, Ti/Ca and Fe/Ca ratios as proxies for terrigenous sediment delivery peaked out in the early Holocene and declined since the mid-Holocene. Benthic foraminifera assemblages do not indicate reef-to-lagoon transport. Alternatively, higher and sustained hydrodynamic energy is probably induced by stronger trade winds and a higher-than-present sea level during the mid-late Holocene. The increase in mid-late Holocene sediment dynamics within the back-reef lagoon is supposed to display sediment-load shedding of sand aprons due to the oversteepening of slopes at sand apron/lagoon edges during their progradation rather than an increase in tropical storm activity during that time.
The influence of sea-level and climate changes on sediment import, composition and distribution in the Bora Bora lagoon during the Holocene is validated. Lagoonal facies succession comprises siderite-rich marly wackestones, foraminifera-siderite wackestones, mollusk-foraminifera marly packstones and mollusk-rich wackestones during the early-mid Holocene, and mudstones since the mid-late Holocene. During the early Holocene, enhanced weathering and iron input from the volcanic island due to wetter climate conditions led to the formation of siderite within the lagoonal sediments. The geochemical composition of these siderites shows that precipitation was driven by microbial activity and iron reduction in the presence of dissolved bicarbonate. Chemical substitutions at grain margins illustrate changes in the oxidation state and probably reflect changes in pore water chemistry due to sea-level rise and climate change (rainfall). In the late Holocene, sediment transport into the lagoon is hampered by motus on the windward side of the lagoon, which led to early submarine lithification within the lagoon.
The Earth’s surface condition we find today is a result of long exposure to metabolism of life forms. Particularly, molecular oxygen in the atmosphere is a feature which developed over time. The first substantial and lasting rise of atmospheric oxygen level happened ≈ 2.5 Ga ago, but localities are reported where transiently elevated oxygen levels appeared before this time-point. To trace the timing and circumstances of the earliest availability of free oxygen in the atmosphere is important to understand the habitats of early microbial life forms on Earth.
This thesis focuses to obtain information of oxygen levels and the related atmospheric cycling of metals in sediments of the 3.5 to 3.2 Ga Barberton Greenstone Belt. First, as iron was a ubiquitous constituent of Archean seawater, I investigated its isotopic composition in minerals of chemical sediments. Hereby, I tried to resolve the changes within the water basin on small scale sedimentary sequence cycles. Second, I focused on the minor constituents of Archean seawater. The Re-Os geochronologic system and the abundance patterns of the platinum-group elements were chosen to integrate information of oxygen promoted weathering of a large source area. To integrate information of a large time interval, the isotopes of uranium were investigated over a large stratigraphic section.
The two key findings of this thesis are:
• Quantitative oxidation of ferrous iron in surface layers of Paleoarchean seawater occurred during the onset and termination of hydrothermal FeIIaq delivery into shallow waters.
• Paleoarchean sedimentary successions of the Barberton Greenstone Belt lack any evidence of transient basin-scale oxygenation.
The Manzimnyama Iron Formation (IF, Fig Tree Group, Barberton Greenstone Belt, South Africa) has been deciphered to exist of cyclic stacks of lithostratigraphic units with varying amounts of iron oxide and carbonate minerals. In-situ femtosecond-Laser-Ablation ICP-MS iron isotope measurements showed that the majority of siderite (γ56Fe ≈ −0.5 ‰) precipitated directly from seawater of γ56Fe ≈ 0 ‰. Ferric iron from the surface layers is preserved in ≤ 1μ m hematite and in magnetite that has been grown within the consolidated sediment. During FeIIaq events, fine-grained hematite (γ56Fe ≈ 2.2 ‰) and magnetite (γ56Fe 0.5 to 0.8 ‰) indicate oxygen levels in surface waters of lower than 0.0002 μM. Upon onset and termination of iron oxide abundance, magnetite with γ56Fe ≈ 0 ‰ indicates that low concentrations of FeIIaq in surface waters were oxidized quantitatively. These observations demonstrate the existence of iron oxidation in Paleoarchean surface waters independent of FeIIaq concentration. This is the first investigation of Paleoarchean IF showing that lithostratigraphic cyclicity can be traced in iron isotopic composition of oxide minerals.
ID-ICP-MS measurement of Re, Ir, Ru, Pt and Pd, trace element (SF-ICP-MS) and ID-MCICP- MS uranium isotope determination have been applied to carbonaceous shale of the Mapepe Fm. (Fig Tree Group) after inverse Aqua Regia leaching and bulk digestion. The sediments reveal a silicified fraction which exhibits a seawater REE signature and a mixture of detrital and meteoritic PGE. Neither enrichment of the redox-sensitive elements Re or Mo nor fractionated uranium isotopes have been found on a stratigraphic interval of several hundred meters. The non-silica fraction shows no depletion of Re which indicates that the detrital material had no contact to oxidizing fluids. ID-TIMS measurements of Re and Os after the CrO3-SO4 Carius Tube method of two sample intervals showed that the Re-Os isotopic systems of the non-silica fractions are identical to two komatiite occurrences. Weltevreden Fm. and Komati Fm. rocks were uplifted, eroded and transported to the deep part of the sedimentary basin without any change to the Re-Os system. Negative fractionated uranium isotopes (γ238U = −0.41 ± 0.01 ‰) associated with detrital Ba-Cr-U occurrences suggest the existence of distal redox-processes that involve uranium species. This study demonstrates that over the time of exposure and deposition of the Mapepe Fm. sedimentation, free oxygen was not available for weathering in the catchment area.
Geoelektrische Methoden sind weit verbreitet und werden häufig zur Erkundung des oberflächennahen Untergrundes eingesetzt. Angewendet werden standardmäßig meist nur linienhafte Anordnungen der Sender- und Empfängerelektroden, die nur wenige Zehner Meter lang sind. Hierdurch haben diese Methoden nur geringe Eindringtiefen. Um größere Eindringtiefen und 3-dimensionale Informationen über den Untergrund zu erhalten, sind in der vorgestellten Studie die Empfänger- und Senderdipole in mehr oder weniger regelmäßigen Abständen über das Untersuchungsgebiet verteilt worden. Mit jeder Empfängerstation sind kontinuierlich die elektrischen Spannungen in bis zu drei Richtungen aufgezeichnet worden. Für die Einspeisungen wurde ein Rechtecksignal verwendet, das sich gut von den Störfrequenzen und den natürlichen Spannungen abhebt. Die Richtungen der Einspeisedipole sind entsprechend den örtlichen Gegebenheiten, jedoch möglichst parallel zu den Messrichtungen, gewählt worden. Zur Auswertung der erhobenen Messdaten wurde ein Programmpaket entwickelt, das eine weitestgehend automatisierte Auswertung der Daten erlaubt. Die Bestimmung der scheinbaren spezifischen Widerstände und ihrer Messfehler wurde an den fouriertransformierten Datenzeitreihen durchgeführt. Hierdurch konnten Störeinflüsse minimiert werden und es wurde möglich selbst stark verrauschte Datensätze auszuwerten. Um die erhobenen Daten interpretieren zu können sind die berechneten scheinbaren spezifischen Widerstände als Grundlage für Inversionen und Modellstudien verwendet worden. Die oben beschriebene Methode wurde im Rahmen der vorliegenden Arbeit in zwei unterschiedlichen Messgebieten angewandt.
Messgebiet im Hohen Vogelsberg
Im Juli 2007 wurde damit begonnen, die Forschungsbohrung Sichenhausen-Eschwald im Hohen Vogelsberg abzuteufen. Ziel war es, Informationen über den strukturellen Aufbau des größten geschlossenen Vulkankomplexes Mitteleuropas zu gewinnen. Die Gesteinsansprache der Tiefbohrung lieferte bereits relativ früh Hinweise darauf, dass ein großer magmatischer Körper aufgeschlossen wurde.Aufgrund der begrenzten räumlichen Aussagekraft der Bohrung und fehlender Geländebefunde war es nicht möglich, den Mechanismus der Platznahme und die Größe des Körpers näher zu beschreiben. Die Kampagne hatte das Ziel diese Lücke zu schließen und ein 3-dimensionales Modell des Untergrundes zu erstellen.In dem annähernd quadratischen Untersuchungsgebiet, das eine Fläche von ca. 25 $km²$ aufweist, wurden 20 Datenlogger zur Aufzeichnung der elektrischen Spannungen aufgebaut. Die Empfängerdipole waren zwischen 20 m und 30 m lang. Insgesamt wurden 36 Stromeinspeisungen mit Stromstärken zwischen 28 A und 40 A an 16 unterschiedlichen Positionen für jeweils 2 bis 3 verschiedene Dipolrichtungen vorgenommen. Die Einspeisedipole waren zwischen 100 m und 300 m lang. Insgesamt konnten 1.439 scheinbare spezifische Widerstände berechnet werden.Die Ergebnisse der Modellierungen und der Inversion der Daten zeigen, dass mit der Forschungsbohrung ein domartiger Körper angebohrt wurde. Anhand der Ergebnisse kann die räumliche Ausdehnung des Körpers eingegrenzt und ein vorher noch nicht kartierter Gang nachgewiesen werden.
Messgebiet im Bereich der Kinzigtalsperre
Das etwa Ost-West verlaufende Kinzigtal bildet die naturräumliche und geologische Grenze zwischen dem vulkanischen Vogelsberg im Norden und dem, in diesem Bereich aus Sedimentgesteinen aufgebauten, Spessart im Süden.Die zwischen Steinau a. d. Str. und Bad Soden-Salmünster befindliche Kinzigtalsperre dient dem Hochwasserschutz und der Regulierung des Pegels der Kinzig bei Dürreperioden. Der aufgestaute See ist relativ flach und weist im Normalstau maximale Tiefen von ca. 6~m auf. Der Stausee ist jedoch über weite Teile etwa 4~m tief. In dieser Kampagne betrug der Abstand zwischen den einzelnen Empfängerstationen etwa 100 m bis 300 m. Es wurde aufgrund der beengten Platzverhältnisse eine Dipollänge von ca. 48 m für die Einspeise- und die Empfängerdipole im Messgebiet gewählt. Insgesamt wurden 14 Empfängerstationen im Messgebiet aufgebaut, von denen sich Neun auf dem Seegrund befanden. Das Messraster orientierte sich am vermuteten Verlauf der Kinzigtalstörung. An 8 Positionen sind in 21 Richtungen elektrische Ströme mit Stärken zwischen 2,2 A und 40 A in den Untergrund eingespeist worden. Es konnten 536 scheinbare spezifische Widerstände berechnet werden. Ziel war es, den Verlauf der Störung näher zu bestimmen und die Tiefe der im Untergrund vorhandenen salinären Grundwässer zu bestimmen. Die Bestimmung des Verlaufs der Kinzigtalstörung sowie die Tiefenbestimmung der salinären Grundwässer war mit den erhobenen Daten jedoch nicht möglich.
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
Klima wird in der industriellen Moderne zum Problem des Wissens. Anders als das Wetter ist das Klima der direkten Wahrnehmung entzogen; ein Umstand, der paradoxerweise immer problematischer zu werden scheint, je mehr Wissen über die komplexen Zusammenhänge von lokalen Wetterereignissen und globalem Klima existiert. Das moderne Klima ist ein wandelbares, globales Phänomen, dessen Erkenntnis doppelt vermittelt ist, insofern es weder sinnlich erfahrbar noch verstandesmäßig erfassbar und somit auf Modellierung angewiesen ist. Diese findet jedoch nicht allein im wissenschaftlichen Kontext statt. Vielmehr ist die Herstellung von Klima als Bedingung von (Lebens-)Wirklichkeit und Zukunft trotz ihres unterschiedlichen epistemischen Status das Ergebnis wechselseitiger Einflussnahme wissenschaftlicher und nicht-wissenschaftlicher Modelle.
This PhD thesis has been carried out within an interdisciplinary cooperational project between the Deutsches Bergbau-Museum Bochum and the Goethe-Universität Frankfurt, which is dedicated to ancient Pb-Ag mining and metal production in the hinterland of the municipium Ulpiana in central Kosovo. Geochemical analysis (OM, XRD, EMP, MC-ICP-MS) of ores, metallurgical (by-) products and metal artefacts allowed to reconstruct the local chaîne opératoire and to decipher significant chronological differences between presumably Roman/late antique and medieval/early modern metallurgical processing. Pb isotope provenance studies documented the relevance of local metal production within the Roman Empire and confirmed the actual existence of a Metalla Dardanica district, which until now solely has been suspected on basis of epigraphy.
The predominant abundance of the by-products matte (Cu, Pb, Fe and Zn sulphides) and speiss (ferrous speiss: Fe-As compounds; base metal speiss: ~(Cu,Ni,Fe,Ag )x(Sb,Sn,As )y ) at smelting sites with a preliminary Roman/late antique dating points to treatment of complex polymetallic ore. Pb isotope analysis demonstrated that the mining district of Shashkoc-Janjevo (partially) supplied six of the ten investigated metallurgical sites. In this mineralisation, parageneses with elevated Cu, As and Sb abundances comprise significant proportions of particularly tennantite-tetrahedrite minerals, chalcopyrite, arsenopyrite and were generated during the early and main stages of ore formation. Later precipitated ore in contrast is marked by a significantly less versatile mineralogy and consists almost exclusively of galena, sphalerite and pyrite/marcasite. Besides increased Cu, As and Sb contents, ore from the main formation stage also exhibits generally higher Ag abundances, which are mainly hosted by fahlore and locally abundant secondary Cu sulphides (chalcocite, digenite and covellite) and oxidised phases (e.g. malachite, azurite). The higher precious metal grades of this ore type, whose geochemical signature (i.e. higher proportions of Cu, As and Sb) is mirrored by the abundance of the metallurgical by-products matte and speiss (almost exclusively found at potentially Roman/late antique smelting sites; see above), presumably were a pivotal factor leading to its preferential exploitation in earlier times. Matte and base metal-rich speiss contain notable amounts of Ag, which are mainly present in Cu-(Fe) sulphides and particularly antimonides ((Cu,Ni)2Sb, Ag3Sb), respectively. While the speiss compounds due to their close association with Pb bullion presumably were cupelled automatically, the metallurgical treatment of matte could not have been proven unambiguously, but overall certainly is highly likely.
The beneficiated ore (i.e. crushed and sorted, potentially also treated by more lavish techniques such as grinding, sieving or wet-mechanical methods) possibly was partially roasted and subsequently together with fluxes and charcoal submitted to the furnaces. The working temperatures approximately ranged between 1100 and 1400 °C. Slags from all presumably Roman/late antique dated and few of their potentially medieval/early modern analogues were produced from smelting of (partially roasted) ore with charcoal and added siliceous material, thus resulting in fayalite-dominant phase assemblages or rarely observed glassy parageneses. Even though several subtypes of fayalite slags have been established on basis of the abundance of Fe-rich oxide phases (i.e. spinel ss and wüstite), late clinopyroxene and the general solidification sequence of the slags, the process conditions (i.e. temperature, fO2, added fluxing agents) must have been widely similar; chemical variations could be explained by varying degrees of interaction of the slag melt with charcoal ash and furnace material. The other investigated metallurgical remains indicate employment of a calcareous flux, which led to formation of Ca-rich olivine-, olivine+clinopyroxene-, clinopyroxene- or melilite-type slags. These types as well as glassy slags were generated at more oxidising conditions outside the fayalite stability field (FMQ buffer equilibrium, cf. Lindsley, 1976) than their olivine-dominant analogues. Conclusions on the furnace construction could be drawn on basis of the typology of the slags, which mostly were tapped into a basin located outside the furnace, but partially (at two presumably medieval/early modern sites) also accumulated in a reservoir within the smelter.
Lead artefacts excavated in Ulpiana could be isotopically related to ores from mineralisations in its vicinity and demonstrate that the resources were at least utilised for local metal production. However, also ship wreck cargo from Israel - including several lead ingots with the inscription 'MET DARD' (Raban, 1999) - and late antique lead-glazed pottery from Serbia and Romania (Walton & Tite, 2010) could be related to a possible Kosovarian/Serbian provenance of the raw material and thus indicate flourishing trade of metal from the Metalla Dardanica district within the Roman Empire.
References:
Lindsley, D. H. (1976). Experimental studies of oxide minerals. In D. Rumble, III (Hrsg.), Oxide minerals (61-88). Reviews in Mineralogy, Volume 3. Washington, DC: Mineralogical Society of America.
Raban, A. (1999). The lead ingots from the wreck site (area K8). Journal of Roman Archaeology, Supplementary Series, 35, 179-188.
Walton, M. S., & Tite, M. S. (2010). Production technology of Roman lead-glazed pottery and its continuance into late antiquity. Archaeometry, 52(5), 733-759.
Ziel der Arbeit war es, die Flugzeitmassenspektrometrie als neue Analysemethode für die instrumentelle Analytik halogenierter Spurengase in der Luft zu etablieren. Die grundle-gende Motivation dafür ist, dass anthropogene Emissionen vieler Vertreter dieser Sub-stanzklasse einen negativen Einfluss auf die Umwelt zeigen: in der Atmosphäre agieren die Substanzen bzw. ihre Abbauprodukte als Katalysatoren für den stratosphärischen Ozonab-bau und verstärken den Strahlungsantrieb der Erde durch Absorption elektromagnetischer Strahlung im sogenannten atmosphärischen Fenster. Um diese Effekte und deren Auswir-kung quantifizieren zu können, ist es notwendig, Konzentrationen und Trends der Substan-zen in der Atmosphäre zu überwachen. Nur so können Gegenmaßnahmen wie Produktions-reglementierungen geplant und bewertet werden. In Kombination mit inverser Modellie-rung können zudem Rückschlüsse über tatsächlich emittierten Mengen gezogen werden. Dies stellt den Anspruch an die Analytik, sehr geringe Mengen dieser Gase sehr präzise quantifizieren zu können, um auch schwache Trends zu erkennen. Zudem muss die Analy-semethode die Möglichkeit zu bieten, mit der wachsenden Anzahl bekannter und zu über-wachender Substanzen Schritt zu halten. Besonders für letzteren Aspekt bietet die Flug-zeitmassenspektrometrie einen entscheidenden Vorteil gegenüber der „konventionellen“ Methode, der Quadrupolmassenspektrometrie: sie zeichnet das gesamten Massenspektrum auf ohne dadurch an Empfindlichkeit einzubüßen. Um das atmosphärische Mischungsver-hältnis von Substanzen im Bereich von pmol mol−1 bis fmol mol−1 bestimmen zu können, muss das Quadrupolmassenspektrometer im Single Ion Monitoring Modus betrieben wer-den – so wird zwar eine hohe Sensitivität erreicht, es wird aber auch nur die Intensität eines bestimmten Masse zu Ladungsverhältnisses (kurz: Masse) zu einem Zeitpunkt aufgezeich-net. Ein Flugzeitmassenspektrometer hingegen extrahiert Ionen mit einer Frequenz im Ki-loherzbereich und zeichnet für jede Extraktion das vollständige Flugzeitspektrum und da-mit Massenspektrum auf.
Aufgabe dieser Arbeit war es, ein Flugzeitmassenspektrometer mit vorgeschalteter Pro-benanreicherungseinheit sowie Gaschromatograph zur Trennung des Subtanzgemisches vor der Detektion aufzubauen und Werkzeuge zur Datenauswertung zu entwickeln. Um einen zukünftigen Feldeinsatz vorzubereiten, sollte der Aufbau möglichst kompakt, mobil und vollständig automatisiert sein. Anschließend sollte Empfindlichkeit, Präzision und dynami-scher Messbereich geprüft, optimiert und die Anwendbarkeit zur Analyse halogenierter Spurengase gezeigt werden. Die Ergebnisse aus der in der vorliegenden Arbeit präsentier-ten Geräteentwicklung finden sich in drei Publikationen wieder, welche in thematischer Reihenfolge die Probenanreicherung (Obersteiner et al., 2016b), den Vergleich von Quadrupol- und Flugzeitmassenspektrometrie (Hoker et al., 2015) sowie Eigenschaften und Anwendung des neuen Aufbaus (Obersteiner et al., 2016a) behandeln. Mit den genannten Aufsätzen ist die Arbeitsgruppe Engel weltweit die erste, welche hochpräzise Analytik ha-logenierter Spurengase routinemäßig mittels Flugzeitmassenspektrometrie durchführt. Der nächste Schritt ist der Übergang von der Laboranwendung zur Feldmessung, z.B. in Form von bodenbasierter in situ Analyse troposphärischer Luftmassen am Taunus Observatorium auf dem Kleinen Feldberg. Da es bisher keine Messstation für die hier beschriebene analy-tische Fragestellung in Deutschland gibt, könnte eine deutliche Verbesserung der Überwa-chung halogenierter Treibhausgase und ozonzerstörender Substanzen in Europa erzielt wer-den. Weiterhin wäre eine Flugzeugapplikation in Zukunft denkbar, welche neben der durch das Flugzeitmassenspektrometer abgedeckten Substanzbandbreite auch von dessen hoher möglicher Spektrenrate profitieren könnte. In Kombination mit Hochgeschwindigkeitsgas-chromatographie könnte eine bisher unerreichte Zeitauflösung der Beprobung der Atmo-sphäre mittels Gaschromatographie-Massenspektrometrie erzielt werden.
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.
Reconstructions of the vegetation of Europe during the Last Glacial Maximum (LGM) are an enigma. Pollen-based analyses have suggested that Europe was largely covered by steppe and tundra, and forests persisted only in small refugia. Climate-vegetation model simulations on the other hand have consistently suggested that broad areas of Europe would have been suitable for forest, even in the depths of the last glaciation. Here we reconcile models with data by demonstrating that the highly mobile groups of hunter-gatherers that inhabited Europe at the LGM could have substantially reduced forest cover through the ignition of wildfires. Similar to hunter-gatherers of the more recent past, Upper Paleolithic humans were masters of the use of fire, and preferred inhabiting semi-open landscapes to facilitate foraging, hunting and travel. Incorporating human agency into a dynamic vegetation-fire model and simulating forest cover shows that even small increases in wildfire frequency over natural background levels resulted in large changes in the forested area of Europe, in part because trees were already stressed by low atmospheric CO2 concentrations and the cold, dry, and highly variable climate. Our results suggest that the impact of humans on the glacial landscape of Europe may be one of the earliest large-scale anthropogenic modifications of the earth system.
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.