550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (673)
- 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 (833) (remove)
Has Fulltext
- yes (833)
Is part of the Bibliography
- no (833)
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 (617)
- 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)
Atmospheric new particle formation is a general phenomenon observed over coniferous forests. So far nucleation is either parameterised as a function of gaseous sulphuric acid concentration only, which is unable to explain the observed seasonality of nucleation events at different measurement sites, or as a function of sulphuric acid and organic molecules. Here we introduce different nucleation parameters based on the interaction of sulphuric acid and terpene oxidation products and elucidate the individual importance. They include basic trace gas and meteorological measurements such as ozone and water vapour concentrations, temperature (for terpene emission) and UV B radiation as a proxy for OH radical formation. We apply these new parameters to field studies conducted at conducted at Finnish and German measurement sites and compare these to nucleation observations on a daily and annual scale. General agreement was found, although the specific compounds responsible for the nucleation process remain speculative. This can be interpreted as follows: During cooler seasons the emission of biogenic terpenes and the OH availability limits the new particle formation while towards warmer seasons the ratio of ozone and water vapour concentration seems to dominate the general behaviour. Therefore, organics seem to support ambient nucleation besides sulphuric acid or an OH-related compound. Using these nucleation parameters to extrapolate the current conditions to prognosed future concentrations of ozone, water vapour and organic concentrations leads to a significant potential increase in the nucleation event number.
In this study we report the set-up of a novel twin chamber technique that uses the comparative method and establishes an appropriate connection of atmospheric and laboratory methods to broaden the tools for investigations. It is designed to study the impact of certain parameters and gases on ambient processes, such as particle formation online, and can be applied in a large variety of conditions. The characterisation of both chambers proved that both chambers operate identically, with a residence time xT (COMPASS1) = 26.5 ± 0.3 min and xT (COMPASS2) = 26.6 ± 0.4 min, at a typical flow rate of 15 L min−1 and a gas leak rate of (1.6 ± 0.8) × 10−5 s−1. Particle loss rates were found to be larger (due to the particles' stickiness to the chamber walls), with an extrapolated maximum of 1.8 × 10−3 s−1 at 1 nm, i.e. a hundredfold of the gas leak rate. This latter value is associated with sticky non-volatile gaseous compounds, too. Comparison measurement showed no significant differences. Therefore operation under atmospheric conditions is trustworthy. To indicate the applicability and the benefit of the system, a set of experiments was conducted under different conditions, i.e. urban and remote, enhanced ozone and terpenes as well as reduced sunlight. In order to do so, an ozone lamp was applied to enhance ozone in one of two chambers; the measurement chamber was protected from radiation by a first-aid cover and volatile organic compounds (VOCs) were added using a small additional flow and a temperature-controlled oven. During the elevated ozone period, ambient particle number and volume increased substantially at urban and remote conditions, but by a different intensity. Protection of solar radiation displayed a clear negative effect on particle number, while terpene addition did cause a distinct daily pattern. E.g. adding β pinene particle number concentration rose by 13% maximum at noontime, while no significant effect was observable during darkness. Therefore, the system is a useful tool for investigating local precursors and the details of ambient particle formation at surface locations as well as potential future feedback processes.
The anthropogenic influence on climate and environment has increased strongly since industrialization about 150 yr ago. The consequences for the atmosphere became more and more apparent and nowadays affect our life quality on Earth progressively. Because of that it is very important to understand the atmospheric processes, on which these effects are based on, in detail. In this study we report the set-up of a novel twin chamber technique that uses the comparative method and establishes an appropriate connection of atmospheric and laboratory methods to broaden the tools for investigations. It is designed to study the impact of certain parameters and gases on ambient processes such as particle formation online and can be applied in a large variety of conditions. The characterisation of both chambers proved that both chambers operate identically with a residence time (xT (COMPASS 1) = 26.5 ± 0.3 min and xT (COMPASS 2) = 26.6 ± 0.4 min) at a typical flow rate of 15 L min−1 and a deposition rate (1.6 ± 0.8) × 10−5 s−1. Comparison measurement showed no significant differences. Therefore operation under atmospheric conditions is trustworthy. To indicate the applicability and the benefit of the system a set of experiments was conducted at different conditions, i.e. urban and remote, enhancing ozone and terpenes as well as reducing sunlight. In the ozone enhanced ambient particle number and volume increased substantially at urban and remote conditions in a different strength. Solar radiation displayed a clear positive effect on particle number as well as terpene addition did at remote conditions. Therefore the system is a useful tool to investigate local precursors, the details of ambient particle formation at surface locations as well as future feedback processes.
Desert dust is one of the most abundant ice nucleating particle types in the atmosphere. Traditionally, clay minerals were assumed to determine the ice nucleation ability of desert dust and constituted the focus of ice nucleation studies over several decades. Recently some feldspar species were identified to be ice active at much higher temperatures than clay minerals, redirecting studies to investigate the contribution of feldspar to ice nucleation on desert dust. However, so far no study has shown the atmospheric relevance of this mineral phase.
For this study four dust samples were collected after airborne transport in the troposphere from the Sahara to different locations (Crete, the Peloponnese, Canary Islands, and the Sinai Peninsula). Additionally, 11 dust samples were collected from the surface from nine of the biggest deserts worldwide. The samples were used to study the ice nucleation behavior specific to different desert dusts. Furthermore, we investigated how representative surface-collected dust is for the atmosphere by comparing to the ice nucleation activity of the airborne samples. We used the IMCA-ZINC setup to form droplets on single aerosol particles which were subsequently exposed to temperatures between 233 and 250 K. Dust particles were collected in parallel on filters for offline cold-stage ice nucleation experiments at 253–263 K. To help the interpretation of the ice nucleation experiments the mineralogical composition of the dusts was investigated. We find that a higher ice nucleation activity in a given sample at 253 K can be attributed to the K-feldspar content present in this sample, whereas at temperatures between 238 and 245 K it is attributed to the sum of feldspar and quartz content present. A high clay content, in contrast, is associated with lower ice nucleation activity. This confirms the importance of feldspar above 250 K and the role of quartz and feldspars determining the ice nucleation activities at lower temperatures as found by earlier studies for monomineral dusts. The airborne samples show on average a lower ice nucleation activity than the surface-collected ones. Furthermore, we find that under certain conditions milling can lead to a decrease in the ice nucleation ability of polymineral samples due to the different hardness and cleavage of individual mineral phases causing an increase of minerals with low ice nucleation ability in the atmospherically relevant size fraction. Comparison of our data set to an existing desert dust parameterization confirms its applicability for climate models. Our results suggest that for an improved prediction of the ice nucleation ability of desert dust in the atmosphere, the modeling of emission and atmospheric transport of the feldspar and quartz mineral phases would be key, while other minerals are only of minor importance.
Desert dust is one of the most abundant ice nucleating particle types in the atmosphere. Tra ditionally, clay minerals were assumed to determine the ice nucleation ability of desert dust and constituted the focus of ice nucleation studies. Only recently some feldspar species were identified to be ice-active at much higher temperatures than clay minerals, redirecting studies to investigate the contribution of feldspar to ice nucleation on desert dust. However, so far no study has shown the atmospheric relevance of this mineral phase.
For this study four dust samples were collected after airborne transport in the troposphere from the Sahara to different locations (Crete, the Peloponnese, Canary Islands and the Sinai Peninsula). Additionally, eleven dust samples were collected from the surface from nine of the biggest deserts worldwide. The samples were used to study the ice nucleation behavior specific to different desert dusts. Furthermore we investigated how representative ice nucleation on surface-collected dust is for that in the atmosphere by comparing to the ice nucleation activity of the airborne samples. We used the IMCA-ZINC set-up to form droplets on single aerosol particles which were subsequently exposed to temperatures between 233 - 250 K. Dust particles were collected in parallel on filters for offline cold stage ice nucleation experiments at 253 - 263 K. To help the interpretation of the results from the ice nucleation experiments the mineralogical composition of the dusts was investigated.We found that a higher ice nucleation activity in a given sample can be attributed at 253 K to the K-feldspar content present in this sample whereas at temperatures between 238 - 245 K it is attributed to the sum of feldspar and quartz content present. A high clay content on the other hand is associated with a lower ice nucleation activity of a sample. This confirms the importance of feldspar at T > 250 K and the role of quartz and feldspars determining the ice nucleation activities at lower T as found by earlier studies for monomineral dust surrogates. Furthermore, we find that milling may lead to a decrease in the ice nucleation ability of polymineral samples due to a change in mineralogical composition in the atmospherically relevant size fraction arising from the different hardness and cleavage of individual mineral phases. Comparison of our comprehensive data set to an existing desert dust parameterization confirms its applicability for climate models. Our results suggest that for an improved prediction of the ice nucleation ability of desert dust in the atmosphere, the modelling of emission and atmospheric transport of the feldspar and quartz mineral phases would be key while other minerals are only of minor importance.
Processes occurring in the tropical upper troposphere and lower stratosphere (UT/LS) are of importance for the global climate, for the stratospheric dynamics and air chemistry, and they influence the global distribution of water vapour, trace gases and aerosols. The mechanisms underlying cloud formation and variability in the UT/LS are of scientific concern as these still are not adequately described and quantified by numerical models. Part of the reasons for this is the scarcity of detailed in-situ measurements in particular from the Tropical Transition Layer (TTL) within the UT/LS. In this contribution we provide measurements of particle number densities and the amounts of non-volatile particles in the submicron size range present in the UT/LS over Southern Brazil, West Africa, and Northern Australia. The data were collected in-situ on board of the Russian high altitude research aircraft M-55 "Geophysica" using the specialised COPAS (COndensation PArticle counting System) instrument during the TROCCINOX (Araçatuba, Brazil, February 2005), the SCOUT-O3 (Darwin, Australia, December 2005), and SCOUT-AMMA (Ouagadougou, Burkina Faso, August 2006) campaigns. The vertical profiles obtained are compared to those from previous measurements from the NASA DC-8 and NASA WB-57F over Costa Rica and other tropical locations between 1999 and 2007. The number density of the submicron particles as function of altitude was found to be remarkably constant (even back to 1987) over the tropical UT/LS altitude band such that a parameterisation suitable for models can be extracted from the measurements. At altitudes corresponding to potential temperatures above 430 K a slight increase of the number densities from 2005/2006 results from the data in comparison to the 1987 to 2007 measurements. The origins of this increase are unknown. By contrast the data from Northern hemispheric mid latitudes do not exhibit such an increase between 1999 and 2006. Vertical profiles of the non-volatile fraction of the submicron particles were also measured by a COPAS channel and are presented here. The resulting profiles of the non-volatile number density fraction show a pronounced maximum of 50% in the tropical TTL over Australia and West Africa. Below and above this fraction is much lower attaining values of 10% and smaller. In the lower stratosphere the fine particles mostly consist of sulphuric acid which is reflected in the low numbers of non-volatile residues measured by COPAS. Without detailed chemical composition measurements the reason for the increase of non-volatile particle fractions cannot yet be given. The long distance transfer flights to Brazil, Australia and West-Africa were executed during a time window of 17 months within a period of relative volcanic quiescence. Thus the data measured during these transfers represent a "snapshot picture" documenting the status of a significant part of the global UT/LS aerosol (with sizes below 1 μm) at low concentration levels 15 years after the last major (i.e., the 1991 Mount Pinatubo) eruption. The corresponding latitudinal distributions of the measured particle number densities are also presented in this paper in order to provide input on the UT/LS background aerosol for modelling purposes.
Processes occurring in the tropical upper troposphere (UT), the Tropical Transition Layer (TTL), and the lower stratosphere (LS) are of importance for the global climate, for stratospheric dynamics and air chemistry, and for their influence on the global distribution of water vapour, trace gases and aerosols. In this contribution we present aerosol and trace gas (in-situ) measurements from the tropical UT/LS over Southern Brazil, Northern Australia, and West Africa. The instruments were operated on board of the Russian high altitude research aircraft M-55 "Geophysica" and the DLR Falcon-20 during the campaigns TROCCINOX (Araçatuba, Brazil, February 2005), SCOUT-O3 (Darwin, Australia, December 2005), and SCOUT-AMMA (Ouagadougou, Burkina Faso, August 2006). The data cover submicron particle number densities and volatility from the COndensation PArticle counting System (COPAS), as well as relevant trace gases like N2O, ozone, and CO. We use these trace gas measurements to place the aerosol data into a broader atmospheric context. Also a juxtaposition of the submicron particle data with previous measurements over Costa Rica and other tropical locations between 1999 and 2007 (NASA DC-8 and NASA WB-57F) is provided. The submicron particle number densities, as a function of altitude, were found to be remarkably constant in the tropical UT/LS altitude band for the two decades after 1987. Thus, a parameterisation suitable for models can be extracted from these measurements. Compared to the average levels in the period between 1987 and 2007 a slight increase of particle abundances was found for 2005/2006 at altitudes with potential temperatures, theta, above 430 K. The origins of this increase are unknown except for increases measured during SCOUT-AMMA. Here the eruption of the Soufrière Hills volcano in the Caribbean caused elevated particle mixing ratios. The vertical profiles from Northern hemispheric mid-latitudes between 1999 and 2006 also are compact enough to derive a parameterisation. The tropical profiles all show a broad maximum of particle mixing ratios (between theta ~ 340 K and 390 K) which extends from below the TTL to above the thermal tropopause. Thus these particles are a "reservoir" for vertical transport into the stratosphere. The ratio of non-volatile particle number density to total particle number density was also measured by COPAS. The vertical profiles of this ratio have a maximum of 50% above 370 K over Australia and West Africa and a pronounced minimum directly below. Without detailed chemical composition measurements a reason for the increase of non-volatile particle fractions cannot yet be given. However, half of the particles from the tropical "reservoir" contain compounds other than sulphuric acid and water. Correlations of the measured aerosol mixing ratios with N2O and ozone exhibit compact relationships for the tropical data from SCOUT-AMMA, TROCCINOX, and SCOUT-O3. Correlations with CO are more scattered probably because of the connection to different pollution source regions. We provide additional data from the long distance transfer flights to the campaign sites in Brazil, Australia, and West-Africa. These were executed during a time window of 17 months within a period of relative volcanic quiescence. Thus the data represent a "snapshot picture" documenting the status of a significant part of the global UT/LS fine aerosol at low concentration levels 15 years after the last major (i.e., the 1991 Mount Pinatubo) eruption. The corresponding latitudinal distributions of the measured particle number densities are presented in this paper to provide data of the UT/LS background aerosol for modelling purposes.
Opportunities and challenges for paleoaltimetry in "small" orogens: insights from the European Alps
(2020)
Many stable isotope paleoaltimetry studies have focused on paleoelevation reconstructions of orogenic plateaus such as the Tibetan or Andean Plateaus. We address the opportunities and challenges of applying stable isotope paleoaltimetry to “smaller” orogens. We do this using a high‐resolution isotope tracking general circulation model (ECHAM5‐wiso) and explore the precipitation δ18O (δ18Op) signal of Cenozoic paleoclimate and topographic change in the European Alps. Results predict a maximum δ18Op change of 4–5‰ (relative to present day) during topographic development of the Alps. This signal of topographic change has the same magnitude as changes in δ18Op values resulting from Pliocene and Last Glacial Maximum global climatic change. Despite the similar magnitude of the isotopic signals resulting from topographic and paleoclimate changes, their spatial patterns across central Europe differ. Our results suggest that an integration of paleoclimate modeling, multiproxy approaches, and low‐elevation reference proxy records distal from an orogen improve topographic reconstructions.
Sesquiterpenes (C15H24) are semi-volatile organic compounds emitted by vegetation and are of interest in atmospheric research because they influence the oxidative capacity of the atmosphere and contribute to the formation of secondary organic aerosols. However, little is known about their emission pattern and no established parameterisation is available for global emission models. The aim of this study is to investigate a Central European spruce forest and its emission response to meteorological and environmental parameters, looking for a parameterisation that incorporates heat and oxidative stress as the main driving forces of the induced emissions. Therefore, a healthy ca. 80 yr old Norway spruce (Picea abies) tree was selected and a dynamical vegetation enclosure technique was applied from April to November 2011. The emissions clearly responded to temperature changes with small variations in the β-factor along the year (βspring = 0.09 ± 0.01, βsummer = 0.12 ± 0.02, βautumn = 0.11 ± 0.02). However, daily calculated values revealed a vast amount of variability in temperature dependencies ((0.02 ± 0.002) < β < (0.27 ± 0.04)) with no distinct seasonality.
By separating the complete dataset in 10 different ozone regimes, we found that in moderately or less polluted atmospheric conditions the main driving force of sesquiterpene emissions is the temperature, but when ambient ozone mixing ratios exceed a critical threshold of (36.6 ± 3.9) ppbv, the emissions become primarily correlated with ozone. Considering the complete dataset, cross correlation analysis resulted in highest correlation with ambient ozone mixing ratios (CCO3 = 0.63 ± 0.01; CCT = 0.47 ± 0.02 at t = 0 h for temperature) with a time shift 2–4 h prior to the emissions. An only temperature dependent algorithm was found to substantially underestimate the induced emissions (20% of the measured; R2 = 0.31). However, the addition of an ozone dependent term improved substantially the fitting between measured and modelled emissions (81% of the modelled emissions could be explained by the measurements; R2 = 0.63), providing confidence about the reliability of the suggested parameterisation for the spruce forest site investigated.
Sesquiterpenes (C15H24) are semi-volatile organic compounds emitted by vegetation and are of interest in atmospheric research because they influence the oxidative capacity of the atmosphere and contribute to the formation of secondary organic aerosols. However, little is known about their emission pattern and no established parameterization is available for global emission models. The aim of this study is to investigate a Central European spruce forest and its emission response to meteorological and environmental parameters, looking for a parameterization that incorporates heat and oxidative stress as the main driving forces of the induced emissions. Therefore, a healthy ca. 80 yr old Norway spruce (Picea abies) tree was selected and a dynamical vegetation enclosure technique was applied from April to November 2011. The emissions clearly responded to temperature changes with small variations in the β-factor along the year (βspring=0.09 ± 0.01, βsummer=0.12 ± 0.02, βautumn=0.11 ± 0.02). However, daily calculated values revealed a vast amount of variability in temperature dependencies ((0.02 ± 0.002)< β<(0.27 ± 0.04)) with no distinct seasonality.
By separating the complete dataset in 10 different ozone regimes, we found that in moderately or less polluted atmospheric conditions the main driving force of sesquiterpene emissions is the temperature, but when ambient ozone mixing ratios exceed a~critical threshold of (36.6 ± 3.9) ppbv, the emissions become primarily correlated with ozone. Considering the complete dataset, cross correlation analysis resulted in highest correlation with ambient ozone mixing ratios (CCO3=0.63 ± 0.01; CCT=0.47 ± 0.02 at t=0 h for temperature) with a time shift 2–4 h prior to the emissions. An only temperature dependent algorithm was found to substantially underestimate the induced emissions (20 % of the measured; R2=0.31). However, the addition of an ozone dependent term improved substantially the fitting between measured and modeled emissions (81 % of the measured; R2=0.63), providing confidence about the reliability of the suggested parameterization for the spruce forest site investigated.
In the original manuscript, Figs. 7–16 included fonts which were not correctly embedded in the file. As such, unless certain propriety software (ArcGIS) is installed on the viewing platform, the figures will appear corrupted. In this Corrigendum, Figs. 7–16 and their captions are reproduced with the fonts correctly embedded. Please find the correct figures below.
The late Miocene palaeorecord provides evidence for a warmer and wetter climate than that of today, and there is uncertainty in the palaeo-CO2 record of at least 200 ppm. We present results from fully coupled atmosphere-ocean-vegetation simulations for the late Miocene that examine the relative roles of palaeogeography (topography and ice sheet geometry) and CO2 concentration in the determination of late Miocene climate through comprehensive terrestrial model-data comparisons. Assuming that these data accurately reflect the late Miocene climate, and that the late Miocene palaeogeographic reconstruction used in the model is robust, then results indicate that:
1. Both palaeogeography and atmospheric CO2 contribute to the proxy-derived precipitation differences between the late Miocene and modern reference climates. However these contributions exibit synergy and so do not add linearly.
2. The vast majority of the proxy-derived temperature differences between the late Miocene and modern reference climates can only be accounted for if we assume a palaeo-CO2 concentration towards the higher end of the range of estimates.
The Late Miocene (∼11.6–5.3 Ma) palaeorecord provides evidence for a warmer and wetter climate than that of today and there is uncertainty in the palaeo-CO2 record of at least 150 ppmv. We present results from fully coupled atmosphere-ocean-vegetation simulations for the Late Miocene that examine the relative roles of palaeogeography (topography and ice sheet geometry) and CO2 concentration in the determination of Late Miocene climate through comprehensive terrestrial model-data comparisons. Assuming that the data accurately reflects the Late Miocene climate, and that the Late Miocene palaeogeographic reconstruction used in the model is robust, then results indicate that the proxy-derived precipitation differences between the Late Miocene and modern can be largely accounted for by the palaeogeographic changes alone. However, the proxy-derived temperatures differences between the Late Miocene and modern can only begin to be accounted for if we assume a palaeo-CO2 concentration towards the higher end of the range of estimates.
Responses of southern ocean seafloor habitats and communities to global and local drivers of change
(2021)
Knowledge of life on the Southern Ocean seafloor has substantially grown since the beginning of this century with increasing ship-based surveys and regular monitoring sites, new technologies and greatly enhanced data sharing. However, seafloor habitats and their communities exhibit high spatial variability and heterogeneity that challenges the way in which we assess the state of the Southern Ocean benthos on larger scales. The Antarctic shelf is rich in diversity compared with deeper water areas, important for storing carbon (“blue carbon”) and provides habitat for commercial fish species. In this paper, we focus on the seafloor habitats of the Antarctic shelf, which are vulnerable to drivers of change including increasing ocean temperatures, iceberg scour, sea ice melt, ocean acidification, fishing pressures, pollution and non-indigenous species. Some of the most vulnerable areas include the West Antarctic Peninsula, which is experiencing rapid regional warming and increased iceberg-scouring, subantarctic islands and tourist destinations where human activities and environmental conditions increase the potential for the establishment of non-indigenous species and active fishing areas around South Georgia, Heard and MacDonald Islands. Vulnerable species include those in areas of regional warming with low thermal tolerance, calcifying species susceptible to increasing ocean acidity as well as slow-growing habitat-forming species that can be damaged by fishing gears e.g., sponges, bryozoan, and coral species. Management regimes can protect seafloor habitats and key species from fishing activities; some areas will need more protection than others, accounting for specific traits that make species vulnerable, slow growing and long-lived species, restricted locations with optimum physiological conditions and available food, and restricted distributions of rare species. Ecosystem-based management practices and long-term, highly protected areas may be the most effective tools in the preservation of vulnerable seafloor habitats. Here, we focus on outlining seafloor responses to drivers of change observed to date and projections for the future. We discuss the need for action to preserve seafloor habitats under climate change, fishing pressures and other anthropogenic impacts.
We have developed and characterized the novel PTR3, a proton transfer reaction-time-of-flight mass spectrometer (PTR-TOF) using a new gas inlet and an innovative reaction chamber design. The reaction chamber consists of a tripole operated with rf voltages generating an electric field only in the radial direction. An elevated electrical field is necessary to reduce clustering of primary hydronium (H3O+) and product ions with water molecules present in the sample gas. The axial movement of the ions is achieved by the sample gas flow only. Therefore, the new design allows a 30-fold longer reaction time and a 40-fold increase in pressure compared to standard PTR-TOF-MS. First calibration tests show sensitivities of up to 18000 counts per second/parts per billion and volume (cps/ppbv) at a mass resolution of >8000 m/Δm (fwhm). The new inlet using center-sampling through a critical orifice reduces wall losses of low volatility compounds. Therefore, the new PTR3 instrument is sensitive to VOC typically present in the ppbv range as well as to semivolatile organic compounds (SVOC) and even highly oxidized organic molecules (HOMs) present in the parts per quadrillion per volume (ppqv) range in the atmosphere.
We conducted measurements of up to the five important short-lived brominated species in the marine boundary layer (MBL) of the mid-latitudes (List/Sylt, North Sea) in June 2009 and of the tropical Western Pacific during the TransBrom ship campaign in October 2009. For the one-week time series in List mean mixing ratios of 2.0, 1.1, 0.2, 0.1 ppt were analysed for CHBr3, CH2Br2, CHBr2Cl and CH2BrCl, with maxima of 5.8 and 1.6 ppt for the two main components CHBr3 and CH2Br2. Along the cruise track in the Western Pacific (between 41° N and 13° S) mean mixing ratios of 1.0, 0.9, 0.2, 0.1 and 0.1 ppt for CHBr3, CH2Br2, CHBrCl2, CHBr2Cl and CH2BrCl were determined. Air samples with coastal influence showed considerably higher mixing ratios than the samples with open ocean origin. Correlation analyses of the two datasets yielded strong linear relationships between the mixing ratios of four of the five species (except for CH2BrCl). Using a combined dataset from the two campaigns, rough estimates of the molar emission ratios between the correlated substances were derived as follows: 9/1/0.3/0.3 for CHBr3/CH2Br2/CHBrCl2/CHBr2Cl. Additional measurements were made in the tropical tropopause layer (TTL) above Teresina (Brazil, 5.07° S, 42.87° W) in June 2008, using balloon-borne cryogenic whole air sampling technique. Near the level of zero clear-sky net radiative heating (LZRH) at 14.8 km about 2.25 ppt organic bromine was bound to the five short-lived species, making up 13 % of total organic bromine (17.82 ppt). CH2Br2 (1.45 ppt) and CHBr3 (0.56 ppt) accounted for 90 % of the budget of short-lived compounds in that region. Near the tropopause (at 17.5 km) organic bromine from short-lived substances was reduced to 1.35 ppt, with 1.07 ppt and 0.12 ppt attributed to CH2Br2 and CHBr3 respectively.
We conducted measurements of the five important short-lived organic bromine species in the marine boundary layer (MBL). Measurements were made in the Northern Hemisphere mid-latitudes (Sylt Island, North Sea) in June 2009 and in the tropical Western Pacific during the TransBrom ship campaign in October 2009. For the one-week time series on Sylt Island, mean mixing ratios of CHBr3, CH2Br2, CHBr2Cl and CH2BrCl were 2.0, 1.1, 0.2, 0.1 ppt, respectively. We found maxima of 5.8 and 1.6 ppt for the two main components CHBr3 and CH2Br2. Along the cruise track in the Western Pacific (between 41° N and 13° S) we measured mean mixing ratios of 0.9, 0.9, 0.2, 0.1 and 0.1 ppt for CHBr3, CH2Br2, CHBrCl2, CHBr2Cl and CH2BrCl. Air samples with coastal influence showed considerably higher mixing ratios than the samples with open ocean origin. Correlation analyses of the two data sets yielded strong linear relationships between the mixing ratios of four of the five species (except for CH2BrCl). Using a combined data set from the two campaigns and a comparison with the results from two former studies, rough estimates of the molar emission ratios between the correlated substances were: 9/1/0.35/0.35 for CHBr3/CH2Br2/CHBrCl2/CHBr2Cl. Additional measurements were made in the tropical tropopause layer (TTL) above Teresina (Brazil, 5° S) in June 2008, using balloon-borne cryogenic whole air sampling technique. Near the level of zero clear-sky net radiative heating (LZRH) at 14.8 km about 2.25 ppt organic bromine was bound to the five short-lived species, making up 13% of total organic bromine (17.82 ppt). CH2Br2 (1.45 ppt) and CHBr3 (0.56 ppt) accounted for 90% of the budget of short-lived compounds in that region. Near the tropopause (at 17.5 km) organic bromine from these substances was reduced to 1.35 ppt, with 1.07 and 0.12 ppt attributed to CH2Br2 and CHBr3, respectively.
A satellite-based climate record of monthly mean surface solar irradiance (SIS) is investigated with regard to possible inhomogeneities in time. The data record is provided by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Climate Monitoring (CM SAF) for the period of 1983 to 2005, covering a disk area between ±70° in latitude and longitude. The Standard Normal Homogeneity Test (SNHT) and two other homogeneity tests are applied with and without the use of reference SIS data (from the Baseline Surface Radiation Network (BSRN) and from the ECMWF (European Centre for Medium-Range Weather Forecasts) ERA -Interim reanalysis. The focus is on the detection of break-like inhomogeneities, which may occur due to satellite or SIS retrieval algorithm changes. In comparison with the few suitable BSRN SIS observation series with limited extension in time (no data before 1992), the CM SAF SIS time series do not show significant inhomogeneities, even though slight discrepancies in the surface measurements appear. The investigation of the full CM SAF SIS domain reveal inhomogeneities related to most of the documented satellite and retrieval changes, but only for relatively small domain fractions (especially in mountainous desert-like areas in Africa). In these regions the retrieval algorithm is not capable of adjusting for the changes of the satellite instruments. For other areas, e.g., Europe, no such breaks in the time series are found. We conclude that the CM SAF SIS data record has to be further assessed and regionally homogenized before climate trend investigations can be conducted.
Lightning climate change projections show large uncertainties caused by limited empirical knowledge and strong assumptions inherent to coarse-grid climate modeling. This study addresses the latter issue by implementing and applying the lightning potential index parameterization (LPI) into a fine-grid convection-permitting regional climate model (CPM). This setup takes advantage of the explicit representation of deep convection in CPMs and allows for process-oriented LPI inputs such as vertical velocity within convective cells and coexistence of microphysical hydrometeor types, which are known to contribute to charge separation mechanisms. The LPI output is compared to output from a simpler flash rate parameterization, namely the CAPE × PREC parameterization, applied in a non-CPM on a coarser grid. The LPI’s implementation into the regional climate model COSMO-CLM successfully reproduces the observed lightning climatology, including its latitudinal gradient, its daily and hourly probability distributions, and its diurnal and annual cycles. Besides, the simulated temperature dependence of lightning reflects the observed dependency. The LPI outperforms the CAPE × PREC parameterization in all applied diagnostics. Based on this satisfactory evaluation, we used the LPI to a climate change projection under the RCP8.5 scenario. For the domain under investigation centered over Germany, the LPI projects a decrease of 4.8% in flash rate by the end of the century, in opposition to a projected increase of 17.4% as projected using the CAPE × PREC parameterization. The future decrease of LPI occurs mostly during the summer afternoons and is related to (i) a change in convection occurrence and (ii) changes in the microphysical mixing. The two parameterizations differ because of different convection occurrences in the CPM and non-CPM and because of changes in the microphysical mixing, which is only represented in the LPI lightning parameterization.
Convection-permitting models (CPMs) have proven their usefulness in representing precipitation on a sub-daily scale. However, investigations on sub-hourly scales are still lacking, even though these are the scales for which showers exhibit the most variability. A Lagrangian approach is implemented here to evaluate the representation of showers in a CPM, using the limited-area climate model COSMO-CLM. This approach consists of tracking 5‑min precipitation fields to retrieve different features of showers (e.g., temporal pattern, horizontal speed, lifetime). In total, 312 cases are simulated at a resolution of 0.01 ° over Central Germany, and among these cases, 78 are evaluated against a radar dataset. The model is able to represent most observed features for different types of convective cells. In addition, the CPM reproduced well the observed relationship between the precipitation characteristics and temperature indicating that the COSMO-CLM model is sophisticated enough to represent the climatological features of showers.
It is common practice to use a 30-year period to derive climatological values, as recommended by the World Meteorological Organization. However this convention relies on important assumptions, of which the validity can be examined by deriving the uncertainty inherent to using a limited time-period for deriving climatological values. In this study a new method, aiming at deriving this uncertainty, has been developed with an application to precipitation for a station in Europe (Westdorpe) and one in Africa (Gulu). The weather generator framework is used to produce synthetic daily precipitation time-series that can also be regarded as alternative climate realizations. The framework consists of an improved Markov model, which shows good performance in reproducing the 5-day precipitation variability. The sub-seasonal, seasonal and the inter-annual signals are introduced in the weather generator framework by including covariates. These covariates are derived from an empirical mode decomposition analysis with an improved stability and significance assessment. Introducing covariates was found to substantially improve the monthly precipitation variability for Gulu. From the weather generator, 1,000 synthetic time-series were produced. The divergence between these time-series demonstrates an uncertainty, inherent to using a 30-year period for mean precipitation, of 11 % for Westdorpe and 15 % for Gulu. The uncertainty for precipitation 10-year return levels was found to be 37 % for both sites.
Convection-permitting climate model are promising tools for improved representation of extremes, but the number of regions for which these models have been evaluated are still rather limited to make robust conclusions. In addition, an integrated interpretation of near-surface characteristics (typically temperature and precipitation) together with cloud properties is limited. The objective of this paper is to comprehensively evaluate the performance of a ‘state-of-the-art’ regional convection-permitting climate model for a mid-latitude coastal region with little orographic forcing. For this purpose, an 11-year integration with the COSMO-CLM model at Convection-Permitting Scale (CPS) using a grid spacing of 2.8 km was compared with in-situ and satellite-based observations of precipitation, temperature, cloud properties and radiation (both at the surface and the top of the atmosphere). CPS clearly improves the representation of precipitation, in especially the diurnal cycle, intensity and spatial distribution of hourly precipitation. Improvements in the representation of temperature are less obvious. In fact the CPS integration overestimates both low and high temperature extremes. The underlying cause for the overestimation of high temperature extremes was attributed to deficiencies in the cloud properties: The modelled cloud fraction is only 46 % whereas a cloud fraction of 65 % was observed. Surprisingly, the effect of this deficiency was less pronounced at the radiation balance at the top of the atmosphere due to a compensating error, in particular an overestimation of the reflectivity of clouds when they are present. Overall, a better representation of convective precipitation and a very good representation of the daily cycle in different cloud types were demonstrated. However, to overcome remaining deficiencies, additional efforts are necessary to improve cloud characteristics in CPS. This will be a challenging task due to compensating deficiencies that currently exist in ‘state-of-the-art’ models, yielding a good representation of average climate conditions. In the light of using the CPS models to study climate change it is necessary that these deficiencies are addressed in future research.
In old and heavily weathered soils, the availability of P might be so small that the primary production of plants is limited. However, plants have evolved several mechanisms to actively take up P from the soil or mine it to overcome this limitation. These mechanisms involve the active uptake of P mediated by mycorrhiza, biotic de-occlusion through root clusters, and the biotic enhancement of weathering through root exudation. The objective of this paper is to investigate how and where these processes contribute to alleviate P limitation on primary productivity. To do so, we propose a process-based model accounting for the major processes of the carbon, water, and P cycle including chemical weathering at the global scale. We use simulation experiments to assess the relative importance of the different uptake mechanisms to alleviate P limitation on biomass production. Implementing P limitation on biomass synthesis allows the assessment of the efficiencies of biomass production across different ecosystems.
We find that active P-uptake is an essential mechanism for sustaining P availability on long time scales, whereas biotic de-occlusion might serve as a buffer on time scales shorter than 10 000 yr. Although active P uptake is essential for reducing P losses by leaching, humid lowland soils reach P limitation after around 100 000 yr of soil evolution. Given the generalized modeling framework, our model results compare reasonably with observed or independently estimated patterns and ranges of P concentrations in soils and vegetation. Furthermore, our simulations suggest that P limitation might be an important driver of biomass production efficiency (the fraction of the gross primary productivity used for biomass growth), and that vegetation on older soils becomes P-limited leading to a smaller biomass production efficiency.
With this study, we provide a theoretical basis for investigating the responses of terrestrial ecosystems to P availability linking geological and ecological time scales under different environmental settings.
In old and heavily weathered soils, the availability of P might be so small that the primary production of plants is limited. However, plants have evolved several mechanisms to actively take up P from the soil or mine it to overcome this limitation. These mechanisms involve the active uptake of P mediated by mycorrhiza, biotic de-occlusion through root clusters, and the biotic enhancement of weathering through root exudation. The objective of this paper is to investigate how and where these processes contribute to alleviate P limitation on primary productivity. To do so, we propose a process-based model accounting for the major processes of the carbon, water, and P cycles including chemical weathering at the global scale. Implementing P limitation on biomass synthesis allows the assessment of the efficiencies of biomass production across different ecosystems. We use simulation experiments to assess the relative importance of the different uptake mechanisms to alleviate P limitation on biomass production. We find that active P uptake is an essential mechanism for sustaining P availability on long timescales, whereas biotic de-occlusion might serve as a buffer on timescales shorter than 10 000 yr. Although active P uptake is essential for reducing P losses by leaching, humid lowland soils reach P limitation after around 100 000 yr of soil evolution. Given the generalized modelling framework, our model results compare reasonably with observed or independently estimated patterns and ranges of P concentrations in soils and vegetation. Furthermore, our simulations suggest that P limitation might be an important driver of biomass production efficiency (the fraction of the gross primary productivity used for biomass growth), and that vegetation on old soils has a smaller biomass production rate when P becomes limiting. With this study, we provide a theoretical basis for investigating the responses of terrestrial ecosystems to P availability linking geological and ecological timescales under different environmental settings.
The timing and duration of leaf deployment strongly regulate earth-atmosphere interactions and biotic processes. Leaf dynamics therefore have major implications for life on earth, including the global energy balance, carbon and water cycles, feedbacks to climate, species extinction risk and agriculture. Evidence of shifts in the timing of leaf deployment and senescence (leaf phenology) as a result of climate change has been accumulating over the past decades, particularly in relation to spring phenology in the northern hemisphere. However, leaf phenological change in other parts of the world has received less attention. This thesis quantifies global phenological change over the past three decades using remotely sensed data. Phenological change was found to be widespread and severe, also in the southern hemisphere. While the detected change testifies of the phenological plasticity of many plant species, it is not clear if the duration of leaf deployment (leaf habit) is equally sensitive to environmental change. Since evergreen and deciduous leaf habits are often distinctly sorted along environmental gradients, ecologists have hypothesised that these patterns result from natural selection for an optimal leaf habit, under a given environmental regime. Such evolutionary convergence can be examined by testing if the physiological niche that is occupied by a particular leaf habit (evergreen or deciduous) is similar among regions with distinct evolutionary histories. Using a process-based model of plant growth and a constructed map of evergreen and deciduous vegetation, the physiological niche of leaf habits was quantified in four global biogeographic realms. Substantial niche overlap was found between the same leaf habit in different realms, suggesting evolutionary convergence of the physiological niche. This implies a sensitivity of leaf habit to environmental change, as environmental variables determine the geographic space where the physiological niche allows a positive carbon balance, and therefore occurrence of the leaf habit. Since the physiological niche consists of the integrated effects of physiological traits and trade-offs, environmental dependencies and leaf habit and phenology, an understanding of the carbon economy of individual plants requires decomposing the physiological niche into its components. Using empirical data on leaf phenology, leaf habit and physiological processes from woody species in a seasonally dry African savanna, a simple carbon balance model was parametrised. Carbon gain varied considerably between species as a result of substantial variation in leaf habit, leaf phenology and physiological traits. The multiple lines of evidence in this thesis therefore suggest that, while convergent selective forces may determine the dominant leaf habit in a particular environment, inter-specific variation is substantial, potentially as a consequence of historical contingencies or competitive interactions.
Soil water potential is crucial to plant transpiration and thus to carbon cycling and biosphere–atmosphere interactions, yet it is difficult to measure in the field. Volumetric and gravimetric water contents are easy and cheap to measure in the field, but can be a poor proxy of plant-available water. Soil water content can be transformed to water potential using soil moisture retention curves. We provide empirically derived soil moisture retention curves for seven soil types in the Kruger National Park, South Africa. Site-specific curves produced excellent estimates of soil water potential from soil water content values. Curves from soils derived from the same geological substrate were similar, potentially allowing for the use of one curve for basalt soils and another for granite soils. It is anticipated that this dataset will help hydrologists and ecophysiologists understand water dynamics, carbon cycling and biosphere–atmosphere interactions under current and changing climatic conditions in the region.
Development of a Bioaerosol single particle detector (BIO IN) for the Fast Ice Nucleus CHamber FINCH
(2010)
In this work we present the setup and first tests of our new BIO IN detector. This detector was constructed to classify atmospheric ice nuclei (IN) for their biological content. It is designed to be coupled to the Fast Ice Nucleus CHamber FINCH. If one particle acts as an ice nucleus, it will be at least partly covered with ice at the end of the development section of the FINCH chamber. The device combines an auto-fluorescence detector and a circular depolarization detector for simultaneous detection of biological material and discrimination between water droplets, ice crystals and non activated large aerosol particles. The excitation of biological material with UV light and analysis of auto-fluorescence is a common principle used for flow cytometry, fluorescence microscopy, spectroscopy and imaging. The detection of auto-fluorescence of airborne single particles demands some more experimental effort. However, expensive commercial sensors are available for special purposes, e.g. size distribution measurements. But these sensors will not fit the specifications needed for the FINCH IN counter (e.g. high sample flow of up 10 LPM). The newly developed -low cost- BIO IN sensor uses a single high-power UV LED for the electronic excitation instead of much more expensive UV lasers. Other key advantages of the new sensor are the low weight, compact size, and the little effect on the aerosol sample, which allows it to be coupled with other instruments for further analysis. The instrument will be flown on one of the first missions of the new German research aircraft "HALO" (High Altitude and LOng range).
In this work we present the setup and first tests of our new BIO IN detector. This detector is designed to classify atmospheric ice nuclei (IN) for their biological content. Biological material is identified via its auto-fluorescence (intrinsic fluorescence) after irradiation with UV radiation. Ice nuclei are key substances for precipitation development via the Bergeron–Findeisen process. The level of scientific knowledge regarding origin and climatology (temporal and spatial distribution) of IN is very low. Some biological material is known to be active as IN even at relatively high temperatures of up to –2°C (e.g. pseudomonas syringae bacteria). These biological IN could have a strong influence on the formation of clouds and precipitation. We have designed the new BIO IN sensor to analyze the abundance of IN of biological origin. The instrument will be flown on one of the first missions of the new German research aircraft ''HALO'' (High Altitude and LOng Range).
This thesis is focusing on the impact of Paratethys and Mediterranean water bodies over the Eurasian climate and the interplay between climate, tectonics and biosphere during the late Miocene. This target was the interval between 12.7 and 7.65 Ma for Paratethys, following the Eastern Paratethys restriction and isolation, and 7.2−6.5 Ma (the early Messinian) in Mediterranean, zooming on the effects of gateway restrictions over the eastern Mediterranean and the new born Aegean domain. In both cases restriction is overlapping with large scale climatic changes and tectonic reconfiguration, leading a sort of symbiotic relationship.
Paratethys was a giant epicontinental sea that covered a large part of Eurasia since Paleogene. Due to the Eurasia-Afro-Arabia collision and formation of the Alpine-Himalayan belt (Rögl, 1999; Popov et al., 2006), the Paratethys was divided during the late Miocene in smaller basins that in time were isolated of each other. The protracted isolation and intense continentalisation of paratethyan realm led to changes in humidity distribution, basin connectivity, sediment sources and salinity. These changes had in turn major consequences over water circulation, water availability, vegetation cover and biota. These changes are more intense after 11.6 Ma, when the Eastern Paratethys lost any sustained marine connection, evolving into an enclosed system with endemic fauna (Harzhauser and Piller, 2007).
Mediterranean Sea is a Mezozoic oceanic relic squeezed between Africa, Europe, Anatolia and Arabia, as Africa continued to subduct beneath the European plate. As opposed to Paratethys, it maintained the open connection with the ocean until Messinian, when the two Atlantic gateways (Betic and Rifian corridors) closed for a short time, isolating the basin. The cut off resulted in a dramatic drop down and onset of evaporitic precipitation in marginal basins, the event receiving the name of Messinian Salinity Crisis (5.97−5.55). The restriction affected all marine ecosystems, due to changes in salinity and stratification of water column.
The main objectives of this thesis were:
(1) build valid paleo-temperatures records for both basins based on biomarkers;
(2) reconstruct the hydrology for the late Miocene time interval;
(3) identify vegetation composition and changes;
(4) identify paleo-fires in the late Miocene sediment records;
(5) identify the biotic response to the overall climate and tectonic changes.
All the above objectives were attained with results published in specific journals (Chapters 5−7).
Based on Panagia section (Taman Peninsula, Russia) the longest Paratethys temperature record was completed (~5 Myr), covering the interval between 12.7 and 7.65 Ma. A comprehensive SST and MAT records was obtained, as well as soil pH and carbon (δ13C) and hydrogen (δ2H) stable isotopic compositions on n-alkanes and alkenones. The main findings are concentrated around three prolonged periods with severe droughts affecting the late Miocene circum-Paratethys region peaking at 9.65, 9.4 and 7.9 Ma, associated with a transition towards open land vegetation, intensification of fire activity and enhanced evaporation and aridity.
The time intervals with dryer conditions recorded in Panagia coincide with periods of mammal turnover and dispersal in Eurasia indicating that major environmental changes occurred in the circum-Paratethys region and Paratethys fragmentation had a great impact on the terrestrial ecosystems, when periods of prolonged droughts generated biotic crises and animal displacements across the Eurasian continent. The δ13CC29n-alkane values and charcoal morphologies from Panagia indicate an increased contribution of C4 plants adapted to drier conditions at 9.66 Ma. Similarly high δ13CC29n-alkane values continue until 9.4 Ma, when in Western Europe increased seasonality accelerated the demise of the evergreen subtropical woodlands and expansion of grasslands from Anatolia and Middle East to Europe.
As a result of basin fragmentation and climatic stress, the Eastern Paratethys sub-basins progressively lost their marine properties and turned into brackish-fresh water bodies fed primarily by riverine input. The shallower areas became in time emerged, obstructing connections and isolating the biota, inducing rapid adjusting or extinctions. Thus, the Paratethys harbored a highly endemic fauna (Rögl, 1999), such as dwarf whales, dolphins, seals (among mammals), as well as fish and other taxa (mollusks, ostracods, diatoms, foraminifera, algae, etc.).
Collectively the data structured and analyzed in chapter five support a model in which the Eastern Paratethys evolved as a largely (en)closed system, registering paleoenvironmental signals that are governed by interbasinal connectivity (or lack of it) and regional climate changes in the basin catchment. Acting as an important source of humidity for Western and Central Asia, the size and areal extent of the Paratethys water body is likely to have had a major impact on hydroclimate patterns in the Eurasian interior, with the cumulative fluctuations in both hydrology and surface temperature enhancing the aridity and seasonality, with different partition of moisture over the year. Our combined data suggests a decoupling of Paratethys from the global system as isolation advanced, dominated by regional tectonics and ultimately the Paratethys volume and areal extent reduction.
In this paper we present evidence that the observed increase in tropical upwelling after the year 2000 may be attributed to a change in the Brewer-Dobson circulation pattern. For this purpose, we use the concept of transit times derived from residual circulation trajectories and different in-situ measurements of ozone and nitrous dioxide. Observations from the Canadian midlatitude ozone profile record, probability density functions of in-situ N2O observations and a shift of the N2O-O3 correlation slopes, taken together, indicate that the increased upwelling in the tropics after the year 2000 appears to have triggered an intensification of tracer transport from the tropics into the extratropics in the lower stratosphere below about 500 K. This finding is corroborated by the fact that transit times along the shallow branch of the residual circulation into the LMS have decreased for the same time period (1993–2003). On a longer time scale (1979–2009), the transit time of the shallow residual circulation branch show a steady decrease of about −1 month/decade over the last 30 years, while the transit times of the deep branch remain unchanged. This highlights the fact that a change in the upwelling across the tropical tropopause is not a direct indicator for changes of the whole Brewer-Dobson circulation.
In this paper we present evidence that the observed increase in tropical upwelling after the year 2000 may be attributed to a change in the Brewer-Dobson circulation pattern. For this purpose, we use the concept of transit times derived from residual circulation trajectories and different in-situ measurements of ozone and nitrous dioxide. Observations from the Canadian midlatitude ozone profile record, probability density functions of in-situ N2O observations and a shift of the N2O-O3 correlation slopes, taken together, indicate that the increased upwelling in the tropics after the year 2000 appears to have triggered an intensification of tracer transport from the tropics into the extratropics in the lower stratosphere below about 500 K. This finding is corroborated by the fact that transit times along the shallow branch of the residual circulation into the LMS have decreased for the same time period (1993–2003). On a longer time scale (1979–2009), the transit time of the shallow residual circulation branch show a steady decrease of about −1 month/decade over the last 30 yr, while the transit times of the deep branch remain unchanged. This highlights that changes in the upwelling across the tropical tropopause are not sufficient as an indicator for changes in the entire Brewer-Dobson circulation.
The seasonality of transport and mixing of air into the lowermost stratosphere (LMS) is studied using distributions of mean age of air and a~mass balance approach, based on in-situ observations of SF6 and CO2 during the SPURT (Spurenstofftransport in der Tropopausenregion, trace gas transport in the tropopause region) aircraft campaigns. Combining the information of the mean age of air and the water vapour distributions we demonstrate that the tropospheric air transported into the LMS above the extratropical tropopause layer (ExTL) originates predominantly from the tropical tropopause layer (TTL). The concept of our mass balance is based on simultaneous measurements of the two passive tracers and the assumption that transport into the LMS can be described by age spectra which are superposition of two different modes. Based on this concept we conclude that the stratospheric influence on LMS composition is strongest in April with tropospheric fractions (α1) below 20% and that the strongest tropospheric signatures are found in October with (α1 greater than 80%. Beyond the fractions, our mass balance concept allows to calculate the associated transit times for transport of tropospheric air from the tropics into the LMS. The shortest transit times (<0.3 years) are derived for the summer, continuously increasing up to 0.8 years by the end of spring. These findings suggest that strong quasi-horizontal mixing across the weak subtropical jet from summer to mid of autumn and the considerably shorter residual transport time-scales within the lower branch of the Brewer-Dobson circulation in summer than in winter dominates the tropospheric influence in the LMS until the beginning of next year's summer.
The seasonality of transport and mixing of air into the lowermost stratosphere (LMS) is studied using distributions of mean age of air and a mass balance approach, based on in-situ observations of SF6 and CO2 during the SPURT (Spurenstofftransport in der Tropopausenregion, trace gas transport in the tropopause region) aircraft campaigns. Combining the information of the mean age of air and the water vapour distributions we demonstrate that the tropospheric air transported into the LMS above the extratropical tropopause layer (ExTL) originates predominantly from the tropical tropopause layer (TTL). The concept of our mass balance is based on simultaneous measurements of the two passive tracers and the assumption that transport into the LMS can be described by age spectra which are superposition of two different modes. Based on this concept we conclude that the stratospheric influence on LMS composition is strongest in April with extreme values of the tropospheric fractions (alpha1) below 20% and that the strongest tropospheric signatures are found in October with alpha1 greater than 80%. Beyond the fractions, our mass balance concept allows us to calculate the associated transit times for transport of tropospheric air from the tropics into the LMS. The shortest transit times (<0.3 years) are derived for the summer, continuously increasing up to 0.8 years by the end of spring. These findings suggest that strong quasi-horizontal mixing across the weak subtropical jet from summer to mid of autumn and the considerably shorter residual transport time-scales within the lower branch of the Brewer-Dobson circulation in summer than in winter dominates the tropospheric influence in the LMS until the beginning of next year's summer.
The complete elastic stiffness tensor of thiourea has been determined from thermal diffuse scattering (TDS) using high-energy photons (100 keV). Comparison with earlier data confirms a very good agreement of the tensor coefficients. In contrast with established methods to obtain elastic stiffness coefficients (e.g. Brillouin spectroscopy, inelastic X-ray or neutron scattering, ultrasound spectroscopy), their determination from TDS is faster, does not require large samples or intricate sample preparation, and is applicable to opaque crystals. Using high-energy photons extends the applicability of the TDS-based approach to organic compounds which would suffer from radiation damage at lower photon energies.
During the APE-THESEO mission in the Indian Ocean the Myasishchev Design Bureau stratospheric research aircraft M55 Geophysica performed a flight over and within the inner core region of tropical cyclone Davina. Measurements of total water, water vapour, temperature, aerosol backscattering, ozone and tracers were made and are discussed here in comparison with the averages of those quantities acquired during the campaign time frame.
Temperature anomalies in the tropical tropopause layer (TTL), warmer than average in the lower part and colder than average in the upper TTL were observed. Ozone was strongly reduced compared to its average value, and thick cirrus decks were present up to the cold point, sometimes topped by a layer of very dry air. Evidence for meridional transport of trace gases in the stratosphere above the cyclone and the perturbed water distribution in the TTL is illustrated and discussed.
During the APE-THESEO mission in the Indian Ocean the Myasishchev Design Bureau stratospheric research aircraft M55 Geophysica performed a flight over and within the inner core region of tropical cyclone Davina. Measurements of total water, water vapour, temperature, aerosol backscattering, ozone and tracers were made and are discussed here in comparison with the averages of those quantities acquired during the campaign time frame. Temperature anomalies in the tropical tropopause layer (TTL), warmer than average in the lower part and colder than average in the upper TTL were observed. Ozone was strongly reduced compared to its average value, and thick cirrus decks were present up to the cold point, sometimes topped by a layer of very dry air. Evidence for meridional transport of trace gases in the stratosphere above the cyclone was observed and perturbed water distribution in the TTL was documented. The paper discuss possible processes of dehydration induced by the cirrus forming above the cyclone, and change in the chemical tracer and water distribution in the lower stratosphere 400–430 K due to meridional transport from the mid-latitudes and link with Davina. Moreover it compares the data prior and after the cyclone passage to discuss its actual impact on the atmospheric chemistry and thermodynamics.
Environmental change impacts on the C- and N-cycle of European forests: a model comparison study
(2013)
Forests are important components of the greenhouse gas balance of Europe. There is considerable uncertainty about how predicted changes to climate and nitrogen deposition will perturb the carbon and nitrogen cycles of European forests and thereby alter forest growth, carbon sequestration and N2O emission. The present study aimed to quantify the carbon and nitrogen balance, including the exchange of greenhouse gases, of European forests over the period 2010–2030, with a particular emphasis on the spatial variability of change. The analysis was carried out for two tree species: European beech and Scots pine. For this purpose, four different dynamic models were used: BASFOR, DailyDayCent, INTEGRATOR and Landscape-DNDC. These models span a range from semi-empirical to complex mechanistic. Comparison of these models allowed assessment of the extent to which model predictions depended on differences in model inputs and structure. We found a European average carbon sink of 0.160 ± 0.020 kgC m−2 yr−1 (pine) and 0.138 ± 0.062 kgC m−2 yr−1 (beech) and N2O source of 0.285 ± 0.125 kgN ha−1 yr−1 (pine) and 0.575 ± 0.105 kgN ha−1 yr−1 (beech). The European average greenhouse gas potential of the carbon sink was 18 (pine) and 8 (beech) times that of the N2O source. Carbon sequestration was larger in the trees than in the soil. Carbon sequestration and forest growth were largest in central Europe and lowest in northern Sweden and Finland, N. Poland and S. Spain. No single driver was found to dominate change across Europe. Forests were found to be most sensitive to change in environmental drivers where the drivers were limiting growth, where changes were particularly large or where changes acted in concert. The models disagreed as to which environmental changes were most significant for the geographical variation in forest growth and as to which tree species showed the largest rate of carbon sequestration. Pine and beech forests were found to have differing sensitivities to environmental change, in particular the response to changes in nitrogen and precipitation, with beech forest more vulnerable to drought. There was considerable uncertainty about the geographical location of N2O emissions. Two of the models BASFOR and LandscapeDNDC had largest emissions in central Europe where nitrogen deposition and soil nitrogen were largest, whereas the two other models identified different regions with large N2O emission. N2O emissions were found to be larger from beech than pine forests and were found to be particularly sensitive to forest growth.
Forests are important components of the greenhouse gas balance of Europe. There is considerable uncertainty about how predicted changes to climate and nitrogen deposition will perturb the carbon and nitrogen cycles of European forests and thereby alter forest growth, carbon sequestration and N2O emission. The present study aimed to quantify the carbon and nitrogen balance, including the exchange of greenhouse gases, of European forests over the period 2010–2030, with a particular emphasis on the spatial variability of change. The analysis was carried out for two tree species: European beech and Scots pine. For this purpose, four different dynamic models were used: BASFOR, DailyDayCent, INTEGRATOR and Landscape-DNDC. These models span a range from semi-empirical to complex mechanistic. Comparison of these models allowed assessment of the extent to which model predictions depended on differences in model inputs and structure. We found a European average carbon sink of 0.160 ± 0.020 kgC m−2 yr−1 (pine) and 0.138 ± 0.062 kgC m−2 yr−1 (beech) and N2O source of 0.285 ± 0.125 kgN ha−1 yr−1 (pine) and 0.575 ± 0.105 kgN ha−1 yr−1 (beech). The European average greenhouse gas potential of the carbon source was 18 (pine) and 8 (beech) times that of the N2O source. Carbon sequestration was larger in the trees than in the soil. Carbon sequestration and forest growth were largest in central Europe and lowest in northern Sweden and Finland, N. Poland and S. Spain. No single driver was found to dominate change across Europe. Forests were found to be most sensitive to change in environmental drivers where the drivers were limiting growth, where changes were particularly large or where changes acted in concert. The models disagreed as to which environmental changes were most significant for the geographical variation in forest growth and as to which tree species showed the largest rate of carbon sequestration. Pine and beech forests were found to have differing sensitivities to environmental change, in particular the response to changes in nitrogen and precipitation, with beech forest more vulnerable to drought. There was considerable uncertainty about the geographical location of N2O emissions. Two of the models BASFOR and LandscapeDNDC had largest emissions in central Europe where nitrogen deposition and soil nitrogen were largest whereas the two other models identified different regions with large N2O emission. N2O emissions were found to be larger from beech than pine forests and were found to be particularly sensitive to forest growth.
Atmospheric particles play an important role in the radiative balance of the Earth, as well as they affect human health and air quality. Hence, the chemical characterization constitutes a crucial task to determinate their properties, sources and fate. Particularly, the analysis of nanoparticles (d<100 nm) represents an analytical challenge, since these particles are abundant in number but have very little mass.
This accumulative thesis focuses on the chemical characterization of nanoparticles, performed in both laboratory and field studies. Here, I present four manuscripts, two of which are my main project as a lead author.
The first manuscript (Caudillo et al., 2021) focuses on the gas and the particle phase originated from biogenic precursor gases (α-pinene and isoprene). The experiments were performed in the CLOUD chamber at CERN to simulate pure biogenic new particle formation. Both gas and particle phases are measured with a nitrate CI-APi-TOF mass spectrometer, while the TD-DMA is coupled to it for particle-phase measurements, this setup allows a direct comparison as both measurements use the identical chemical ionization and detector. This study demonstrates the suitability of the TD-DMA for measuring newly formed nanoparticles and it confirms that isoprene suppresses new particle formation but contributes to the growth of newly formed particles.
The second manuscript (Caudillo et al., 2022) presents an intercomparison of four different techniques (including the TD-DMA) for measuring the chemical composition of SOA nanoparticles. The measurements were conducted in the CLOUD chamber. The intercomparison was done by contrasting the observed chemical composition, the calculated volatility, and the thermal desorption behavior (for the thermal desorption techniques). The methods generally agreed on the most important compounds that are found in the nanoparticles. However, they did see different parts of the organic spectrum. Potential explanations for these differences are suggested.
The third manuscript (Ungeheuer al., 2022) presents both laboratory and ambient measurements to investigate the ability of lubricant oil to form new particles. These new particles are an important source of ultrafine particles in the areas nearby large airports. The ambient measurements were performed downwind of Frankfurt International Airport, and it was found that the fraction of lubricant oil is largest in the smallest particles. In the laboratory, the main finding was that evaporated lubricant oil nucleates and forms new particles rapidly. The results suggest that nucleation of lubricant oil and subsequent particle growth can occur in the cooling exhaust plumes of aircraft-turbofans.
The fourth manuscript (Wang et al., 2022) is a new particle formation study in the CLOUD chamber at CERN. This study shows that nitric acid, sulfuric acid, and ammonia interact synergistically and rapidly form particles under upper free tropospheric conditions. These particles can grow by condensation (driven by the availability of ammonia) up to CCN sizes and INP particles. The ability of these particles to act as a CCN and INP was also investigated and it was found to be as efficient as for desert dust. This mechanism constitutes an important finding and it can account for previous observations of high concentrations of ammonia and ammonium nitrate over the Asia monsoon region.
Biogenic organic precursors play an important role in atmospheric new particle formation (NPF). One of the major precursor species is α-pinene, which upon oxidation can form a suite of products covering a wide range of volatilities. Highly oxygenated organic molecules (HOMs) comprise a fraction of the oxidation products formed. While it is known that HOMs contribute to secondary organic aerosol (SOA) formation, including NPF, they have not been well studied in newly formed particles due to their very low mass concentrations. Here we present gas- and particle-phase chemical composition data from experimental studies of α-pinene oxidation, including in the presence of isoprene, at temperatures (−50 and −30 ∘C) and relative humidities (20 % and 60 %) relevant in the upper free troposphere. The measurements took place at the CERN Cosmics Leaving Outdoor Droplets (CLOUD) chamber. The particle chemical composition was analyzed by a thermal desorption differential mobility analyzer (TD-DMA) coupled to a nitrate chemical ionization–atmospheric pressure interface–time-of-flight (CI-APi-TOF) mass spectrometer. CI-APi-TOF was used for particle- and gas-phase measurements, applying the same ionization and detection scheme. Our measurements revealed the presence of C8−10 monomers and C18−20 dimers as the major compounds in the particles (diameter up to ∼ 100 nm). Particularly, for the system with isoprene added, C5 (C5H10O5−7) and C15 compounds (C15H24O5−10) were detected. This observation is consistent with the previously observed formation of such compounds in the gas phase. However, although the C5 and C15 compounds do not easily nucleate, our measurements indicate that they can still contribute to the particle growth at free tropospheric conditions. For the experiments reported here, most likely isoprene oxidation products enhance the growth of particles larger than 15 nm. Additionally, we report on the nucleation rates measured at 1.7 nm (J1.7 nm) and compared with previous studies, we found lower J1.7 nm values, very likely due to the higher α-pinene and ozone mixing ratios used in the present study.
Biogenic organic precursors play an important role in atmospheric new particle formation (NPF). One of the major precursor species is α-pinene, which upon oxidation can form a suite of products covering a wide range of volatilities. Highly oxygenated organic molecules (HOMs) comprise a fraction of the oxidation products formed. While it is known that HOMs contribute to secondary organic aerosol (SOA) formation, including NPF, they have not been well studied in newly formed particles due to their very low mass concentrations. Here we present gas- and particle-phase chemical composition data from experimental studies of α-pinene oxidation, including in the presence of isoprene, at temperatures (−50 and −30 ∘C) and relative humidities (20 % and 60 %) relevant in the upper free troposphere. The measurements took place at the CERN Cosmics Leaving Outdoor Droplets (CLOUD) chamber. The particle chemical composition was analyzed by a thermal desorption differential mobility analyzer (TD-DMA) coupled to a nitrate chemical ionization–atmospheric pressure interface–time-of-flight (CI-APi-TOF) mass spectrometer. CI-APi-TOF was used for particle- and gas-phase measurements, applying the same ionization and detection scheme. Our measurements revealed the presence of C8−10 monomers and C18−20 dimers as the major compounds in the particles (diameter up to ∼ 100 nm). Particularly, for the system with isoprene added, C5 (C5H10O5−7) and C15 compounds (C15H24O5−10) were detected. This observation is consistent with the previously observed formation of such compounds in the gas phase. However, although the C5 and C15 compounds do not easily nucleate, our measurements indicate that they can still contribute to the particle growth at free tropospheric conditions. For the experiments reported here, most likely isoprene oxidation products enhance the growth of particles larger than 15 nm. Additionally, we report on the nucleation rates measured at 1.7 nm (J1.7 nm) and compared with previous studies, we found lower J1.7 nm values, very likely due to the higher α-pinene and ozone mixing ratios used in the present study.
The Weissert Event ~133 million years ago marked a profound global cooling that punctuated the Early Cretaceous greenhouse. We present modelling, high-resolution bulk organic carbon isotopes and chronostratigraphically calibrated sea surface temperature (SSTs) based on an organic paleothermometer (the TEX86 proxy), which capture the Weissert Event in the semi-enclosed Weddell Sea basin, offshore Antarctica (paleolatitude ~54 °S; paleowater depth ~500 meters). We document a ~3–4 °C drop in SST coinciding with the Weissert cold end, and converge the Weddell Sea data, climate simulations and available worldwide multi-proxy based temperature data towards one unifying solution providing a best-fit between all lines of evidence. The outcome confirms a 3.0 °C ( ±1.7 °C) global mean surface cooling across the Weissert Event, which translates into a ~40% drop in atmospheric pCO2 over a period of ~700 thousand years. Consistent with geologic evidence, this pCO2 drop favoured the potential build-up of local polar ice.
PolarCAP – A deep learning approach for first motion polarity classification of earthquake waveforms
(2022)
Highlights
• We present PolarCAP, a deep learning model that can classify the polarity of a waveform with a 98% accuracy.
• The first-motion polarity of seismograms is a useful parameter, but its manual determination can be laborious and imprecise.
• We demonstrate that in several cases the model can assign trace polar-ity more accurately than a human analyst.
Abstract
The polarity of first P-wave arrivals plays a significant role in the effective determination of focal mechanisms specially for smaller earthquakes. Manual estimation of polarities is not only time-consuming but also prone to human errors. This warrants a need for an automated algorithm for first motion polarity determination. We present a deep learning model - PolarCAP that uses an autoencoder architecture to identify first-motion polarities of earth-quake waveforms. PolarCAP is trained in a supervised fashion using more than 130,000 labelled traces from the Italian seismic dataset (INSTANCE) and is cross-validated on 22,000 traces to choose the most optimal set of hyperparameters. We obtain an accuracy of 0.98 on a completely unseen test dataset of almost 33,000 traces. Furthermore, we check the model generalizability by testing it on the datasets provided by previous works and show that our model achieves a higher recall on both positive and negative polarities.
Current atmospheric models do not include secondary organic aerosol (SOA) production from gas-phase reactions of polycyclic aromatic hydrocarbons (PAHs). Recent studies have shown that primary semivolatile emissions, previously assumed to be inert, undergo oxidation in the gas phase, leading to SOA formation. This opens the possibility that low-volatility gas-phase precursors are a potentially large source of SOA. In this work, SOA formation from gas-phase photooxidation of naphthalene, 1-methylnaphthalene (1-MN), 2-methylnaphthalene (2-MN), and 1,2-dimethylnaphthalene (1,2-DMN) is studied in the Caltech dual 28-m3 chambers. Under high-NOx conditions and aerosol mass loadings between 10 and 40 microg m-3, the SOA yields (mass of SOA per mass of hydrocarbon reacted) ranged from 0.19 to 0.30 for naphthalene, 0.19 to 0.39 for 1-MN, 0.26 to 0.45 for 2-MN, and constant at 0.31 for 1,2-DMN. Under low-NOx conditions, the SOA yields were measured to be 0.73, 0.68, and 0.58, for naphthalene, 1-MN, and 2-MN, respectively. The SOA was observed to be semivolatile under high-NOx conditions and essentially nonvolatile under low-NOx conditions, owing to the higher fraction of ring-retaining products formed under low-NOx conditions. When applying these measured yields to estimate SOA formation from primary emissions of diesel engines and wood burning, PAHs are estimated to yield 3–5 times more SOA than light aromatic compounds. PAHs can also account for up to 54% of the total SOA from oxidation of diesel emissions, representing a potentially large source of urban SOA.
Current atmospheric models do not include secondary organic aerosol (SOA) production from gas-phase reactions of polycyclic aromatic hydrocarbons (PAHs). Recent studies have shown that primary semivolatile emissions, previously assumed to be inert, undergo oxidation in the gas phase, leading to SOA formation. This opens the possibility that low-volatility gas-phase precursors are a potentially large source of SOA. In this work, SOA formation from gas-phase photooxidation of naphthalene, 1-methylnaphthalene (1-MN), 2-methylnaphthalene (2-MN), and 1,2-dimethylnaphthalene (1,2-DMN) is studied in the Caltech dual 28-m3 chambers. Under high-NOx conditions and aerosol mass loadings between 10 and 40 μg m, the SOA yields (mass of SOA per mass of hydrocarbon reacted) ranged from 0.19 to 0.30 for naphthalene, 0.19 to 0.39 for 1-MN, 0.26 to 0.45 for 2-MN, and constant at 0.31 for 1,2-DMN. Under low-NOx conditions, the SOA yields were measured to be 0.73, 0.68, and 0.58, for naphthalene, 1-MN, and 2-MN, respectively. The SOA was observed to be semivolatile under high-NOx conditions and essentially nonvolatile under low-NOx conditions, owing to the higher fraction of ring-retaining products formed under low-NOx conditions. When applying these measured yields to estimate SOA formation from primary emissions of diesel engines and wood burning, PAHs are estimated to yield 3–5 times more SOA than light aromatic compounds. PAHs can also account for up to 54% of the total SOA from oxidation of diesel emissions, representing a potentially large source of urban SOA.
Highlights
• Constrictional structures range from dome-and-basin folds to coeval folds and boudins.
• Under bulk constriction, the competent layer rotates slower than a passive plane.
• Extension-parallel and –perpendicular folds grow simultaneously.
• Extension-perpendicular folds affect previous boudins.
Abstract
We conducted scaled analogue modelling to show the influence of varying single layer initial orientation on the geometry of folds and boudins in a bulk constrictional strain field. The initial angle between the plane of shortening and the competent layer (θZ(i)) was incrementally increased from 0° to 90° by multiples of 11.25°. While the amount of layer thickening decreased with increasing θZ(i), the deformation structures produced range from pure dome-and-basin folds to coeval folds and boudins. Based on the attitude of fold axes, there are extension-parallel (FEPR) and extension-perpendicular (FEPP) folds, with axes subparallel and subperpendicular to the principal stretching axis (X), respectively. Coeval growth of FEPR folds and boudins occurred when θZ(i) > ca. 25°. The FEPP folds can be subdivided into a first type which affect the entire layer (if θZ(i) ranges between 11.25 and 78.75°) and a second type, referred to as FBEPP folds, which are affecting pre-existing boudins if θZ(i) > 45°. The interlimb angle of all types of folds increases with increasing θZ(i). Folds and boudins similar to the ones produced in this study can be found in salt domes and in tectonites of subduction zones.
We report on HCFC-22 data acquired by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) in the reduced spectral resolution nominal observation mode. The data cover the period from January 2005 to April 2012 and the altitude range from the upper troposphere (above cloud top altitude) to about 50 km. The profile retrieval was performed by constrained nonlinear least squares fitting of modelled spectra to the measured limb spectral radiances. The spectral ν4-band at 816.5 ± 13 cm−1 was used for the retrieval. A Tikhonov-type smoothing constraint was applied to stabilise the retrieval. In the lower stratosphere, we find a global volume mixing ratio of HCFC-22 of about 185 pptv in January 2005. The rate of linear growth in the lower latitudes lower stratosphere was about 6 to 7 pptv year−1 in the period 2005–2012. The profiles obtained were compared with ACE-FTS satellite data v3.5, as well as with MkIV balloon profiles and cryosampler balloon measurements. Between 13 and 22 km, average agreement within −3 to +5 pptv (MIPAS – ACE) with ACE-FTS v3.5 profiles is demonstrated. Agreement with MkIV solar occultation balloon-borne measurements is within 10–20 pptv below 30 km and worse above, while in situ cryosampler balloon measurements are systematically lower over their full altitude range by 15–50 pptv below 24 km and less than 10 pptv above 28 km. MIPAS HCFC-22 time series below 10 km altitude are shown to agree mostly well to corresponding time series of near-surface abundances from the NOAA/ESRL and AGAGE networks, although a more pronounced seasonal cycle is obvious in the satellite data. This is attributed to tropopause altitude fluctuations and subsidence of polar winter stratospheric air into the troposphere. A parametric model consisting of constant, linear, quasi-biennial oscillation (QBO) and several sine and cosine terms with different periods has been fitted to the temporal variation of stratospheric HCFC-22 for all 10°-latitude/1-to-2-km-altitude bins. The relative linear variation was always positive, with relative increases of 40–70 % decade−1 in the tropics and global lower stratosphere, and up to 120 % decade−1 in the upper stratosphere of the northern polar region and the southern extratropical hemisphere. Asian HCFC-22 emissions have become the major source of global upper tropospheric HCFC-22. In the upper troposphere, monsoon air, rich in HCFC-22, is instantaneously mixed into the tropics. In the middle stratosphere, between 20 and 30 km, the observed trend is inconsistent with the trend at the surface (corrected for the age of stratospheric air), hinting at circulation changes. There exists a stronger positive trend in HCFC-22 in the Southern Hemisphere and a more muted positive trend in the Northern Hemisphere, implying a potential change in the stratospheric circulation over the observation period.
We report on HCFC-22 data acquired by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) in reduced spectral resolution nominal mode in the period from January 2005 to April 2012 from version 5.02 level-1b spectral data and covering an altitude range from the upper troposphere (above cloud top altitude) to about 50 km. The profile retrieval was performed by constrained nonlinear least squares fitting of measured limb spectral radiances to modelled spectra. The spectral v4-band at 816.5 ± 13 cm-1 was used for the retrieval. A Tikhonov-type smoothing constraint was applied to stabilise the retrieval. In the lower stratosphere, we find a global volume mixing ratio of HCFC-22 of about 185 pptv in January 2005. The linear growth rate in the lower latitudes lower stratosphere was about 6 to 7 pptv yr-1 in the period 2005–2012. The obtained profiles were compared with ACE-FTS satellite data v3.5, as well as with MkIV balloon profiles and in situ cryosampler balloon measurements. Between 13 and 22 km, average agreement within -3 to +5 pptv (MIPAS–ACE) with ACE-FTS v3.5 pro files is demonstrated. Agreement with MkIV solar occultation balloon-borne measurements is within 10–20 pptv below 30 km and worse above, while in situ cryosampler balloon measurements are systematically lower over their full altitude range by 15– 50 pptv below 24 km and less than 10 pptv above 28 km. Obtained MIPAS HCFC-22 time series below 10 km altitude are shown to agree mostly well to corresponding time series of near-surface abundances from NOAA/ESRL and AGAGE networks, although a more pronounced seasonal cycle is obvious in the satellite data, probably due to tropopause altitude fluctuations and subsidence of polar winter stratospheric air into the troposphere. A parametric model consisting of constant, linear, quasi-biennial oscillation (QBO) and several sine and cosine terms with different periods has been fitted to the temporal variation of stratospheric HCFC-22 for all 10° latitude/1 to 2 km altitude bins. The relative linear variation was always positive, with relative increases of 40–70%decade-1 in the tropics and global lower stratosphere, and up to 120%decade-1 in the upper stratosphere of the northern polar region and the southern extratropical hemisphere. In the middle stratosphere between 20 and 30 km, the observed trend is not consistent with the age of stratospheric air-corrected trend at ground, but stronger positive at the Southern Hemisphere and less strong increasing in the Northern Hemisphere, hinting towards changes in the stratospheric circulation over the observation period.
First-principles modeling techniques offer the ability to simulate a wide range of systems under different physical conditions, such as temperature, pressure, and composition, without relying on empirical knowledge. Density functional theory (DFT), a quantum mechanical method, has become an exceptionally successful framework for materials science modeling. Employing DFT makes it possible to gain valuable insights into the fundamental state of a system, enabling the reliable determination of equilibrium crystal structures. Over time, DFT has become an essential tool that can be incorporated into various schemes for predicting the properties of a material related to its structure, insulating/metallic behavior, magnetism, and optics. DFT is regularly applied in numerous fields, spanning from fundamental subjects in condensed matter physics to the study of large-scale phenomena in geosciences. In the latter, the effectiveness of DFT stems from its ability to simulate the properties found on the Earth, other planets, and meteorites, which may pose challenges for their direct study or laboratory investigation.
In this thesis, a comprehensive examination of a family of monosulfides and a perovskite heterostructure was conducted. These materials are relevant for their potential applications in technology, energy harvesting, and in the case of monosulfides, their speculated abundance on the planet Mercury.
Firstly, a DFT approach was used to analyze two non-magnetic monosulfides, CaS and MgS. We determined their structural properties and then focused on the modeling of their reflectivity in the infrared region. The calculation of the reflectivity considered both harmonic and anharmonic contributions. In the harmonic limit, the non-analytic correction was employed to accurately determine the LO/TO splitting, which is necessary to delimit the retstrahlend band, that is, the maximum of the reflectivity. The anharmonic effects given by up to three-phonon and isotopic scatterings, which were included using perturbation theory, primarily smeared the reflectivity spectra edges in the high-wave region.
Secondly, four polymorphs of MnS were studied using a combination of first-principles methods to simulate their antiferromagnetic (AFM) and paramagnetic (PM) states. The integration of DFT+$U$ with special quasirandom structures (SQS) supercells, and occupation matrix control techniques was crucial for achieving convergence, structural optimization accuracy, and obtaining finite energy band gaps and local magnetic moments in the PM phases. The addition of the Hubbard $U$ correction was necessary to treat the highly-correlated Mn $d$-electrons. The success of our approach was clear based on our electronic structure predictions for the PM rock-salt B1-MnS polymorph. Experimentally this phase has been observed to be an insulator, but multiple \emph{ab initio} works resulted previously in metallic behavior. Our computations, on the other hand, predicted insulating and magnetic properties that compare well with available measurements. Additionally, the pressure-field stability of the four MnS polymorphs was studied. In the case of the PM phases, B1-MnS was identified to be the most stable up to about 21 GPa, then transforming into the B31-MnS polymorph. This finding was in close agreement with high-pressure experiments reporting a similar phase transformation. The optical properties of B1-, B4-, and B31-MnS were also simulated. The SQS technique was used to obtain soft-mode-free phonon band structures within the harmonic approximation. Then, the anharmonic effects were included, and the reflectivity was calculated for B1-MnS and B4-MnS. In both cases, a good agreement for the LO/TO splitting was achieved in comparison to experimental results.
Lastly, the oxygen-deficient heterostructure of LaAlO$_{3-\delta}$ /SrTiO$_{3-\delta}$ was investigated also employing DFT+$U$, with a particular emphasis on the potential impact of vacancy clustering at the interface. Six distinct configurations of pairs of vacancies were studied and their energies were compared to find the most stable one. The orbital reconstruction of Ti orbitals was also examined based on their location with respect to the vacancies and the local magnetic moments were calculated. The final results showed that linearly arranged vacancies located opposite to Ti ions give the most energetically stable configuration.
Small-scale thermal upwellings under the northern East African Rift from S travel time tomography
(2016)
There is a long-standing debate over how many and what types of plumes underlie the East African Rift and whether they do or do not drive its extension and consequent magmatism and seismicity. Here we present a new tomographic study of relative teleseismic S and SKS residuals that expands the resolution from previous regional studies below the northern East African Rift to image structure from the surface to the base of the transition zone. The images reveal two low-velocity clusters, below Afar and west of the Main Ethiopian Rift, that extend throughout the upper mantle and comprise several smaller-scale (about 100 km diameter), low-velocity features. These structures support those of our recent P tomographic study below the region. The relative magnitude of S to P residuals is around 3.5, which is consistent with a predominantly thermal nature of the anomalies. The S and P velocity anomalies in the low-velocity clusters can be explained by similar excess temperatures in the range of 100–200°C, consistent with temperatures inferred from other seismic, geochemical, and petrological studies. Somewhat stronger VS anomalies below Afar than west of the Main Ethiopian Rift may include an expression of volatiles and/or melt in this region. These results, together with a comparison with previous larger-scale tomographic models, indicate that these structures are likely small-scale upwellings with mild excess temperatures, rising from a regional thermal boundary layer at the base of the upper mantle.
Metasomatic reaction zones between mafic and ultramafic rocks exhumed from subduction zones provide a window into mass-transfer processes at high pressure. However, accurate interpretation of the rock record requires distinguishing high-pressure metasomatic processes from inherited oceanic signatures prior to subduction. We integrated constraints from bulk-rock geochemical compositions and petrophysical properties, mineral chemistry, and thermodynamic modeling to understand the formation of reaction zones between juxtaposed metagabbro and serpentinite as exemplified by the Voltri Massif (Ligurian Alps, Italy). Distinct zones of variably metasomatized metagabbro are dominated by chlorite, amphibole, clinopyroxene, epidote, rutile, ilmenite, and titanite between serpentinite and eclogitic metagabbro. Whereas the precursor serpentinite and oxide gabbro formed and were likely already in contact in an oceanic setting, the reaction zones formed by diffusional Mg-metasomatism between the two rocks from prograde to peak, to retrograde conditions in a subduction zone. Metasomatism of mafic rocks by Mg-rich fluids that previously equilibrated with serpentinite could be widespread along the subduction interface, within the subducted slab, and the mantle wedge. Furthermore, the models predict that talc formation by Si-metasomatism of serpentinite in subduction zones is limited by pressure-dependent increase in the silica activity buffered by the serpentine-talc equilibrium. Elevated activities of aqueous Ca and Al species would also favor the formation of chlorite and garnet. Accordingly, unusual conditions or processes would be required to stabilize abundant talc at high P-T conditions. Alternatively, a different set of mineral assemblages, such as serpentine- or chlorite-rich rocks, may be controlling the coupling-decoupling transition of the plate interface.
Talc formation via silica-metasomatism of ultramafic rocks is believed to play key roles in subduction zone processes. Yet, the conditions of talc formation remain poorly constrained. We used thermodynamic reaction-path models to assess the formation of talc at the slab-mantle interface and show that it is restricted to a limited set of pressure–temperature conditions, protolith, and fluid compositions. In contrast, our models predict that chlorite formation is ubiquitous at conditions relevant to the slab-mantle interface of subduction zones. The scarcity of talc and abundance of chlorite is evident in the rock record of exhumed subduction zone terranes. Talc formation during Si-metasomatism may thus play a more limited role in volatile cycling, strain localization, and in controlling the decoupling-coupling transition of the plate interface. Conversely, the observed and predicted ubiquity of chlorite corroborates its prominent role in slab-mantle interface processes that previous studies attributed to talc.
Key Points:
Limited talc formation by Si-metasomatism of ultramafic rocks in subduction zones
Chlorite formation is likely pervasive at the slab-mantle interface
Preferential formation of chlorite has wide-ranging chemical and physical implications for subduction zone processes
Plain Language Summary: In subduction zones, talc can form during chemical reactions of mantle rocks with silica-enriched fluids at the interface between descending oceanic plates and the overriding mantle. Its formation and distribution in subduction zones are believed to affect the volatile budget, rheological properties, and the down-dip limit of the decoupling of the slab-mantle interface. Therefore, illuminating the conditions that facilitate talc formation at high pressure-temperature conditions is key in assessing its roles in fundamental subduction zone processes. Using thermodynamic reaction-path models, we show that the formation of talc at the slab-mantle interface is restricted to a limited set of environmental conditions, because its formation is highly sensitive to the compositions of the mantle rocks and reactant fluids. Contrary to common belief, talc is unlikely to form in high abundance in ultramafic rocks metasomatized by Si-rich slab-derived fluids. Rather, our models predict the ubiquitous formation of chlorite along with other silicate minerals during Si-metasomatism due to the competing effects from other dissolved components that favor their formation over talc. This study calls into question the importance of talc during Si-metasomatism in subduction zones but highlights the more predominant role of chlorite.
The mechanisms of transfer of crustal material from the subducting slab to the overlying mantle wedge are still debated. Mélange rocks, formed by mixing of sediments, oceanic crust, and ultramafics along the slab-mantle interface, are predicted to ascend as diapirs from the slab-top and transfer their compositional signatures to the source region of arc magmas. However, the compositions of melts that result from the interaction of mélanges with a peridotite wedge remain unknown. Here we present experimental evidence that melting of peridotite hybridized by mélanges produces melts that carry the major and trace element abundances observed in natural arc magmas. We propose that differences in nature and relative contributions of mélanges hybridizing the mantle produce a range of primary arc magmas, from tholeiitic to calc-alkaline. Thus, assimilation of mélanges into the wedge may play a key role in transferring subduction signatures from the slab to the source of arc magmas.
Terrestrial climate and ecosystem evolution during ‘Greenhouse Earth’ phases of the early Paleogene remain incompletely known. Particularly, paleobotanical records from high southern latitudes are giving only limited insights into the Paleocene and early Eocene vegetation of the region. Hence, data from continuous well-calibrated sequences are required to make progress with the reconstruction of terrestrial climate and ecosystem dynamics from the southern latitudes during the early Paleogene.
In order to elucidate the terrestrial conditions from the high southern latitudes during the early Paleogene, terrestrial palynology was applied in the present study to two well-dated deep-marine sediment cores located at the Australo-Antarctic region: (i) IODP Site U1356 (Wilkes Land margin, East Antarctica) and (ii) ODP Site 1172 (East Tasman Plateau, southwest Pacific Ocean). The studied sequence from IODP Site U1356 comprises mid-shelfal sediments from the early to middle Eocene (53.9 – 46 million years ago [Ma]). For the ODP Site 1172, the studied succession is characterized by sediments deposited in shallow marine environments of the middle Paleocene to the early Eocene (60.7 – 54.2 Ma).
Based on the obtained pollen and spores (sporomorphs) results from the studied sequences of Site U1356 and Site 1172, this study aims to: (1) decipher the terrestrial climate conditions along the Australo-Antarctic region from the middle Paleocene to the middle Eocene; (2) evaluate the structure, diversity and compositional patterns of forests that throve in the Australo-Antarctic region during the early Paleogene; (3) understand the response of forests from the high southern latitudes to the climate dynamics from the early Paleogene; (4) establish a connection between the generated terrestrial palynomorph data and published Sea Surface Temperatures (SSTs) from the same cores.
To decipher the terrestrial climatic conditions on the Australo-Antarctic region, this study relies on the nearest living relative (NLR) concept that assumes that fossil taxa have similar climate requirements as their modern counterparts. This approach was applied to the sporomorph results of Site U1356 and Site 1172, following mainly the bioclimatic analysis. With regard to the structure and diversity patterns of the vegetation from the same region, the present study presents combined qualitative (i.e., reconstruction of the vegetation based mainly on the habitats of the known living relatives) and quantitative (i.e., application of ordination techniques, rarefaction and diversity indices) analyses of the fossil sporomorphs results.
The overall results from the paleoclimatic and vegetation reconstruction approaches applied in the present study, indicate that temperate and paratropical forests during the early Paleogene throve under different climatic conditions on the Wilkes Land margin and on Tasmania, at paleolatitudes of ∼70°S and ∼65°S, respectively.
Specifically, the sporomorph results from Site U1356, suggest that a highly diverse forest similar to present-day forests from New Caledonia was thriving on Antarctica during the early Eocene (53.9 – 51.9 Ma). These forests were characterized by the presence of termophilous taxa that are restricted today to tropical and subtropical settings, notably Bombacoideae, Strasburgeria, Beauprea, Spathiphyllum, Anacolosa and Lygodium. In combination with MBT/CBT paleotemperature results, they provide strong evidence for near-tropical warmth at least in the coastal lowlands along the Wilkes Land margin. The coeval presence of frost tolerant taxa such as Nothofagus, Araucariaceae and Podocarpaceae during the early Eocene on the same record suggests that paratropical forests were thriving along the Wilkes Land margin. Due to the presence of this kind of vegetation, it is possible to suggest that forests in this region were subject to a climatic gradient related to differences in elevation and/or the proximity to the coastline.
By the middle Eocene, the paratropical forests that characterized the vegetation of the early Eocene on the Wilkes Land margin were replaced by low diversity temperate forests dominated by Nothofagus, and similar to present-day cool-temperate forests from New Zealand. The dominance of these forests and the absence of thermophilous elements together with the lower temperatures suggested by the MBT/CBT and the sporomorph-based temperatures indicate consistently cooler conditions during this time interval.
With regard to the sporomorph results of Site 1172, this study suggests that three vegetation types were thriving on Tasmania from the middle Paleocene to the early Eocene under different climatic conditions. During the middle to late Paleocene, warm-temperate forests dominated by Podocarpaceae and Araucariaceae were the prevailing vegetation on Tasmania. The dominance of these forests was interrupted by the transient predominance of cool-temperate forests dominated by Nothofagus and Araucariaceae across the middle/late Paleocene transition interval (~59.5 to ~59.0 Ma). This cool-temperate forest was characterized by a lack of frost-sensitive elements (i.e., palms and cycads) indicating cooler conditions with harsher winters on Tasmania during this time interval. By the early Eocene, and linked with the Paleocene Eocene Thermal Maximum (PETM), Paleocene temperate forests dominated by gymnosperms were replaced by paratropical rainforests with the remarkable presence of the tropical mangrove palm Nypa during the PETM and the earliest Eocene. The overall results from Site U1356 and Site 1172, provide a new assessment of the terrestrial climatic conditions in the Australo-Antarctic region for validating climate models and understanding the response of high-latitude terrestrial ecosystems to the climate dynamics of the early Paleogene on southern latitudes.
The climatic conditions in the higher latitudes during the early Paleogene were further unravelled by comparing the obtained terrestrial and marine results. The integration of the obtained sporomorph data with previously published TEX86-based SSTs from Site 1172 documents that the vegetation dynamics were closely linked with the temperature evolution from the Australo-Antarctic region. Moreover, the comparison of TEX86-based SSTs and sporomorph-based climatic estimations from Site 1172 suggests a warm-season bias of both calibrations of TEX86 (i.e., TEX86Hand TEX86H), when this proxy is applied to high southern latitudes records of the early Paleogene.
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.
Reconstructing the early Paleogene climate dynamics of terrestrial settings in the high southern latitudes is important to assess the role of high-latitude physical and biogeochemical processes in the global climate system. However, whereas a number of high-quality Paleogene climate records has become available for the marine realm of the high southern latitudes over the recent past, the long-term evolution of coeval terrestrial climates and ecosystems is yet poorly known. We here explore the climate and vegetation dynamics on Tasmania from the middle Paleocene to the early Eocene (60.7–54.2 Ma) based on a sporomorph record from Ocean Drilling Program (ODP) Site 1172 on the East Tasman Plateau. Our results show that three distinctly different vegetation types thrived on Tasmania under a high-precipitation regime during the middle Paleocene to early Eocene, with each type representing different temperature conditions: (i) warm-temperate forests dominated by gymnosperms that were dominant during the middle and late Paleocene; (ii) cool-temperate forests dominated by southern beech (Nothofagus) and araucarians across the middle/late Paleocene transition interval (~59.5 to ~59.0 Ma); and (iii) paratropical forests rich in ferns that were established during and in the wake of the Paleocene–Eocene Thermal Maximum (PETM). The transient establishment of cool-temperate forests lacking any frost-sensitive elements (i.e., palms and cycads) across the middle/late Paleocene transition interval indicates markedly cooler conditions, with the occurrence of frosts in winter, on Tasmania during that time. The integration of our sporomorph data with previously published TEX86-based sea-surface temperatures from ODP Site 1172 documents that the vegetation dynamics on Tasmania were closely linked with the temperature evolution in the Tasman sector of the Southwest Pacific region. Moreover, the comparison of our season-specific climate estimates for the sporomorph assemblages from ODP Site 1172 with the TEX86L- and TEX86H-based temperature data suggests a warm-season bias of both calibrations for the early Paleogene of the high southern latitudes.
A new, two-channel instrument for simultaneous NO3 and N2O5 monitoring was used to make the first comprehensive set of nocturnal NOx measurements (NO, NO2, NO3 and N2O5) at the Taunus Observatory, a rural mountain site (Kleiner Feldberg) in South-western Germany. In May 2008, NO3 and N2O5 mixing ratios were well above the instrumental detection limit (a few ppt) on all nights of the campaign and were characterised by large variability resulting from inhomogeneously distributed sinks. The concentrations of NO3, N2O5 and NO2 were consistent with the equilibrium constant, K2, defining the rates of formation and thermal dissociation of N2O5. A steady-state lifetime analysis showed that nocturnal NOx losses were generally dominated by reaction of NO3 with volatile organic compounds in this forested region, with N2O5 uptake to aerosols of secondary importance. Analysis of a limited dataset obtained at high relative humidity indicated that the loss of N2O5 by reaction with water vapour is less efficient (> factor 3) than derived using laboratory kinetic data. The fraction of NOx present as NO3 and N2O5 reached ≈20% on some nights, with night-time losses of NOx competing with daytime losses.
A new, two-channel instrument for simultaneous NO3 and N2O5 monitoring was used to make the first comprehensive set of nocturnal NOx measurements (NO, NO2, NO3 and N2O5) at the Taunus Observatory, a rural mountain site (Kleiner Feldberg) in South-western Germany. In May 2008, NO3 and N2O5 mixing ratios were well above the instrumental detection limit (a few ppt) on all nights of the campaign and were characterised by large variability. The concentrations of NO3, N2O5 and NO2 were consistent with the equilibrium constant, K2, defining the rates of formation and thermal dissociation of N2O5. A steady-state lifetime analysis is consistent with the loss of nocturnal NOx being dominated by the reaction of NO3 with volatile organic compounds in this forested region, with N2O5 uptake to aerosols of secondary importance. Analysis of a limited dataset obtained at high relative humidity indicated that the loss of N2O5 by reaction with water vapour is less efficient (>factor 3) than derived using laboratory kinetic data. The fraction of NOx present as NO3 and N2O5 reached ~20% on some nights, with night-time losses of NOx competing with daytime losses.
We combine structural analysis of fractures with 22 U–Pb dates measured in fracture-filling carbonate cements from bed-parallel fibrous calcite veins (beef), conjugated veins and faults within the Vaca Muerta Formation along the Andean fold and thrust belt in the Neuquén Basin. The measured ages constrain accurately the relationships between overpressures caused by hydrocarbon generation and Andean compression as mechanisms for natural fracturing and vein formation. Two generations of fibres have been identified in beef. The first one, consists of dark fibres from the inner zones, which are perpendicular to bedding and contain abundant cone-in-cone structures and hydrocarbon inclusions. U–Pb dating of these fibres yielded Early to Late Cretaceous ages from 116.7 ± 17.7 to 78.8 ± 10.2 Ma. The second generation of fibres corresponds to the outer zones and consists of white fibres oblique to bedding, indicating growth during layer-parallel shortening. Bed-perpendicular veins cutting beef yielded Late Cretaceous-late Palaeocene dates from 72.8 ± 22.4 to 60.9 ± 10.4 Ma. Eocene ages from 52.0 ± 2.9 to 42.2 ± 18.9 Ma were measured in bed-parallel slip surfaces and reverse and strike-slip faults, whereas Miocene dates from 13.9 ± 2.6 to 6.2 ± 1.1 Ma were measured in E-W calcite veins. U–Pb dating of veins, structural analysis of fractures and subsidence curves, indicate that beef inner zones formed in the oil window during burial of the Neuquén basin, and that tectonic stresses could enhance their formation. Beef outer zones and bed-perpendicular veins formed during E-W Late Cretaceous-late Palaeocene layer-parallel shortening. Contrarily, late Palaeocene-late Eocene bed-parallel slip surfaces and faults and Miocene E-W veins formed during NE-SW and E-W syn-to post-folding deformation, respectively. In both cases, syn-to post-folding compression occurred synchronously with forelandward migration of magmatic activity attributed to flat subduction of the Pacific slab beneath the Andes.
Non-forest ecosystems, dominated by shrubs, grasses and herbaceous plants, provide ecosystem services including carbon sequestration and forage for grazing, and are highly sensitive to climatic changes. Yet these ecosystems are poorly represented in remotely sensed biomass products and are undersampled by in situ monitoring. Current global change threats emphasize the need for new tools to capture biomass change in non-forest ecosystems at appropriate scales. Here we developed and deployed a new protocol for photogrammetric height using unoccupied aerial vehicle (UAV) images to test its capability for delivering standardized measurements of biomass across a globally distributed field experiment. We assessed whether canopy height inferred from UAV photogrammetry allows the prediction of aboveground biomass (AGB) across low-stature plant species by conducting 38 photogrammetric surveys over 741 harvested plots to sample 50 species. We found mean canopy height was strongly predictive of AGB across species, with a median adjusted R2 of 0.87 (ranging from 0.46 to 0.99) and median prediction error from leave-one-out cross-validation of 3.9%. Biomass per-unit-of-height was similar within but different among, plant functional types. We found that photogrammetric reconstructions of canopy height were sensitive to wind speed but not sun elevation during surveys. We demonstrated that our photogrammetric approach produced generalizable measurements across growth forms and environmental settings and yielded accuracies as good as those obtained from in situ approaches. We demonstrate that using a standardized approach for UAV photogrammetry can deliver accurate AGB estimates across a wide range of dynamic and heterogeneous ecosystems. Many academic and land management institutions have the technical capacity to deploy these approaches over extents of 1–10 ha−1. Photogrammetric approaches could provide much-needed information required to calibrate and validate the vegetation models and satellite-derived biomass products that are essential to understand vulnerable and understudied non-forested ecosystems around the globe.
Two types of particles exist in the atmosphere, primary and secondary particles. While primary particles such as soot, mineral dust, sea salt particles or pollen are introduced directly as particles into the atmosphere, secondary particles are formed in the atmosphere by condensation of gases. The formation of such new aerosol particles takes place frequently and at a broad variety of atmospheric conditions and geographic locations. A considerable fraction of the atmospheric particles is formed by such nucleation processes. The newly formed particles may grow by condensation to sizes where they are large enough to act as cloud condensation nuclei and therefore may affect cloud properties. The fundamental processes of aerosol nucleation are described and typical atmospheric observations are discussed. Two recent studies are introduced that potentially change our current understanding of atmospheric nucleation substantially.
Number concentrations of total and non-volatile aerosol particles with size diameters >0.01 µm as well as particle size distributions (0.4–23 µm diameter) were measured in situ in the Arctic lower stratosphere (10–20.5 km altitude). The measurements were obtained during the campaigns European Polar Stratospheric Cloud and Lee Wave Experiment (EUPLEX) and Envisat-Arctic-Validation (EAV). The campaigns were based in Kiruna, Sweden, and took place from January to March 2003. Measurements were conducted onboard the Russian high-altitude research aircraft Geophysica using the low-pressure Condensation Nucleus Counter COPAS (COndensation PArticle Counter System) and a modified FSSP 300 (Forward Scattering Spectrometer Probe). Around 18–20 km altitude typical total particle number concentrations nt range at 10–20 cm−3 (ambient conditions). Correlations with the trace gases nitrous oxide (N2O) and trichlorofluoromethane (CFC-11) are discussed. Inside the polar vortex the total number of particles >0.01 µm increases with potential temperature while N2O is decreasing which indicates a source of particles in the above polar stratosphere or mesosphere. A separate channel of the COPAS instrument measures the fraction of aerosol particles non-volatile at 250°C. Inside the polar vortex a much higher fraction of particles contained non-volatile residues than outside the vortex (~24% outside vortex). This is most likely due to a strongly increased fraction of meteoritic material in the particles which is transported downward from the mesosphere inside the polar vortex. The high fraction of non-volatile residual particles gives therefore experimental evidence for downward transport of mesospheric air inside the polar vortex. It is also shown that the fraction of non-volatile residual particles serves directly as a suitable experimental vortex tracer. Nanometer-sized meteoritic smoke particles may also serve as nuclei for the condensation of gaseous sulfuric acid and water in the polar vortex and these additional particles may be responsible for the increase in the observed particle concentration at low N2O. The number concentrations of particles >0.4 µm measured with the FSSP decrease markedly inside the polar vortex with increasing potential temperature, also a consequence of subsidence of air from higher altitudes inside the vortex. Another focus of the analysis was put on the particle measurements in the lowermost stratosphere. For the total particle density relatively high number concentrations of several hundred particles per cm3 at altitudes below ~14 km were observed in several flights. To investigate the origin of these high number concentrations we conducted air mass trajectory calculations and compared the particle measurements with other trace gas observations. The high number concentrations of total particles in the lowermost stratosphere are probably caused by transport of originally tropospheric air from lower latitudes and are potentially influenced by recent particle nucleation.
Number concentrations of total and non-volatile aerosol particles with size diameters >0.01 μm as well as particle size distributions (0.4–23 μm diameter) were measured in situ in the Arctic lower stratosphere (10–20.5 km altitude). The measurements were obtained during the campaigns European Polar Stratospheric Cloud and Lee Wave Experiment (EUPLEX) and Envisat-Arctic-Validation (EAV). The campaigns were based in Kiruna, Sweden, and took place from January to March 2003. Measurements were conducted onboard the Russian high-altitude research aircraft Geophysica using the low-pressure Condensation Nucleus Counter COPAS (COndensation PArticle Counter System) and a modified FSSP 300 (Forward Scattering Spectrometer Probe). Around 18–20 km altitude typical total particle number concentrations nt range at 10–20 cm−3 (ambient conditions). Correlations with the trace gases nitrous oxide (N2O) and trichlorofluoromethane (CFC-11) are discussed. Inside the polar vortex the total number of particles >0.01 μm increases with potential temperature while N2O is decreasing which indicates a source of particles in the above polar stratosphere or mesosphere. A separate channel of the COPAS instrument measures the fraction of aerosol particles non-volatile at 250°C. Inside the polar vortex a much higher fraction of particles contained non-volatile residues than outside the vortex (~67% inside vortex, ~24% outside vortex). This is most likely due to a strongly increased fraction of meteoric material in the particles which is transported downward from the mesosphere inside the polar vortex. The high fraction of non-volatile residual particles gives therefore experimental evidence for downward transport of mesospheric air inside the polar vortex. It is also shown that the fraction of non-volatile residual particles serves directly as a suitable experimental vortex tracer. Nanometer-sized meteoric smoke particles may also serve as nuclei for the condensation of gaseous sulfuric acid and water in the polar vortex and these additional particles may be responsible for the increase in the observed particle concentration at low N2O. The number concentrations of particles >0.4 μm measured with the FSSP decrease markedly inside the polar vortex with increasing potential temperature, also a consequence of subsidence of air from higher altitudes inside the vortex. Another focus of the analysis was put on the particle measurements in the lowermost stratosphere. For the total particle density relatively high number concentrations of several hundred particles per cm3 at altitudes below ~14 km were observed in several flights. To investigate the origin of these high number concentrations we conducted air mass trajectory calculations and compared the particle measurements with other trace gas observations. The high number concentrations of total particles in the lowermost stratosphere are probably caused by transport of originally tropospheric air from lower latitudes and are potentially influenced by recent particle nucleation.
This study reports and discusses new radiometric ages, petrographical and volcanological observations and whole rock geochemical data of the rocks of the Rudnik Mts. volcano-intrusive complex. The complex hosts a Pb-Zn-Ag deposit and belongs to the Serbo-Macedonian metallogenetic belt. Two distinct igneous events are distignuished. The first occurred >30 Ma and was mainly characterized by extrusive and shallow intrusive dacites and andesites and was unrelated to mineralization. The second igneous event occurred <23 Ma and was highly heterogeneous in terms of volcanic products and petrographic varieties, but with predominance of quartzlatites. The dacite-andesites (first event) and the quartzlatites (second event) are geochemically similar and display a calc-alkaline affinity and highly incompatible element enriched patterns on spider diagrams, but the younger quartzlatites are richer in K2O, Rb and Ba and poorer in Sr. This is taken as evidence that mixing between an ultrapotassic lamprophyre/lamproite magma and an acid calc-alkaline (dacite-like) magma was essential petrogenetic processes during the second event. The proposed simplified volcanological model suggests that this mixing was responsible for triggering strongly explosive volcanic activity as well as for providing conditions for active hydrothermal and mineralization processes. The observed link between a specific magmatic phase and ore deposit formation can be a general phennomenon in the Balkans, and must be addressed by further and more advanced studies.
This article presents an environmental remote sensing application using a UAV that is specifically aimed at reducing the data gap between field scale and satellite scale in soil erosion monitoring in Morocco. A fixed-wing aircraft type Sirius I (MAVinci, Germany) equipped with a digital system camera (Panasonic) is employed. UAV surveys are conducted over different study sites with varying extents and flying heights in order to provide both very high resolution site-specific data and lower-resolution overviews, thus fully exploiting the large potential of the chosen UAV for multi-scale mapping purposes. Depending on the scale and area coverage, two different approaches for georeferencing are used, based on high-precision GCPs or the UAV’s log file with exterior orientation values respectively. The photogrammetric image processing enables the creation of Digital Terrain Models (DTMs) and ortho-image mosaics with very high resolution on a sub-decimetre level. The created data products were used for quantifying gully and badland erosion in 2D and 3D as well as for the analysis of the surrounding areas and landscape development for larger extents.
This article presents a multiscale approach for detecting and monitoring soil erosion phenomena (i.e. gully erosion) in the agro-industrial area around the city of Taroudannt, Souss basin, Morocco. The study area is characterized as semi-arid with an annual average precipitation of 200 mm. Water scarcity, high population dynamics and changing land use towards huge areas of irrigation farming present numerous threats to sustainability. The agro-industry produces citrus fruits and vegetables in monocropping, mainly for the European market. Badland areas strongly affected by gully erosion border the agricultural areas as well as residential areas. To counteract the significant loss of land, land-leveling measures are attempted to create space for plantations and greenhouses. In order to develop sustainable approaches to limit gully growth the detection and monitoring of gully systems is fundamental. Specific gully sites are monitored with unmanned aerial vehicle (UAV) taking small-format aerial photographs (SFAP). This enables extremely high-resolution analysis (SFAP resolution: 2-10 cm) of the actual size of the gully channels as well as a detailed continued surveillance of their growth. Transferring the methodology on a larger scale using Quickbird satellite data (resolution: 60 cm) leads to the possibility of a large-scale analysis of the whole area around the city of Taroudannt (Area extent: ca. 350 km²). The results will then reveal possible relationships of gully growth and agro-industrial management and may even illustrate further interdependencies. The main objective is the identification of areas with high gully-erosion risk due to non-sustainable land use and the development of mitigation strategies for the study area.
Climatology of morphology and cloud-radiative properties of marine low-level mixed-phase clouds
(2023)
Marine stratocumuli cover about 40 - 60% of the ocean surface. They self-organize into different morphological regimes. The two organized cellular regimes are called open and closed mesoscale-cellular convective (MCC) clouds. In mid-to-high latitudes, open and closed cells are the two most frequent types of MCC clouds. In particular, many MCC clouds consist of a mixture of vapor, liquid droplets, and ice particles, referred to as mixed-phase clouds (MPCs). Even for the same cloud fraction, the albedo of open cells is, on average, lower than that of closed MCC clouds. Cloud phase and morphology individually influence the cloud radiative effect. Thus, this thesis investigates the relationships between the cloud phase, MCC organization, cell size, and differences regarding the cloud-radiative effect.
This thesis focuses on space-borne retrievals to achieve extensive temporal and spatial coverage. The liDAR-raDAR (DARDAR) version 2 product collocates two active and one passive satellite: CloudSat, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), and Moderate Resolution Imaging Spectroradiometer (MODIS). The cloud phase of DARDAR is vertically integrated to establish a single cloud phase at each data point. The MCC classification data set based on the liquid water path (LWP) of MODIS scenes is collocated with the DARDAR product to determine the MCC organization. Cell-size statistics of both MCC clouds are obtained using a marker-based image segmentation method on MODIS reflectance scenes. In addition, based on MODIS reflectance scenes, a convolutional neural network (CNN) is developed to classify open and closed MCC scenes to avoid missing mature MPCs with a low LWP.
The first part of this thesis explores the relationships between cloud phase, morphology, and cloud albedo in the Southern Ocean (SO). At a given cloud-top temperature (CTT), seasonal changes in the mixed-phase fraction, defined as the number of MPCs divided by the sum of MPC and supercooled liquid cloud (SLC) pixels, are stronger than the morphological changes. Therefore, external factors seem to influence these changes instead of morphology. The dependence of cloud phase on cloud-top height (CTH) is more substantial than on CTT in clouds with CTHs below 2.5 km. The previously observed acceleration of closed-to-open transition in MPCs, known as preconditioning, is not the primary driver of climatological cloud morphology statistics in the SO. The morphological differences in cloud albedo are more pronounced in SLCs than in MPCs. This change in albedo alters the cloud radiative effect in the SO by 21Wm−2 to 39Wm−2 depending onseason and cloud phase.
Open and closed MCC clouds exhibit larger equivalent cell diameters in the MPCs than in SLCs in austral summer, whereas, in austral winter, the SLCs are larger. The cell’s aspect ratio accounts for varying CTHs. Closed cells have smaller aspect ratios than open cells, so their cell diameter is smaller, independent of CTH. While the seasonal differences in closed cells are due to changes in CTH, the seasonal aspect ratio differences in open cells are mainly caused by MPCs. With increasing aspect ratios, the cloud albedo decreases in both open and closed MCC clouds, with the most substantial decrease in open MPCs clouds. This leads to cloud-radiative changes of 60 - 75Wm−2 in the SO, depending on cloud phase and aspect ratio.
The established CNN exhibits a good accuracy of 80.6%, with even higher accuracies in the Open (85.5%) and Closed (87.3%) categories. The global MCC climatology based on the CNN generally agrees well with previous MCC distributions. The most notable difference occurs in the Northern Hemisphere (NH) in boreal winter, with a higher occurrence frequency of closed and open MCC clouds. This might indicate missing MPCs in previous studies based on the LWP and some restricted to warm cloud scenes. Thus, the developed CNN seems to better represent the different morphologies in MPCs than in previous classifications.
In conclusion, this thesis shows that understanding the dependencies of cloud phase, cloud morphology, and cell size is important to enhance predictions of the cloud-radiative effect and thus, it is important to evaluate how cloud phase, cloud morphology, and cellsize change in a warming climate.
Community trait assembly in highly diverse tropical rainforests is still poorly understood. Based on more than a decade of field measurements in a biodiversity hotspot of southern Ecuador, we implemented plant trait variation and improved soil organic matter dynamics in a widely used dynamic vegetation model (the Lund-Potsdam-Jena General Ecosystem Simulator, LPJ-GUESS) to explore the main drivers of community assembly along an elevational gradient. In the model used here (LPJ-GUESS-NTD, where NTD stands for nutrient-trait dynamics), each plant individual can possess different trait combinations, and the community trait composition emerges via ecological sorting. Further model developments include plant growth limitation by phosphorous (P) and mycorrhizal nutrient uptake. The new model version reproduced the main observed community trait shift and related vegetation processes along the elevational gradient, but only if nutrient limitations to plant growth were activated. In turn, when traits were fixed, low productivity communities emerged due to reduced nutrient-use efficiency. Mycorrhizal nutrient uptake, when deactivated, reduced net primary production (NPP) by 61–72% along the gradient. Our results strongly suggest that the elevational temperature gradient drives community assembly and ecosystem functioning indirectly through its effect on soil nutrient dynamics and vegetation traits. This illustrates the importance of considering these processes to yield realistic model predictions.
Projections of future changes in runoff can have important implications for water resources and flooding. In this study, runoff projections from ISI-MIP (Inter-sectoral Impact Model Intercomparison Project) simulations forced with HadGEM2-ES bias-corrected climate data under the Representative Concentration Pathway 8.5 have been analysed. Projections of change from the baseline period (1981–2010) to the future (2070–2099) from a number of different ecosystems and hydrological models were studied. The differences between projections from the two types of model were looked at globally and regionally. Typically, across different regions the ecosystem models tended to project larger increases and smaller decreases in runoff than the hydrological models. However, the differences varied both regionally and seasonally. Sensitivity experiments were also used to investigate the contributions of varying CO2 and allowing vegetation distribution to evolve on projected changes in runoff. In two out of four models which had data available from CO2 sensitivity experiments, allowing CO2 to vary was found to increase runoff more than keeping CO2 constant, while in two models runoff decreased. This suggests more uncertainty in runoff responses to elevated CO2 than previously considered. As CO2 effects on evapotranspiration via stomatal conductance and leaf-area index are more commonly included in ecosystems models than in hydrological models, this may partially explain some of the difference between model types. Keeping the vegetation distribution static in JULES runs had much less effect on runoff projections than varying CO2, but this may be more pronounced if looked at over a longer timescale as vegetation changes may take longer to reach a new state.
Comparing projections of future changes in runoff from hydrological and biome models in ISI-MIP
(2013)
Future changes in runoff can have important implications for water resources and flooding. In this study, runoff projections from ISI-MIP (Inter-sectoral Impact Model Intercomparison Project) simulations forced with HadGEM2-ES bias-corrected climate data under the Representative Concentration Pathway 8.5 have been analysed for differences between impact models. Projections of change from a baseline period (1981–2010) to the future (2070–2099) from 12 impacts models which contributed to the hydrological and biomes sectors of ISI-MIP were studied. The biome models differed from the hydrological models by the inclusion of CO2 impacts and most also included a dynamic vegetation distribution. The biome and hydrological models agreed on the sign of runoff change for most regions of the world. However, in West Africa, the hydrological models projected drying, and the biome models a moistening. The biome models tended to produce larger increases and smaller decreases in regionally averaged runoff than the hydrological models, although there is large inter-model spread. The timing of runoff change was similar, but there were differences in magnitude, particularly at peak runoff. The impact of vegetation distribution change was much smaller than the projected change over time, while elevated CO2 had an effect as large as the magnitude of change over time projected by some models in some regions. The effect of CO2 on runoff was not consistent across the models, with two models showing increases and two decreases. There was also more spread in projections from the runs with elevated CO2 than with constant CO2. The biome models which gave increased runoff from elevated CO2 were also those which differed most from the hydrological models. Spatially, regions with most difference between model types tended to be projected to have most effect from elevated CO2, and seasonal differences were also similar, so elevated CO2 can partly explain the differences between hydrological and biome model runoff change projections. Therefore, this shows that a range of impact models should be considered to give the full range of uncertainty in impacts studies.
The objective of the present doctoral thesis was to investigate the occurrence, distribution, and behaviour of six hydrophilic ethers: ethyl tert-butyl ether (ETBE), 1,4-dioxane, ethylene glycol dimethyl ether (monoglyme), diethylene glycol dimethyl ether (diglyme), triethylene glycol dimethyl ether (triglyme), and tetraethylene glycol dimethyl ether (tetraglyme) in surface-, waste-, ground- and drinking water samples. Solid phase extraction and gas chromatography/mass spectrometry were used to analyze the six hydrophilic ethers. Altogether more than 150 surface water samples, almost 100 of each groundwater and wastewater samples, and 10 raw and drinking water samples were analyzed during the research project.
Initially, the method was validated in order to simultaneously determine the analytes of interest in various aquatic environments. A solid phase extraction method that uses coconut charcoal (Resprep® activated coconut charcoal, Restek) or carbon molecular sieve material (SupelcleanTM Envi-CarbTM Plus, Supelco) for analyte absorption were found suitable for determination of ETBE, 1,4-dioxane, and glymes in surface-, drinking-, ground- and wastewater samples. Precision and accuracy of both methods was demonstrated for all analytes of interest. The recovery of target compounds from the ultrapure water spiked at 1.0 µg L−1 was between 86.8 % and 98.2 %, with relative standard deviation below 6 %. The samples spiked at 10.0 µg L−1 gave slightly higher recovery of 90.6 % to 112.2 % with a relative standard deviation below 3.4 % for each analyte. Detection and quantification limits in ultrapure water and surface waters were furthermore established. The limit of quantitation (LOQ) in ultrapure water ranged between 0.024 µg L−1 to 0.057 µg L−1 using Restek cartridges, and 0.030 µg L−1 to 0.069 µg L−1 using Supelco cartridges. In the surface water samples the calculated LOQ was 0.032 µg L−1 to 0.067µg L−1 using coconut charcoal material and 0.032 µg L−1 to 0.052 µg L−1 using the carbon molecular sieve material. Moreover, stability of the unpreserved and preserved water samples as well as the extracts was determined. Preservation of samples with sodium bisulfate (at 1 gram per Liter) resulted in much better stability of the ethers in water samples. Subsequently, 27 samples obtained from seven surface water bodies in Germany (Rivers Rhine, Lippe, Main, Oder, Rur, Schwarzbach and Wesel-Datteln Canal) were analyzed for the six hydrophilic ethers. ETBE was present in only two surface waters (Rhine River and Wesel-Datteln Canal) with concentrations close to the LOQ (up to 0.065 µg L−1). 1,4-Dioxane was detected in all of the water samples at concentrations reaching 1.93 µg L–1. Monoglyme was identified only in the Main and Rhine Rivers at the maximum concentration of 0.114 µg L–1 and 0.427 µg L–1, respectively. Very high concentrations (up to 1.73 µg L−1) of diglyme, triglyme, and tetraglyme were detected in the samples from the Oder River. These glymes were also detected in the Rhine River; however the concentrations did not exceed 0.200 µg L–1. Furthermore, tetraglyme was detected in the Main River at an average concentration of 0.409 µg L–1 (n = 6) and in one sample from the Rur River at 0.192 µg L–1.
Four sampling campaigns were conducted at the Oderbruch polder between October 2009 and May 2012, in order to study the behavior of the hydrophilic ethers and organophosphates during riverbank filtration and in the anoxic aquifer. Moreover the suitability of these target compounds was assessed for their use as groundwater organic tracers. At the time of each sampling campaign, concentrations of triglyme and tetraglyme in the Oder River were between 20–185 ng L–1 (n = 4) and 273¬–1576 ng L–1 (n = 4). Monoglyme, diglyme, and 1,4-dioxane were analyzed only during the two last sampling campaigns. At that time, the concentration of diglyme in Oder River was 65¬–94 ng L-1 (n = 2) and 1,4-dioxane 1610¬–3290 ng L–1 (n = 2). In the drainage ditch, following bank filtration, concentrations of ethers ranged between 1090 ng L–1 and 1467 ng L–1 for 1,4-dioxane, 23¬ng L–1 and 41 ng L–1 for diglyme, 37 ng L–1 and 149 ng L–1 for triglyme, and 496 ng L–1 and 1403 ng L–1 for tetraglyme. In the anoxic aquifer, 1,4-dioxane showed the greatest persistence during the groundwater passage. At the distance of 1150 m from the river and an estimated groundwater age of 41.9 years, a concentration above 200 ng L−1 was detected. A positive correlation was found for the inorganic tracer chloride (Cl−) with 1,4-dioxane and tetraglyme. Similarities in the behavior of Cl− and the organic compound suggested that 1,4-dioxane and tetraglyme are controlled by the same hydraulic process and therefore can be used as additional tracers to study the dynamics of the groundwater system. These results show that high concentrations of ethers are present in the surface water and are not removed during bank filtration processes. Moreover, the hydrophilic ethers persist in the anoxic aquifer and little or no degradation is expected, supporting, their possible application as organic tracers.
A separate sampling project was conducted for 1,4-dioxane that focused primarily on its fate in the aquatic environment. This study provided missing information on the extent of water pollution with 1,4-dioxane is Germany. Numerous waste-, surface-, ground- and drinking water samples were collected in order to determine the persistence of 1,4-dioxane in the aquatic environment. The occurrence of 1,4-dioxane was determined in wastewater samples from four municipal sewage treatment plants (STP). The influent and effluent samples were collected during weekly campaigns. The average influent concentrations in all four plants ranged from 262 ± 32 ng L−1 to 834 ± 480 ng L−1, whereas the average effluents concentrations were between 267 ± 35 ng L−1 and 62,260 ± 36,000 ng L−1. The source of increased 1,4-dioxane concentrations in one of the effluents was identified to originate from impurities in the methanol used in the postanoxic denitrification process. Spatial and temporal distribution of 1,4-dioxane in the river Main, Rhine, and Oder was also examined. Concentrations reaching 2,200 ng L−1 in the Oder River, and 860 ng L−1 in both Main and Rhine River were detected. The average load during the sampling was estimated to be 6.5 kg d−1 in the Main, 34.1 kg d−1 in the Oder, and 134.5 kg d−1 in the Rhine River. In all of the sampled rivers, concentrations of 1,4-dioxane increased with distance from the mouth of the river and were found to negatively correlate with the discharge of the river. In order to determine if 1,4-dioxane can reach drinking water supplies, samples from a Rhine River bank filtration site and potable water from two drinking water production facilities were analyzed for the presence of 1,4-dioxane in the raw water and finished potable water. The raw water (following bank filtration) contained 650 ng L−1 to 670 ng L−1 of 1,4-dioxane, whereas the concentration in the finished drinking water fell only to 600 ng L−1 and 490 ng L−1, respectively.
During the final project, investigations of the source identification of high glyme concentrations in the Oder River were carried out. During four sampling campaigns between January, 2012 and April, 2013, 50 samples from the Oder River in the Oderbruch region and Poland were collected. During the first two samplings in the Oderbruch polder, glymes were detected at concentration reaching 0.07 µg L-1 (diglyme), 0.54 µg L−1 (triglyme) and 1.73 µg L−1 (tetraglyme) in the Oder River. The extensive sampling campaign of the Oder River (about 500 km) in Poland helped to identify the area of possible glyme entry into the river. During that sampling the maximum concentrations of triglyme and tetraglyme were 0.46 µg L−1 and 2.21 µg L−1, respectively. A closer investigation of the identified area of pollution, helped to determine the possible sources of glymes in the Oder River. Hence, the final sampling focused on the Kaczawa River, a left tributary of the Oder River and Czarna Woda, a left tributary of Kaczawa River. Moreover, samples from an industrial wastewater treatment plant were collected. Samples from Czarna Woda stream and Kaczawa River contained even higher concentrations of diglyme, triglyme, and tetraglyme, reaching 5.18 µg L−1, 12.87 µg L−1 and 80.81 µg L−1, respectively. Finally, three water samples from a wastewater treatment plant receiving influents from a copper smelter were analyzed. Diglyme, triglyme, and tetraglyme were present at an average concentration of 569 µg L−1, 4300 µg L−1, and 65900 µg L−1, respectively in the wastewater. Further research helped to identify the source of the glymes in the wastewater. The gas desulfurization process – Solinox implemented in the nearby copper smelter uses glymes as physical absorption medium for sulfur dioxide.
Results of this doctoral research provide important information about the occurrence, distribution, and behavior of hydrophilic ethers: 1,4-dioxane, monoglyme, diglyme, triglyme, and tetraglyme in the aquatic environment. A method capable of analyzing a wide range of ether compounds: from a volatile ETBE to a high molecular weight tetraglyme was validated. 1,4-Dioxane and tetraglyme were found to be applicable as organic tracers, since they are not easily attenuated during bank filtration and the anoxic groundwater passage. The extent of water pollution with 1,4-dioxane was shown in waste-, surface-, ground-, and drinking waters. One source of extremely high concentrations of 1,4-dioxane in a municipal sewage treatment plant applying postanoxic denitrification was identified, however more information is needed on the entry of 1,4-dioxane into surface waters. Moreover, 1,4-dioxane was present in drinking water samples from river bank filtration, which demonstrates its persistence in the aquatic environment and its low degradation potential during bank filtration and subsequent water treatment. Furthermore, this was the first study that focused primarily on identifying sources of glymes in surface waters. Glymes find a widespread use in industrial sectors, hence establishing their origin in the surface water is difficult (as with 1,4-dioxane). In this work, a gas desulphurization process was identified to be a dominating source of glyme pollution in the Oder River.
Remote sensing data are essential for monitoring the Earth’s surface waters, especially since the amount of publicly available in-situ data is declining. Satellite altimetry provides valuable information on the water levels and variations of lakes, reservoirs and rivers. In combination with satellite imagery, the derived time series allow the monitoring of lake storage changes and river discharge. However, satellite altimetry is limited in terms of its spatial resolution due to its measurement geometry, only providing information in the nadir direction beneath the satellite’s orbit. In a case study in the Mississippi River Basin (MRB), this study investigates the potential and limitations of past and current satellite missions for the monitoring of basin-wide storage changes. For that purpose, an automated target detection is developed and the extracted lake surfaces are merged with the satellites’ tracks. This reveals that the current altimeter configuration misses about 80% of all lakes larger than 0.1 km2 in the MRB and 20% of lakes larger than 10 km2, corresponding to 30% and 7% of the total water area, respectively. Past altimetry configurations perform even more poorly. From the larger water bodies represented by a global hydrology model, at least 91% of targets and 98% of storage changes are captured by the current altimeter configuration. This will improve significantly with the launch of the planned Surface Water and Ocean Topography (SWOT) mission.
In this study, we construct a new monthly zonal mean carbon dioxide (CO2) distribution from the upper troposphere to the stratosphere over the 2000–2010 time period. This reconstructed CO2 product is based on a Lagrangian backward trajectory model driven by ERA-Interim reanalysis meteorology and tropospheric CO2 measurements. Comparisons of our CO2 product to extratropical in situ measurements from aircraft transects and balloon profiles show remarkably good agreement. The main features of the CO2 distribution include (1) relatively large mixing ratios in the tropical stratosphere; (2) seasonal variability in the extratropics, with relatively high mixing ratios in the summer and autumn hemisphere in the 15–20 km altitude layer; and (3) decreasing mixing ratios with increasing altitude from the upper troposphere to the middle stratosphere ( ∼ 35 km). These features are consistent with expected variability due to the transport of long-lived trace gases by the stratospheric Brewer–Dobson circulation. The method used here to construct this CO2 product is unique from other modelling efforts and should be useful for model and satellite validation in the upper troposphere and stratosphere as a prior for inversion modelling and to analyse features of stratosphere–troposphere exchange as well as the stratospheric circulation and its variability.
In this study, we aim to reconstruct a relevant and new database of monthly zonal mean distribution of carbon dioxide (CO2) at global scale extending from the upper-troposphere (UT) to stratosphere (S). This product can be used for model and satellite validation in the UT/S, as a prior for inversion modelling and mainly to analyse a plausible feature of the stratospherictropospheric exchange as well as the stratospheric circulation and its variability. To do so, we investigate the ability of a Lagrangian trajectory model guided by ERA-Interim reanalysis to construct the CO2 abundance in the UT/S. From 10 year backward trajectories and tropospheric observations of CO2, we reconstruct upper-tropospheric and stratospheric CO2 over the period 2000–2010. The inter-comparisons of the reconstructed CO2 with mid-latitude vertical profiles measured by balloon samples as well as quasi-horizontal air samples from ER-2 aircraft during SOLVE and CONTRAIL campaigns exhibit a remarkable agreement. That demonstrates the potential of Lagrangian model to reconstruct CO2 in the UT/S. The zonal mean distribution exhibits relatively large CO2 in the tropical stratosphere due to the seasonal variation of the tropical upwelling of Brewer-Dobson circulation. During winter and spring, the tropical pipe is relatively isolated but is less confined during summer and autumn so that high CO2 values are more readily transported out of the tropics to the mid- and high latitude stratosphere. The shape of the vertical profiles suggests that relatively high CO2 above 20 km altitude mainly enter the stratosphere through tropical upwelling. CO2 mixing ratio is relatively low in the polar and tropical regions above 25 km. On average the CO2 mixing ratio decreases with altitude by 6-8 ppmv from the UT to stratosphere (e.g. up to 35 km) and is nearly constant with altitude.
The CLOUD (Cosmics Leaving OUtdoor Droplets) experiment at CERN is studying the nucleation and growth of aerosol particles under atmospheric conditions, and their activation into cloud droplets. A key feature of the CLOUD experiment is precise control of the experimental parameters. Temperature uniformity and stability in the chamber are important since many of the processes under study are sensitive to temperature and also to contaminants that can be released from the stainless steel walls by upward temperature fluctuations. The air enclosed within the 3 m CLOUD chamber is equipped with several arrays (strings) of high precision, fast-response thermometers to measure its temperature. Here we present a study of the air temperature uniformity inside the CLOUD chamber under various experimental conditions. Measurements were performed under calibration conditions and run conditions, which are distinguished by the flow rate of fresh air and trace gases entering the chamber: 20 l/min and up to 210 l/min, respectively. During steady-state calibration runs between −70 °C and +20 °C, the air temperature uniformity is better than +/−0.06 °C in the radial direction and +/−0.1 °C in the vertical direction. Larger non-uniformities are present during experimental runs, depending on the temperature control of the make-up air and trace gases (since some trace gases require elevated temperatures until injection into the chamber). The temperature stability is a few times 0.01 °C over periods of several hours during either calibration or steady-state run conditions. During rapid adiabatic expansions to activate cloud droplets and ice particles, the chamber walls are up to 10 °C warmer than the enclosed air. This results in larger non-uniformities while the air returns to its equilibrium temperature with time constant of about 200 s.
Projections from coarse-grid global circulation models are not suitable for regional estimates of water balance or trends of extreme precipitation and temperature, especially not in complex terrain. Thus, downscaling of global to regionally resolved projections is necessary to provide input to integrated water resources management approaches for river basins like the Upper Danube River Basin (UDRB) and the Upper Brahmaputra River Basin (UBRB).
This paper discusses the application of the regional climate model COSMO-CLM as a dynamical downscaling tool. To provide accurate data the COSMO-CLM model output was post-processed by statistical means. This downscaling chain performs well in the baseline period 1971 to 2000. However, COSMO-CLM performs better in the UDRB than in the UBRB because of a longer application experience and a less complex climate in Europe.
Different climate change scenarios were downscaled for the time period 1960–2100. The projections show an increase of temperature in both basins and for all seasons. The values are generally higher in the UBRB with the highest values occurring in the region of the Tibetan Plateau. Annual precipitation shows no substantial change. However, seasonal amounts show clear trends, for instance an increasing amount of spring precipitation in the UDRB. Again, the largest trends for different precipitation statistics are projected in the region of the Tibetan Plateau. Here, the projections show up to 50% longer dry periods in the months June to September with a simultaneous increase of about 10% for the maximum amount of precipitation on five consecutive days. For the Assam region in India, the projections also show an increase of 25% in the number of consecutive dry days during the monsoon season leading to prolonged monsoon breaks.
Melting inside earth is a common phenomenon and can be observed in many different regions where melt travels through the mantle and crust to eventually reach the surface where it crystallizes to build large volcanic provinces, whole stratigraphic layers of flood basalts, or even the oceanic crust. Often, melt reaching the surface is a good source of information. It can be used to achieve a better understanding about processes taking place in deeper regions inside the mantle and it is therefore essential to fundamentally understand melting and melt percolation processes. In order to achieve a deeper understanding, the aim of this thesis is to investigate processes that are connected to melting by using numerical models.
The physical model used is a so called two-phase flow model which describes the ability of melt to percolate through a viscously deforming, partially molten matrix. A famous feature of two-phase flow are solitary porosity waves, which are waves of locally higher porosity ascending through a partially molten background, keeping its shape constant, driven by decompaction and compaction of the solid matrix in front and behind the wave.
The viscosity law for shear- and volume viscosity was strongly simplified in most previous studies that modeled solitary waves. Often the porosity dependency is underestimated or its influence on the volume viscosity is even neglected, leading to too high viscosities. In this work more realistic laws are used that strongly decrease for small melt fractions. Those laws are incorporated into a 2D Finite Difference mantle convection code with two-phase flow to study the ascent of solitary porosity waves.
The model results show that an initial Gaussian-shaped wave rapidly evolves into a solitary wave with a certain amplitude, traveling upwards with constant velocity. Even though strongly weaker viscosities are used, the effect on dispersion curves and wave shape are only minor as long as the background porosity is rather small. The results are still in agreement to semi-analytical solutions which neglect shear stresses in the melt segregation equation. Higher background porosities and wave amplitudes lead to significant decrease in phase velocity and wave width, as the viscosity is strongly effected. However, the models show that solitary waves are still a possible mechanism for more realistic matrix viscosities.
While the ascending of porosity waves are mostly described by the movement of fluid melt, partially molten regions inside Earth trigger upwelling of both, solid and fluid phases, which can be called diapirism. While diapirs can have a wide variety of wavelengths, porosity waves are restricted to a few times the compaction length. The size of a melt perturbation in terms of compaction length therefore describes whether material is transported by diapirism or porosity waves. In this thesis we study the transition from diapiric rise to solitary porosity waves by systematically changing the size of a porosity perturbation from 1.8 to 120 times the compaction length. In case of a perturbation of the size of a few times the compaction length a single porosity wave will emerge, either with a positive or negative vertical matrix flux and if melt is not allowed to move relative to the matrix a diapir will emerge. In between these physical end members a regime can be observed where the partially molten perturbation will split up into numerous solitary waves, whose phase velocity is low compared to the Stokes velocity and the swarm of solitary waves will ascend jointly as a diapir, slowly elongating due to a higher amplitude main solitary wave.
Solitary waves will always emerge from a melt perturbation as long as two-phase flow is enabled, but the time for a solitary wave to emerge increases non-linearly with the perturbation radius in terms of compaction length. In nature, in many cases this time might be too long for solitary waves to emerge.
Another important feature when it comes to two-phase flow is the transport of trace elements in melt. Incompatible elements prefer to go into the melt, which eventually enriches the area where it crystallizes again. In order to model this redistribution, the code FDCON was extended to allow for fully consistent transport of elements in melt, including melting, freezing and re-equilibration with time. A 2D model, a simple representation of a volcanic back arc, is set up to investigate the behavior of trace elements. The influence of retention number and re-equilibration time is examined. Lava-lamp like convection can be observed in the lower part of the model, producing melt, that eventually leads to enrichment in trace elements in the upper high-viscous layer. The total enrichment in this layer approaches an asymptotic value and a 0D model is introduced to recreate this behavior.
In partially molten regions inside the earth melt buoyancy may trigger upwelling of both solid and fluid phases, i.e. diapirism. If the melt is allowed to move separately with respect to the matrix, melt perturbations may evolve into solitary porosity waves. While diapirs may form on a wide range of scales, porosity waves are restricted to sizes of a few times the compaction length. Thus, the size of a partially molten perturbation controls whether a diapir or a porosity wave will emerge. We study the transition from diapiric rise to solitary porosity waves by solving the two-phase flow equations of conservation of mass and momentum in 2D with porosity dependent matrix viscosity. We systematically vary the initial size of a porosity perturbation from 1 to 100 times the compaction length. If the perturbation is much larger than a regular solitary wave, its Stokes velocity is large and therefore faster than the segregating melt. Consequently, the fluid is not able to form a porosity wave and a diapir emerges. For small perturbations solitary waves emerge, either with a positive or negative vertical matrix velocity inside. In between the diapir and solitary wave regimes we observe a third regime of solitary wave induced focusing of melt. In these cases, diapirism is dominant but the fluid is still fast enough to locally build up small solitary waves which rise slightly faster than the diapir and form finger like structures at the front of the diapir. In our numerical simulations the width of these fingers is controlled by the compaction length or the grid size, whichever is larger. In cases where the compaction length becomes similar to or smaller than the grid size the finger-like leading solitary porosity waves are no more properly resolved, and too big and too fast waves may be the result. Therefore, one should be careful in large scale two-phase flow modelling with melt focusing especially when compaction length and grid size are of similar order.
In partially molten regions inside the Earth, melt buoyancy may trigger upwelling of both solid and fluid phases, i.e., diapirism. If the melt is allowed to move separately with respect to the matrix, melt perturbations may evolve into solitary porosity waves. While diapirs may form on a wide range of scales, porosity waves are restricted to sizes of a few times the compaction length. Thus, the size of a partially molten perturbation in terms of compaction length controls whether material is dominantly transported by porosity waves or by diapirism. We study the transition from diapiric rise to solitary porosity waves by solving the two-phase flow equations of conservation of mass and momentum in 2D with porosity-dependent matrix viscosity. We systematically vary the initial size of a porosity perturbation from 1.8 to 120 times the compaction length. If the perturbation is of the order of a few compaction lengths, a single solitary wave will emerge, either with a positive or negative vertical matrix flux. If melt is not allowed to move separately to the matrix a diapir will emerge. In between these end members we observe a regime where the partially molten perturbation will split up into numerous solitary waves, whose phase velocity is so low compared to the Stokes velocity that the whole swarm of waves will ascend jointly as a diapir, just slowly elongating due to a higher amplitude main solitary wave. Only if the melt is not allowed to move separately to the matrix will no solitary waves build up, but as soon as two-phase flow is enabled solitary waves will eventually emerge. The required time to build them up increases nonlinearly with the perturbation radius in terms of compaction length and might be too long to allow for them in nature in many cases.
Motions on planetary spatial scales in the atmosphere are governed by the planetary geostrophic equations. However, little attention has been paid to the interaction between the baroclinic and barotropic flows within the planetary geostrophic scaling. This is the focus of the present study, which utilizes planetary geostrophic equations for a Boussinesq fluid supplemented by a novel evolution equation for the barotropic flow. The latter is affected by meridional momentum flux due to baroclinic flow and drag by the surface wind. The barotropic wind, on the other hand, affects the baroclinic flow through buoyancy advection. Via a relaxation towards a prescribed buoyancy profile the model produces realistic major features of the zonally symmetric wind and temperature fields. We show that there is considerable cancellation between the barotropic and the baroclinic surface zonal mean zonal winds. Linear and nonlinear model responses to steady diabatic zonally asymmetric forcing are investigated, and the arising stationary waves are interpreted in terms of analytical solutions. We also study the problem of baroclinic instability on the sphere within the present model.
Balloon-borne stratospheric BrO measurements : comparison with Envisat/SCIAMACHY BrO limb profiles
(2005)
For the first time, results of all four existing stratospheric BrO profiling instruments, are presented and compared with reference to the SLIMCAT 3-dimensional chemical transport model (3-D CTM). Model calculations are used to infer a BrO profile validation set, measured by 3 different balloon sensors, for the new Envisat/SCIAMACHY (ENVIronment SATellite/SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY) satellite instrument. The balloon observations include (a) balloon-borne in situ resonance fluorescence detection of BrO, (b) balloon-borne solar occultation DOAS measurements (Differential Optical Absorption Spectroscopy) of BrO in the UV, and (c) BrO profiling from the solar occultation SAOZ (Systeme d'Analyse par Observation Zenithale) balloon instrument. Since stratospheric BrO is subject to considerable diurnal variation and none of the measurements are performed close enough in time and space for a direct comparison, all balloon observations are considered with reference to outputs from the 3-D CTM. The referencing is performed by forward and backward air mass trajectory calculations to match the balloon with the satellite observations. The diurnal variation of BrO is considered by 1-D photochemical model calculation along the trajectories. The 1-D photochemical model is initialised with output data of the 3-D model with additional constraints on the vertical transport, the total amount and photochemistry of stratospheric bromine as given by the various balloon observations. Total [Bry]=(20.1±2.8)pptv obtained from DOAS BrO observations at mid-latitudes in 2003, serves as an upper limit of the comparison. Most of the balloon observations agree with the photochemical model predictions within their given error estimates. First retrieval exercises of BrO limb profiling from the SCIAMACHY satellite instrument agree to <±50% with the photochemically-corrected balloon observations, and tend to show less agreement below 20 km.
Balloon-borne stratospheric BrO measurements: comparison with Envisat/SCIAMACHY BrO limb profiles
(2006)
For the first time, results of all four existing stratospheric BrO profiling instruments, are presented and compared with reference to the SLIMCAT 3-dimensional chemical transport model (3-D CTM). Model calculations are used to infer a BrO profile validation set, measured by 3 different balloon sensors, for the new Envisat/SCIAMACHY (ENVIronment SATellite/SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY) satellite instrument. The balloon observations include (a) balloon-borne in situ resonance fluorescence detection of BrO, (b) balloon-borne solar occultation DOAS measurements (Differential Optical Absorption Spectroscopy) of BrO in the UV, and (c) BrO profiling from the solar occultation SAOZ (Systeme d'Analyse par Observation Zenithale) balloon instrument. Since stratospheric BrO is subject to considerable diurnal variation and none of the measurements are performed close enough in time and space for a direct comparison, all balloon observations are considered with reference to outputs from the 3-D CTM. The referencing is performed by forward and backward air mass trajectory calculations to match the balloon with the satellite observations. The diurnal variation of BrO is considered by 1-D photochemical model calculation along the trajectories. The 1-D photochemical model is initialised with output data of the 3-D model with additional constraints on the vertical transport, the total amount and photochemistry of stratospheric bromine as given by the various balloon observations. Total [Bry]=(20.1±2.8)pptv obtained from DOAS BrO observations at mid-latitudes in 2003, serves as an upper limit of the comparison. Most of the balloon observations agree with the photochemical model predictions within their given error estimates. First retrieval exercises of BrO limb profiling from the SCIAMACHY satellite instrument agree to <±50% with the photochemically-corrected balloon observations, and tend to show less agreement below 20 km.
Pollen-based climate reconstructions were performed on two high-resolution pollen marines cores from the Alboran and Aegean Seas in order to unravel the climatic variability in the coastal settings of the Mediterranean region between 15 000 and 4000 years BP (the Lateglacial, and early to mid-Holocene). The quantitative climate reconstructions for the Alboran and Aegean Sea records focus mainly on the reconstruction of the seasonality changes (temperatures and precipitation), a crucial parameter in the Mediterranean region. This study is based on a multi-method approach comprising 3 methods: the Modern Analogues Technique (MAT), the recent Non-Metric Multidimensional Scaling/Generalized Additive Model method (NMDS/GAM) and Partial Least Squares regression (PLS). The climate signal inferred from this comparative approach confirms that cold and dry conditions prevailed in the Mediterranean region during the Oldest and Younger Dryas periods, while temperate conditions prevailed during the Bølling/Allerød and the Holocene. Our records suggest a West/East gradient of decreasing precipitation across the Mediterranean region during the cooler Late-glacial and early Holocene periods, similar to present-day conditions. Winter precipitation was highest during warm intervals and lowest during cooling phases. Several short-lived cool intervals (i.e. Older Dryas, another oscillation after this one (GI-1c2), Gerzensee/Preboreal Oscillations, 8.2 ka event, Bond events) connected to the North Atlantic climate system are documented in the Alboran and Aegean Sea records indicating that the climate oscillations associated with the successive steps of the deglaciation in the North Atlantic area occurred in both the western and eastern Mediterranean regions. This observation confirms the presence of strong climatic linkages between the North Atlantic and Mediterranean regions.
Pollen-based climate reconstructions were performed on two high-resolution pollen – marines cores from the Alboran and Aegean Seas in order to unravel the climatic variability in the coastal settings of the Mediterranean region between 15 000 and 4000 cal yrs BP (the Lateglacial, and early to mid-Holocene). The quantitative climate reconstructions for the Alboran and Aegean Sea records focus mainly on the reconstruction of the seasonality changes (temperatures and precipitation), a crucial parameter in the Mediterranean region. This study is based on a multi-method approach comprising 3 methods: the Modern Analogues Technique (MAT), the recent Non-Metric Multidimensional Scaling/Generalized Additive Model method (NMDS/GAM) and Partial Least Squares regression (PLS). The climate signal inferred from this comparative approach confirms that cold and dry conditions prevailed in the Mediterranean region during the Heinrich event 1 and Younger Dryas periods, while temperate conditions prevailed during the Bølling/Allerød and the Holocene. Our records suggest a West/East gradient of decreasing precipitation across the Mediterranean region during the cooler Late-glacial and early Holocene periods, similar to present-day conditions. Winter precipitation was highest during warm intervals and lowest during cooling phases. Several short-lived cool intervals (i.e., Older Dryas, another oscillation after this one (GI-1c2), Gerzensee/Preboreal Oscillations, 8.2 ka event, Bond events) connected to the North Atlantic climate system are documented in the Alboran and Aegean Sea records indicating that the climate oscillations associated with the successive steps of the deglaciation in the North Atlantic area occurred in both the western and eastern Mediterranean regions. This observation confirms the presence of strong climatic linkages between the North Atlantic and Mediterranean regions.
With the Earth system being about to leave Holocene conditions and thus the known safe operating space for humanity, frameworks such as the Planetary Boundaries (PBs) and the Sustainable Development Goals (SDGs) provide quantitative metrics to guide sustainability transformations. In order to strive, not only for compliance with the PBs but also for societal well-being, some approaches attempt to combine both PBs and SDGs within a single assessment. We focus on two prominent examples, the “Doughnut” by Kate Raworth and the #SDGinPB project of the 2018 report to the Club of Rome, which are not only aimed at public outreach, but also at guiding policy-making. To meet these objectives, the approaches should possess a certain accuracy in determining the progress in achieving the SDGs and in complying with the PBs. We evaluate, by using a multi-indicator approach for comparison, whether both approaches’ limited set of indicators can still represent the SDGs’ complexity. This comparative approach estimates the progress in achieving SDGs, especially in the Global North, to be significantly lower. Based on these results and against the approaches’ purposes, we discuss their simplifications and at which point the results are no longer reliable. We conclude that global assessments can be an important factor in initiating transformative processes by stimulating public discourse, but that the actual implementation of these would require approaches with greater recognition of local particularities.
In this study we investigate the scaling of precipitation extremes with temperature in the Mediterranean region by assessing against observations the present day and future regional climate simulations performed in the frame of the HyMeX and MED-CORDEX programs. Over the 1979–2008 period, despite differences in quantitative precipitation simulation across the various models, the change in precipitation extremes with respect to temperature is robust and consistent. The spatial variability of the temperature–precipitation extremes relationship displays a hook shape across the Mediterranean, with negative slope at high temperatures and a slope following Clausius–Clapeyron (CC)-scaling at low temperatures. The temperature at which the slope of the temperature–precipitation extreme relation sharply changes (or temperature break), ranges from about 20 °C in the western Mediterranean to <10 °C in Greece. In addition, this slope is always negative in the arid regions of the Mediterranean. The scaling of the simulated precipitation extremes is insensitive to ocean–atmosphere coupling, while it depends very weakly on the resolution at high temperatures for short precipitation accumulation times. In future climate scenario simulations covering the 2070–2100 period, the temperature break shifts to higher temperatures by a value which is on average the mean regional temperature change due to global warming. The slope of the simulated future temperature–precipitation extremes relationship is close to CC-scaling at temperatures below the temperature break, while at high temperatures, the negative slope is close, but somewhat flatter or steeper, than in the current climate depending on the model. Overall, models predict more intense precipitation extremes in the future. Adjusting the temperature–precipitation extremes relationship in the present climate using the CC law and the temperature shift in the future allows the recovery of the temperature–precipitation extremes relationship in the future climate. This implies negligible regional changes of relative humidity in the future despite the large warming and drying over the Mediterranean. This suggests that the Mediterranean Sea is the primary source of moisture which counteracts the drying and warming impacts on relative humidity in parts of the Mediterranean region.
Extensive black shale deposits formed in the Early Cretaceous South Atlantic, supporting the notion that this emerging ocean basin was a globally important site of organic carbon burial. The magnitude of organic carbon burial in marine basins is known to be controlled by various tectonic, oceanographic, hydrological, and climatic processes acting on different temporal and spatial scales, the nature and relative importance of which are poorly understood for the young South Atlantic. Here we present new bulk and molecular geochemical data from an Aptian–Albian sediment record recovered from the deep Cape Basin at Deep Sea Drilling Project (DSDP) Site 361, which we combine with general circulation model results to identify driving mechanisms of organic carbon burial. A multimillion-year decrease (i.e., Early Aptian–Albian) in organic carbon burial, reflected in a lithological succession of black shale, gray shale, and red beds, was caused by increasing bottom water oxygenation due to abating hydrographic restriction via South Atlantic–Southern Ocean gateways. These results emphasize basin evolution and ocean gateway development as a decisive primary control on enhanced organic carbon preservation in the Cape Basin at geological timescales (> 1 Myr). The Early Aptian black shale sequence comprises alternations of shales with high (> 6 %) and relatively low (∼ 3.5 %) organic carbon content of marine sources, the former being deposited during the global Oceanic Anoxic Event (OAE) 1a, as well as during repetitive intervals before and after OAE 1a. In all cases, these short-term intervals of enhanced organic carbon burial coincided with strong influxes of sediments derived from the proximal African continent, indicating closely coupled climate–land–ocean interactions. Supported by our model results, we show that fluctuations in weathering-derived nutrient input from the southern African continent, linked to changes in orbitally driven humidity and aridity, were the underlying drivers of repetitive episodes of enhanced organic carbon burial in the deep Cape Basin. These results suggest that deep marine environments of emerging ocean basins responded sensitively and directly to short-term fluctuations in riverine nutrient fluxes. We explain this relationship using the lack of wide and mature continental shelf seas that could have acted as a barrier or filter for nutrient transfer from the continent into the deep ocean.
Extensive black shale deposits formed in the Early Cretaceous South Atlantic, supporting the notion that this emerging ocean basin was a globally important site of organic carbon burial. The magnitude of organic carbon burial in marine basins is known to be controlled by various tectonic, oceanographic, hydrological, and climatic processes acting on different temporal and spatial scales, the nature and relative importance of which are poorly understood for the young South Atlantic. Here we present new bulk and molecular geochemical data from an Aptian–Albian sediment record recovered from the deep Cape Basin at Deep Sea Drilling Project (DSDP) Site 361, which we combine with general circulation model results to identify driving mechanisms of organic carbon burial. A multi-million year decrease (i.e. Early Aptian–Albian) in organic carbon burial, reflected in a lithological succession of black shale, gray shale, and red beds, was caused by increasing bottom water oxygenation due to abating tectonic restriction via South Atlantic-Southern Ocean gateways. These results emphasize basin evolution and ocean gateway development as a decisive primary control on enhanced organic carbon preservation in the Cape Basin at geological time scales (>1 Myr). The Early Aptian black shale sequence comprises alternations of shales with high (>5%) and relatively low (~3%) organic carbon content of marine sources, the former being deposited during the global Oceanic Anoxic Event (OAE) 1a, as well as during repetitive events before and after OAE 1a. In all cases, these short-term events of enhanced organic carbon burial coincided with strong influxes of sediments derived from the proximal African continent, indicating closely coupled climate–land–ocean interactions. Supported by our model results, we propose that fluctuations in weathering-derived nutrient input from the southern African continent, linked to fluctuations in pCO2 and/or orbitally driven humidity/aridity, were the underlying drivers of short-term organic carbon burial in the deep Cape Basin. These results suggest that deep marine environments of emerging ocean basins responded sensitively and directly to short term fluctuations in riverine nutrient fluxes. We explain this relationship by the lack of wide and mature continental shelf seas that could have acted as a barrier or filter for nutrient transfer from the continent into the deep ocean.
A hygroscopicity tandem differential mobility analyzer (HTDMA) was used to measure the water uptake (hygroscopicity) of secondary organic aerosol (SOA) formed during the chemical and photochemical oxidation of several organic precursors in a smog chamber. Electron ionization mass spectra of the non-refractory submicron aerosol were simultaneously determined with an aerosol mass spectrometer (AMS), and correlations between the two different signals were investigated. SOA hygroscopicity was found to strongly correlate with the relative abundance of the ion signal m/z 44 expressed as a fraction of total organic signal (f44). m/z 44 is due mostly to the ion fragment CO2+ for all types of SOA systems studied, and has been previously shown to strongly correlate with organic O/C for ambient and chamber OA. The analysis was also performed on ambient OA from two field experiments at the remote site Jungfraujoch, and the megacity Mexico City, where similar results were found. A simple empirical linear relation between the hygroscopicity of OA at subsaturated RH, as given by the hygroscopic growth factor (GF) or "ϰorg" parameter, and f44 was determined and is given by ϰorg = 2.2 × f44 − 0.13. This approximation can be further verified and refined as the database for AMS and HTDMA measurements is constantly being expanded around the world. The use of this approximation could introduce an important simplification in the parameterization of hygroscopicity of OA in atmospheric models, since f44 is correlated with the photochemical age of an air mass.
A hygroscopicity tandem differential mobility analyzer (HTDMA) was used to measure the water uptake (hygroscopicity) of secondary organic aerosol (SOA) formed during the chemical and photochemical oxidation of several organic precursors in a smog chamber. Electron ionization mass spectra of the non-refractory submicron aerosol were simultaneously determined with an aerosol mass spectrometer (AMS), and correlations between the two different signals were investigated. SOA hygroscopicity was found to strongly correlate with the relative abundance of the ion signal m/z 44 expressed as a fraction of total organic signal (f44). m/z 44 is due mostly to the ion fragment CO2+ for all types of SOA systems studied, and has been previously shown to strongly correlate with organic O/C for ambient and chamber OA. The analysis was also performed on ambient OA from two field experiments at the remote site Jungfraujoch, and the megacity Mexico City, where similar results were found. A simple empirical linear relation between the hygroscopicity of OA at subsaturated RH, as given by the hygroscopic growth factor (GF) or "κorg" parameter, and f44 was determined and is given by κorg=2.2×f44−0.13. This approximation can be further verified and refined as the database for AMS and HTDMA measurements is constantly being expanded around the world. The use of this approximation could introduce an important simplification in the parameterization of hygroscopicity of OA in atmospheric models, since f44 is correlated with the photochemical age of an air mass.
During a 4-week run in October–November 2006, a pilot experiment was performed at the CERN Proton Synchrotron in preparation for the CLOUD1 experiment, whose aim is to study the possible influence of cosmic rays on clouds. The purpose of the pilot experiment was firstly to carry out exploratory measurements of the effect of ionising particle radiation on aerosol formation from trace H2SO4 vapour and secondly to provide technical input for the CLOUD design. A total of 44 nucleation bursts were produced and recorded, with formation rates of particles above the 3 nm detection threshold of between 0.1 and 100 cm−3s−1, and growth rates between 2 and 37 nm h−1. The corresponding H2SO4 concentrations were typically around 106 cm−3 or less. The experimentally-measured formation rates and H2SO4 concentrations are comparable to those found in the atmosphere, supporting the idea that sulphuric acid is involved in the nucleation of atmospheric aerosols. However, sulphuric acid alone is not able to explain the observed rapid growth rates, which suggests the presence of additional trace vapours in the aerosol chamber, whose identity is unknown. By analysing the charged fraction, a few of the aerosol bursts appear to have a contribution from ion-induced nucleation and ion-ion recombination to form neutral clusters. Some indications were also found for the accelerator beam timing and intensity to influence the aerosol particle formation rate at the highest experimental SO2 concentrations of 6 ppb, although none was found at lower concentrations. Overall, the exploratory measurements provide suggestive evidence for ion-induced nucleation or ion-ion recombination as sources of aerosol particles. However in order to quantify the conditions under which ion processes become significant, improvements are needed in controlling the experimental variables and in the reproducibility of the experiments. Finally, concerning technical aspects, the most important lessons for the CLOUD design include the stringent requirement of internal cleanliness of the aerosol chamber, as well as maintenance of extremely stable temperatures (variations below 0.1°C).
During a 4-week run in October–November 2006, a pilot experiment was performed at the CERN Proton Synchrotron in preparation for the Cosmics Leaving OUtdoor Droplets (CLOUD) experiment, whose aim is to study the possible influence of cosmic rays on clouds. The purpose of the pilot experiment was firstly to carry out exploratory measurements of the effect of ionising particle radiation on aerosol formation from trace H2SO4 vapour and secondly to provide technical input for the CLOUD design. A total of 44 nucleation bursts were produced and recorded, with formation rates of particles above the 3 nm detection threshold of between 0.1 and 100 cm -3 s -1, and growth rates between 2 and 37 nm h -1. The corresponding H2O concentrations were typically around 106 cm -3 or less. The experimentally-measured formation rates and htwosofour concentrations are comparable to those found in the atmosphere, supporting the idea that sulphuric acid is involved in the nucleation of atmospheric aerosols. However, sulphuric acid alone is not able to explain the observed rapid growth rates, which suggests the presence of additional trace vapours in the aerosol chamber, whose identity is unknown. By analysing the charged fraction, a few of the aerosol bursts appear to have a contribution from ion-induced nucleation and ion-ion recombination to form neutral clusters. Some indications were also found for the accelerator beam timing and intensity to influence the aerosol particle formation rate at the highest experimental SO2 concentrations of 6 ppb, although none was found at lower concentrations. Overall, the exploratory measurements provide suggestive evidence for ion-induced nucleation or ion-ion recombination as sources of aerosol particles. However in order to quantify the conditions under which ion processes become significant, improvements are needed in controlling the experimental variables and in the reproducibility of the experiments. Finally, concerning technical aspects, the most important lessons for the CLOUD design include the stringent requirement of internal cleanliness of the aerosol chamber, as well as maintenance of extremely stable temperatures (variations below 0.1 °C)
The present work was devised to address the systematic analysis of samples from a range of Roman non-ferrous metal artefacts from different archaeological contexts and sites in the Roman provinces of Germania Superior. One of the focal points of this study is the provenancing of different lead objects from five important Roman settlements between 15 BC and the beginning of fourth century AD. For this purpose, measurements were made on lead and copper ore samples from the Siegerland, Eifel, Hunsrück and Lahn-Dill area in Germany and supplemented with data from the literature to create a data bank of lead isotope ratios of European deposits. Compositional analysis of lead objects by Electron Microprobe analysis showed that Romans were able to purify lead from ore up to 99%. Multi-Collector Inductively Coupled Plasma Mass-Spectrometry was used to determine the source of lead, which played an important role in nearly all aspects of Roman life. Lead isotope ratios were measured for ore samples from German deposits from the eastern side of the Rhine (Siegerland, Lahn-Dill, Ems) and the western side of the Rhine (Eifel, Hunsrück), which contained enough ore reserves to answer the increasing local demand and are believed to have been mined during the Roman period. This data together with those from Mediterranean ore deposits from the literature was used to establish a data bank. The Mediterranean ore deposits range from Cambrian (high 207Pb/206Pb) to tertiary (lower 207Pb/206Pb) values. In particular, the Cypriot deposits are younger, while the Spanish deposits fall either with the younger Sardic ores or close to the older Cypriot ores. The lead isotope ratios of most German ore deposits fall in between the 208Pb/206Pb vs. 207Pb/206Pb ratios of Sardinia and Cyprus, where the lead isotope signature of ore deposits from France and Britain are also found. Over 240 lead objects were measured from Wallendorf (second century BC to first century AD) Dangstetten (15-8 BC), Waldgirmes (AD 1-10), Mainz (AD 1-300), Martberg (first to fourth centuries AD) & Trier (third to fourth centuries AD). Comparing the lead isotope ratios of lead objects and those from German ores shows that the source of over 85 percent of objects are Eifel ore deposits, but the Roman’s had also imported lead from the Southern Massif Central and from Great Britain. A further topic of this work was the systematic study of the variation of copper isotope ratios in different copper minerals and the mechanisms, which controls copper isotope fractionation in ores deposits. For this purpose, copper isotope analyses were made by Multi-Collector Inductively Coupled Plasma Mass-Spectrometry from a series of hydrothermal copper sulphides and their alteration products. Copper and lead isotope ratios were measured in coexisting phases of chalcopyrite and malachite and also coexisting malachite and azurite. No significant fractionation was observed in malachite-azurite phases, but in chalcopyrite-malachite coexisting phases, malachite always shows a positive fractionation to heavier isotope values. Zhu et al. and Larson et al. showed that isotopic variations in copper principally reflect mass fractionation in response to low temperature processes rather than source heterogeneity. The low temperature ore formation processes are mostly represented by weathering of primary sulphide ores to produce secondary carbonate phases and therefore are usually observed on the surface of ore deposits, which were probably removed during the early Bronze Age. Using this concept, copper isotope ratios were measured in some Early Bronze Age copper alloys and Roman copper alloys. However, no large copper isotope fractionation has been observed. Lead and copper isotope ratios were measured on samples from the Kupferschiefer. Two profiles were investigated; 1) Sangerhausen, which was not directly influenced by the oxidizing brines of Rote Fäule and 2) Oberkatz, where both Rote Fäule-controlled and structure-controlled mineralization were observed. Results from maturation studies of organic matter suggest the maximum temperature affecting the Kupferschiefer did not exceed 130°C. delta-65-Cu ranges between -0.78-+0.58‰, shows a positive correlation with copper concentration. Maximum temperature in the Kupferschiefer profile from Oberkatz is supposed to be around 150°C. delta-65Cu in this profile ranges between -0.71-+0.68‰. The pattern of copper isotope fractionation and copper concentration is same as the for profile of Sangerhausen. Origina lead isotope ratios are strongly overprinted by high concentrations of uranium in bottom of both profiles causing more radiogenic lead.
Large-scale hydrological modelling has become increasingly wide-spread during the last decade. An annual workshop series on large-scale hydrological modelling has provided, since 1997, a forum to the German-speaking community for discussing recent developments and achievements in this research area. In this paper we present the findings from the 2007 workshop which focused on advances and visions in large-scale hydrological modelling. We identify the state of the art, difficulties and research perspectives with respect to the themes "sensitivity of model results", "integrated modelling" and "coupling of processes in hydrosphere, atmosphere and biosphere". Some achievements in large-scale hydrological modelling during the last ten years are presented together with a selection of remaining challenges for the future.
Quantification of spatially and temporally resolved water flows and water storage variations for all land areas of the globe is required to assess water resources, water scarcity and flood hazards, and to understand the Earth system. This quantification is done with the help of global hydrological models (GHMs). What are the challenges and prospects in the development and application of GHMs? Seven important challenges are presented. (1) Data scarcity makes quantification of human water use difficult even though significant progress has been achieved in the last decade. (2) Uncertainty of meteorological input data strongly affects model outputs. (3) The reaction of vegetation to changing climate and CO2 concentrations is uncertain and not taken into account in most GHMs that serve to estimate climate change impacts. (4) Reasons for discrepant responses of GHMs to changing climate have yet to be identified. (5) More accurate estimates of monthly time series of water availability and use are needed to provide good indicators of water scarcity. (6) Integration of gradient-based groundwater modelling into GHMs is necessary for a better simulation of groundwater–surface water interactions and capillary rise. (7) Detection and attribution of human interference with freshwater systems by using GHMs are constrained by data of insufficient quality but also GHM uncertainty itself. Regarding prospects for progress, we propose to decrease the uncertainty of GHM output by making better use of in situ and remotely sensed observations of output variables such as river discharge or total water storage variations by multi-criteria validation, calibration or data assimilation. Finally, we present an initiative that works towards the vision of hyperresolution global hydrological modelling where GHM outputs would be provided at a 1-km resolution with reasonable accuracy.
Long-term average groundwater recharge, which is equivalent to renewable groundwater resources, is the major limiting factor for the sustainable use of groundwater. Compared to surface water resources, groundwater resources are more protected from pollution, and their use is less restricted by seasonal and inter-annual flow variations. To support water management in a globalized world, it is necessary to estimate groundwater recharge at the global scale. Here, we present a best estimate of global-scale long-term average diffuse groundwater recharge (i.e. renewable groundwater resources) that has been calculated by the most recent version of the WaterGAP Global Hydrology Model WGHM (spatial resolution of 0.5° by 0.5°, daily time steps). The estimate was obtained using two state-of-the-art global data sets of gridded observed precipitation that we corrected for measurement errors, which also allowed to quantify the uncertainty due to these equally uncertain data sets. The standard WGHM groundwater recharge algorithm was modified for semi-arid and arid regions, based on independent estimates of diffuse groundwater recharge, which lead to an unbiased estimation of groundwater recharge in these regions. WGHM was tuned against observed long-term average river discharge at 1235 gauging stations by adjusting, individually for each basin, the partitioning of precipitation into evapotranspiration and total runoff. We estimate that global groundwater recharge was 12 666 km3/yr for the climate normal 1961–1990, i.e. 32% of total renewable water resources. In semi-arid and arid regions, mountainous regions, permafrost regions and in the Asian Monsoon region, groundwater recharge accounts for a lower fraction of total runoff, which makes these regions particularly vulnerable to seasonal and inter-annual precipitation variability and water pollution. Average per-capita renewable groundwater resources of countries vary between 8 m3/(capita yr) for Egypt to more than 1 million m3/(capita yr) for the Falkland Islands, the global average in the year 2000 being 2091 m3/(capita yr). Regarding the uncertainty of estimated groundwater resources due to the two precipitation data sets, deviation from the mean is 1.1% for the global value, and less than 1% for 50 out of the 165 countries considered, between 1 and 5% for 62, between 5 and 20% for 43 and between 20 and 80% for 10 countries. Deviations at the grid scale can be much larger, ranging between 0 and 186 mm/yr.