550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (672)
- Doctoral Thesis (116)
- Conference Proceeding (14)
- Working Paper (11)
- Book (6)
- Part of a Book (4)
- Master's Thesis (3)
- Preprint (2)
- Review (2)
- Bachelor Thesis (1)
Language
- English (832) (remove)
Has Fulltext
- yes (832)
Is part of the Bibliography
- no (832)
Keywords
- climate change (11)
- Climate change (7)
- COSMO-CLM (6)
- precipitation (6)
- Atmospheric chemistry (5)
- Palaeoclimate (5)
- Biogeochemistry (4)
- Geochemistry (4)
- Palaeoceanography (4)
- Bewässerung (3)
Institute
- Geowissenschaften (616)
- Geowissenschaften / Geographie (126)
- Biodiversität und Klima Forschungszentrum (BiK-F) (61)
- Senckenbergische Naturforschende Gesellschaft (45)
- Geographie (22)
- Biowissenschaften (21)
- Institut für Ökologie, Evolution und Diversität (6)
- Institut für sozial-ökologische Forschung (ISOE) (5)
- Frankfurt Institute for Advanced Studies (FIAS) (4)
- Extern (3)
Holocene climate was characterised by variability on multi-centennial to multi-decadal time scales. In central Europe, these fluctuations were most pronounced during winter. Here we present a new record of past winter climate variability for the last 10.8 ka based on four speleothems from Bunker Cave, Western Germany. Due to its central European location, the cave site is particularly well suited to record changes in precipitation and temperature in response to changes in the North Atlantic realm. We present high resolution records of δ18O, δ13C values and Mg/Ca ratios. We attribute changes in the Mg/Ca ratio to variations in the meteoric precipitation. The stable C isotope composition of the speleothems most likely reflects changes in vegetation and precipitation and variations in the δ18O signal are interpreted as variations in meteoric precipitation and temperature. We found cold and dry periods between 9 and 7 ka, 6.5 and 5.5 ka, 4 and 3 ka as well as between 0.7 to 0.2 ka. The proxy signals in our stalagmites compare well with other isotope records and, thus, seem representative for central European Holocene climate variability. The prominent 8.2 ka event and the Little Ice Age cold events are both recorded in the Bunker cave record. However, these events show a contrasting relationship between climate and δ18O, which is explained by different causes underlying the two climate anomalies. Whereas the Little Ice Age is attributed to a pronounced negative phase of the North Atlantic Oscillation, the 8.2 ka event was triggered by cooler conditions in the North Atlantic due to a slowdown of the Thermohaline Circulation.
The Izu–Bonin–Mariana volcanic arc is situated at a convergent plate margin where subduction initiation triggered the formation of MORB-like forearc basalts as a result of decompression melting and near-trench spreading. International Ocean Discovery Program (IODP) Expedition 352 recovered samples within the forearc basalt stratigraphy that contained unusual macroscopic globular textures hosted in andesitic glass (Unit 6, Hole 1440B). It is unclear how these andesites, which are unique in a stratigraphic sequence dominated by forearc basalts, and the globular textures therein may have formed. Here, we present detailed textural evidence, major and trace element analysis, as well as B and Sr isotope compositions, to investigate the genesis of these globular andesites. Samples consist of K2O-rich basaltic globules set in a glassy groundmass of andesitic composition. Between these two textural domains a likely hydrated interface of devitrified glass occurs, which, based on textural evidence, seems to be genetically linked to the formation of the globules. The andesitic groundmass is Cl rich (ca. 3000 µg/g ), whereas globules and the interface are Cl poor (ca. 300 µg/g ). Concentrations of fluid-mobile trace elements also appear to be fractionated in that globules and show enrichments in B, K, Rb, Cs, and Tl, but not in Ba and W relative to the andesitic groundmass, whereas the interface shows depletions in the latter, but is enriched in the former. Interestingly, globules and andesitic groundmass have identical Sr isotopic composition within analytical uncertainty ( 87Sr∕86Sr of 0.70580 ± 10 ), indicating that they likely formed from the same source. However, globules show high δ11 B (ca. + 7 ‰ ), whereas their host andesites are isotopically lighter (ca. – 1 ‰ ), potentially indicating that whatever process led to their formation either introduced heavier B isotopes to the globules, or induced stable isotope fractionation of B between globules and their groundmass. Based on the bulk of the textural information and geochemical data obtained from these samples, we conclude that these andesites likely formed as a result of the assimilation of shallowly altered oceanic crust (AOC) during forearc basaltic magmatism. Assimilation likely introduced radiogenic Sr, as well as heavier B isotopes to comparatively unradiogenic and low δ11B forearc basalt parental magmas (average 87Sr∕86Sr of 0.703284). Moreover, the globular textures are consistent with their formation being the result of fluid-melt immiscibility that was potentially induced by the rapid release of water from assimilated AOC whose escape likely formed the interface. If the globular textures present in these samples are indeed the result of fluid-melt immiscibility, then this process led to significant trace element and stable isotope fractionation. The textures and chemical compositions of the globules highlight the need for future experimental studies aimed at investigating the exsolution process with respect to potential trace element and isotopic fractionation in arc magmas that have perhaps not been previously considered.
In this thesis, laboratory investigations have been conducted to investigate several processes occurring during the melt segregation (crystal settling and compaction processes), as well as during emplacement of plutons. With the help of three different sets of centrifuge experiments rates of these three magmatic processes have been evaluated. In the first series of the centrifuge experiments, the diapiric ascent of buoyant material from two source layers at different depths was studied. Through five models, the hypothesis of ascending diapirs was tested and it was demonstrated whether a rising diapir ascends straight upward or if its ascent might be deviated by another buoyant, softer – and consequently easier to travel through – layer which is located within the overburden strata. We were interested under which conditions they can be formed. For this purpose we placed perturbations on top of both the buoyant layers; either with a set-off of both the protrusions (for three of these experiments), or with both protrusion sitting directly on top of each other (for one of the experiments). In the first experiment, we omitted the perturbations, to test which pathways diapirs take which grow from natural Rayleigh-Taylor instabilities. Three others experiments differed in the viscosity contrast between the overburden and the buoyant material. Through the experimental runs, the effects of different overburden viscosities and perturbation positions on the number of the diapirs were observed. The modeling results show that two diapirs rising from the offset perturbations do not take the same pathway through the overburden layer. Rather, each diapir takes a different pathway, with the deeper diapir piercing through its overburden while rising, regardless if it was a buoyant layer or denser overburden layers. However, when the two perturbations were situated directly above each other in the different PDMS layers, this resulted in the formation of one big diapir rather than several smaller ones, and the overburden layer was less deformed than with offset perturbations. Diapiric structures as those derived from the models without perturbation and where the perturbation are offset occur within Great Kavir Basin (Iran), where numerous salt diapirs grew from several salt horizons, which show a similar spatial distribution. The resulting structure observed in the model where the two perturbations situated directly above each other, is close to what is observed in composite batholiths such as the Flasergranitoid Zone within the Bergsträßer Odenwald Crystalline Complex (Germany). The second series of models were aimed to study crystal settling within a magma. For this purpose experiments with an artificial magma of 30 vol% olivine in 70 vol% basaltic melt were conducted to elucidate the formation mechanisms and time scales of gravitational cumulates. Through the experiments, two physical processes have been observed: (i) purely mechanical compaction, and (ii) chemical compaction induced by dissolution and re-precipitation of settled crystals. The results reveals that the mechanical settling of the dense olivine suspension occurs at about 1/6 the speed of simple Stokes settling, and a sedimentation exponent n of 4.1 is found. Evidences of chemical compaction induced by dissolution and re-precipitation of settled crystals have been highlighted by a detailed analysis of the fine structure of olivine grain boundaries. This last has revealed (1) the presence of Ca, which is characteristic only for MORB-melt, at the interface of two adjacent Ol-grains even when no melt is present; (2) a not fully crystallized boundary layer between two adjacent olivine grains. The crystal size distribution curves and the grain size growth exponent n ~3.6 indicate that diffusion controlled Ostwald ripening is the dominant crystal growth mechanism in concentrated magmatic suspensions. Finally, the formation times in natural olivine adcumulates have been calculated. The last series of centrifuge experiments deals with the crystal-melt settling-floating mechanism in a system composed of natural two pyroxene gabbro. The results have revealed a vertical evolution of the major and trace elements in the melt phase. Then, a numerical modelling of the sedimentation process of the crystals has been made in order to describe the compaction evolution with time. In comparing the numerical simulation with the centrifuge modelling, the stratification of the compacted layer in the runs is reproduced in numerical models. Moreover, on the base of the numerical and centrifuge modelling, a sedimentation exponent describing a deviation of settling in concentrated suspensions from Stokes sedimentation has been evaluated. Finally, the numerical simulation is applied to the Muskox intrusion to estimate the formation time and the melt fraction evolution in using the hindered sedimentation model calculations.
Bayesian Networks are computer-based environmental models that are frequently used to support decision-making under uncertainty. Under data scarce conditions, Bayesian Networks can be developed, parameterized, and run based on expert knowledge only. However, the efficiency of expert-based Bayesian Network modeling is limited by the difficulty in deriving model inputs in the time available during expert workshops. This thesis therefore aimed at developing a simple and robust method for deriving conditional probability tables from expert estimates in a time-efficient way. The design and application of this new elicitation and conversion method is demonstrated using a case study in Xinjiang, Northwest China. The key characteristics of this method are its time-efficiency and the approach to use different conversion tables based on varying levels of confidence. Although the method has its limitations, e.g. it can only be applied for variables with one conditioning variable; it provides the opportunity to support the parameterization of Bayesian Networks which would otherwise remain half-finished due to time constraints. In addition, a case study in the Murray-Darling Basin, Australia, is used to compare Bayesian Network types and software to improve the presentation clarity of large Bayesian Networks. Both case studies aimed at gaining insights on how to improve the applicability of Bayesian Networks to support environmental management.
Historically, the expansion of soy plantations has been a major driver of land-use/cover change (LUCC) in Brazil. While a series of recent public actions and supply-chain commitments reportedly curbed the replacement of forests by soy, the expansion of the agricultural commodity still poses a considerable threat to the Amazonian and Cerrado biomes. Identification of areas under high risk of soy expansion is thus paramount to assist conservation efforts in the region. We mapped the areas suitable for undergoing transition to soy plantations in the Legal Amazon with a machine-learning approach adopted from the ecological modeling literature. Simulated soy expansion for the year 2014 exhibited favorable validation scores compared to other LUCC models. We then used our model to simulate how potential future infrastructure improvements would affect the 2014 probabilities of soy occurrence in the region. In addition to the 2.3 Mha of planted soy in the Legal Amazon in 2014, our model identified another 14.7 Mha with high probability of soy conversion in the region given the infrastructure conditions at that time. Out of those, pastures and forests represented 9.8 and 0.4 Mha, respectively. Under the new infrastructure scenarios simulated, the Legal Amazonian area under high risk of soy conversion increased by up to 2.1 Mha (14.6%). These changes led to up to 11.4 and 51.4% increases in the high-risk of conversion areas of pastures and forests, respectively. If conversion occurs in the identified high-risk areas, at least 4.8 Pg of CO2 could be released into the atmosphere, a value that represents 10 times the total CO2 emissions of Brazil in 2014. Our results highlight the importance of targeting conservation policies and enforcement actions, including the Soy Moratorium, to mitigate future forest cover loss associated with infrastructure improvements in the region.
In situ measurements of ice crystal size distributions in tropical upper troposphere/lower stratosphere (UT/LS) clouds were performed during the SCOUT-AMMA campaign over West Africa in August 2006. The cloud properties were measured with a Forward Scattering Spectrometer Probe (FSSP-100) and a Cloud Imaging Probe (CIP) operated aboard the Russian high altitude research aircraft M-55 Geophysica with the mission base in Ouagadougou, Burkina Faso. A total of 117 ice particle size distributions were obtained from the measurements in the vicinity of Mesoscale Convective Systems (MCS). Two to four modal lognormal size distributions were fitted to the average size distributions for different potential temperature bins. The measurements showed proportionately more large ice particles compared to former measurements above maritime regions. With the help of trace gas measurements of NO, NOy, CO2, CO, and O3 and satellite images, clouds in young and aged MCS outflow were identified. These events were observed at altitudes of 11.0 km to 14.2 km corresponding to potential temperature levels of 346 K to 356 K. In a young outflow from a developing MCS ice crystal number concentrations of up to (8.3 ± 1.6) cm−3 and rimed ice particles with maximum dimensions exceeding 1.5 mm were found. A maximum ice water content of 0.05 g m−3 was observed and an effective radius of about 90 μm. In contrast the aged outflow events were more diluted and showed a maximum number concentration of 0.03 cm−3, an ice water content of 2.3 × 10−4 g m−3, an effective radius of about 18 μm, while the largest particles had a maximum dimension of 61 μm.
Close to the tropopause subvisual cirrus were encountered four times at altitudes of 15 km to 16.4 km. The mean ice particle number concentration of these encounters was 0.01 cm−3 with maximum particle sizes of 130 μm, and the mean ice water content was about 1.4 × 10−4 g m−3. All known in situ measurements of subvisual tropopause cirrus are compared and an exponential fit on the size distributions is established for modelling purposes.
A comparison of aerosol to ice crystal number concentrations, in order to obtain an estimate on how many ice particles may result from activation of the present aerosol, yielded low ratios for the subvisual cirrus cases of roughly one cloud particle per 30 000 aerosol particles, while for the MCS outflow cases this resulted in a high ratio of one cloud particle per 300 aerosol particles.
In-situ measurements of ice crystal size distributions in tropical upper troposphere/lower stratosphere (UT/LS) clouds were performed during the SCOUT-AMMA campaign over West Africa in August 2006. The cloud properties were measured with a Forward Scattering Spectrometer Probe (FSSP-100) and a Cloud Imaging Probe (CIP) operated aboard the Russian high altitude research aircraft M-55 ''Geophysica'' with the mission base in Ouagadougou, Burkina Faso. A total of 117 ice particle size distributions were obtained from the measurements in the vicinity of Mesoscale Convective Systems (MCS). Two or three modal lognormal size distributions were fitted to the average size distributions for different potential temperature bins. The measurements showed proportionate more large ice particles compared to former measurements above maritime regions. With the help of trace gas measurements of NO, NOy, CO2, CO, and O3, and satellite images clouds in young and aged MCS outflow were identified. These events were observed at altitudes of 11.0 km to 14.2 km corresponding to potential temperature levels of 346 K to 356 K. In a young outflow (developing MCS) ice crystal number concentrations of up to 8.3 cm−3 and rimed ice particles with maximum dimensions exceeding 1.5 mm were found. A maximum ice water content of 0.05 g m−3 was observed and an effective radius of about 90 μm. In contrast the aged outflow events were more diluted and showed a maximum number concentration of 0.03 cm−3, an ice water content of 2.3 × 10−4 g m−3, an effective radius of about 18 μm, while the largest particles had a maximum dimension of 61 μm.
Close to the tropopause subvisual cirrus were encountered four times at altitudes of 15 km to 16.4 km. The mean ice particle number concentration of these encounters was 0.01 cm−3 with maximum particle sizes of 130 μm, and the mean ice water content was about 1.4 × 10−4 g m−3. All known in-situ measurements of subvisual tropopause cirrus are compared and an exponential fit on the size distributions is established in order to give a parameterisation for modelling.
A comparison of aerosol to ice crystal number concentrations, in order to obtain an estimate on how many ice particles result from activation of the present aerosol, yielded low activation ratios for the subvisual cirrus cases of roughly one cloud particle per 30 000 aerosol particles, while for the MCS outflow cases this resulted in a high ratio of one cloud particle per 300 aerosol particles.
Climate change and its impacts already pose considerable challenges for societies that will further increase with global warming (IPCC, 2014a, b). Uncertainties of the climatic response to greenhouse gas emissions include the potential passing of large-scale tipping points (e.g. Lenton et al., 2008; Levermann et al., 2012; Schellnhuber, 2010) and changes in extreme meteorological events (Field et al., 2012) with complex impacts on societies (Hallegatte et al., 2013). Thus climate change mitigation is considered a necessary societal response for avoiding uncontrollable impacts (Conference of the Parties, 2010). On the other hand, large-scale climate change mitigation itself implies fundamental changes in, for example, the global energy system. The associated challenges come on top of others that derive from equally important ethical imperatives like the fulfilment of increasing food demand that may draw on the same resources. For example, ensuring food security for a growing population may require an expansion of cropland, thereby reducing natural carbon sinks or the area available for bio-energy production. So far, available studies addressing this problem have relied on individual impact models, ignoring uncertainty in crop model and biome model projections. Here, we propose a probabilistic decision framework that allows for an evaluation of agricultural management and mitigation options in a multi-impact-model setting. Based on simulations generated within the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP), we outline how cross-sectorally consistent multi-model impact simulations could be used to generate the information required for robust decision making.
Using an illustrative future land use pattern, we discuss the trade-off between potential gains in crop production and associated losses in natural carbon sinks in the new multiple crop- and biome-model setting. In addition, crop and water model simulations are combined to explore irrigation increases as one possible measure of agricultural intensification that could limit the expansion of cropland required in response to climate change and growing food demand. This example shows that current impact model uncertainties pose an important challenge to long-term mitigation planning and must not be ignored in long-term strategic decision making.
In order to achieve climate change mitigation, long-term decisions are required that must be reconciled with other societal goals that draw on the same resources. For example, ensuring food security for a growing population may require an expansion of crop land, thereby reducing natural carbon sinks or the area available for bio-energy production. Here, we show that current impact-model uncertainties pose an important challenge to long-term mitigation planning and propose a new risk-assessment and decision framework that accounts for competing interests.
Based on cross-sectorally consistent simulations generated within the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP) we discuss potential gains and limitations of additional irrigation and trade-offs of the expansion of agricultural land as two possible response measures to climate change and growing food demand. We describe an illustrative example in which the combination of both measures may close the supply demand gap while leading to a loss of approximately half of all natural carbon sinks.
We highlight current limitations of available simulations and additional steps required for a comprehensive risk assessment.
Both, gas and particle scavenging contribute to the transport of organic compounds by ice crystals in the troposphere. To simulate these processes an experimental setup was developed to form airborne ice crystals under atmospheric conditions. Experiments were performed in a wall independent reactor (WIR) installed in a walk-in cold chamber maintained constantly at -20°C. Aerosol particles were added to the carrier gas of ambient air by an aerosol generator to allow heterogeneous ice formation. Temperature variations and hydrodynamic conditions of the WIR were investigated to determine the conditions for ice crystal formation and crystal growth by vapour deposition. In detail, the dependence of temperature variations from flow rate and temperature of the physical wall as well as temperature variations with an increasing reactor depth were studied. The conditions to provide a stable aerosol concentration in the carrier gas flow were also studied. The temperature distribution inside the reactor was strongly dependent on flow rate and physical wall temperature. At an inlet temperature of -20°C, a flow rate of 30 L•min exp -1 and a physical wall temperature of +5°C turned out to provide ideal conditions for ice formation. At these conditions a sharp and stable laminar down draft "jet stream" of cold air in the centre of the reactor was produced. Temperatures measured at the chamber outlet were kept well below the freezing point in the whole reactor depth of 1.0 m. Thus, melting did not affect ice formation and crystal growth. The maximum residence time for airborne ice crystals was calculated to at 40 s. Ice crystal growth rates increased also with increasing reactor depth. The maximum ice crystal growth rate was calculated at 2.82 mg• exp -1. Further, the removal efficiency of the cleaning device for aerosol particles was 99.8% after 10 min. A reliable particle supply was attained after a preliminary lead time of 15 min. Thus, the minimum lead time was determined at 25 min. Several test runs revealed that the WIR is suitable to perform experiments with airborne ice crystals.
The continental expression of global cooling during the Miocene Climate Transition in Central Asia is poorly documented, as the tectonically active setting complicates the correlation of Neogene regional and global climatic developments. This study presents new geochemical data (CaSO4 content, carbonate δ13C and δ18O) from the endorheic alluvial‐lacustrine Aktau succession (Ili Basin, south‐east Kazakhstan) combined with findings from the previously published facies evolution. Time series analysis revealed long‐eccentricity forcing of the paleohydrology throughout the entire succession, split into several facies‐dependent segments. Orbital tuning, constrained by new laser ablation U‐Pb dates and a preexisting magnetostratigraphy, places the succession in a 5.0 Ma long interval in the middle to late Miocene (15.6 to 10.6 Ma). The long‐term water accumulation in the Ili Basin followed the timing of the Miocene Climate Transition, suggesting increased precipitation in the catchment area in response to climate cooling and stronger westerly winds. This was paced by minima of the 2.4 Ma eccentricity cycle, which favored the establishment of a discharge playa (~14.3 Ma) and a perennial lake (12.6 to 11.8 Ma). Furthermore, low obliquity amplitudes (nodes) caused a transient weakening of the westerlies at ~13.7 to 13.5 Ma and at ~12.7 Ma, resulting in negative hydrological budgets and salinization. Flooding of the windward Ili Basin coeval with aridification in the leeward basins suggests that the Tian Shan was a climate boundary already in the middle Miocene. Our results emphasize the impact of climate fluctuations on the westerlies' strength and thus on Central Asian hydrology.
Micromorphology is a suitable method to study the contents and stratigraphic relationships of pit fills. Within the ramparts of Corneşti-Iarcuri, fill layers of a pit were sampled. Th e pit fill was macroscopically divided into primary and secondary fill due to striking differences. These differences could be verified and concretized micromorphologically.
Strain localization in the lithosphere and the formation, evolution, and maintenance of resulting plate boundaries play a crucial role in plate tectonics and thermo‐chemical mantle convection. Previously activated lithospheric deformation zones often appear to maintain a “memory” of weakening, leading to tectonic inheritance within plate reorganizations including the Wilson cycle. Different mechanisms have been proposed to explain such strain localization, but it remains unclear which operates on what spatio‐temporal scales, and how to best incorporate them in large‐scale mantle convection models. Here, we analyze two candidates, (1), grain‐size sensitive rheology and, (2), damage‐style parameterizations of yield, stress which are sometimes used to approximate the former. Grain‐size reduction due to dynamic recrystallization can drive localization in the ductile domain, and grain growth provides a time‐dependent rheological hardening component potentially enabling the preservation of rheological heterogeneities. We compare the dynamic weakening and hardening effects as well as the timescales of strength evolution for a composite rheology including grain‐size dynamics with a pseudo‐plastic rheology including damage‐ (or “strain”‐) dependent weakening. We explore the implications of different proposed grain‐size evolution laws, and test to which extent strain‐dependent rheologies can mimic the weakening and hardening effects of the more complex micro‐physical behavior. Such an analysis helps to better understand the parallels and differences between various strain‐localization modeling approaches used in different tectonics and geodynamics communities. More importantly, our results contribute to efforts to identify the key ingredients of strain‐localization and damage hysteresis within plate tectonics and how to represent those in planetary‐scale modeling.
Diamonds growing in the Earth’s mantle often trap inclusions of fluids that are highly saline in composition. These fluids are thought to emerge from deep in subduction zones and may also be involved in the generation of some of the kimberlite magmas. However, the source of these fluids and the mechanism of their transport into the mantle lithosphere are unresolved. Here, we present experimental results showing that alkali chlorides are stable solid phases in the mantle lithosphere below 110 km. These alkali chlorides are formed by the reaction of subducted marine sediments with peridotite and show identical K/Na ratios to fluid inclusions in diamond. At temperatures >1100°C and low pressures, the chlorides are unstable; here, potassium is accommodated in mica and melt. The reaction of subducted sediments with peridotite explains the occurrence of Mg carbonates and the highly saline fluids found in diamonds and in chlorine-enriched kimberlite magmas.
The crude oil constituents benzene, toluene, ethylbenzene, and the three xylene isomers (BTEX) are the dominating groundwater contaminants originating from surface spill accidents by oil production facilities and with gasoline and jet fuel. Thereby BTEX posing a threat to the world´s scarce drinking water resources due to their water solubility and toxicity. An active remediation cleanup involving a BTEX event proves not only to be very expensive but almost impossible when it comes to the complete removal of contaminants from the subsurface. A favoured and common practice is combining an active remediation process focussing on the source of contamination coupled together with the monitoring of the residual contamination in the subsurface (monitored natural attenuation; MNA). MNA include all naturally occuring biological, chemical and physical processes in the subsurface. The general goal of this work was to improve the knowledge of biodegradation of aromatic hydrocarbons under anaerobic conditions in groundwater. For this groundwater and soil at the former military underground storage tank (UST) site Schäferhof – Süd near Nienburg/Weser (Niedersachsen, Germany) were sampled and analysed. The investigations were done in collaboration of the Umweltbundesamt, the universitys of Frankfurt and Bremen and the alphacon GmbH Ganderkesee. To investigate the extent of groundwater contamination, the terminal electron acceptor processes (TEAPs) and the metabolites of BTEX degradation in groundwater, six observation wells were sampled at regular intervals between January 2002 and September 2004. The wells were positioned in order to cover the upstream, the source area and the downstream of the presumed contamination source. Additionally, vertical sediment profiles were sampled and investigated with respect to spreading and concentration of BTEX in the subsurface. A large residual contamination involving BTEX is present in soil and groundwater at the studied locality. Maximum BTEX concentration values of 17 mg/kg were recorded in analysing sediment in the unsaturated zone. In the capillary fringe, values of 450 mg/kg were recorded (October 2004) and in the saturated zone maximum values of 6.7 mg/kg BTEX were detected. The groundwater samples indicate increasing BTEX concentrations in the groundwater flow direction (from 532 µg/l up to 3300 µg/l (mean values)). Biodegradation of aromatic hydrocarbons under anaerobic conditions in the sub surface at contaminated sites is characterised by generation of metabolites. From the monoaromatic hydrocarbons BTEX metabolites such as benzoic acid (BA) and the methylated homologs and C1-and C2-benzyl-succinic acids (BSA) are generated as intermediates. A solid-phase extraction method based on octadecyl-bonded silica sorbent has been developed to concentrate such metabolite compounds from water samples followed by derivatization and gas chromatography/mass spectrometry (GC/MS) of the extracts. The recovery rate range between 75 and 97%. The method detection limit was 0.8 µg/l. Organic acids were identified as metabolic by-products of biodegradation. Benzoic acid, C1-, C2- and C3-benzoic acid were determined in all contaminated wells with considerable concentrations. Furthermore, the depletion of the dominant terminal electron acceptors (TEAs) oxygen, nitrate, and sulphate and the production of dissolved ferrous iron and methane in groundwater indicate biological mediated processes in the plume evidently proving the occurrence of NA. A large overlap of different redox zones at the studied part of the plume has been observed. A important finding in this study is the strong influence of groundwater level fluctuations on the BTEX concentration in groundwater. A very dry summer in 2003 was recorded during the monitoring period, resulting on site in a drop of the groundwater level to 1.7 m and a concomitant increase of BTEX concentrations from 240 µg/l to 1300 µg/l. The groundwater level fluctuations, natural degradation and retention processes essentially influence BTEX concentrations in the groundwater. Groundwater level fluctuations have by far a stronger influence than the influence of biological degradation. Increasing BTEX concentrations are hence not a consequence of limited biological degradation. Another part of the study was to observe the isotopic fractionation of the electron acceptor Fe(III), due to biologically mediated reduction of Fe(III) to the watersoluble Fe(II) at the site and first field data are presented. Both groundwater and sediment samples were analysed with respect to their Fe isotopic compositions using high mass resolution Multi Collector-Inductively Coupled Plasma-Mass Spectrometry (MC-ICP-MS). The delta56Fe -values of groundwater samples taken from observation wells located downstream of the source area were isotopically lighter than delta56Fe -values obtained from groundwater in the uncontaminated well. The Fe isotopic composition of most parts of the sediment profile was similar to the Fe isotopic composition of uncontaminated groundwater. Thus, a significant iron isotope fractionation can be observed between sediment and groundwater downstream of the BTEX contamination.
This work analyses several granitic bodies of the Variscan Orogen of Central and Western Europe in order to improve our knowledge about different aspects of their evolution, regarding their ascent and emplacement mechanisms, as well as their deformation history. In the Iberian Massif two granitoid bodies, namely the La Bazana pluton and the Nisa-Alburquerque batholith, were studied in order to decipher their ascent and emplacement history. The La Bazana pluton is a small, sub-circular body in map view that intruded into rocks of the Ossa-Morena Zone in the core of a late upright antiform. Its three-dimensional drop-pipe shape, its internal dome foliation pattern and the structure of the host rock suggest that the magma ascended and emplaced diapirically. The Nisa-Alburquerque batholith is a large body that intruded into rocks of the Central Iberian Zone, the Central Unit, and the Ossa-Morena Zone. Its cartographic shape is elongate and parallel to the NW—SE to WNW—ESE Variscan structures. In the light of the available structural data and the gravimetric models, the intrusion is viewed as a continuous lateral magma flow from the eastern root guided towards the west through the southern limb of a kilometre-scale antiform. As mass-transfer mechanisms, a combination of rigid translation of the country rocks, stoping, and possibly ballooning is proposed. In the Bohemian Massif several small granitoid bodies showing a strong solid-state deformation were studied in order to integrate their tectonometamorphic history in the geotectonic framework of the south-western Bohemian Massif, focusing principally on the deformation phase referred to as D3. Four ductile deformation phases are proposed for the study area. D1 produced high-temperature fabrics under upper amphibolite to granulite facies conditions. Its kinematics is unknown. D2 occurred under amphibolite to upper greenschist facies conditions under N—S to NNW—SSE compression. It is responsible for a subvertical NW—SE striking foliation in migmatites developed under dextral simple shear and for the deformation at the Bayerischer Pfahl shear-zone system at its earlier stages. Many granitoid dykes and stocks were found to be affected by sinistral shear along subvertical planes trending ENE to ESE. Since this deformation, which is called D3 in the present work, is not compatible with a N—S to NNW—SSE compression, it is proposed that these sinistral shear zones in granites do not belong to the Bayerischer Pfahl shear-zone system and constitute themselves a separated one, which is called “D3 shear-zone system”. D3 took place under upper greenschist to lower amphibolite facies conditions (~480-550°C). Both the intrusion and the deformation of the granites affected by D3 occurred at deep to intermediate levels of the crust, whereas the deformation took place under NE—SW compression. Datings on two of the deformed granites yielded 324.4 ± 0.8 Ma and 315.0 ± 1.0 Ma: Thus, the age of D3 is most probably ~315 Ma. The intrusion of most of the sheared granitoids was pre-kinematic with respect to D3. After D3 the N—S to NNW—SSE compression which governed D2 was restored, giving way to the next deformation phase D4, which was linked to further deformation at and next to the principal shears of the Bayerischer Pfahl shear-zone system under greenschist facies conditions. The causes for the change of the stress field leading to a NE—SW compression during D3 might be related to (1) global changes in the dynamics of the tectonic plates in late Variscan times, (2) orogenic collapse leading to the sinking of the Teplá-Barrandian and lateral extrusion of the surrounding Moldanubian rocks, (3) distortion of the regional stress field by local intrusion of large stocks, such as the Saldenburg granite of the Fürstenstein Massif, or (4) distortion of the regional stress field due to the existence of ephemeral releasing bends in the Bayerischer Pfahl shear zone during its early evolution.
The calcareous substrate of spring-fed fens makes them unique islands of biodiversity, hosting endangered, vulnerable, and protected vascular plants. Hence, spring-fed fens ecosystems require special conservation attention because many of them are destroyed (e.g. drained, forested) and it is extremely difficult or even impossible to restore the unique hydrogeological and geochemical conditions enabling their function. The long-term perspective of paleoecological studies allows indication of former wetland ecosystem states and provides understanding of their development over millennia. To examine the late Holocene dynamics of a calcareous spring-fed fen (Raganu Mire) ecosystem on the Baltic Sea coast (Latvia) in relation to environmental changes, substrate and human activity, we have undertaken high-resolution analyses of plant macrofossils, pollen, mollusc, stable carbon (δ13C) and oxygen (δ18O) isotopes combined with radiocarbon dating (AMS) in three coring locations. Our study revealed that peat deposits began accumulating ca. 7000 cal. yr BP and calcareous deposits (tufa) from 1450 cal. yr BP, coinciding with regional hydrological changes. Several fire events occurred between 4000 and 1600 cal. yr BP, which appeared to have had a limited effect on local vegetation. The most significant changes in the forest and peatland ecosystems were at 3200 cal. yr BP associated with a dry climate stage and high fire activity, and then between 1400 and 500 cal. yr BP potentially associated with temperature changes during the Medieval Climate Anomaly (MCA) and Little Ice Age. Hydrological disturbances in the peatland catchment from 1400 cal. yr BP were most likely strengthened by human activity (deforestation) in this region. The relationship between the development of this peatland and changes in its catchment area, such as land cover changes or fluctuations in groundwater levels, suggest that protection and restoration of spring-fed fen ecosystems should also include the surrounding catchment. The presence of calcareous sediments, as well as appropriate temperature and local hydrological conditions appear to be the most crucial factors controlling Cladium marisus populations in our site - currently at the eastern limit of its distribution in Europe.
The future physiology of marine phytoplankton will be impacted by a range of changes in global ocean conditions, including salinity regimes that vary spatially and on a range of short- to geological timescales. Coccolithophores have global ecological and biogeochemical significance as the most important calcifying marine phytoplankton group. Previous research has shown that the morphology of their exoskeletal calcified plates (coccoliths) responds to changing salinity in the most abundant coccolithophore species, Emiliania huxleyi. However, the extent to which these responses may be strain-specific is not well established. Here we investigated the growth response of six strains of E. huxleyi under low (ca. 25) and high (ca. 45) salinity batch culture conditions and found substantial variability in the magnitude and direction of response to salinity change across strains. Growth rates declined under low and high salinity conditions in four of the six strains but increased under both low and high salinity in strain RCC1232 and were higher under low salinity and lower under high salinity in strain PLYB11. When detailed changes in coccolith and coccosphere size were quantified in two of these strains that were isolated from contrasting salinity regimes (coastal Norwegian low salinity of ca. 30 and Mediterranean high salinity of ca. 37), the Norwegian strain showed an average 26% larger mean coccolith size at high salinities compared to low salinities. In contrast, coccolith size in the Mediterranean strain showed a smaller size trend (11% increase) but severely impeded coccolith formation in the low salinity treatment. Coccosphere size similarly increased with salinity in the Norwegian strain but this trend was not observed in the Mediterranean strain. Coccolith size changes with salinity compiled for other strains also show variability, strongly suggesting that the effect of salinity change on coccolithophore morphology is likely to be strain specific. We propose that physiological adaptation to local conditions, in particular strategies for plasticity under stress, has an important role in determining ecotype responses to salinity.
Spatial variations of nitrogen trace gas emissions from tropical mountain forests in Nyungwe, Rwanda
(2012)
Globally, tropical forest soils represent the second largest source of N2O and NO. However, there is still considerable uncertainty on the spatial variability and soil properties controlling N trace gas emission. Therefore, we carried out an incubation experiment with soils from 31 locations in the Nyungwe tropical mountain forest in southwestern Rwanda. All soils were incubated at three different moisture levels (50, 70 and 90 % water filled pore space (WFPS)) at 17 °C. Nitrous oxide emission varied between 4.5 and 400 μg N m−2 h−1, while NO emission varied from 6.6 to 265 μg N m−2 h−1. Mean N2O emission at different moisture levels was 46.5 ± 11.1 (50 %WFPS), 71.7 ± 11.5 (70 %WFPS) and 98.8 ± 16.4 (90 %WFPS) μg N m−2 h−1, while mean NO emission was 69.3 ± 9.3 (50 %WFPS), 47.1 ± 5.8 (70 %WFPS) and 36.1 ± 4.2 (90 %WFPS) μg N m−2 h−1. The latter suggests that climate (i.e. dry vs. wet season) controls N2O and NO emissions. Positive correlations with soil carbon and nitrogen indicate a biological control over N2O and NO production. But interestingly N2O and NO emissions also showed a positive correlation with free iron and a negative correlation with soil pH (only N2O). The latter suggest that chemo-denitrification might, at least for N2O, be an important production pathway. In conclusion improved understanding and process based modeling of N trace gas emission from tropical forests will benefit from spatially explicit trace gas emission estimates linked to basic soil property data and differentiating between biological and chemical pathways for N trace gas formation.
Globally, tropical forest soils represent the second largest source of N2O and NO. However, there is still considerable uncertainty on the spatial variability and soil properties controlling N trace gas emission. To investigate how soil properties affect N2O and NO emission, we carried out an incubation experiment with soils from 31 locations in the Nyungwe tropical mountain forest in southwestern Rwanda. All soils were incubated at three different moisture levels (50, 70 and 90% water filled pore space (WFPS)) at 17 °C. Nitrous oxide emission varied between 4.5 and 400 μg N m−2 h−1, while NO emission varied from 6.6 to 265 μg N m−2 h−1. Mean N2O emission at different moisture levels was 46.5 ± 11.1 (50% WFPS), 71.7 ± 11.5 (70% WFPS) and 98.8 ± 16.4 (90% WFPS) μg N m−2 h−1, while mean NO emission was 69.3 ± 9.3 (50% WFPS), 47.1 ± 5.8 (70% WFPS) and 36.1 ± 4.2 (90% WFPS) μg N m−2 h−1. The latter suggests that climate (i.e. dry vs. wet season) controls N2O and NO emissions. Positive correlations with soil carbon and nitrogen indicate a biological control over N2O and NO production. But interestingly N2O and NO emissions also showed a negative correlation (only N2O) with soil pH and a positive correlation with free iron. The latter suggest that chemo-denitrification might, at least for N2O, be an important production pathway. In conclusion improved understanding and process based modeling of N trace gas emission from tropical forests will not only benefit from better spatial explicit trace gas emission and basic soil property monitoring, but also by differentiating between biological and chemical pathways for N trace gas formation.
Tropical forest soils are a significant source for the greenhouse gas N2O as well as for NO, a precursor of tropospheric ozone. However, current estimates are uncertain due to the limited number of field measurements. Furthermore, there is considerable spatial and temporal variability of N2O and NO emissions due to the variation of environmental conditions such as soil properties, vegetation characteristics and meteorology. In this study we used a process-based model (ForestDNDC-tropica) to estimate N2O and NO emissions from tropical highland forest (Nyungwe) soils in southwestern Rwanda. To extend the model inputs to regional scale, ForestDNDC-tropica was linked to an exceptionally large legacy soil dataset. There was agreement between N2O and NO measurements and the model predictions though the ForestDNDC-tropica resulted in considerable lower emissions for few sites. Low similarity was specifically found for acidic soil with high clay content and reduced metals, indicating that chemo-denitrification processes on acidic soils might be under-represented in the current ForestDNDC-tropica model. The results showed that soil bulk density and pH are the most influential factors driving spatial variations in soil N2O and NO emissions for tropical forest soils. The area investigated (1113 km2) was estimated to emit ca. 439 ± 50 t N2O-N yr−1 (2.8–5.5 kg N2O-N ha−1 yr−1) and 244 ± 16 t NO-N yr−1 (0.8–5.1 kg N ha−1 yr−1). Consistent with less detailed studies, we confirm that tropical highland rainforest soils are a major source of atmospheric N2O and NO.
The accessories perovskite, pyrochlore, zirconolite, calzirtite and melanite from carbonatites and carbonate-rich foidites from the Kaiserstuhl are variously suited for the in situ determination of their U–Pb ages and Sr, Nd- and Hf-isotope ratios by LA-ICP-MS. The 143Nd/144Nd ratios may be determined precisely in all five phases, the 176Hf/177Hf ratios only in calzirtite and the 87Sr/86Sr ratios in perovskites and pyrochlores. The carbonatites and carbonate-rich foidites belong to one of the three magmatic groups that Schleicher et al. (1990) distinguished in the Kaiserstuhl on the basis of their Sr, Nd and Pb isotope ratios. Tephrites, phonolites and essexites (nepheline monzogabbros) form the second and limburgites (nepheline basanites) and olivine nephelinites the third. Our 87Sr/86Sr isotope data from the accessories overlap with the carbonatite and olivine nephelinite fields defined by Schleicher et al. (1990) but exhibit a much narrower range. These and the εNd and εHf values plot along the mantle array in the field of oceanic island basalts relatively close to mid-ocean ridge basalts. Previously reported K–Ar, Ar–Ar and fission track ages for the Kaiserstuhl lie between 16.2 and 17.8 Ma. They stem entirely from the geologically older tephrites, phonolites and essexites. No ages existed so far for the geologically younger carbonatites and carbonate-rich foidites except for one apatite fission track age (15.8 Ma). We obtained precise U–Pb ages for zirconolites and calzirtites of 15.66, respectively 15.5 Ma (± 0.1 2σ) and for pyrochlores of 15.35 ± 0.24 Ma. Only the perovskites from the Badberg soevite yielded a U–P concordia age of 14.56 ± 0.86 Ma while the perovskites from bergalites (haüyne melilitites) only gave 206Pb/238U and 208Pb/232Th ages of 15.26 ± 0.21, respectively, 15.28 ± 0.48 Ma. The main Kaiserstuhl rock types were emplaced over a time span of 1.6 Ma almost 1 million years before the carbonatites and carbonate-rich foidites. These were emplaced within only 0.32 Ma.
A massive occurrence of microbial carbonates, including abundant sponge remains, within the Devonian Elbingerode Reef Complex was likely deposited in a former cavity of the fore-reef slope during the early Frasnian. It is suggested that the formation of microbial carbonate was to a large part favored by the activity of heterotrophic, i.e., sulfate-reducing bacteria, in analogy to Quaternary coral reef microbialites. The Elbingerode Reef Complex is an example of an oceanic or Darwinian barrier reef system. In modern barrier reef settings, microbialite formation is commonly further facilitated by weathering products from the central volcanic islands. The Devonian microbialites of the Elbingerode Reef Complex occur in the form of reticulate and laminated frameworks. Reticulate framework is rich in hexactinellid glass sponges, the tissue decay of which led to the formation of abundant micrite as well as peloidal and stromatactis textures. Supposed calcimicrobes such as Angusticellularia (formerly Angulocellularia) and Frutexites, also known from cryptic habitats, were part of the microbial association. The microbial degradation of sponge tissue likely also contributed to the laminated framework accretion as evidenced by the occurrence of remains of so-called “keratose” demosponges. Further typical textures in the microbialite of the Elbingerode Reef Complex include zebra limestone, i.e., the more or less regular intercalation of microbial carbonate and cement. Elevated concentrations of magnesium in the microbialite as compared to the surrounding metazoan (stromatoporoid-coral) reef limestone suggests that the microbialite of the Elbingerode Reef Complex was initially rich in high-magnesium calcite, which would be yet another parallel to modern, cryptic coral reef microbial carbonates. Deposition and accretion of the microbialite largely occurred in oxygenated seawater with suboxic episodes as indicated by the trace element (REE + Y) data.
Accretion rates of Holocene tropical coral reefs in three areas in the Atlantic, Pacific and Indian Oceans have been quantified in 79 dated core sections in 34 reef cores from Belize, the Maldives and French Polynesia. Holocene vertical reef accretion rate averages 5.05 m/kyr and has decreased during the past 10 kyr. Accretion rates in branched and massive coral facies are statistically similar. Reef accretion rate is positively correlated with the rate of sea‐level rise, that is the degree of creation of accommodation space, and with climate as expressed in a Holocene sea surface temperature anomaly. Accommodation space is also created by subsidence, but at a rate one to two orders of magnitude lower than that created by glacio‐eustasy (0.04 to 0.16 m/kyr). Lagoonal background sedimentation in adjacent reef lagoons averages 0.89 m/kyr as measured in 72 dated core sections in 28 cores. Lagoonal carbonate sedimentation on top of underlying mangrove peat usually starts after a considerable hiatus of ca 3 kyr on average. The lagoonal background sedimentation rate increased during the Holocene, probably due to deepening. The differences between vertical reef accretion and lagoonal background sedimentation rates are a major factor in the production of the widely known saucer shapes typical of tropical reefs and carbonate platforms, that is the creation of unfilled accommodation space. Reef core recovery, used as a proxy for reef consolidation, and core depth exhibit a statistically negative correlation based on data from 326 core barrels. Recovery and marine cement abundance (average volume 8.6%) also decrease from windward to leeward core positions. These observations are presumably a result of both a decrease in the rate of sea‐level rise that is the increase in time available for submarine cementation during the Holocene and the amount of flushing of reef interstices by marine waters.
Constructive waterfalls
(1911)
The excavation of valleys by waterfalls is one of the best known and most effective processes by which rivers cut down the surface of the earth. The influence of waterfalls is usually regarded as solely destructive, and as always helping to lower the land. They undermine and cut backward the rock faces over which they fall : by this recession they excavate deep gorges ; and the existence of these gorges enables the adjacent country to be lowered to the level of the valIey floors. The waterfalls, moreover, empty any lakes they rnay reach in their retreat, while the ravines below the falls may drain the springs and thus desiccate the neighbouring hihlands. Observations in various countries had suggested to me that waterfalls may sometimes be constructive in stead of destructive, and that they may reserse their usual procedure, advancing instead of retreating, filling valleys instead of excavating them, and forrning alluvial plains and lakes instead of destroying them. The best illustrations I have seen of such advancing, constructive waterfalls are on some rivers of Dalmatia and Bosnia, where they occur in various stages of development. ...
The assumption that mankind is able to have an in uence on global or regional climate, respectively, due to the emission of greenhouse gases, is often discussed. This assumption is both very important and very obscure. In consequence, it is necessary to clarify definitively which meteorological elements (climate parameters) are in uencend by the anthropogenic climate impact, and to which extent in which regions of the world. In addition, to be able to interprete such an information properly, it is also necessary to know the magnitude of the different climate signals due to natural variability (for example due to volcanic or solar activity) and the magnitide of stochastic climate noise. The usual tool of climatologists, general circulation models (GCM) suffer from the problem that they are at least quantitatively uncertain with regard to the regional patterns of the behaviour of climate elements and from the lack of accurate information about long-term (decadal and centennial) forcing. In contrast to that, statistical methods as used in this study have the advantage to test hypotheses directly based on observational data. So, we focus to the very reality of climate variability as it has occurred in the past. We apply two strategies of time series analyis with regard to the observed climate variables under consideration. First, each time series is splitted into its variation components. This procedure is called 'structure-oriented time series separation'. The second strategy called 'cause-oriented time series separation' matches various time series representing various forcing mechanisms with those representing the climate behaviour (climate elements). In this way it can be assessed which part of observed climate variability can be explained by this (combined) forcing and which part remains unexplained.
We present simulations with the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the Arctic winter 2002/2003. We integrated a Lagrangian denitrification scheme into the three-dimensional version of CLaMS that calculates the growth and sedimentation of nitric acid trihydrate (NAT) particles along individual particle trajectories. From those, we derive the HNO3 downward flux resulting from different particle nucleation assumptions. The simulation results show a clear vertical redistribution of total inorganic nitrogen (NOy), with a maximum vortex average permanent NOy removal of over 5 ppb in late December between 500 and 550 K and a corresponding increase of NOy of over 2 ppb below about 450 K. The simulated vertical redistribution of NOy is compared with balloon observations by MkIV and in-situ observations from the high altitude aircraft Geophysica. Assuming a globally uniform NAT particle nucleation rate of 3.4·10−6 cm−3 h−1 in the model, the observed denitrification is well reproduced. In the investigated winter 2002/2003, the denitrification has only moderate impact (<=10%) on the simulated vortex average ozone loss of about 1.1 ppm near the 460 K level. At higher altitudes, above 600 K potential temperature, the simulations show significant ozone depletion through NOx-catalytic cycles due to the unusual early exposure of vortex air to sunlight.
We present simulations with the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the Arctic winter 2002/2003. We integrated a Lagrangian denitrification scheme into the three-dimensional version of CLaMS that calculates the growth and sedimentation of nitric acid trihydrate (NAT) particles along individual particle trajectories. From those, we derive the HNO3 downward flux resulting from different particle nucleation assumptions. The simulation results show a clear vertical redistribution of total inorganic nitrogen ( ), with a maximum vortex average permanent removal of over 5ppb in late December between 500 and 550K and a corresponding increase of of over 2ppb below about 450K. The simulated vertical redistribution of is compared with balloon observations by MkIV and in-situ observations from the high altitude aircraft Geophysica. Assuming a globally uniform NAT particle nucleation rate of 7.8x10-6cm-3h-1 in the model, the observed denitrification is well reproduced.
In the investigated winter 2002/2003, the denitrification has only moderate impact (≤14%) on the simulated vortex average ozone loss of about 1.1ppm near the 460K level. At higher altitudes, above 600K potential temperature, the simulations show significant ozone depletion through -catalytic cycles due to the unusual early exposure of vortex air to sunlight.
The Match method for the quantification of polar chemical ozone loss is investigated mainly with respect to the impact of the transport of air masses across the vortex edge. For the winter 2002/03, we show that significant transport across the vortex edge occurred and was simulated by the Chemical Lagrangian Model of the Stratosphere. In-situ observations of inert tracers and ozone from HAGAR on the Geophysica aircraft and balloon-borne sondes, and remote observations from MIPAS on the ENVISAT satellite were reproduced well by CLaMS. The model even reproduced a small vortex remnant that remained a distinct feature until June 2003 and was also observed in-situ by a balloon-borne whole air sampler. We use this CLaMS simulation to quantify the impact of transport across the vortex edge on ozone loss estimates from the Match method. We show that a time integration of the determined vortex average ozone loss rates, as performed in Match, results in a larger ozone loss than the polar vortex average ozone loss in CLaMS. The determination of the Match ozone loss rates is also influenced by the transport of air across the vortex edge. We use the model to investigate how the sampling of the ozone sondes on which Match is based represents the vortex average ozone loss rate. Both the time integration of ozone loss and the determination of ozone loss rates for Match are evaluated using the winter 2002/2003 CLaMS simulation. These impacts can explain the majority of the differences between CLaMS and Match column ozone loss. While the investigated effects somewhat reduce the apparent discrepancy in January ozone loss rates reported earlier, a distinct discrepancy between simulations and Match remains. However, its contribution to the accumulated ozone loss over the winter is not large.
The Match method for quantification of polar chemical ozone loss is investigated mainly with respect to the impact of mixing across the vortex edge onto this estimate. We show for the winter 2002/03 that significant mixing across the vortex edge occurred and was accurately modeled by the Chemical Lagrangian Model of the Stratosphere. Observations of inert tracers and ozone in-situ from HAGAR on the Geophysica aircraft and sondes and also remote from MIPAS on ENVISAT were reproduced well. The model even reproduced a small vortex remnant that was isolated until June 2003 and was observed in-situ by a balloon-borne whole air sampler. We use this CLaMS simulation to quantify the impact of cross vortex edge mixing on the results of the Match method. It is shown that a time integration of the determined vortex average ozone loss rates as performed in Match results in larger ozone loss than the polar vortex average ozone loss in CLaMS. Also, the determination of the Match ozone loss rates can be influenced by mixing. This is especially important below 430 K, where ozone outside the vortex is lower than inside and the vortex boundary is not a strong transport barrier. This effect and further sampling effects cause an offset between vortex average ozone loss rates derived from Match and deduced from CLaMS with an even sampling for the entire vortex. Both, the time-integration of ozone loss and the determination of ozone loss rates for Match are evaluated using the winter 2002/03 CLaMS simulation. These impacts can explain the differences between CLaMS and Match column ozone loss. While the investigated effects somewhat reduce the apparent discrepancy in January ozone loss rates, a discrepancy between simulations and Match remains. However, its contribution to the accumulated ozone loss over the winter is not large.
We present a new experimental dataset on the impact of the heavy halogens chlorine, bromine and iodine on the Raman water bands concerning pressure and their concentration at room temperature. These experiments were conducted at ambient temperature, with variations in halogen concentration and pressure ranging from 0 to 1.4 GPa.
The strength of the Raman water band shift change increases with the ionic size from chlorine, over bromine, to iodine. Our experiments further demonstrate that increased pressure diminishes the impact of the halogen shift change to a varying extent for each of the three halogens. This finding can have significant implications for the salinity calculation of fluid inclusions in minerals such as quartz or olivine. Particularly in the low salinity range, the concentration can be markedly underestimated if the pressure effect is neglected. For experiments in diamond anvil cells involving halogens dissolved in water, the change in Raman water band shifts can serve either as a new tool to monitor pressure, or to monitor the salinity.
Highlights
• We find DBrfluid/melt = 1.19 to 3.92 for experimental Br degassing from basalt magma into aqueous fluids.
• D <1 under almost dry conditions propose only little Br degassing for dry intra-plate volcanism relative to volcanic arcs.
• An annual global Br flux of 23.5–72.9 × 109 g/y into the atmosphere was calculated.
Abstract
We present the first in-situ partitioning data for bromine between a natural basaltic melt and a coexisting fluid. For this study hydrothermal diamond anvil cell experiments at pressures up to 1.7 GPa were conducted. We combined laser heating to melt the basalt glass with external heating to lower the temperature gradient in the cell and to initiate circulation for the aqueous fluid. Bromine concentrations were measured in-situ with X-ray fluorescence in the basaltic melts, glasses, and in the fluid. From the results we calculated partition coefficients of DBrfluid/melt = 1.19 to 3.92 in the range of 0.4 to 1 GPa for aqueous fluids. Experiments with neon as the surrounding fluid (DBrfluid/melt = 0.38 ± 0.01 at 1.1 GPa) suggest that Br-release from a basalt into volatiles that have no bonding affinity with Br is weak. This should be the case for dry intra-plate volcanic eruptions. From the experimentally gained partition coefficients and from global Br concentration values in melt inclusions of arc magmas, we calculated an annual global Br flux of 23.5–72.9 × 109 g/y.
Occurrence and sources of 2,4,7,9-tetramethyl-5-decyne-4,7-diol (TMDD) in the aquatic environment
(2011)
The aim of the present study was to identify the sources of 2,4,7,9-tetramethyl-5-decyne-4,7-diol (TMDD) into the aquatic environment and to investigate its occurrence in rivers and wastewater treatment plants (WWTPs). Therefore, TMDD was analyzed in 441 wastewater samples from influents and effluents of 27 municipal WWTPs, in 6 sludge samples, in 52 wastewater samples from 3 sewage systems of municipal WWTPs, in 489 surface samples from 24 rivers, in 9 wastewater samples of 3 paper-recycling industries and in 65 groundwater samples. TMDD was also analyzed in household paper products, in 23 samples of toilet
papers, in 5 types of paper towels and in 12 types of paper tissues. The samples were collected between 2007 and 2011. The water samples were extracted with solid phase extraction (SPE) and the household paper samples with Soxhlet extraction. Gas chromatography-mass spectrometry (GC-MS) was used for quantification purposes. Between November 2007 and January 2008, TMDD was detected in the river Rhine at Worms with permanent high concentrations (up to 1330 ng/L). The results showed that TMDD is uniformly distributed across the river at Worms. An increase of the mean TMDD concentration from approximately 500 ng/L to 1000 ng/L was registered in January 2008. Due to the minor fluctuations of the TMDD concentration during the sampling period it is expected that the input of TMDD into the river is continuous. Therefore, TMDD might rather originate from effluents of municipal WWTPs than from temporal sources. The mean TMDD load based on the analysis of 147 water samples collected in the River Rhine was 62.8 kg/d which is equivalent to 23 t/a suggesting that TMDD must be used and/or produced in high quantities in order to be found in those high concentrations. To determine if TMDD is discharged by effluents of municipal WWTPs into the rivers, 24 hours influent and effluent samples of four municipal WWTPs in the Frankfurt/Rhine-Main metropolitan region were collected during November 2008 and February 2010 and analyzed for TMDD. The TMDD influent concentrations varied between 134 ng/L and 5846 ng/L and the effluent concentrations between <LOQ (limit of quantitation) and 3539 ng/L. The TMDD elimination rates in the four WWTPs varied between 33% and 68%. The results showed that effluents of municipal WWTPs are an important source of TMDD in the aquatic environment because TMDD is not completely removed from the sewage during the wastewater treatment. Weekly and daily variations of the TMDD concentration in the influents of two municipal WWTPs indicated that both private households and indirect industrial dischargers contribute to the introduction of TMDD into the municipal sewage systems. A more detailed study of the TMDD elimination rate in the different wastewater treatment stages was carried out in the WWTP Niederrad/Griesheim in Frankfurt am Main. The results showed that the removal of TMDD is mainly carried out during the aerobic biological treatments, where the elimination rate was 46%. In contrast, during the anoxic treatment the removal efficiency was only 1.4% and during the mechanical treatment the elimination rate was 19%. To determine the sources of TMDD in the sewage, household paper products (paper tissues, toilet papers and paper towels) were analyzed for TMDD using Soxhlet extraction. TMDD was detected in 83% of the samples (n=40). The highest mean TMDD concentrations were found in recycled toilet paper (0.20 μg/g) and in paper towels (0.11 μg/g). In paper tissues and non-recycled toilet paper the mean TMDD concentrations were lower 0.080 μg/g and 0.025 μg/g respectively. According to these results the high TMDD influent concentrations found previously in municipal WWTPs (mean 1.20 μg/L) cannot be explained due to migration of TMDD from the household paper products into the sewage. Thus indirect industrial dischargers are the cause of the high influent TMDD concentrations. Effluents of municipal WWTPs with different indirect industrial dischargers (textile-, metal processing-, food processing-, electroplating-, paper-recycling- and printing ink factories) were analyzed. The highest mean TMDD concentrations were found in the effluents of municipal WWTPs that have paper-recycling (71.3 μg/L) and printing ink factories (138 μg/L) as indirect industrial dischargers. These results were confirmed by analyzing process wastewater of three paper-recycling factories located in Germany. High TMDD concentrations were detected and fluctuated between 1.83 μg/L and 113 μg/L. TMDD was also analyzed in the wastewater of a non-recycling-paper factory but its concentration was much lower (0.066 μg/L) indicating that TMDD is introduced into the processing water during the papermaking process due to the use of waste paper. Analyses of wastewater samples from different parts of the sewage pipes of a municipal WWTP in Hesse, which receives the wastewater from a printing ink factory, were carried out. The TMDD concentration in the wastewater sample from the sewage pipe of the printing ink factory was much higher (3,300 μg/L) than the TMDD concentration detected in the other wastewater samples from the sewage system (0.030 μg/L – 0.89 g/L). These results confirm the printing ink production as one of the principal sources of TMDD in the sewage. Analysis of surface water samples of the River Modau downstream from the effluent of the WWTP Nieder-Ramstadt showed TMDD concentrations of up to 28.0 μg/L. These high TMDD concentrations might be caused by the indirect wastewater discharges of a paint factory connected to the municipal sewage system. These results indicate that TMDD is introduced into the municipal WWTPs principally by indirect industrial dischargers and they are mainly paint and printing ink factories. The paper-recycling factories also represent an important source of TMDD in municipal WWTPs but indirectly. According to statements given by the representatives of two paper recycling factories neither TMDD or any other TMDD containing product is used or added during the papermaking process. Therefore, TMDD is washed out from the printing inks of the coloured waste paper and concentrated in the process wastewater in the closed water circuits of paper-recycling factories reaching rivers and municipal WWTPs. The occurrence and distribution of TMDD in surface waters in Germany was also studied. The results showed that TMDD is widely distributed across different rivers systems in the federal states of Hesse, North-Rhine-Westphalia, Bavaria, Baden-Wuerttemberg and Rhineland-Palatinate. In Hesse, TMDD was detected in the some of main rivers with mean concentrations of 812 ng/L (Schwarzbach, Hessian Ried), 374 ng/L (Kinzig), 393 ng/L (Main, at Frankfurt), 539 ng/L (Werra), 326 ng/L (Fulda), 151 ng/L (Emsbach) and 161 ng/L (Nidda). In small rivers (creeks) the mean TMDD concentrations varied between <LOQ (Diemel, Urselbach) and 1890 ng/L (Darmbach). The results showed that the TMDD concentrations in creeks are highly influenced by both effluents of WWTPs and by the distance between the sampling point and the nearest WWTP. Surface samples from sampling locations downstream from WWTPs dischargers showed higher TMDD concentrations (mean 518 ng/L) than sampling locations upstream from WWTPs dischargers (mean 35.1 ng/L). The behavior of TMDD during bank filtration was investigated at two locations, at a water utility company at the Lower River Rhine (urban area) and at the Oderbruch polder (rural area). The results indicated that TMDD is removed from the surface water by bank filtration at both sampling locations. The removal process is probably carried out in the first meters of the aquifer (hyporheic zone) by biodegradation processes, since TMDD does not tend to be absorbed by sediments and it was not found in the groundwater of monitoring wells. In groundwater samples from the Hessian Ried (n=23) TMDD was found only in five samples and the highest TMDD concentration was 135 ng/L. According to these results, TMDD does not represent a concern for drinking water in Germany, since it does not reach the groundwater with high concentrations and it has a low toxicity potential. The input of TMDD into the North Sea was estimated to be 60.7 t/a by considering the mean transported loads of TMDD by the River Rhine at Wesel (58.3 t/a) and Meuse in the Netherlands (2.40 t/a). The estimated discharge of TMDD by German municipal WWTPs (8.19 t/a) and paper-recycling factories (9.24 t/a) into rivers seems to be too low considering that the mean TMDD load in the River Rhine downstream from Wesel is 58.3 t/a. However, due to the high density of population and industries at the Lower Rhine it is expected that more relevant sources of TMDD are located along the Rhine River increasing the transported load. According to the results of this PhD project TMDD is a non-ionic surfactant contained in products, which are applied on surfaces (printing inks and paints) and has the potential to reach the aquatic environment. Therefore, TMDD should fulfill the requirement of a biodegradability of 80% established by the “Law on the Environmental Impact of Detergents and Cleaning Products” in Germany. However, due to the partial elimination rates of TMDD obtained in municipal WWTPs (between 33% and 68%) and to the absence of information about the execution of the biodegradation test on TMDD, it is unknown if TMDD is in accordance with this law. Otherwise, its use as surfactant in such products is questionable.
The East African Rift System (EARS) was initiated in the Eocene epoch between 50 and 21 Ma probably due to the influence of mantle plumes that caused volcanism, flood basalts and rifting extensions in Ethiopa and the Afar region. As a result of magmatic intrusions and adiabatic decompression melting within the lithosphere caused by the impact of the Kenya plume, there was a southward propagation of the EARS of about 30 – 15 Ma from Ethiopia to Kenya, which coincide with the occurrence of volcanism. The EARS developed towards the south along the margins of the Tanzania Craton between 15 and 8 Ma. Previous findings of low-velocity anomalies within the upper mantle and the mantle transition zone indicate an upwelling of hot mantle material in the vicinity of the Afar region and the East African Rift. This study includes the analysis of P- and S-receiver functions in order to determine further impacts on the lithosphere from below. The aim was to determine the topographic undulations of further boundary layers and to identify their variability owing to the rifting processes and the formation of the EARS. The study area included the Tanzania Craton and the surrounding rift branches of the East African Rift System.
The region of the Rwenzori Mountains can be analysed in detail because of the large dataset of the RiftLink project. The use of the P-receiver function technique and the H-K stacking method enabled to determine different vP /vS ratios depending on the tectonic setting in the Rwenzori region: Rift shoulders (vP /vS =1.74), Albert Rift segment (vP /vS =1.80), Edward Rift segment (vP /vS =1.87) and Rwenzori Mountains (vP /vS =1.86). To determine the topography of the Moho, it is necessary to take into account the thickness of the sedimentary layer, the surface topography, the azimuthal variations in crustal thickness and the impact of local anomalies. After correcting these effects on the Moho depths, significant variations in Moho topography could be determined. The Moho depths range from 29 to 39 km beneath the rift shoulders of the Albertine Rift. Within the rift valley, the crustal thickness varies between 25 – 31 km in the Edward Rift segment and 22 – 30 km in the Albert Rift segment. An averaged crustal thickness of about 26 km within the rift valley indicates the lack of the crustal root beneath the Rwenzoris. Similar variations in crustal thickness were determined by using an automatic procedure for analysing S-receiver functions that was developed in this study.
The S-receiver functions are created by applying a rotation criterion in order to rotate the Z, N and E components into the L, Q and T components. It is necessary to perform trial rotations using different incident and azimuth angles to determine the correct rotation angles. The latter are identified by the use of the rotation criterion, including the amplitude ratio of the converted Moho signal to the direct S/SKS-wave signal. The L component is rotated correctly in the direction of the incident shear wave in the case of the maximum amplitude ratio. After analysing the frequency content of the receiver functions in order to sort out harmonic and long-periodic traces, the individual Moho signals are checked for consistency in order to remove atypic signals. To increase the signal-to-noise ratios on the traces, the S-receiver functions are stacked. For this purpose, the signals of the direct shear waves must originate from similar epicenters. On the basis of similar ray paths, the receiver functions show comparable waveforms and converted signals. To perform the stacking procedure, it is necessary to merge the datasets of the adjacent stations in order to obtain a sufficient number of receiver functions. This analysis is based on the assumption that the incident seismic waves arriving at the adjacent stations penetrate to some extent the same underground structures in the case of similar wave propagation paths. This approach accounts for the fact that the converted signals do not result exclusively from the piercing points at the boundary layers. Further signals originate from the conversions at the boundary layer within the Fresnel Zone. The piercing points are derived from the significant signals in the receiver functions. Depending on the order of arrival of the converted phases on the traces, the signals are attributed to the theoretical discontinuities DIS1, DIS2, DIS3 and DIS4. However, partly due to the low signal-to-noise ratios on the traces, it is difficult to identify the real conversions on the traces and to ensure that the converted signals are attributed to the correct boundary layers. For this reason, it is necessary to check the consistency of the conversion depths among each other. In the case of inconsistent conversion depths, the corresponding signals are either adjusted to another seismic boundary layer or removed from the dataset. To verify the functionality of the automatic procedure and to determine the resolvability with respect to two boundary layers, several models are tested including horizontal and dipping discontinuities. To resolve distinct discontinuities, their depths must differ by at least 60 km, otherwise, due to similar depth ranges of the different boundary layers, the converted signals cannot be separated from each other. As a consequence, the converted signals that originate from different discontinuities are attributed to a single one. Further tests including break-off edges of seismic discontinuities are performed to check the attributions of the converted signals to the discontinuities. Owing to the varying number of boundary layers, the converted signals cannot be attributed to the discontinuities according to the order of their arrivals on the traces. It is necessary to correct their attributions to the seismic discontinuities in order to resolve the boundary layers.
The crust-mantle boundary and further discontinuities within the lithospheric mantle are investigated by applying this automatic procedure. Depending on the tectonic setting, the conversion depths of the Moho range from about 30 – 45 km beneath the western rift shoulder to 20 – 35 km within the rift valley up to 30 – 40 km beneath the eastern rift shoulder. The long wavelengths of the shear waves hamper the correct identification of the converted phases in the S-receiver functions. With respect to the relative differences in conversion depth, the topographic undulations of the crust-mantle boundary are consistent with the Moho depths derived from P-receiver functions. In contrast to the Rwenzori region, it is difficult to resolve completely the trend of the Moho in the remaining area of the East African Rift due to the small dataset provided by IRIS. The results exibit an increase in crustal thickness to up to 45 km in the region of the Cenozoic volcanics such as Virunga, Kivu, Rungwe and Kenya. The greatest Moho depths of more than 50 km are located near Mount Kilimanjaro. In addition to the Moho, the analysis of the S-receiver functions revealed two further boundary layers at depths of 60 – 140 km and 110 – 260 km, which are associated with a mid-lithospheric discontinuity and the lithosphere-asthenosphere boundary, respectively. The shallowest conversion depths of the LAB are focussed to small-scale regions within the rift branches, namely the northern Albertine Rift, the Chyulu Hills and the Mozambique Belt, which are located around the Tanzania Craton. The larger thickness of the lithosphere beneath the cratonic terrain indicates that the Tanzania Craton is not significantly eroded. However, there are indications that the lithosphere beneath the craton and the rift branches is penetrated by ascending asthenospheric melts to depths of up to 140 and 60 km, respectively. The top of the ascending melts is associated with the occurrence of the mid-lithospheric discontinuity. The shallowest conversion depths of this boundary layer (60 – 90 km) are related to the rifted areas of the EARS and the Cenozoic volcanic provinces, which are located along the Albertine Rift, the Kenya Rift and the Rukwa-Malawi rift zones. The deepest conversion depths of up to 140 km are related to the Rwenzori Belt, the Ugandan Basement Complex and the interior of the Tanzania Craton.
The Early Permian coal is of great value in the Tengxian Coalfield, Shandon Province, Eastern China. This work deals with the new data focusing on mineralogical characteristics in the Early Permian Shanxi Formation No. 3 coal from the Jinyuan Mine. The Jinyuan coal is a low ash and highly volatile A bituminous coal. Minerals in the No. 3 coal mainly comprise of kaolinite, ankerite, illite, calcite, siderite, and quartz, with varying compositions of trace amounts of pyrite, jarosite, bassanite, anatase, and rutile. According to mineral assemblage in the coal plies, three Types (A to C) can be identified in the No. 3 coal. The dominant minerals in Type A are poorly-ordered kaolinite, illite, quartz, pyrite, and jarosite. Type B is mainly composed of well-ordered kaolinite, illite, siderite, ankerite, and calcite. Type C, with just one sample (JY-3-7c), which contains high proportions of calcite (54%) and ankerite (34%). Terrigenous minerals are elevated in coal plies that typically have relatively high contents of ash yield. The formation of syngenetic pyrite was generally due to seawater, while the sulphate minerals (jarosite and coquimbite) were derived from the oxidation of pyrite. Epigenetic vein-like or fracture-fillings carbonate minerals (ankerite, calcite, and siderite), kaolinite, and pyrite, as well as authigenic quartz were derived from the influx of hydrothermal fluids during different periods, from the authigenic to epigenetic. The paragonite in the coal may have been formed by the precipitated from Na-rich hydrothermal fluids. No effects of magmatic intrusion on mineralogy were investigated in this research.
Determining how the elevation of the Northern Andes has evolved over time is of paramount importance for understanding the response of the Northern Andes to deformational and geodynamic processes and its role as an orographic barrier for atmospheric vapor transport over geologic time. However, a fundamental requirement when using stable isotope data for paleoaltimetry reconstructions is knowledge about the present-day changes of δ18O and δD with elevation (isotopic lapse rate). This study defines the present-day river isotopic lapse rate near the Equator (∼3°S) based on analysis of δ18O and δD of surface waters collected from streams across the Western Cordillera and the Inter-Andean depression in Southern Ecuador. The results for the two domains show a decrease of δ18O with elevation which fits a linear regression with a slope of −0.18‰/100 m (R2 = 0.73, n = 83). However, we establish a present-day lapse rate of −0.15‰/100 m for δ18O (R2 = 0.88, n = 19) and -1.4‰/100 m for δD (R2 = 0.93, n = 19) from water samples collected along the west facing slopes of the Western Ecuadorian Cordillera which is mainly subject to moisture transport from the Pacific. We argue that this empirical relationship, consistent with those obtained in different tropical areas of the world, can inform stable isotope paleoaltimetry reconstructions in tropical latitudes.
Strong perturbations of the Arctic stratosphere during the winter 2002/2003 by planetary waves led to enhanced stretching and folding of the vortex. On two occasions the vortex in the lower stratosphere split into two secondary vortices that re-merged after some days. As a result of these strong disturbances the role of transport in and out of the vortex was stronger than usual. An advection and mixing simulation with the Chemical Lagrangian Model of the Stratosphere (CLaMS) utilising a suite of inert tracers tagging the original position of the air masses has been carried out. The results show a variety of synoptic and small scale features in the vicinity of the vortex boundary, especially long filaments peeling off the vortex edge and being slowly mixed into the mid latitude environment. The vortex folding events, followed by re-merging of different parts of the vortex led to strong filamentation of the vortex interior. During January, February, and March 2003 flights of the Russian high-altitude aircraft Geophysica were performed in order to probe the vortex, filaments and in one case the merging zone between the secondary vortices. Comparisons between CLaMS results and observations obtained from the Geophysica flights show in general good agreement. Several areas affected by both transport and strong mixing could be identified, allowing explanation of many of the structures observed during the flights. Furthermore, the CLaMS simulations allow for a quantification of the air mass exchange between mid latitudes and the vortex interior. The simulation suggests that after the formation of the vortex was completed, its interior remaind relatively undisturbed. Only during the two re-merging events were substantial amounts of extra-vortex air transported into the polar vortex. When in March the vortex starts weakening additional influence from lower latitudes becomes apparent in the model results. In the lower stratosphere export of vortex air leads only to a fraction of about 5% polar air in mid latitudes by the end of March. An upper limit for the contribution of ozone depleted vortex air on mid-latitude ozone loss is derived, indicating that the maximum final impact of dilution is on the order of 50%.
Strong perturbations of the Arctic stratosphere during the winter 2002/2003 by planetary waves led to enhanced stretching and folding of the vortex. On two occasions the vortex in the lower stratosphere split into two secondary vortices that re-merged after some days. As a result of these strong disturbances the role of transport in and out of the vortex was stronger than usual. An advection and mixing simulation with the Chemical Lagrangian Model of the Stratosphere (CLaMS) utilising a suite of inert tracers tagging the original position of the air masses has been carried out. The results show a variety of synoptic and small scale features in the vicinity of the vortex boundary, especially long filaments peeling off the vortex edge and being slowly mixed into the mid latitude environment. The vortex folding events, followed by re-merging of different parts of the vortex led to strong filamentation of the vortex interior. During January, February, and March 2003 flights of the Russian high-altitude aircraft Geophysica were performed in order to probe the vortex, filaments and in one case the merging zone between the secondary vortices. Comparisons between CLaMS results and observations obtained from the Geophysica flights show in general good agreement.
The International Halocarbons in Air Comparison Experiment (IHALACE) was conducted to document relationships between calibration scales among various laboratories that measure atmospheric greenhouse and ozone depleting gases. Six stainless steel cylinders containing natural and modified natural air samples were circulated among 19 laboratories. Results from this experiment reveal relatively good agreement among commonly used calibration scales for a number of trace gases present in the unpolluted atmosphere at pmol mol−1 (parts per trillion) levels, such as chlorofluorocarbons (CFCs), hydrochlorofluorocarbons (HCFCs), and hydrofluorocarbons (HFCs). Some scale relationships were found to be consistent with those derived from bi-lateral experiments or from analysis of atmospheric data, while others revealed discrepancies. The transfer of calibration scales among laboratories was found to be problematic in many cases, meaning that measurements tied to a common scale may not, in fact, be compatible. These results reveal substantial improvements in calibration over previous comparisons. However there is room for improvement in communication and coordination of calibration activities with respect to the measurement of halogenated and related trace gases.
Die vorliegende Arbeit liefert einen Beitrag zum Verständnis der Rolle des RO x bei der troposphärischen Ozonbildung. Troposphärisches Ozon (O 3 ) spielt eine wichtige Rolle bei der Selbstreinigung der Atmosphäre. Andererseits führen erhöhte Ozonkonzentrationen zu gesundheitlichen Beeinträchtigungen beim Menschen und Schäden an Pflanzen und Umwelt. Die Anwesenheit von flüchtigen organischen Verbindungen (VOCs) führt zur Bildung von Peroxyradikalen (RO x ), die das normale photochemische Gleichgewicht zwischen Ozon und Stickoxiden zu Gunsten erhöhter OzonKonzentrationen verschieben. Im Rahmen der Arbeit wurde ein chemischer Verstärker zur Messung der GesamtPeroxyradikalkonzentration gebaut. RO x reagiert im Einlass des Gerätes mit hinzugefügtem NO und CO in einer Kettenreaktion und bildet dabei NO 2 . Dieses wird mit einem Luminoldetektor nachgewiesen. Der Detektor wird alle 2 Stunden kalibriert. Die Kettenlänge wird durch eine Kalibrierung des Gerätes mit HO 2 Radikalen bestimmt, die durch die Photolyse von H 2 O gebildet werden. Der Verstärkungsfaktor wurde in Bezug auf eine Querempfindlichkeit gegen Wasserdampf korrigiert. Die Messgenauigkeit ist etwa 70% bei 60% relativer Feuchte. Messungen am Taunus Observatorium auf dem Kleinen Feldberg in den Sommermonaten der beiden Jahre 1998 und 1999 werden vorgestellt. Die Ozon und RO x Konzentrationen sind gut miteinander korreliert. Allerdings ist die Tagestemperatur die für die Ozon und RO x Konzentrationen bei weitem wichtigste Einflussgröße und ist daher der beste Parameter zur statistischen Beschreibung von photochemischen Vorgängen. Auf der Grundlage der Messungen am Kleinen Feldberg wurde ein einfaches statistisches Modell zur Vorhersage des Ozonmaximums erstellt. Mit den Parametern Temperatur und Ozonkonzentration am Vortag konnte das statistische Modell bereits 80% der Variation der Ozonkonzentration erklären. Durch die Berücksichtigung der RO x Messungen am Vormittag konnte lediglich eine Verbesserung der erklärten Varianz um 0.5% erzielt werden. Um einen Hinweis auf den Einfluss anthropogener Emissionen zu bekommen, wurde der Wochengang von Ozon, RO x und NO x ebenfalls untersucht. Die Zunahme des Ozonmischungsverhältnisses am Wochenende bei gleichzeitigem Rückgang des Mischungsverhältnisses der Stickoxide wird damit erklärt, dass am Kleinen Feldberg eine VOClimitierte Situation vorgefunden wurde. Die Ozonbildungsrate auf Basis der Reaktion zwischen RO x und NO wurde für Tage mit einem Maximum der Globalstrahlung über 600 W m tdatensatz niedrig (r = 0,46). Die beobachtete Änderung des Ozonmischungsverhältnisses wurde mit dem berechneten mittleren Tagesgang der Ozonbildungsrate verglichen. Die Ozonbildungsrate lag um die Mittagszeit bei etwa 5 ppbv h Verlustprozesse zu erklären. Am Abend werden etwa 2 ppbv O 3 pro Stunde abgebaut. Im Rahmen einer Messkampagne im Juni/Juli 2000 am Meteorologischen Observatorium Hohenpeißenberg fanden Messungen der Konzentrationen von RO x , OH, einer Reihe von VOCs, und anderen relevanten Spurengasen statt. Die Messdaten werden mit Hilfe eines auf der Annahme des lokalen photostationären Gleichgewichts der Radikale basierenden Modells interpretiert. Die Modellergebnisse stimmten sehr gut mit den Messungen überein. Die Überschätzung der Konzentration an 2 Tagen wurde durch den Einfluss sauerstoffhaltiger VOCs erklärt. Das '' Recycling" der HO 2 Radikale (die Reaktion zwischen HO 2 und NO) ist die wichtigste Quelle für OH und die wichtigste Senke für RO x . Durch die erhöhte NOKonzentration am Vormittag wird HO 2 sehr schnell in OH umgewandelt, das wiederum für die VOCOxidation und RO x Bildung verantwortlich ist. Die wichtigste OHSenke und RO x Quelle ist die Oxidation von Isopren und den Terpenen. Um die Rolle der photochemischen Ozonbildung auf regionaler Skala zu untersuchen, wurden Ozonmessungen aus ganz Deutschland auf unterschiedlichen zeitlichen und räumlichen Skalen statistisch untersucht. Die Netto Änderungsrate der Ozonkonzentration war tagsüber an 3 nahe zusammenliegenden Stationen sehr ähnlich. Die OzonMessdaten von 277 deutschen Messstationen wurden mit den an einer Waldmessstelle nahe Königstein gemessenen Ozonwerten korreliert. Die Ozonmessungen in Königstein erklären 50% der Varianz der sommerlichen Ozonmessungen zwischen 11:00 und 16:00 MEZ an Stationen, die in einem Umkreis von etwa 250 km von Königstein liegen. Auf das ganze Jahr bezogen, liegt diese ''charakteristische Entfernung" bei etwa 350 km. Diese Ergebnisse deuten darauf hin, dass die Prozesse, die einen wichtigen Einfluss auf die Ozonkonzentration ausüben, auf regionalen Skalen von einigen hundert Kilometern aktiv sind. Zusammenfassend lässt sich sagen, dass die gemessenen RO x Konzentrationen mit den aufgrund der Oxidation der VOCs durch OH berechneten Konzentrationen konsistent sind. Obwohl die RO x Konzentationen für die chemische Modellierung von Bedeutung sind, tragen RO x Messungen nur wenig zu einer Verbesserung der Qualität von kurzfristigen statistischen Ozonprognosen bei. Keywords: Ozone, Troposphere, Peroxy Radicals, Free Radicals, Photochemistry, Chemical Amplifier
Measurements of OH, the sum of peroxy radicals (ROx), non-methane hydrocarbons (NMHCs) and various other trace gases were made at the Meteorological Observatory Hohenpeissenberg in June 2000. The data from an intensive measurement period characterised by high solar insolation (18-21 June) are analysed. The maximum midday OH concentration ranged between 4.5 x 106 molecules cm-3 and 7.4 x 106 molecules cm-3. The maximum total ROx mixing ratio increased from about 55 pptv on 18 June to nearly 70 pptv on 20 and 21 June. A total of 64 NMHCs, including isoprene and monoterpenes, were measured every 1 to 6 hours. The oxidation rate of the NMHCs by OH was calculated and reached a total of over 14 x 106 molecules cm-3 s-1 on two days. A simple photostationary state balance model was used to simulate the ambient OH and ROx concentrations with the measured data as input. The model was able to reproduce the main features of the diurnal profiles of both OH and ROx. The model results proved to be most sensitive to assumptions about the mixing ratio of formaldehyde (HCHO), which was included as a proxy for carbonyl compounds, and about the partitioning between HO2 and RO2. The measured OH concentration and ROx mixing ratios were reproduced well by assuming the presence of 3 ppbv HCHO and a ratio HO2/RO2 between 1:1 and 1:2. The most important source of OH, and conversely the greatest sink for ROx, was the recycling of HO2 radicals to OH. This reaction was responsible for the recycling of more than 45 x 106 molecules cm-3 s-1 on two days. The most important sink for OH, and the largest source of ROx, was the oxidation of NMHCs, in particular, of isoprene and the monoterpenes.
Measurements of OH, total peroxy radicals, non-methane hydrocarbons (NMHCs) and various other trace gases were made at the Meteorological Observatory Hohenpeissenberg in June 2000. The data from an intensive measurement period characterised by high solar insolation (18-21 June) are analysed. The maximum midday OH concentration ranged between 4.5x106 molecules cm-3 and 7.4x106 molecules cm-3. The maximum total ROx (ROx =OH+RO+HO2+RO2) mixing ratio increased from about 55 pptv on 18 June to nearly 70 pptv on 20 and 21 June. A total of 64 NMHCs, including isoprene and monoterpenes, were measured every 1 to 6 hours. The oxidation rate of the NMHCs by OH was calculated and reached a total of over 14x106 molecules cm-3 s-1 on two days. A simple photostationary state balance model was used to simulate the ambient OH and peroxy radical concentrations with the measured data as input. This approach was able to reproduce the main features of the diurnal profiles of both OH and peroxy radicals. The balance equations were used to test the effect of the assumptions made in this model. The results proved to be most sensitive to assumptions about the impact of unmeasured volatile organic compounds (VOC), e.g. formaldehyde (HCHO), and about the partitioning between HO2 and RO2. The measured OH concentration and peroxy radical mixing ratios were reproduced well by assuming the presence of 3 ppbv HCHO as a proxy for oxygenated hydrocarbons, and a HO2/ RO2 ratio between 1:1 and 1:2. The most important source of OH, and conversely the greatest sink for peroxy radicals, was the recycling of HO2 radicals to OH. This reaction was responsible for the recycling of more than 45x106 molecules cm-3 s-1 on two days. The most important sink for OH, and the largest source of peroxy radicals, was the oxidation of NMHCs, in particular, of isoprene and the monoterpenes.
The biomarker record in two different lakes in central Europe, Lake Albano and Lake Constance, is used to reflect environmental changes and lake system response during the Late Glacial and Holocene. Extractable organic compounds in lake sediments, which can be assigned to their biological source (biomarkers) function as fingerprints of past aquatic or land plant organisms. Using gas chromatography coupled with mass spectrometry, 21 different biomarkers (predominantly steroids and triterpenoids) as well as a variety of n-alkanes, nalkanols, and n-alkanoic acids could be identified in the sediment records of Lake Albano and Lake Constance. In the Holocene sediments of Lake Albano, the distribution of biomarkers such as dinosterol (dinoflagellates), isoarborinol, and diplopterol (aquatic organisms) indicate three biomarker zones: The period between 0-3,800 years BP (zone 3) is characterized by high concentrations of these biomarkers and others such as tetrahymanol and diploptene. Conversely, zone 2 (3,800-6,500 years BP) shows very low concentrations of all autochthonous biomarkers. In zone 1 (6,500–11,480 years BP), dinosterol, isoarborinol, and diplopterol range on a relatively high level, whereas diploptene and tetrahymanol display comparatively low concentrations. The results suggest at least two distinct changes in the predominance of primary producers during the Holocene, which are related to changes in the lake system such as lake mixing and water column stratification. This interpretation is consistent with previous investigations of Lake Albano sediments including pigment and hydrogen index data (Ariztegui et al., 1996b; Guilizzoni et al., 2002). Allochthonous biomarkers such as long-chain n-alkanes, amyrenones and friedelin indicate a development from forest to a more open landscape from 6,000 and 5.000 years BP, respectively. After a period of high concentrations during the first half of the Holocene, all biomarkers derived from deciduous trees exhibit relatively low values until around 1,000 years BP. Again, this is consistent with results from previous pollen investigations (Ariztegui et al., 2000). The sediment core from Upper Lake Constance comprises the Late Glacial and Holocene. It was analysed for biomarkers and inorganic tracers in order to compare the biomarker results with other proxy data from the same core. Magnetic susceptibility (MS) was measured to get a high-resolution stratigraphic framework of the core and to obtain further information about changes of the proportions of allochthonous and autochthonous input. Enhanced concentrations and accumulation rates of dinosterol (biomarker for dinoflagellates) and biogenic calcite give evidence of increasing lake productivity at the beginning of the Holocene followed by a decrease in bioproductivity after around 7,000 years BP. Younger Dryas sediments are characterized by low amounts of both dinosterol and biogenic calcite indicating a low productivity. The comparison of the concentrations and accumulation rates of b-sitosterol and stigmastanol with parameters reflecting lake productivity suggests that both steroids in Lake Constance sediments are mainly derived from terrigenous sources. Biomarkers as well as concentrations and accumulation rates of allochthonous inorganic compounds such as titanium, magnesium and strontium indicate a slightly enhanced allochthonous input after 8,500 years BP. Significant increase of erosive matter input from enhanced soil erosion is not observed before 4,000 years BP. This can be attributed to the combined effects of precipitation increase as a result of climatic deterioration and anthropogenic deforestation which is consistent with observations from other lakes in Central Europe. The MS record of Lake Constance confirms these results by tracing the climatically induced shifts of more intense bioproduction (low MS caused by increased calcite deposition) during the ‘climatic optimum’. This is followed by increasing input of terrigenous sediment compounds during colder and wetter periods which lead to higher MS values in the lake sediments. The occurrence of tetrahymanol in Lake Constance sediments questions the unambiguous use of tetrahymanol as an indicator for water column stratification. Anaerobic organic macroaggregates within the oxygenated, photic zone of the water column have to be considered as a possible living space for anaerobic microorganisms containing tetrahymanol. The direct comparison of two very different lakes Albano and Constance with respect to biomarkers indicating climate or environmental change provides a contribution to the recent biomarker research for a better understanding of biomarkers in lacustrine sediments.
An eclogite barometer has profound importance in the study of upper mantle processes and potential application to diamond prospecting. Studies on the partitioning of Li between clinopyroxene (cpx) and garnet (grt) in natural samples have shown that this particular element is very sensitive to changes in pressure and could be calibrated as the barometer demanded for bimineralic eclogites. Experiments were performed from 4 to 13 GPa and 1100 to 1400°C in the CMAS (CaO, MgO, Al2O3, SiO2) system with Li added as Li3PO4 to quantify this pressure dependence into a barometer expressed in the following equation: P= (0.00255*T-lnKd)/0.2351 where P is in GPa, T is in °C and Kd is defined as the partition coefficient of Li (in ppm) between clinopyroxene and garnet. The experimental pressures are reproduced to ± 0.38 GPa (1σ) by this equation. This barometer is strictly applicable only to CMAS. Experiments at 1300°C, 8-12 GPa showed that Henry’s Law is fulfilled for Li partitioning between cpx and grt in the concentration range of approximately 0.01 – 1 wt% Li. Direct application of the equation to experiments in natural systems performed at 1300°C from 4 GPa to 13 GPa consistently overestimates pressures by approximately 2 GPa. Our previous experiments in the system CaO-MgO-Al2O3-SiO2 + Li3PO4 showed that the partitioning of Li between garnet and clinopyroxene is pressure dependent in eclogitic bulk compositions. This supports experimentally the hypothesis of Seitz et al. (2003), based on the analysis of Li in eclogitic xenoliths and inclusions in diamond, that the partitioning of this particular element between clinopyroxene and garnet is very sensitive to changes in pressure and could be calibrated as a barometer for bimineralic eclogites. In order to calibrate this pressure dependence into a barometer, experiments were performed in natural systems using starting materials sourced from a well preserved eclogitic xenolith from the Roberts Victor kimberlite pipe (South Africa) to extrapolate our findings in CMAS to natural systems. Sixteen multianvil experiments were performed from 4-13 GPa and 1100-1500°C. Our findings reinforced the general trend we observed in the CMAS system, that KdLi cpx-grt for Li decreases with increasing P, and that at P ≥ 12 GPa, garnet is able to incorporate more Li than clinopyroxene. Multiple linear regression was applied to our experimental results to create the barometer: P = (0.000963*T – ln KdLi cpx-grt + 1.581) / 0.252 Where P is pressure in GPa, T is temperature in °C and KdLi cpx-grt is defined as the partitioning coefficient of Li obtained by dividing the concentration of Li in cpx by the concentration of Li in garnet. This barometer reproduces the experimental conditions to ± 0.2 GPa. It is applicable to eclogitic xenoliths, to garnet pyroxenites and to peridotitic and eclogitic inclusions in diamond. Application of the barometer to diamond bearing xenoliths results in pressures in the diamond stability field. Clinopyroxene is easily corrupted in xenoliths and also preferentially takes in Li during short lived metasomatic processes. Care must be taken therefore to analyse primary, unaltered clinopyroxene. Our preliminary application to natural samples shows that the barometer can be applied beyond the experimental range to pressures down to 3 GPa. Seventeen eclogitic xenoliths were chosen from a sample set of greater than 200 for their fresh microscopic and macroscopic appearance and were analyzed for Li content in coexisting garnet (grt) and clinopyroxene (cpx). These samples can be subdivided into two groups on the basis of Mg in cpx (cpfu: cations per formula unit, based on 6 oxygens): Group 1 with Mg > 0.75, and Group 2 with Mg < 0.75. Group 1 xenoliths show lower Li contents in both grt and cpx compared to Group 2. The Li barom ter calibrated in Hanrahan et al. (2009b)/Chapter 3 was applied to these samples as well as available literature data to obtain pressures of provenance - Group 2 xenoliths often provide pressures that appear unrealistic for eclogitic xenoliths. In light of observed crystal chemical relations in the natural samples, a new fitting procedure was applied to the experimental data presented in Chapter 3. This new fit appears to be more realistic than the previous fit, although a strong relationship with Mg# remains present, suggesting that Li-barometry is, at present, only applicable to Mg-rich eclogites. Inclusions in diamond, with the exception of eclogitic inclusions of coexisting majorite and cpx, often yield pressures that are inconsistent with the pressures required for diamond formation. Although an interesting observation when comparing all of the data is that inclusions in diamond have significantly higher average Li concentrations compared to xenoliths, which suggests that Li is highly present in the fluids from which diamonds form in the mantle, an observation which was previously made for the deep mantle as a result of high Li in ferropericlase inclusions in diamond (Seitz et al. 2003).
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.
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.
Oceanic islands only comprise a small amount of the Earth’s land area but harbour a disproportionate amount of global biodiversity. This vast diversity is not only reflected in the taxonomic uniqueness of island biota but also in the remarkable evolution of functional traits. Functional traits, i.e. measurable characteristics that strongly influence the fitness of species, determine how a species responds to its environment and can help to gain more insights into the biogeographical, ecological and evolutionary processes that have shaped island biodiversity. However, research in island biogeography has primarily focused on species richness, and knowledge of functional trait patterns on oceanic islands is scarce. Hence, in this dissertation, I have explored how trait-based approaches can increase our understanding of how biodiversity on oceanic islands assembles and how it is driven by the environment. The Canary Islands (Spain) are a particularly suitable model system to investigate patterns and drivers of biodiversity. The archipelago is characterised by a high variation in environmental heterogeneity and inhabits a unique and well-described native flora. Therefore, I have investigated five principal research questions using the flora (Spermatophytes) of the Canary Islands as a study object. First, I have analysed how climate and biogeography shape the assembly of the Canary Islands flora using a novel trait-based approach. Second, the question of whether rare climates link to functional trait distinctiveness in the native Canary Islands flora was addressed. Third, I have examined how intraspecific trait variation is represented in the native flora of oceanic islands focusing on the succulent scrub of La Palma (Canary Islands). Fourth, this dissertation investigated whether scientific floras can be reliable sources for trait data of plants native to oceanic islands. Finally, I have explored how climate change may impact the native Canary Islands flora by analysing possible climate change-induced shifts in plant species distribution and plant traits.
The results of my dissertation expand the understanding of the importance of biogeography and the environment in determining the functional composition of island floras. I have assessed that traits of endemic plant species did not expand the functional trait space of the Canary Islands but were packed with the ones of non-endemic species. This result hints at a trait convergence in endemic species, possibly driven by non-adaptive speciation processes. Moreover, I have evidenced that humidity is a critical driver of functional diversity in native plant assemblages and particularly leads to a high trait convergence in arid environments via environmental filtering. In contrast, alien species have expanded the Canary Islands flora’s functional trait space. I further have shown that in contrast to native species assemblages, alien species assemblages are characterised by an increasing functional diversity with increasing aridity. This contrasting pattern of functional diversity could pose a potential risk to the native flora of the Canary Islands as a low functional diversity is expected to reduce the resilience of species assemblages to the establishment of more functionally diverse alien plant species. However, in this dissertation, I also have revealed that endemic plant species on the Canary Islands show a high intraspecific variation in arid environments, possibly as an adaptation to environmental stress. Intraspecific variation could help endemic plant species have a competitive advantage over alien species and be more resilient to environmental changes. Furthermore, in this dissertation, I have shown that scientific floras and taxonomic monographs could be used to gain information on quantitative functional traits of plants native to oceanic islands. This finding is particularly relevant for advances in trait-based research, as coverage of trait data for oceanic island floras is extremely poor in global trait databases. Hence, for some of the studies included in this dissertation, trait data were retrieved from scientific floras and taxonomic monographs and used to answer novel scientific research questions. Thus, I have used trait data from the literature to analyse the effect of climate change on the range size of plants native to the Canary Islands. Identifying plant species of particular conservation concern is critical on oceanic islands as many island species have limited distributions and small population sizes, and their niche tracking is impeded by insularity. I have revealed that single-island endemic plants gain less and lose more climatically suitable areas than archipelago endemic and non-endemic native plants due to a climate change-induced decrease in precipitation until 2100...
Questions: Both species turnover and intraspecific trait variation can affect plant assemblage dynamics along environmental gradients. Here, we asked how community assemblage patterns in relation to species turnover and intraspecific variation differ between endemic and non-endemic species. We hypothesized that endemic species show lower intraspecific variation than non-endemic species because they tend to have high rates of in situ speciation, whereas non-endemic species are expected to have a larger gene pool and higher phenotypic plasticity.
Location: La Palma, Canary Islands.
Methods: We established 44 sampling sites along a directional gradient of precipitation, heat load, soil nitrogen, phosphorus and pH. Along this gradient, we estimated species abundances and measured three traits (plant height, leaf area and leaf thickness) on perennial endemic and non-endemic plant species. In total, we recorded traits for 1,223 plant individuals of 43 species. Subsequently, we calculated community-weighted mean traits to measure the relative contribution of species turnover, intraspecific variation and their covariation along the analysed gradient.
Results: The contribution of intraspecific variation to total variation was similar in endemic and non-endemic assemblages. For plant height, intraspecific variation explained roughly as much variation as species turnover. For leaf area and leaf thickness, intraspecific variation explained almost no variation. Species turnover effects mainly drove trait responses along the environmental gradient, but intraspecific variation was important for responses in leaf area to precipitation.
Conclusions: Despite their distinct evolutionary history, endemic and non-endemic plant assemblages show similar patterns in species turnover and intraspecific variation. Our results indicate that species turnover is the main component of trait variation in the underlying study system. However, intraspecific variation can increase individual species’ fitness in response to precipitation. Overall, our study challenges the theory that intraspecific trait variation is more important for the establishment of non-endemic species compared with endemic species.
Climatic and biogeographical drivers of functional diversity in the flora of the Canary Islands
(2022)
Aim: Functional traits can help us to elucidate biogeographical and ecological processes driving assemblage structure. We analysed the functional diversity of plant species of different evolutionary origins across an island archipelago, along environmental gradients and across geological age, to assess functional aspects of island biogeographical theory.
Location: Canary Islands, Spain.
Major taxa studied: Spermatophytes.
Time period: Present day.
Methods: We collected data for four traits (plant height, leaf length, flower length and fruit length) associated with resource acquisition, competitive ability, reproduction and dispersal ability of 893 endemic, non-endemic native and alien plant species (c. 43% of the Canary Island flora) from the literature. Linking these traits to species occurrences and composition across a 500 m × 500 m grid, we calculated functional diversity for endemic, non-endemic native and alien assemblages using multidimensional functional hypervolumes and related the resulting patterns to climatic (humidity) and island biogeographical (geographical isolation, topographic complexity and geological age) gradients.
Results: Trait space of endemic and non-endemic native species overlapped considerably, and alien species added novel trait combinations, expanding the overall functional space of the Canary Islands. We found that functional diversity of endemic plant assemblages was highest in geographically isolated and humid grid cells. Functional diversity of non-endemic native assemblages was highest in less isolated and humid grid cells. In contrast, functional diversity of alien assemblages was highest in arid ecosystems. Topographic complexity and geological age had only a subordinate effect on functional diversity across floristic groups.
Main conclusions: We found that endemic and non-endemic native island species possess similar traits, whereas alien species tend to expand functional space in ecosystems where they have been introduced. The spatial distribution of the functional diversity of floristic groups is very distinct across environmental gradients, indicating that species assemblages of different evolutionary origins thrive functionally in dissimilar habitats.
Sedimentary charcoal records are widely used to reconstruct regional changes in fire regimes through time in the geological past. Existing global compilations are not geographically comprehensive and do not provide consistent metadata for all sites. Furthermore, the age models provided for these records are not harmonised and many are based on older calibrations of the radiocarbon ages. These issues limit the use of existing compilations for research into past fire regimes. Here, we present an expanded database of charcoal records, accompanied by new age models based on recalibration of radiocarbon ages using IntCal20 and Bayesian age-modelling software. We document the structure and contents of the database, the construction of the age models, and the quality control measures applied. We also record the expansion of geographical coverage relative to previous charcoal compilations and the expansion of metadata that can be used to inform analyses. This first version of the Reading Palaeofire Database contains 1676 records (entities) from 1480 sites worldwide. The database (RPDv1b – Harrison et al., 2021) is available at https://doi.org/10.17864/1947.000345.
Deriving stratospheric age of air spectra using an idealized set of chemically active trace gases
(2019)
Analysis of stratospheric transport from an observational point of view is frequently realized by evaluation of the mean age of air values from long-lived trace gases. However, this provides more insight into general transport strength and less into its mechanism. Deriving complete transit time distributions (age spectra) is desirable, but their deduction from direct measurements is difficult. It is so far primarily based on model work. This paper introduces a modified version of an inverse method to infer age spectra from mixing ratios of short-lived trace gases and investigates its basic principle in an idealized model simulation. For a full description of transport seasonality the method includes an imposed seasonal cycle to gain multimodal spectra. An ECHAM/MESSy Atmospheric Chemistry (EMAC) model simulation is utilized for a general proof of concept of the method and features an idealized dataset of 40 radioactive trace gases with different chemical lifetimes as well as 40 chemically inert pulsed trace gases to calculate pulse age spectra. It is assessed whether the modified inverse method in combination with the seasonal cycle can provide matching age spectra when chemistry is well-known. Annual and seasonal mean inverse spectra are compared to pulse spectra including first and second moments as well as the ratio between them to assess the performance on these timescales. Results indicate that the modified inverse age spectra match the annual and seasonal pulse age spectra well on global scale beyond 1.5 years of mean age of air. The imposed seasonal cycle emerges as a reliable tool to include transport seasonality in the age spectra. Below 1.5 years of mean age of air, tropospheric influence intensifies and breaks the assumption of single entry through the tropical tropopause, leading to inaccurate spectra, in particular in the Northern Hemisphere. The imposed seasonal cycle wrongly prescribes seasonal entry in this lower region and does not lead to a better agreement between inverse and pulse age spectra without further improvement. Tests with a focus on future application to observational data imply that subsets of trace gases with 5 to 10 species are sufficient for deriving well-matching age spectra. These subsets can also compensate for an average uncertainty of up to ±20 % in the knowledge of chemical lifetime if a deviation of circa ±10 % in modal age and amplitude of the resulting spectra is tolerated.
Analysis of stratospheric transport from an observational point of view is frequently realized by evaluation of mean age of air values from long-lived trace gases. However, this provides more insight into general transport strength and less into its mechanism. Deriving complete transit time distributions (age spectra) is desirable, but their deduction from direct measurements is difficult and so far primarily achieved by assumptions about dynamics and spectra themselves. This paper introduces a modified version of an inverse method to infer age spectra from mixing ratios of short-lived trace gases. For a full description of transport seasonality the formulation includes an imposed seasonal cycle to gain multimodal spectra. The EMAC model simulation used for a proof of concept features an idealized dataset of 40 radioactive trace gases with different chemical lifetimes as well as 40 chemically inert pulsed trace gases to calculate pulse age spectra. Annual and seasonal mean inverse spectra are compared to pulse spectra including first and second moments as well as the ratio between them to assess the performance on these time scales. Results indicate that the modified inverse age spectra match the annual and seasonal pulse age spectra well on global scale beyond 1.5 years mean age of air. The imposed seasonal cycle emerges as a reliable tool to include transport seasonality in the age spectra. Below 1.5 years mean age of air, tropospheric influence intensifies and breaks the assumption of single entry through the tropical tropopause, leading to inaccurate spectra in particular in the northern hemisphere. The imposed seasonal cycle wrongly prescribes seasonal entry in this lower region and does not lead to a better agreement between inverse and pulse age spectra without further improvement. As the inverse method aims for future implementation on in situ observational data, possible critical factors for this purpose are delineated finally.
Many natural minerals exist in the form of a solid solution. The systematic changes in structural and physical properties of oxide solid solutions are of geological importance and allow for wide applications. In order to understand the composition-structureproperty relations, substitutional solid solutions of CuxZn2−xTiO4, ZnxMg1−xTi2O5 and CuxMg1−xTi2O5 have been synthesised by mechanochemical activation assisted solid state synthesis. Self-propagating high-temperature synthesis has been employed to achieve the interstitial solid solutions of Ti5Si3Zx (Z refers to the element boron or oxygen).
The changes in the crystal structure and physical properties due to the formation of solid solutions are investigated by employing X-ray diffraction, neutron diffraction, Raman spectroscopy, low-temperature heat capacity, thermal expansion, scanning electron microscopy, UV-vis spectroscopy, plane-wave ultrasound spectroscopy and density functional theory calculations.
This paper compares measurements of gaseous and particulate emissions from a wide range of biomass-burning plumes intercepted by the NASA DC-8 research aircraft during the three phases of the ARCTAS-2008 experiment: ARCTAS-A, based out of Fairbanks, Alaska USA (3 April to 19 April 2008); ARCTAS-B based out of Cold Lake, Alberta, Canada (29 June to 13 July 2008); and ARCTAS-CARB, based out of Palmdale, California, USA (18 June to 24 June 2008). Extensive investigations of boreal fire plume evolution were undertaken during ARCTAS-B, where four distinct fire plumes that were intercepted by the aircraft over a range of down-wind distances (0.1 to 16 hr transport times) were studied in detail. Based on these analyses, there was no evidence for ozone production and a box model simulation of the data confirmed that net ozone production was slow (on average 1 ppbv h−1 in the first 3 h and much lower afterwards) due to limited NOx. Peroxyacetyl nitrate concentrations (PAN) increased with plume age and the box model estimated an average production rate of ~80 pptv h−1 in the first 3 h. Like ozone, there was also no evidence for net secondary inorganic or organic aerosol formation. There was no apparent increase in aerosol mass concentrations in the boreal fire plumes due to secondary organic aerosol (SOA) formation; however, there were indications of chemical processing of the organic aerosols. In addition to the detailed studies of boreal fire plume evolution, about 500 smoke plumes intercepted by the NASA DC-8 aircraft were segregated by fire source region. The normalized excess mixing ratios (i.e. ΔX/ΔCO) of gaseous (carbon dioxide, acetonitrile, hydrogen cyanide, toluene, benzene, methane, oxides of nitrogen (NOx), ozone, PAN) and fine aerosol particulate components (nitrate, sulfate, ammonium, chloride, organic aerosols and water soluble organic carbon) of these plumes were compared.
This paper compares measurements of gaseous and particulate emissions from a wide range of biomass-burning plumes intercepted by the NASA DC-8 research aircraft during the three phases of the ARCTAS-2008 experiment: ARCTAS-A, based out of Fairbanks, Alaska, USA (3 April to 19 April 2008); ARCTAS-B based out of Cold Lake, Alberta, Canada (29 June to 13 July 2008); and ARCTAS-CARB, based out of Palmdale, California, USA (18 June to 24 June 2008). Approximately 500 smoke plumes from biomass burning emissions that varied in age from minutes to days were segregated by fire source region and urban emission influences. The normalized excess mixing ratios (NEMR) of gaseous (carbon dioxide, acetonitrile, hydrogen cyanide, toluene, benzene, methane, oxides of nitrogen and ozone) and fine aerosol particulate components (nitrate, sulfate, ammonium, chloride, organic aerosols and water soluble organic carbon) of these plumes were compared. A detailed statistical analysis of the different plume categories for different gaseous and aerosol species is presented in this paper.
The comparison of NEMR values showed that CH4 concentrations were higher in air-masses that were influenced by urban emissions. Fresh biomass burning plumes mixed with urban emissions showed a higher degree of oxidative processing in comparison with fresh biomass burning only plumes. This was evident in higher concentrations of inorganic aerosol components such as sulfate, nitrate and ammonium, but not reflected in the organic components. Lower NOx NEMRs combined with high sulfate, nitrate and ammonium NEMRs in aerosols of plumes subject to long-range transport, when comparing all plume categories, provided evidence of advanced processing of these plumes.
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 self-consistent 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.
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.
Enhanced aridification of Central Asia driven by the combined effects of orogenic surface uplift, Paratethys retreat, changes in atmospheric moisture transport and global cooling is one of the most prominent Cenozoic climate change events of the Northern Hemisphere. Deciphering regional long-term patterns of Central Asian hydrology is, therefore, a key element in understanding the role of Northern Hemisphere mid-latitude drying in the global hydrological system. Pedogenic carbonates record information of the paleoecosystem, the paleogeography, hydrology, tectonic and climatic conditions as well as the sedimentary regime during the time in which they formed. The calcrete-bearing paleosols in continental sedimentary basins yield the possibility for understanding these processes, mechanisms and controlling factors. This study characterizes long-term paleoenvironmental conditions between the late Oligocene and early Miocene in SE Kazakhstan based on stable isotopes, elemental geochemistry and laser ablation U-Pb geochronology from alluvial, fluvial and pedogenic deposits. The main topics addressed in this thesis are paleosol and calcrete formation processes in the light of geochemistry, tectonic and climatic influences on both, as well as the paleoenvironmental and hydrological conditions which led to these processes.
The sedimentological succession of the 14 km-long Kendyrlisai Valley in the Ili Basin, SE Kazakhstan in Central Asia, yields exceptional well-developed paleosols which provide an insight into the relationship between tectonics and sedimentation as well as soil formation processes. The valley accommodates more than 600 m of siliciclastic sediments deposited in a fluvial to alluvial environment in the lower part which grade into lacustrine to fluvial deposition with minor gypsum accumulation in the upper part. The yet undated sediments of Kendyrlisai Valley are compared with a biostratigraphically dated and well described succession in an adjacent location. The correlation of the two successions revealed a proximal-distal relation with the Kendyrlisai Valley succession represents the basin margin position with coarse alluvial-fluvial sediments. The calcrete nodules in Kendyrlisai Valley paleosols give the opportunity of U-Pb dating with LA-ICP-MS which uses the variation in both U and Pb to directly date the calcite. The U-Pb dating revealed a late Oligocene-early Miocene age for the investigated section. The analysis of facies and geochemistry of the paleosol profiles in Kendyrlisai Valley gives insight into calcretization stages, pedogenic processes and landscape stability. Kendyrlisai Valley paleosols show a low to moderate weathering intensity with the highest weathering intensity between 0.5 and 0.8 m depth in the paleosol profile. The comparison between acid leachable and non-acid leachable fraction indicate an incorporation of Ca and Ba in incipient calcrete calcite, whereas clay minerals adsorb Mg and Sr. During the evolution from early to more mature calcrete stages, i.e. calcrete nodules, Ba is lost by dissolution and subsequently adsorbed on clay minerals. The nodules consist almost exclusively of pure calcium carbonate with rare substitution of Ca by Mg indicating a constant supply of Ca by weathering and solution during calcrete formation. The occurrence of Mg-bearing clay minerals indicates high evaporative conditions with alkaline waters (pH >9).
Sampled calcrete nodule cross sections reveal more than one weathering and calcrete formation phase with a higher supply of Ca, Mg and Fe during early phases. The organic and inorganic carbon isotope composition of the calcrete nodules reflects C3 metabolism under occasionally moisture stress, resulting in higher δ13C values and lower respiration rates. This study also presents calculated atmospheric pCO2 values from the calcrete nodules with the equation from Cerling et al. (1999). The average calculated CO2 concentration for Kendyrlisai Valley paleosols is 313 ± 110 ppmv. The comparison with data from other studies during the late Oligocene–early Miocene reveal lower values for Kendyrlisai Valley paleosols, which may stem from an underestimation of the plant-respired CO2.
The knowledge of the variability within the paleosols and calcrete nodules enables a reliable paleoclimatic interpretation for the Kendyrlisai Valley and beyond the Ili Basin. Sedimentary facies and geochemical weathering indices suggest an increased surface and groundwater discharge fed by orographically-enhanced precipitation in the Tien Shan hinterland. In contrast, pedogenic stable isotope data and elevated rates of Mg fixation in clay minerals mirror enhanced rates of evaporation in the vadose zone due to protracted aridification. This study posits that pronounced surface uplift of the Tien Shan Mountains during the Oligocene–Miocene transition promoted regionally increased orographic precipitation and the development of fluvial discharge systems.
The comparison with other studies in adjacent basins creates a precipitation pattern for Central Asia during the late Oligocene to early Miocene. The westerlies supplied Central Asia with atmospheric moisture from the West, possibly from the Eastern Paratethys. The uplifting Tien Shan Mountain ranges captured this moisture on their luv-side, resulting in a pronounced orographic rainfall in the adjacent Ili and Issyk Kul Basins. The Tarim Basin and the Valley of Lakes experienced drier conditions due to the rain shadow effect on the lee-side of the Tien Shan and Altai Mountains. In course of this findings, the thesis highlights the crucial role the Tien Shan Mountains play and had been playing in former times as an orographic barrier for the distribution of atmospheric moisture across Central Asia.
Global warming, changes in the hydrological cycle and enhanced marine primary productivity all have been invoked to have contributed to the occurrence of widespread ocean anoxia during the Cenomanian-Turonian Oceanic Anoxic Event (OAE2; ~ 94 Ma), but disentangling these factors on a regional scale has remained problematic. We generated palynological and organic geochemical records that allow the separation of these forcing factors in a core spanning the OAE2 from Wunstorf, Lower Saxony Basin (LSB; North Gemany), which exhibits cyclic black shale–marl alternations related to the orbital precession cycle.
Despite the widely varying depositional conditions complicating the interpretation of the obtained records, TEX86H indicates that sea-surface temperature (SST) evolution in the LSB during OAE2 resembles that of previously studied sites throughout the proto-North Atlantic. Cooling during the so-called Plenus Cold Event interrupted black shale deposition during the early stages of OAE2. However, TEX86 does not vary significantly across marl–black shale alternations, suggesting that temperature variations did not force the formation of the cyclic black shale horizons. Relative (i.e., with respect to marine palynomorphs) and absolute abundances of pollen and spores are elevated during phases of black shale deposition, indicative of enhanced precipitation and run-off. High abundances of cysts from inferred heterotrophic and euryhaline dinoflagellates supports high run-off, which likely introduced additional nutrients to the epicontinental shelf resulting in elevated marine primary productivity.
We conclude that orbitally-forced enhanced precipitation and run-off, in tandem with elevated marine primary productivity, were critical in cyclic black shale formation on the northwest European epicontinental shelf and potentially for other OAE2 sections in the proto-Atlantic and Western Interior Seaway at similar latitudes as well.
Global warming, changes in the hydrological cycle and enhanced marine primary productivity all have been invoked as having contributed to the occurrence of widespread ocean anoxia during the Cenomanian–Turonian oceanic anoxic event (OAE2; ~94 Ma), but disentangling these factors on a regional scale has remained problematic. In an attempt to separate these forcing factors, we generated palynological and organic geochemical records using a core spanning the OAE2 from Wunstorf, Lower Saxony Basin (LSB; northern Germany), which exhibits cyclic black shale–marl alternations related to the orbital precession cycle.
Despite the widely varying depositional conditions complicating the interpretation of the obtained records, TEX86H indicates that sea-surface temperature (SST) evolution in the LSB during OAE2 resembles that of previously studied sites throughout the proto-North Atlantic. Cooling during the so-called Plenus Cold Event interrupted black shale deposition during the early stages of OAE2. However, TEX86 does not vary significantly across black shale–marl alternations, suggesting that temperature variations did not force the formation of the cyclic black shale horizons. Relative (i.e., with respect to marine palynomorphs) and absolute abundances of pollen and spores are elevated during phases of black shale deposition, indicative of enhanced precipitation and run-off. High abundances of cysts from inferred heterotrophic and euryhaline dinoflagellates supports high run-off, which likely introduced additional nutrients to the epicontinental shelf resulting in elevated marine primary productivity.
We conclude that orbitally forced enhanced precipitation and run-off, in tandem with elevated marine primary productivity, were critical in cyclic black shale formation on the northern European epicontinental shelf and potentially for other OAE2 sections in the proto-Atlantic and Western Interior Seaway at similar latitudes as well.
The complex magnetotelluric (MT) apparent resistivity tensor can be decomposed into two real tensors, the apparent resistivity and the resistivity phase tensors, which represent relationships between the observed electric field at a point on the Earth's surface and an associated apparent current density. We explain the differences between these tensors and conventional estimates of apparent resistivity and phase for simple resistivity environments and demonstrate, using canonical models in 1‐D and 2‐D environments, that both tensors are more sensitive to vertical and horizontal resistivity gradients than their conventional counterparts. The properties of the new tensors are explained using electromagnetic induction theory and the effects of associated charges at resistivity boundaries. We introduce a new way to plot tensor ellipses, which brings significant improvements to the interpretation of MT data, using appropriate visualization software. The apparent resistivity tensor gives information about the magnitude and direction of apparent resistivity subsurface structures and has a strong response to vertical resistivity contrasts. The resistivity phase tensor is highly sensitive to vertical boundaries and the associated fields in the TM mode. It is also free from static distortions under the same conditions implied for the conventional phase tensor. These findings have prompted a study in the potential of the new tensors for 3‐D inversions. The results from a 3‐D inversion of a canonical oblique conductor straddling two quarter spaces show distinct improvements in resolving the boundaries of the conductor and open a promising field for future studies.
In November 2016, magnetotelluric (MT) data were collected at the Ceboruco Volcano in cooperation with the Centro de Sismología y Volcanología de Occidente (SisVoc, Universidad de Guadalajara, Mexico). The Ceboruco is a 2280 m high stratovolcano, located in Nayarit State, Mexico. It is placed in the central part of the Tepic-Zacoalco Rift (TZR), which constitutes the north-western end of the Trans-Mexican Volcanic Belt. Together with Chapala and Colima (in the Jalisco Block), they form the triple rift system developed as a consequence of the ongoing subduction of the Rivera and Cocos oceanic plates beneath the North American continental crust. Although its last eruption occurred in 1870, it is the most active volcano in the area, showing volcanic-earthquake activity together with ongoing vapor emissions. The survey was part of a geothermal project (CeMIEGeo-P24) and focused on the determination of electrical conductivity properties to characterize the deep structure and the geothermal potential of the Volcano. Frequency dependent magnetotelluric response functions were calculated from 25 broadband MT stations, which covered an area of 10 x 10 km2 including its crater, calderas and foreland. The results were interpreted using anisotropic 3-D forward modelling and isotropic 3-D inversion approaches, considering strong topographical effects. The final resistivity model implies a highly conductive layer, reaching from near-surface to approximately 2 km depth, which might be related to a hydrothermal system. Here, mineralized fluids and clay minerals can cause high conductivities around 1 S/m. For longer periods, the principal axes of the MT response tensors (phase tensor, apparent resistivity tensor) are in good agreement with the strike direction of the underlying rift system. However, they are not rendered by the isotropic inversion. Thus the data suggest an anisotropic electrical conductivity at greater depth with its principal axis determined by the response tensors.
The Ceboruco is a 2280 m high stratovolcano located in Nayarit State, Mexico. Despite its last eruption which occurred in 1870, it is the most active volcano in the area, showing volcanicearthquake activity together with ongoing vapor emissions. The magnetotelluric survey was carried out in November 2016. It was part of a geothermal project (CeMIEGeo-P24) and focused on the determination of the electrical conductivity distribution in the subsurface of the volcano.
The Magnetotelluric Apparent Resistivity Tensor, as introduced by Brown (2016), can be decomposed into an amplitude and a phase tensor. The fundamental physics behind those new tensors were presented in Hering et al. (2019), using canonical models in 1-D (isotropic and anisotropic) and 2-D resistivity environments. Here, the tensors are introduced for a high-quality data set, where their interpretational benefits become very obvious. Additionally, results from an isotropic 3-D inversion are presented and compared to an alternative 3-D anisotropic forward model.
During the Holocene, North American ice sheet collapse and rapid sea-level rise reconnected the Black Sea with the global ocean. Rapid meltwater releases into the North Atlantic and associated climate change arguably slowed the pace of Neolithisation across southeastern Europe, originally hypothesized as a catastrophic flooding that fueled culturally-widespread deluge myths. However, we currently lack an independent record linking the timing of meltwater events, sea-level rise and environmental change with the timing of Neolithisation in southeastern Europe. Here, we present a sea surface salinity record from the Northern Aegean Sea indicative of two meltwater events at ~8.4 and ~7.6 kiloyears that can be directly linked to rapid declines in the establishment of Neolithic sites in southeast Europe. The meltwater events point to an increased outflow of low salinity water from the Black Sea driven by rapid sea level rise >1.4 m following freshwater outbursts from Lake Agassiz and the final decay of the Laurentide ice sheet. Our results shed new light on the link between catastrophic sea-level rise and the Neolithisation of southeastern Europe, and present a historical example of how coastal populations could have been impacted by future rapid sea-level rise.
Immersion freezing is the most relevant heterogeneous ice nucleation mechanism through which ice crystals are formed in mixed-phase clouds. In recent years, an increasing number of laboratory experiments utilizing a variety of instruments have examined immersion freezing activity of atmospherically relevant ice-nucleating particles. However, an intercomparison of these laboratory results is a difficult task because investigators have used different ice nucleation (IN) measurement methods to produce these results. A remaining challenge is to explore the sensitivity and accuracy of these techniques and to understand how the IN results are potentially influenced or biased by experimental parameters associated with these techniques.
Within the framework of INUIT (Ice Nuclei Research Unit), we distributed an illite-rich sample (illite NX) as a representative surrogate for atmospheric mineral dust particles to investigators to perform immersion freezing experiments using different IN measurement methods and to obtain IN data as a function of particle concentration, temperature (T), cooling rate and nucleation time. A total of 17 measurement methods were involved in the data intercomparison. Experiments with seven instruments started with the test sample pre-suspended in water before cooling, while 10 other instruments employed water vapor condensation onto dry-dispersed particles followed by immersion freezing. The resulting comprehensive immersion freezing data set was evaluated using the ice nucleation active surface-site density, ns, to develop a representative ns(T) spectrum that spans a wide temperature range (−37 °C < T < −11 °C) and covers 9 orders of magnitude in ns.
In general, the 17 immersion freezing measurement techniques deviate, within a range of about 8 °C in terms of temperature, by 3 orders of magnitude with respect to ns. In addition, we show evidence that the immersion freezing efficiency expressed in ns of illite NX particles is relatively independent of droplet size, particle mass in suspension, particle size and cooling rate during freezing. A strong temperature dependence and weak time and size dependence of the immersion freezing efficiency of illite-rich clay mineral particles enabled the ns parameterization solely as a function of temperature. We also characterized the ns(T) spectra and identified a section with a steep slope between −20 and −27 °C, where a large fraction of active sites of our test dust may trigger immersion freezing. This slope was followed by a region with a gentler slope at temperatures below −27 °C. While the agreement between different instruments was reasonable below ~ −27 °C, there seemed to be a different trend in the temperature-dependent ice nucleation activity from the suspension and dry-dispersed particle measurements for this mineral dust, in particular at higher temperatures. For instance, the ice nucleation activity expressed in ns was smaller for the average of the wet suspended samples and higher for the average of the dry-dispersed aerosol samples between about −27 and −18 °C. Only instruments making measurements with wet suspended samples were able to measure ice nucleation above −18 °C. A possible explanation for the deviation between −27 and −18 °C is discussed. Multiple exponential distribution fits in both linear and log space for both specific surface area-based ns(T) and geometric surface area-based ns(T) are provided. These new fits, constrained by using identical reference samples, will help to compare IN measurement methods that are not included in the present study and IN data from future IN instruments.
Immersion freezing is the most relevant heterogeneous ice nucleation mechanism through which ice crystals are formed in mixed-phase clouds. In recent years, an increasing number of laboratory experiments utilizing a variety of instruments have examined immersion freezing activity of atmospherically relevant ice nucleating particles (INPs). However, an inter-comparison of these laboratory results is a difficult task because investigators have used different ice nucleation (IN) measurement methods to produce these results. A remaining challenge is to explore the sensitivity and accuracy of these techniques and to understand how the IN results are potentially influenced or biased by experimental parameters associated with these techniques.
Within the framework of INUIT (Ice Nucleation research UnIT), we distributed an illite rich sample (illite NX) as a representative surrogate for atmospheric mineral dust particles to investigators to perform immersion freezing experiments using different IN measurement methods and to obtain IN data as a function of particle concentration, temperature (T), cooling rate and nucleation time. Seventeen measurement methods were involved in the data inter-comparison. Experiments with seven instruments started with the test sample pre-suspended in water before cooling, while ten other instruments employed water vapor condensation onto dry-dispersed particles followed by immersion freezing. The resulting comprehensive immersion freezing dataset was evaluated using the ice nucleation active surface-site density (ns) to develop a representative ns(T) spectrum that spans a wide temperature range (−37 °C < T < −11 °C) and covers nine orders of magnitude in ns.
Our inter-comparison results revealed a discrepancy between suspension and dry-dispersed particle measurements for this mineral dust. While the agreement was good below ~ −26 °C, the ice nucleation activity, expressed in ns, was smaller for the wet suspended samples and higher for the dry-dispersed aerosol samples between about −26 and −18 °C. Only instruments making measurement techniques with wet suspended samples were able to measure ice nucleation above −18 °C. A possible explanation for the deviation between −26 and −18 °C is discussed. In general, the seventeen immersion freezing measurement techniques deviate, within the range of about 7 °C in terms of temperature, by three orders of magnitude with respect to ns. In addition, we show evidence that the immersion freezing efficiency (i.e., ns) of illite NX particles is relatively independent on droplet size, particle mass in suspension, particle size and cooling rate during freezing. A strong temperature-dependence and weak time- and size-dependence of immersion freezing efficiency of illite-rich clay mineral particles enabled the ns parameterization solely as a function of temperature. We also characterized the ns (T) spectra, and identified a section with a steep slope between −20 and −27 °C, where a large fraction of active sites of our test dust may trigger immersion freezing. This slope was followed by a region with a gentler slope at temperatures below −27 °C. A multiple exponential distribution fit is expressed as ns(T) = exp(23.82 × exp(−exp(0.16 × (T + 17.49))) + 1.39) based on the specific surface area and ns(T) = exp(25.75 × exp(−exp(0.13 × (T + 17.17))) + 3.34) based on the geometric area (ns and T in m−2 and °C, respectively). These new fits, constrained by using an identical reference samples, will help to compare IN measurement methods that are not included in the present study and, thereby, IN data from future IN instruments.
Droplets produced in a cloud condensation nuclei chamber (CCNC) as a function of supersaturation have been separated from unactivated aerosol particles using counterflow virtual impaction. Residual material after droplets were evaporated was chemically analyzed with an Aerodyne Aerosol Mass Spectrometer (AMS) and the Particle Analysis by Laser Mass Spectrometry (PALMS) instrument. Experiments were initially conducted to verify activation conditions for monodisperse ammonium sulfate particles and to determine the resulting droplet size distribution as a function of supersaturation. Based on the observed droplet size, the counterflow virtual impactor cut-size was set to differentiate droplets from unactivated interstitial particles. Validation experiments were then performed to verify that only droplets with sufficient size passed through the counterflow virtual impactor for subsequent analysis. A two-component external mixture of monodisperse particles was also exposed to a supersaturation which would activate one of the types (hygroscopic salts) but not the other (polystyrene latex spheres or adipic acid). The mass spectrum observed after separation indicated only the former, validating separation of droplets from unactivated particles. Results from ambient measurements using this technique and AMS analysis were inconclusive, showing little chemical differentiation between ambient aerosol and activated droplet residuals, largely due to low signal levels. When employing as single particle mass spectrometer for compositional analysis, however, we observed enhancement of sulfate in droplet residuals.
Droplets produced in a cloud condensation nucleus chamber as a function of supersaturation have been separated from unactivated aerosol particles using counterflow virtual impaction. Residual material after droplets were evaporated was chemically analyzed with an Aerodyne Aerosol Mass Spectrometer and the Particle Analysis by Laser Mass Spectrometry instrument. Experiments were initially conducted to verify activation conditions for monodisperse ammonium sulfate particles and to determine the resulting droplet size distribution as a function of supersaturation. Based on the observed droplet size, the counterflow virtual impactor cut-size was set to differentiate droplets from unactivated interstitial particles. Validation experiments were then performed to verify that only droplets with sufficient size passed through the counterflow virtual impactor for subsequent analysis. A two-component external mixture of monodisperse particles was also exposed to a supersaturation which would activate one of the types (ammonium sulfate) but not the other (polystyrene latex spheres). The mass spectrum observed after separation indicated only the former, validating separation of droplets from unactivated particles. Results from atmospheric measurements using this technique indicate that aerosol particles often activate predominantly as a function of particle size. Chemical composition is not irrelevant, however, and we observed enhancement of sulfate in droplet residuals using single particle analysis.
The international research project RECONCILE has addressed central questions regarding polar ozone depletion, with the objective to quantify some of the most relevant yet still uncertain physical and chemical processes and thereby improve prognostic modelling capabilities to realistically predict the response of the ozone layer to climate change. This overview paper outlines the scope and the general approach of RECONCILE, and it provides a summary of observations and modelling in 2010 and 2011 that have generated an in many respects unprecedented dataset to study processes in the Arctic winter stratosphere. Principally, it summarises important outcomes of RECONCILE including (i) better constraints and enhanced consistency on the set of parameters governing catalytic ozone destruction cycles, (ii) a better understanding of the role of cold binary aerosols in heterogeneous chlorine activation, (iii) an improved scheme of polar stratospheric cloud (PSC) processes that includes heterogeneous nucleation of nitric acid trihydrate (NAT) and ice on non-volatile background aerosol leading to better model parameterisations with respect to denitrification, and (iv) long transient simulations with a chemistry-climate model (CCM) updated based on the results of RECONCILE that better reproduce past ozone trends in Antarctica and are deemed to produce more reliable predictions of future ozone trends. The process studies and the global simulations conducted in RECONCILE show that in the Arctic, ozone depletion uncertainties in the chemical and microphysical processes are now clearly smaller than the sensitivity to dynamic variability.
Significant reductions in stratospheric ozone occur inside the polar vortices each spring when chlorine radicals produced by heterogeneous reactions on cold particle surfaces in winter destroy ozone mainly in two catalytic cycles, the ClO dimer cycle and the ClO/BrO cycle. Chlorofluorocarbons (CFCs), which are responsible for most of the chlorine currently present in the stratosphere, have been banned by the Montreal Protocol and its amendments, and the ozone layer is predicted to recover to 1980 levels within the next few decades. During the same period, however, climate change is expected to alter the temperature, circulation patterns and chemical composition in the stratosphere, and possible geo-engineering ventures to mitigate climate change may lead to additional changes. To realistically predict the response of the ozone layer to such influences requires the correct representation of all relevant processes. The European project RECONCILE has comprehensively addressed remaining questions in the context of polar ozone depletion, with the objective to quantify the rates of some of the most relevant, yet still uncertain physical and chemical processes. To this end RECONCILE used a broad approach of laboratory experiments, two field missions in the Arctic winter 2009/10 employing the high altitude research aircraft M55-Geophysica and an extensive match ozone sonde campaign, as well as microphysical and chemical transport modelling and data assimilation. Some of the main outcomes of RECONCILE are as follows: (1) vortex meteorology: the 2009/10 Arctic winter was unusually cold at stratospheric levels during the six-week period from mid-December 2009 until the end of January 2010, with reduced transport and mixing across the polar vortex edge; polar vortex stability and how it is influenced by dynamic processes in the troposphere has led to unprecedented, synoptic-scale stratospheric regions with temperatures below the frost point; in these regions stratospheric ice clouds have been observed, extending over >106km2 during more than 3 weeks. (2) Particle microphysics: heterogeneous nucleation of nitric acid trihydrate (NAT) particles in the absence of ice has been unambiguously demonstrated; conversely, the synoptic scale ice clouds also appear to nucleate heterogeneously; a variety of possible heterogeneous nuclei has been characterised by chemical analysis of the non-volatile fraction of the background aerosol; substantial formation of solid particles and denitrification via their sedimentation has been observed and model parameterizations have been improved. (3) Chemistry: strong evidence has been found for significant chlorine activation not only on polar stratospheric clouds (PSCs) but also on cold binary aerosol; laboratory experiments and field data on the ClOOCl photolysis rate and other kinetic parameters have been shown to be consistent with an adequate degree of certainty; no evidence has been found that would support the existence of yet unknown chemical mechanisms making a significant contribution to polar ozone loss. (4) Global modelling: results from process studies have been implemented in a prognostic chemistry climate model (CCM); simulations with improved parameterisations of processes relevant for polar ozone depletion are evaluated against satellite data and other long term records using data assimilation and detrended fluctuation analysis. Finally, measurements and process studies within RECONCILE were also applied to the winter 2010/11, when special meteorological conditions led to the highest chemical ozone loss ever observed in the Arctic. In addition to quantifying the 2010/11 ozone loss and to understand its causes including possible connections to climate change, its impacts were addressed, such as changes in surface ultraviolet (UV) radiation in the densely populated northern mid-latitudes.
Hydro-climatic causes of widespread floods in central Europe : on rain-on-snow and Vb-cyclone events
(2021)
The presented work investigates the hydro-meteorological and hydro-climatological drivers of widespread floods in Central Europe during the past century. Due to the strong seasonality of the detected flood drivers, the thesis is divided into two parts: the first part focuses on widespread winter floods and the second one on extreme summer floods. For analysing past flood events, we profited from the dynamically downscaled centennial ERA-20C reanalysis (continuously from 1901—2010). The downscaling was performed over Europe with a coupled regional atmosphere-ocean model (COSMO-CLM+NEMO) to represent the water cycle more realistic. These high resolution atmospheric data allowed us to study the four-dimensional atmospheric state during selected floods during the early decades of the 20th century for the first time with such a high temporal and spatial resolution.
During the winter half-year, the observed floods were particularly widespread. High peak discharges were recorded simultaneously in the Rhine, Elbe, and Danube catchments. Most of these trans-basin floods were compound events caused by rainfall during extensive snowmelt (i.e., rain-on-snow events). Interestingly, the winter flood time series exhibited a remarkable high flood frequency during the 1940s and 1980s, while other decades were flood-poor. We detected a synchronization of the inter-annual flood frequency with the superposition of the North Atlantic Oscillation (NAO) and the Scandinavian pattern (SCA). The negative NAO phase is often associated with large snowfall and cyclone tracks over southern Europe, while the negative SCA pattern correlates with total precipitation in the affected river catchments.
During the summer half-year, most extreme floods in Central Europe were caused by so-called Vb-cyclones propagating from the Mediterranean Sea north-eastward to Central Europe. So far in the literature, only a few Vb-events, which occurred during the past two decades, have been analysed. We extended the previous case studies by several past Vb-cyclone floods since 1900. We investigated the processes that intensify Vb-cyclone precipitation with Lagrangian moisture-source diagnostics and the parametric transfer entropy measure TE-linear. Overall, an enhanced and dynamically driven moisture uptake over the Mediterranean Sea was found to be characteristic for Vb-events with heavy precipitation. This is supported by high information exchange from evaporation over the western basin of the Mediterranean Sea towards heavy precipitation in the Odra catchment. The dominating moisture uptake regions during the investigated events were, however, the European continent and the North Sea. A possible cause could be the pre-moistening of non-saturated continental moisture sources upstream of the affected river catchments as indicated by significant information exchange from land surface evaporation and soil moisture content along the Vb-cyclone pathway. Besides, evaporation over the Mediterranean Sea might contribute to Vb-cyclone intensification in the early stages of their development through latent heat release. On the catchment scale, orographic rainfall and convective precipitation further enhance the flood triggering rainfall. As expected, the Vb-cyclones mainly trigger precipitation along west-east orientated mountain ranges such as the Alps or Ore mountains due to their meridional pathway. Remarkably, during summer, we detected a convective fraction of up to 90% during the afternoons of individual days and up to 23% on average (based on convective cell tracking and convection-permitting simulations of selected flood events since 1900).
The presented analyses deepened the knowledge on atmospheric and hydroclimatic drivers of widespread floods in Central Europe. This will serve as a basis for future studies on the predictability of floods induced by rain-on-snow and Vb-cyclone precipitation events in the context of a changing climate.
Several past summer floods in Central Europe were associated with so-called Vb‑cyclones propagating from the Mediterranean Sea north-eastward to Central Europe. This study illustrates the usefulness of the parametric transfer entropy measure TE‑linear in investigating heavy Vb‑cyclone precipitation events in the Odra catchment (Poland). With the application of the TE‑linear approach, we confirm the impact of the Mediterranean Sea on precipitation intensification. Moreover, we also detect significant information exchange to Vb‑cyclone precipitation from evaporation over the European continent along the typical Vb‑cyclone pathway. Thus, the Mediterranean Sea could enhance the Vb‑cyclone precipitation by pre-moistening continental moisture source regions that contribute to precipitation downstream in the investigated catchments. Overall, the transfer entropy approach with the measure TE‑linear proved to be computationally effective and complementary to traditional methods such as Lagrangian and Eulerian diagnostics.
Moisture sources of heavy precipitation in Central Europe in synoptic situations with Vb-cyclones
(2022)
During the past century, several extreme summer floods in Central Europe were associated with so-called Vb-cyclones propagating from the Mediterranean Sea north-eastward to Central Europe. The processes intensifying the precipitation in synoptic situations with Vb-cyclones in the Danube, Elbe, and Odra catchments are only partially understood. Our study aims to investigate these processes with Lagrangian moisture-source diagnostics for 16 selected Vb-events. Moreover, we analyse the characteristics of typical moisture source regions during 1107 Vb-events from 1901 to 2010 based on ERA-20C reanalysis dynamically downscaled with COSMO-CLM+NEMO. We observe moisture contributions by various source regions highlighting the complex dynamical interplay of different air masses leading to moisture convergence in synoptic situations with Vb-cyclones. Overall, up to 80% of the precipitation originates from the European continent, indicating the importance of continental moisture recycling, especially within the respective river catchment. Other major moisture uptake regions are the North Sea, the Baltic Sea, the North Atlantic, and for a few events the Black Sea. Remarkably, anomalies in these oceanic source regions show no connection to precipitation amounts in synoptic situations with Vb-cyclones. In contrast, the Vb-cyclones with the highest precipitation are associated with anomalously high evaporation in the Mediterranean Sea, even though the Mediterranean Sea is only a minor moisture source region on average. Interestingly, the evaporation anomalies are not connected with sea-surface temperature but with wind-speed anomalies (Spearman’s rank correlation coefficient R≈0.7, significant with p<0.01) indicating mainly dynamically driven evaporation. The particular role of the Mediterranean Sea hints towards possible importance of Mediterranean moisture for the early-stage intensification of Vb-cyclones and the pre-moistening of the continental uptake regions upstream of the target catchments.
Water footprints have been proposed as sustainability indicators, relating the consumption of goods like food to the amount of water necessary for their production and the impacts of that water use in the source regions. We have further developed the existing water footprint methodology by globally resolving virtual water flows and import and source regions at 5 arc minutes spatial resolution, and by assessing local impacts of export production. Applying this method to three exemplary cities, Berlin, Delhi and Lagos, we find major differences in amounts, composition, and origin of green and blue virtual water imports, due to differences in diets, trade integration and crop water productivities in the source regions. While almost all of Delhi's and Lagos' virtual water imports are of domestic origin, Berlin on average imports from more than 4000 km distance, in particular soy (livestock feed), coffee and cocoa. While 42% of Delhi's virtual water imports are blue water based, the fractions for Berlin and Lagos are 2% and 0.5%, respectively, roughly equal to local drinking water abstractions of these cities. Some of the external source regions of Berlin's virtual water imports appear to be critically water scarce and/or food insecure. However for deriving recommendations on sustainable consumption and trade, further analysis of context-specific costs and benefits associated with export production will be required.
Water footprints have been proposed as sustainability indicators, relating the consumption of goods like food to the amount of water necessary for their production and the impacts of that water use in the source regions. We further developed the existing water footprint methodology, by globally resolving virtual water flows from production to consumption regions for major food crops at 5 arcmin spatial resolution. We distinguished domestic and international flows, and assessed local impacts of export production. Applying this method to three exemplary cities, Berlin, Delhi and Lagos, we find major differences in amounts, composition, and origin of green and blue virtual water imports, due to differences in diets, trade integration and crop water productivities in the source regions. While almost all of Delhi's and Lagos' virtual water imports are of domestic origin, Berlin on average imports from more than 4000 km distance, in particular soy (livestock feed), coffee and cocoa. While 42% of Delhi's virtual water imports are blue water based, the fractions for Berlin and Lagos are 2 and 0.5%, respectively, roughly equal to the water volumes abstracted in these two cities for domestic water use. Some of the external source regions of Berlin's virtual water imports appear to be critically water scarce and/or food insecure. However, for deriving recommendations on sustainable consumption and trade, further analysis of context-specific costs and benefits associated with export production will be required.
From July 2002 to March 2004 the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) aboard the European Space Agency´s Environmental Satellite (Envisat) measured nearly continuously mid infrared limb radiance spectra. These measurements are utilised to retrieve the global distribution of the chlorofluorocarbon CFC-11 by applying a new fast forward model for Envisat MIPAS and an accompanying optimal estimation retrieval processor. A detailed analysis shows that the total retrieval errors of the individual CFC-11 volume mixing ratios are typically below 10% in the altitude range 10 to 25 km and that the systematic components dominate. Contribution of a priori information to the retrieval results are less than 5 to 10% and the vertical resolution of the observations is about 3 to 4 km in the same vertical range. The data are successfully validated by comparison with several other space experiments, an air-borne in-situ instrument, measurements from ground-based networks, and independent Envisat MIPAS analyses. The retrieval results from 425 000 Envisat MIPAS limb scans are compiled to provide a new climatological data set of CFC-11. The climatology shows significantly lower CFC-11 abundances in the lower stratosphere compared with the Reference Atmospheres for MIPAS (RAMstan V3.1) climatology. Depending on the atmospheric conditions the differences between the climatologies are up to 30 to 110 ppt (45 to 150%) at 19 to 27 km altitude. Additionally, time series of CFC-11 mean abundance and variability for five latitudinal bands are presented. The observed CFC-11 distributions can be explained by the residual mean circulation and large-scale eddy-transports in the upper troposphere and lower stratosphere. The new CFC-11 data set is well suited for further scientific studies.
We present the application of Time-of-Flight Mass Spectrometry (TOF MS) for the analysis of halocarbons in the atmosphere, after cryogenic sample preconcentration and gas chromatographic separation. For the described field of application, the Quadrupole Mass Spectrometer (QP MS) is the state-of-the-art detector. This work aims at comparing two commercially available instruments, a QP MS and a TOF MS with respect to mass resolution, mass accuracy, sensitivity, measurement precision and detector linearity. Both mass spectrometers are operated on the same gas chromatographic system by splitting the column effluent to both detectors. The QP MS had to be operated in optimised Single Ion Monitoring (SIM) mode to achieve a sensitivity which could compete with the TOF MS. The TOF MS provided full mass range information in any acquired mass spectrum without losing sensitivity. Whilst the QP MS showed the performance already achieved in earlier tests, the sensitivity of the TOF MS was on average higher than that of the QP MS in the "operational" SIM mode by a factor of up to 3 reaching detection limits of less than 0.2 pg. Measurement precision determined for the whole analytical system was up to 0.2% depending on substance and sampled volume. The TOF MS instrument used for this study displayed significant non-linearities of up to 10% for two third of all analysed substances.
We present the application of time-of-flight mass spectrometry (TOF MS) for the analysis of halocarbons in the atmosphere after cryogenic sample preconcentration and gas chromatographic separation. For the described field of application, the quadrupole mass spectrometer (QP MS) is a state-of-the-art detector. This work aims at comparing two commercially available instruments, a QP MS and a TOF MS, with respect to mass resolution, mass accuracy, stability of the mass axis and instrument sensitivity, detector sensitivity, measurement precision and detector linearity. Both mass spectrometers are operated on the same gas chromatographic system by splitting the column effluent to both detectors. The QP MS had to be operated in optimised single ion monitoring (SIM) mode to achieve a sensitivity which could compete with the TOF MS. The TOF MS provided full mass range information in any acquired mass spectrum without losing sensitivity. Whilst the QP MS showed the performance already achieved in earlier tests, the sensitivity of the TOF MS was on average higher than that of the QP MS in the "operational" SIM mode by a factor of up to 3, reaching detection limits of less than 0.2 pg. Measurement precision determined for the whole analytical system was up to 0.2% depending on substance and sampled volume. The TOF MS instrument used for this study displayed significant non-linearities of up to 10% for two-thirds of all analysed substances.
The TTL is the transition layer between the tropical troposphere and stratosphere, and is the main region where tropospheric air enters the stratosphere. In this thesis different transport processes are studied by using in situ measurements of tracers. Long-lived tracers were measured with the High Altitude Gas Analyzer (HAGAR) on board the M55 Geophysica aircraft. The instrument was developed by the University of Frankfurt and measures the long-lived tracers CO2, N2O, CFC-12, CFC-11, H-1211, SF6, CH4 and H2 with two gas chromatographic channels and a CO2 sensor (LICOR). The measurements are supported by CO and O3 measurements of other instruments. Two campaigns were conducted to obtain measurements in the TTL: SCOUT-O3 (November/December 2005 in Darwin, Australia) and AMMA-SCOUT-O3 (August 2006 in Ouagadougou, Burkina Faso). After a general introduction of the thesis in chapters one and two, the third chapter describes the findings during this last campaign. Five local flights are analyzed to study the different transport processes that occur in the tropical tropopause layer above West-Africa: deep convection up to the level of main convective outflow, vertical mixing after overshooting of air in deep convection, horizontal inmixing from the extratropical lower stratosphere, and horizontal transport across the subtropical barrier. Main findings are that the TTL over West-Africa is mostly influenced by remote convection. The subtropical barrier is not a strong barrier but more a region of transition between the extratropical and the tropical stratosphere. Chapter 4 presents the results obtained during the SCOUT-O3 campaign. From the eight local flights the last four flights (051129, 051130a, 051130b, 051205) show enhanced values of ozone, CO and CO2 between 355 and 380 K potential temperature in comparison with the first four flights (051116, 051119, 051123, 051125). Horizontal inmixing from the extra-tropical stratosphere and influence of the local convective system Hector cannot explain the enhanced values of the two flights on 30 November Therefore, other possible explanations for these enhanced CO, CO2 and ozone levels are proposed. The first explanation is vertical mixing in the vicinity of the jet stream. However, the jet cannot explain the differences between the flights on 30 November and the flights on 29 November and 5 December. Another possible explanation is influence of polluted boundary layer air masses from the Indonesian region. Especially air sampled during the flights on November 30 crossed large parts of northern Indonesia between 8 and 10 days before the measurements. Convective uplift of biomass burning and other pollution plumes can transport CO and ozone precursors into the upper troposphere, where they can significantly enhance the ozone production. The last chapter deals with the vertical ascent rate in the TTL and uses measurements of both the SCOUT-O3 and AMMA-SCOUT-O3 campaign as well as data from previous aircraft campaigns (TROCCINOX and APE-THESEO). Time scales and residence times for mean vertical transport in the background TTL are estimated for different seasons and over different geographic regions using in situ observations of CO2 and long-lived tracers. The vertical transport time scales are constrained using the seasonal variation of CO2 in the tropical troposphere as a “tracer clock” for vertical ascent. Two methods are applied to calculate the residence time in the layer between 360 and 390 K potential temperature. The first method uses the slope of the CO2 index, the second method fits the CO2 index directly to the measurements assuming a constant ascent rate. The first method yields residence times for Australia,West Africa, and Brazil of the same order, 35-45 days to 380 K and 50 days to 390 K (where no value can be derived for Australia as the slope is changing approximately one month before the campaign). For APE-THESEO, the method does not yield reasonable results. The best estimates using the second method show moderate residence times between 360 and 390 K of 60±25 days SCOUT-O3 (NH autumn) and 43±8 days for AMMA/SCOUT-O3 (NH summer). These results agree well with the results calculated using the first method. For APE-THESEO and TROCCINOX the best fits yield shorter residence times of 23±7 and 40±10 days, respectively, both during winter. These results correspond well to the expectations based on the seasonal variation of the Brewer-Dobson circulation.
Tracer measurements in the tropical tropopause layer during the AMMA/SCOUT-O3 aircraft campaign
(2009)
We present airborne in situ measurements made during the AMMA (African Monsoon Multidisciplinary Analysis)/SCOUT-O3 campaign between 31 July and 17 August 2006 on board the M55 Geophysica aircraft, based in Ouagadougou, Burkina Faso. CO2 and N2O were measured with the High Altitude Gas Analyzer (HAGAR), CO was measured with the Cryogenically Operated Laser Diode (COLD) instrument, and O3 with the Fast Ozone ANalyzer (FOZAN). We analyze the data obtained during five local flights to study the dominant transport processes controlling the tropical tropopause layer (TTL) above West-Africa: deep convection up to the level of main convective outflow, overshooting of deep convection, horizontal inmixing across the subtropical tropopause, and horizontal transport across the subtropical barrier. Except for the flight of 13 August, distinct minima in CO2 indicate convective outflow of boundary layer air in the TTL. The CO2 profiles show that the level of main convective outflow was mostly located between 350 and 360 K, and for 11 August reached up to 370 K. While the CO2 minima indicate quite significant convective influence, the O3 profiles suggest that the observed convective signatures were mostly not fresh, but of older origin. When compared with the mean O3 profile measured during a previous campaign over Darwin in November 2005, the O3 minimum at the main convective outflow level was less pronounced over Ouagadougou. Furthermore O3 mixing ratios were much higher throughout the whole TTL and, unlike over Darwin, rarely showed low values observed in the regional boundary layer. Signatures of irreversible mixing following overshooting of convective air were scarce in the tracer data. Some small signatures indicative of this process were found in CO2 profiles between 390 and 410 K during the flights of 4 and 8 August, and in CO data at 410 K on 7 August. However, the absence of expected corresponding signatures in other tracer data makes this evidence inconclusive, and overall there is little indication from the observations that overshooting convection has a profound impact on TTL composition during AMMA. We find the amount of photochemically aged air isentropically mixed into the TTL across the subtropical tropopause to be not significant. Using the N2O observations we estimate the fraction of aged extratropical stratospheric air in the TTL to be 0.0±0.1 up to 370 K during the local flights, increasing above this level to 0.2±0.15 at 390 K. The subtropical barrier, as indicated by the slope of the correlation between N2O and O3 between 415 and 490 K, does not appear as a sharp border between the tropics and extratropics, but rather as a gradual transition region between 10 and 25° N latitude where isentropic mixing between these two regions may occur.
Tracer measurements in the tropical tropopause layer during the AMMA/SCOUT-O3 aircraft campaign
(2010)
We present airborne in situ measurements made during the AMMA (African Monsoon Multidisciplinary Analysis)/SCOUT-O3 campaign between 31 July and 17 August 2006 on board the M55 Geophysica aircraft, based in Ouagadougou, Burkina Faso. CO<sub>2</sub> and N<sub>2</sub>O were measured with the High Altitude Gas Analyzer (HAGAR), CO was measured with the Cryogenically Operated Laser Diode (COLD) instrument, and O<sub>3</sub> with the Fast Ozone ANalyzer (FOZAN). We analyze the data obtained during five local flights to study the dominant transport processes controlling the tropical tropopause layer (TTL) above West-Africa: deep convection up to the level of main convective outflow, overshooting of deep convection, horizontal inmixing across the subtropical tropopause, and horizontal transport across the subtropical barrier. Except for the flight of 13 August, distinct minima in CO<sub>2</sub> indicate convective outflow of boundary layer air in the TTL. The CO<sub>2</sub> profiles show that the level of main convective outflow was mostly located between 350 and 360 K, and for 11 August reached up to 370 K. While the CO<sub>2</sub> minima indicate quite significant convective influence, the O<sub>3</sub> profiles suggest that the observed convective signatures were mostly not fresh, but of older origin. When compared with the mean O<sub>3</sub> profile measured during a previous campaign over Darwin in November 2005, the O<sub>3</sub> minimum at the main convective outflow level was less pronounced over Ouagadougou. Furthermore O<sub>3</sub> mixing ratios were much higher throughout the whole TTL and, unlike over Darwin, rarely showed low values observed in the regional boundary layer. Signatures of irreversible mixing following overshooting of convective air were scarce in the tracer data. Some small signatures indicative of this process were found in CO<sub>2</sub> profiles between 390 and 410 K during the flights of 4 and 8 August, and in CO data at 410 K on 7 August. However, the absence of expected corresponding signatures in other tracer data makes this evidence inconclusive, and overall there is little indication from the observations that overshooting convection has a profound impact on TTL composition during AMMA. We find the amount of photochemically aged air isentropically mixed into the TTL across the subtropical tropopause to be not significant. Using the N<sub>2</sub>O observations we estimate the fraction of aged extratropical stratospheric air in the TTL to be 0.0±0.1 up to 370 K during the local flights, increasing above this level to 0.2±0.15 at 390 K. The subtropical barrier, as indicated by the slope of the correlation between N<sub>2</sub>O and O<sub>3</sub> between 415 and 490 K, does not appear as a sharp border between the tropics and extratropics, but rather as a gradual transition region between 10 and 25° N latitude where isentropic mixing between these two regions may occur.
A new global synthesis and biomization of long (> 40 kyr) pollen-data records is presented, and used with simulations from the HadCM3 and FAMOUS climate models to analyse the dynamics of the global terrestrial biosphere and carbon storage over the last glacial–interglacial cycle. Global modelled (BIOME4) biome distributions over time generally agree well with those inferred from pollen data. The two climate models show good agreement in global net primary productivity (NPP). NPP is strongly influenced by atmospheric carbon dioxide (CO2) concentrations through CO2 fertilization. The combined effects of modelled changes in vegetation and (via a simple model) soil carbon result in a global terrestrial carbon storage at the Last Glacial Maximum that is 210–470 Pg C less than in pre-industrial time. Without the contribution from exposed glacial continental shelves the reduction would be larger, 330–960 Pg C. Other intervals of low terrestrial carbon storage include stadial intervals at 108 and 85 kaBP, and between 60 and 65 kaBP during Marine Isotope Stage 4. Terrestrial carbon storage, determined by the balance of global NPP and decomposition, influences the stable carbon isotope composition (δ 13C) of seawater because terrestrial organic carbon is depleted in 13C. Using a simple carbon-isotope mass balance equation we find agreement in trends between modelled ocean δ 13C based on modelled land carbon storage, and palaeo-archives of ocean δ 13C, confirming that terrestrial carbon storage variations may be important drivers of ocean δ 13 C changes.
A new global synthesis and biomization of long (> 40 kyr) pollen-data records is presented and used with simulations from the HadCM3 and FAMOUS climate models and the BIOME4 vegetation model to analyse the dynamics of the global terrestrial biosphere and carbon storage over the last glacial–interglacial cycle. Simulated biome distributions using BIOME4 driven by HadCM3 and FAMOUS at the global scale over time generally agree well with those inferred from pollen data. Global average areas of grassland and dry shrubland, desert, and tundra biomes show large-scale increases during the Last Glacial Maximum, between ca. 64 and 74 ka BP and cool substages of Marine Isotope Stage 5, at the expense of the tropical forest, warm-temperate forest, and temperate forest biomes. These changes are reflected in BIOME4 simulations of global net primary productivity, showing good agreement between the two models. Such changes are likely to affect terrestrial carbon storage, which in turn influences the stable carbon isotopic composition of seawater as terrestrial carbon is depleted in 13C.
A comprehensive evaluation of seasonal backward trajectories initialized in the Northern Hemisphere lowermost stratosphere (LMS) has been performed to investigate the origin of air parcels and the main mechanisms determining characteristic structures in H2O and CO within the LMS. In particular we explain the fundamental role of the transit time since last tropopause crossing (tTST) for the chemical structure of the LMS as well as the feature of the extra-tropical tropopause transition layer (ExTL) as identified from CO profiles. The distribution of H2O in the background LMS above Θ=320 K and 340 K in northern winter and summer, respectively, is found to be governed mainly by the saturation mixing ratio, which in turn is determined by the Lagrangian Cold Point (LCP) encountered by each trajectory. Most of the backward trajectories from this region in the LMS experienced their LCP in the tropics and sub-tropics. The transit time since crossing the tropopause from the troposphere to the stratosphere (tTST) is independent of the H2O value of the air parcel. TST often occurs 20 days after trajectories have encountered their LCP. CO, on the other hand, depends strongly on tTST due to its finite lifetime. The ExTL as identified from CO measurements is then explained as a layer of air just above the tropopause, which on average encountered TST fairly recently.
A comprehensive evaluation of seasonal backward trajectories initialized in the northern hemisphere lowermost stratosphere (LMS) has been performed to investigate the factors that determine the temporal and spatial structure of troposphere-to-stratosphere-transport (TST) and it’s impact on the LMS. In particular we explain the fundamental role of the transit time since last TST (tTST) for the chemical composition of the LMS. According to our results the structure of the LMS can be characterized by a layer with tTST<40 days forming a narrow band around the local tropopause. This layer extends about 30K above the local dynamical tropopause, corresponding to the extratropical tropopause transition layer (ExTL) as identified by CO. The LMS beyond this layer shows a relatively well defined separation as marked by an aprupt transition to longer tTST indicating less frequent mixing and a smaller fraction of tropospheric air. Thus the LMS constitutes a region of two well defined regimes of tropospheric influence. These can be characterized mainly by different transport times from the troposphere and different fractions of tropospheric air. Carbon monoxide (CO) mirrors this structure of tTST due to it’s finite lifetime on the order of three months. Water vapour isopleths, on the other hand, do not uniquely indicate TST and are independent of tTST, but are determined by the Lagrangian Cold Point (LCP) of air parcels. Most of the backward trajectories from the LMS experienced their LCP in the tropics and sub-tropics, and TST often occurs 20 days after trajectories have encountered their LCP. Therefore, ExTL properties deduced from CO and H2O provide totally different informations on transport and particular TST for the LMS.
Knowledge about the initial tectonic and depositional dynamics, as well as the influence of early rifting on climate and environmental evolution remains speculative to a large extent, because sediments are usually deeply buried. Within the East African Rift System, inversion tectonics uplifted a few of these successions to the surface hence presenting rare windows into the pre-rift depositional history. One such example, an exceptional 700 m long and up to 60 m high fresh road cut provided the opportunity to study in detail initial rift successions of the southern Albertine Rift (Western Uganda). This focusses on the basal and poorly known Middle to Late Miocene in order to unravel the climatic, environmental, hydrological and tectonic evolution of the initial Albertine Rift. A large and robust multi-proxy dataset was gathered comprising 169 m of stratigraphic thickness, which spans from 14.5 to 4.9 Ma according to a revised lithostratigraphic model. Fieldwork comprised logging of the sedimentary record, spectral gamma ray, magnetic susceptibility and 2D wall mapping with photomosaics. Additionally, the sections were sampled for bulk mineral and clay mineral analysis. The succession exposes a suite of lithofacies and architectural elements detailing the evolution of a fluvio-lacustrine system. Five depositional environments were identified which show an overall back-stepping trend from an alluvial plain to a delta plain and finally palustrine/shallow lacustrine conditions. Mesoscale base-level cycles, preservation potential of architectural elements, and stacking pattern exhibit limited accommodation space. However, it increases over time. This overall trend indicates increasing tectonic subsidence, which can be explained by flexural downwarp within the pre-rift phase and in the upper part grading into fault-controlled crustal extension of the syn-rift phase, which more and more disrupted a large-scale river system. From the Middle Miocene up to the early Pliocene, this study revealed that palaeoclimate trends become marked by increasing and more fluctuating Th concentrations, loss of feldspar, intercalated lenses of hydroxosulphate minerals, and a shift from smectite-dominated to kaolinite-dominated clays. These signals are all interpreted as detrital except for the hydroxosulphates, and they mirror the increasing intensity of chemical weathering and stripping of soils in the catchment. A trend towards increasing humidity is supported by an increase in lacustrine sediment facies and a lake-level rise. Nevertheless, intercalation of hydroxosulphate, ferricretes and pedogenised horizons prove ongoing seasonality and dry intervals. Finally, based on a revised stratigraphic model a sequence stratigraphic correlation of the outcrop's depositional cycles with basin-scale cycles is presented. According to these cycles, transition from the pre-rift to the syn-rift stage is marked by an unconformity and a tectonic pulse in the latest Miocene. However, the response of fluvial supply, the depositional system as well as climate conditions are less punctuated and characterised by gradual trends and temporal delays. The long pre-rift phase (ca 10 Myr) and the gradual transition to the syn-rift phase is in accordance with the active rifting model, which is based on thermal thinning of the lithosphere by asthenospheric upwelling.
Emissions of halogenated very short-lived substances (VSLS) are poorly constrained. However, their inclusion in global models is required to simulate a realistic inorganic bromine (Bry) loading in both the troposphere, where bromine chemistry perturbs global oxidizing capacity, and in the stratosphere, where it is a major sink for ozone (O3). We have performed simulations using a 3-D chemical transport model (CTM) including three top-down and a single bottom-up derived emission inventory of the major brominated VSLS bromoform (CHBr3) and dibromomethane (CH2Br2). We perform the first concerted evaluation of these inventories, comparing both the magnitude and spatial distribution of emissions. For a quantitative evaluation of each inventory, model output is compared with independent long-term observations at National Oceanic and Atmospheric Administration (NOAA) ground-based stations and with aircraft observations made during the NSF HIAPER Pole-to-Pole Observations (HIPPO) project. For CHBr3, the mean absolute deviation between model and surface observation ranges from 0.22 (38%) to 0.78 (115%) parts per trillion (ppt) in the tropics, depending on emission inventory. For CH2Br2, the range is 0.17 (24%) to 1.25 (167%) ppt. We also use aircraft observations made during the 2011 "Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere" (SHIVA) campaign, in the tropical West Pacific. Here, the performance of the various inventories also varies significantly, but overall the CTM is able to reproduce observed CHBr3 well in the free troposphere using an inventory based on observed sea-to-air fluxes. Finally, we identify the range of uncertainty associated with these VSLS emission inventories on stratospheric bromine loading due to VSLS (BryVSLS). Our simulations show BryVSLS ranges from ~ 4.0 to 8.0 ppt depending on the inventory. We report an optimised estimate at the lower end of this range (~ 4 ppt) based on combining the CHBr3 and CH2Br2 inventories which give best agreement with the compilation of observations in the tropics.
Emissions of halogenated very short-lived substances (VSLS) are poorly constrained. However, their inclusion in global models is required to simulate a realistic inorganic bromine (Bry) loading in both the troposphere, where bromine chemistry perturbs global oxidising capacity, and in the stratosphere, where it is a major sink for ozone (O3). We have performed simulations using a 3-D chemical transport model (CTM) including three top-down and a single bottom-up derived emission inventory of the major brominated VSLS bromoform (CHBr3) and dibromomethane (CH2Br2). We perform the first concerted evaluation of these inventories, comparing both the magnitude and spatial distribution of emissions. For a quantitative evaluation of each inventory, model output is compared with independent long-term observations at National Oceanic and Atmospheric Administration (NOAA) ground-based stations and with aircraft observations made during the NSF (National Science Foundation) HIAPER Pole-to-Pole Observations (HIPPO) project. For CHBr3, the mean absolute deviation between model and surface observation ranges from 0.22 (38%) to 0.78 (115%) parts per trillion (ppt) in the tropics, depending on emission inventory. For CH2Br2, the range is 0.17 (24%) to 1.25 (167%) ppt. We also use aircraft observations made during the 2011 Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere (SHIVA) campaign, in the tropical western Pacific. Here, the performance of the various inventories also varies significantly, but overall the CTM is able to reproduce observed CHBr3 well in the free troposphere using an inventory based on observed sea-to-air fluxes. Finally, we identify the range of uncertainty associated with these VSLS emission inventories on stratospheric bromine loading due to VSLS (BryVSLS). Our simulations show BryVSLS ranges from ~4.0 to 8.0 ppt depending on the inventory. We report an optimised estimate at the lower end of this range (~4 ppt) based on combining the CHBr3 and CH2Br2 inventories which give best agreement with the compilation of observations in the tropics.
The first concerted multi-model intercomparison of halogenated very short-lived substances (VSLS) has been performed, within the framework of the ongoing Atmospheric Tracer Transport Model Intercomparison Project (TransCom). Eleven global models or model variants participated, simulating the major natural bromine VSLS, bromoform (CHBr3) and dibromomethane (CH2Br2), over a 20-year period (1993-2012). The overarching goal of TransCom-VSLS was to provide a reconciled model estimate of the stratospheric source gas injection (SGI) of bromine from these gases, to constrain the current measurement-derived range, and to investigate inter-model differences
due to emissions and transport processes. Models ran with standardised idealised chemistry, to isolate differences due to transport, and we investigated the sensitivity of results to a range of VSLS emission inventories. Models were tested in their ability to reproduce the observed seasonal and spatial distribution of VSLS at the surface, using measurements from NOAA’s long-term global monitoring network, and in the tropical troposphere, using recent aircraft measurements - including high altitude observations from the NASA Global Hawk platform.
The models generally capture the seasonal cycle of surface CHBr3 and CH2Br2 well, with a strong model measurement correlation (r ≥0.7) and a low sensitivity to the choice of emission inventory, at most sites. In a given model, the absolute model-measurement agreement is highly sensitive to the choice of emissions and inter-model differences are also apparent, even when using the same inventory, highlighting the challenges faced in evaluating such inventories at the global scale. Across the ensemble, most consistency is found within the tropics where most of the models (8 out of 11) achieve optimal agreement to surface CHBr3 observations using the lowest of the three CHBr3 emission inventories tested (similarly, 8 out of 11 models for CH2Br2). In general, the models are able to reproduce well observations of CHBr3 and CH2Br2 obtained in the tropical tropopause layer (TTL) at various locations throughout the Pacific. Zonal variability in VSLS loading in the TTL is generally consistent among models, with CHBr3 (and to a lesser extent CH2Br2) most elevated over the tropical West Pacific during boreal winter. The models also indicate the Asian Monsoon during boreal summer to be an important pathway for VSLS reaching the stratosphere, though the strength of this signal varies considerably among models.
We derive an ensemble climatological mean estimate of the stratospheric bromine SGI from CHBr3 and CH2Br2 of 2.0 (1.2-2.5) ppt, ∼57% larger than the best estimate from the most re- cent World Meteorological Organization (WMO) Ozone Assessment Report. We find no evidence for a long-term, transport-driven trend in the stratospheric SGI of bromine over the simulation period. However, transport-driven inter-annual variability in the annual mean bromine SGI is of the order of a ±5%, with SGI exhibiting a strong positive correlation with ENSO in the East Pacific
The first concerted multi-model intercomparison of halogenated very short-lived substances (VSLS) has been performed, within the framework of the ongoing Atmospheric Tracer Transport Model Intercomparison Project (TransCom). Eleven global models or model variants participated (nine chemical transport models and two chemistry–climate models) by simulating the major natural bromine VSLS, bromoform (CHBr3) and dibromomethane (CH2Br2), over a 20-year period (1993–2012). Except for three model simulations, all others were driven offline by (or nudged to) reanalysed meteorology. The overarching goal of TransCom-VSLS was to provide a reconciled model estimate of the stratospheric source gas injection (SGI) of bromine from these gases, to constrain the current measurement-derived range, and to investigate inter-model differences due to emissions and transport processes. Models ran with standardised idealised chemistry, to isolate differences due to transport, and we investigated the sensitivity of results to a range of VSLS emission inventories. Models were tested in their ability to reproduce the observed seasonal and spatial distribution of VSLS at the surface, using measurements from NOAA's long-term global monitoring network, and in the tropical troposphere, using recent aircraft measurements – including high-altitude observations from the NASA Global Hawk platform.
The models generally capture the observed seasonal cycle of surface CHBr3 and CH2Br2 well, with a strong model–measurement correlation (r ≥ 0.7) at most sites. In a given model, the absolute model–measurement agreement at the surface is highly sensitive to the choice of emissions. Large inter-model differences are apparent when using the same emission inventory, highlighting the challenges faced in evaluating such inventories at the global scale. Across the ensemble, most consistency is found within the tropics where most of the models (8 out of 11) achieve best agreement to surface CHBr3 observations using the lowest of the three CHBr3 emission inventories tested (similarly, 8 out of 11 models for CH2Br2). In general, the models reproduce observations of CHBr3 and CH2Br2 obtained in the tropical tropopause layer (TTL) at various locations throughout the Pacific well. Zonal variability in VSLS loading in the TTL is generally consistent among models, with CHBr3 (and to a lesser extent CH2Br2) most elevated over the tropical western Pacific during boreal winter. The models also indicate the Asian monsoon during boreal summer to be an important pathway for VSLS reaching the stratosphere, though the strength of this signal varies considerably among models.
We derive an ensemble climatological mean estimate of the stratospheric bromine SGI from CHBr3 and CH2Br2 of 2.0 (1.2–2.5) ppt, ∼ 57 % larger than the best estimate from the most recent World Meteorological Organization (WMO) Ozone Assessment Report. We find no evidence for a long-term, transport-driven trend in the stratospheric SGI of bromine over the simulation period. The transport-driven interannual variability in the annual mean bromine SGI is of the order of ±5 %, with SGI exhibiting a strong positive correlation with the El Niño–Southern Oscillation (ENSO) in the eastern Pacific. Overall, our results do not show systematic differences between models specific to the choice of reanalysis meteorology, rather clear differences are seen related to differences in the implementation of transport processes in the models.
Precipitation extremes with devastating socioeconomic consequences within the South American Monsoon System (SAMS) are expected to become more frequent in the near future. The complexity in SAMS behavior, however, poses severe challenges for reliable future projections. Thus, robust paleomonsoon records are needed to constrain the high spatiotemporal variability in the response of SAMS rainfall to different climatic drivers. This study uses Ti/Ca ratios from X‐ray fluorescence scanning of a sediment core retrieved off eastern Brazilian to trace precipitation changes over the past 322 Kyr. The results indicate that despite the spatiotemporal complexity of the SAMS, insolation forcing is the primary pacemaker of variations in the monsoonal system. Additional modulation by atmospheric pCO2 suggests that SAMS intensity over eastern Brazil will be suppressed by rising CO2 emissions in the future. Lastly, our record reveals an unprecedented strong and persistent wet period during Marine Isotope Stage 6 driven by anomalously strong trade winds.
Video and image data are regularly used in the field of benthic ecology to document biodiversity. However, their use is subject to a number of challenges, principally the identification of taxa within the images without associated physical specimens. The challenge of applying traditional taxonomic keys to the identification of fauna from images has led to the development of personal, group, or institution level reference image catalogues of operational taxonomic units (OTUs) or morphospecies. Lack of standardisation among these reference catalogues has led to problems with observer bias and the inability to combine datasets across studies. In addition, lack of a common reference standard is stifling efforts in the application of artificial intelligence to taxon identification. Using the North Atlantic deep sea as a case study, we propose a database structure to facilitate standardisation of morphospecies image catalogues between research groups and support future use in multiple front-end applications. We also propose a framework for coordination of international efforts to develop reference guides for the identification of marine species from images. The proposed structure maps to the Darwin Core standard to allow integration with existing databases. We suggest a management framework where high-level taxonomic groups are curated by a regional team, consisting of both end users and taxonomic experts. We identify a mechanism by which overall quality of data within a common reference guide could be raised over the next decade. Finally, we discuss the role of a common reference standard in advancing marine ecology and supporting sustainable use of this ecosystem.
The growth of aerosol due to the aqueous phase oxidation of sulfur dioxide by ozone was measured in laboratory-generated clouds created in the Cosmics Leaving OUtdoor Droplets (CLOUD) chamber at the European Organization for Nuclear Research (CERN). Experiments were performed at 10 and −10 °C, on acidic (sulfuric acid) and on partially to fully neutralised (ammonium sulfate) seed aerosol. Clouds were generated by performing an adiabatic expansion – pressurising the chamber to 220 hPa above atmospheric pressure, and then rapidly releasing the excess pressure, resulting in a cooling, condensation of water on the aerosol and a cloud lifetime of approximately 6 min. A model was developed to compare the observed aerosol growth with that predicted using oxidation rate constants previously measured in bulk solutions. The model captured the measured aerosol growth very well for experiments performed at 10 and −10 °C, indicating that, in contrast to some previous studies, the oxidation rates of SO2 in a dispersed aqueous system can be well represented by using accepted rate constants, based on bulk measurements. To the best of our knowledge, these are the first laboratory-based measurements of aqueous phase oxidation in a dispersed, super-cooled population of droplets. The measurements are therefore important in confirming that the extrapolation of currently accepted reaction rate constants to temperatures below 0 °C is correct.