550 Geowissenschaften
Refine
Year of publication
- 2021 (92)
- 2016 (78)
- 2020 (58)
- 2022 (53)
- 2019 (48)
- 2011 (44)
- 2015 (44)
- 2010 (42)
- 2013 (42)
- 2014 (40)
- 2009 (39)
- 2017 (38)
- 2008 (35)
- 2012 (31)
- 2018 (31)
- 2005 (20)
- 2023 (20)
- 2007 (15)
- 2006 (11)
- 2003 (10)
- 2024 (10)
- 2002 (9)
- 2004 (5)
- 2001 (4)
- 1996 (2)
- 1838 (1)
- 1871 (1)
- 1885 (1)
- 1904 (1)
- 1905 (1)
- 1908 (1)
- 1910 (1)
- 1911 (1)
- 1921 (1)
- 1925 (1)
- 1929 (1)
- 1958 (1)
- 1961 (1)
- 1962 (1)
- 1964 (1)
- 1969 (1)
- 1973 (1)
- 1974 (1)
- 1976 (1)
- 1981 (1)
- 1993 (1)
- 1998 (1)
- 2000 (1)
Document Type
- Article (685)
- 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 (845) (remove)
Has Fulltext
- yes (845)
Is part of the Bibliography
- no (845)
Keywords
- climate change (11)
- Climate change (7)
- COSMO-CLM (6)
- Palaeoclimate (6)
- precipitation (6)
- Atmospheric chemistry (5)
- Biogeochemistry (5)
- Palaeoceanography (5)
- Geochemistry (4)
- Bewässerung (3)
Institute
- Geowissenschaften (619)
- Geowissenschaften / Geographie (134)
- Biodiversität und Klima Forschungszentrum (BiK-F) (62)
- Senckenbergische Naturforschende Gesellschaft (47)
- Geographie (22)
- Biowissenschaften (21)
- Institut für Ökologie, Evolution und Diversität (8)
- Institut für sozial-ökologische Forschung (ISOE) (5)
- Frankfurt Institute for Advanced Studies (FIAS) (4)
- Medizin (4)
Chemical ozone loss in winter 1991–1992 is recalculated based on observations of the HALOE satellite instrument, ER-2 aircraft measurements and balloon data. HALOE satellite observations are shown to be reliable in the lower stratosphere below 400 K, at altitudes where profiles are most likely disturbed by the enhanced sulfate aerosols, as a result of the Mt. Pinatubo eruption in June 1991. Very large chemical ozone loss was observed below 400 K from Kiruna balloon observations between December and March 1992. Additionally, for the two winters after the Mt. Pinatubo eruption, HALOE satellite observations show a stronger extent of chemical ozone loss at lower altitudes compared to other Arctic winter between 1991 and 2003. In stipe of already occurring deactivation of chlorine in March 1992, Mipas-B and LPMA balloon observations indicate still chlorine activation at lower altitudes, consistent with observed chemical ozone loss occurring between February and March and April. Enhanced chemical ozone loss in the Arctic winter 1991–1992 as calculated in earlier studies is corroborated here.
Reliable identification of chondrules, calcium-aluminum-rich inclusions (CAIs), carbonate grains, and Ca-phosphate grains at depth within untouched, unprepared chondritic samples by a nondestructive analytical method, such as synchrotron X-ray fluorescence (SXRF) computed tomography (CT), is an essential first step before intrusive analytical and sample preparation methods are performed. The detection of a local Ca-enrichment could indicate the presence of such a component, all of which contain Ca as major element and/or Ca-bearing minerals, allowing it to be precisely located at depth within a sample. However, the depth limitation from which Ca-K fluorescence can travel through a chondrite sample (e.g., ∼115 µm through material of 1.5 g cm−3) to XRF detectors leaves many Ca-bearing components undetected at deeper depths. In comparison, Sr-K lines travel much greater distances (∼1700 µm) through the same sample density and are, thus, detected from much greater depths. Here, we demonstrate a clear, positive, and preferential correlation between Ca and Sr and conclude that Sr-detection can be used as proxy for the presence of Ca (and, thus, Ca-bearing components) throughout mm-sized samples of carbonaceous chondritic material. This has valuable implications, especially for sample return missions from carbonaceous C-type asteroids, such as Ryugu or Bennu. Reliable localization, identification, and targeted analysis by SXRF of Ca-bearing chondrules, CAIs, and carbonates at depth within untouched, unprepared samples in the initial stages of a multianalysis investigation insures the valuable information they hold of pre- and post-accretion processes in the early solar system is neither corrupted nor destroyed in subsequent processing and analyses.
Ecophysiological studies on Antarctic cryptophytes to assess whether climatic changes such as ocean acidification and enhanced stratification affect their growth in Antarctic coastal waters in the future are lacking so far. This is the first study that investigated the combined effects of increasing availability of pCO2 (400 and 1000 µatm) and irradiance (20, 200 and 500 μmol photons m−2 s −1) on growth, elemental composition and photophysiology of the Antarctic cryptophyte Geminigera cryophila. Under ambient pCO2, this species was characterized by a pronounced sensitivity to increasing irradiance with complete growth inhibition at the highest light intensity. Interestingly, when grown under high pCO2 this negative light effect vanished and it reached highest rates of growth and particulate organic carbon production at the highest irradiance compared to the other tested experimental conditions. Our results for G. cryophila reveal beneficial effects of ocean acidification in conjunction with enhanced irradiance on growth and photosynthesis. Hence, cryptophytes such as G. cryophila may be potential winners of climate change, potentially thriving better in more stratified and acidic coastal waters and contributing in higher abundance to future phytoplankton assemblages of coastal Antarctic waters.
About half of present-day cloud condensation nuclei originate from atmospheric nucleation, frequently appearing as a burst of new particles near midday1. Atmospheric observations show that the growth rate of new particles often accelerates when the diameter of the particles is between one and ten nanometres2,3. In this critical size range, new particles are most likely to be lost by coagulation with pre-existing particles4, thereby failing to form new cloud condensation nuclei that are typically 50 to 100 nanometres across. Sulfuric acid vapour is often involved in nucleation but is too scarce to explain most subsequent growth5,6, leaving organic vapours as the most plausible alternative, at least in the planetary boundary layer7,8,9,10. Although recent studies11,12,13 predict that low-volatility organic vapours contribute during initial growth, direct evidence has been lacking. The accelerating growth may result from increased photolytic production of condensable organic species in the afternoon2, and the presence of a possible Kelvin (curvature) effect, which inhibits organic vapour condensation on the smallest particles (the nano-Köhler theory)2,14, has so far remained ambiguous. Here we present experiments performed in a large chamber under atmospheric conditions that investigate the role of organic vapours in the initial growth of nucleated organic particles in the absence of inorganic acids and bases such as sulfuric acid or ammonia and amines, respectively. Using data from the same set of experiments, it has been shown15 that organic vapours alone can drive nucleation. We focus on the growth of nucleated particles and find that the organic vapours that drive initial growth have extremely low volatilities (saturation concentration less than 10−4.5 micrograms per cubic metre). As the particles increase in size and the Kelvin barrier falls, subsequent growth is primarily due to more abundant organic vapours of slightly higher volatility (saturation concentrations of 10−4.5 to 10−0.5 micrograms per cubic metre). We present a particle growth model that quantitatively reproduces our measurements. Furthermore, we implement a parameterization of the first steps of growth in a global aerosol model and find that concentrations of atmospheric cloud concentration nuclei can change substantially in response, that is, by up to 50 per cent in comparison with previously assumed growth rate parameterizations.
Derivation and characterization of a new filter for nonlinear high-dimensional data assimilation
(2015)
Data assimilation (DA) combines model forecasts with real-world observations to achieve an optimal estimate of the state of a dynamical system. The quality of predictions in nonlinear and chaotic systems such as atmospheric or oceanic circulation is strongly sensitive to the initial conditions. Therefore, beyond the consistent reconstruction of past states, a primary relevance of advanced DA methods concerns the proper model initialization. The ensemble Kalman filter (EnKF) and its deterministic variants, mostly square root filters such as the ensemble transform Kalman filter (ETKF), represent a popular alternative to variational DA schemes. They are applied in a wide range of research and operations. Their forecast step employs an ensemble integration that fully respects the nonlinear nature of the analyzed system. In the analysis step, they implicitly assume the prior state and observation errors to be Gaussian. Consequently, in nonlinear systems, the mean and covariance of the analysis ensemble are biased and these filters remain suboptimal. In contrast, the fully nonlinear, non-Gaussian particle filter (PF) relies on Bayes' theorem without further assumptions, which guarantees an exact asymptotic behavior. However, it is exposed to weight collapse, particularly in higher-dimensional settings, known as the curse of dimensionality.
This work presents a new method to obtain an analysis ensemble with mean and covariance that exactly match the corresponding Bayesian estimates. This is achieved by a deterministic matrix square root transformation of the forecast ensemble, and subsequently a suitable random rotation that significantly contributes to filter stability while preserving the required second-order statistics. The forecast step remains as in the ETKF. The algorithm, which is fairly easy to implement and computationally efficient, is referred to as the nonlinear ensemble transform filter (NETF). The limitation with respect to fully-nonlinear filtering is that the NETF only considers the mean and covariance of the Bayesian analysis density, neglecting higher-order moments.
The properties and performance of the proposed algorithm are investigated via a set of experiments. The results indicate that such a filter formulation can increase the analysis quality, even for relatively small ensemble sizes, compared to other ensemble filters in nonlinear, non-Gaussian scenarios. They also confirm that localization enhances the applicability of this PF-inspired scheme in larger-dimensional systems. Finally, the novel filter is coupled to a large-scale ocean general circulation model with a realistic observation scenario. The NETF remains stable with a small ensemble size and shows a consistent behavior. Additionally, its analyses exhibit low estimation errors, as revealed by a comparison with a free ensemble integration and the ETKF. The results confirm that, in principle, the filter can be applied successfully and as simple as the ETKF in high-dimensional problems. No further modifications are needed, even though the algorithm is only based on the particle weights. Thus, it is able to overcome the curse of dimensionality, even in deterministic systems. This proves that the NETF constitutes a promising and user-friendly method for nonlinear high-dimensional DA.
Long-term average groundwater recharge representing the sustainable groundwater resources is modeled as a 0.5° by 0.5° grid on global scale by the WaterGAP Global Hydrology Model. Due to uncertainties of estimating groundwater recharge, especially in semiarid and arid regions, independent estimates are used for calibrating the model. This work compiled a new set of independent groundwater recharge estimates based on a work of Scanlon et al. (2006). The 59 independent estimates, together with an already existing independent estimates compilation, are used for the evaluation of two WGHM variants; one variant is modeling with an improved more realistically distributed daily precipitation dataset.
The objective of this thesis is the evaluation of the modeled data of the WaterGAP Global Hydrology Model (WGHM). The analysis of the impact of the new Watch Forcing Data (WFD) precipitation dataset on the modeled groundwater recharge tends to result in lower values in humid and higher values in (semi-)arid regions compared to the WGHM standard variant. Comparing both WGHM variants to the independent estimates compilations, representing (semi-)arid regions, the WGHM variant shows over- and underestimations especially of the low values and the WGHM WFD variant shows a bias toward overestimation especially for values below 4 mm/yr. The analysis of texture, hydrogeology and vegetation/ land cover could not give satisfying explanations for the discrepancies, but derived from the groundwater recharge measurement methods analysis indirect/ localized recharge seems to be a significant factor causing underestimations, as resulted in the comparison of the independent estimates based on Scanlon et al. 2006 with the WGHM variants.
Over the last several decades, spinel-structured minerals with the chemical formula AB2O4 (where A and B stand for divalent and trivalent cations, respectively) have attracted more and more attention, particularly with regards to their breakdown at high pressures and temperatures and the nature of the so-called "post-spinel" phases. Spinel-structured phases with different endmember compositions, like magnetite (Fe3O4), hercynite (FeAl2O4) or spinel (MgAl2O4), are known to breakdown differently at high pressure-temperature conditions (e.g., Akaogi et al. 1999; Schollenbruch et al. 2010; Woodland et al. 2012). Such phases are of particular interest when they incorporate ferric (Fe3+) and ferrous (Fe2+) cations as this makes their stability sensitive to redox conditions. Since magnetite and magnesioferrite (MgFe3+ 2O4) have been found as inclusions in diamond (e.g., Stachel et al. 1998; Harte et al. 1999; Wirth et al. 2014; Palot et al. 2016; Jacob et al. 2016), understanding their phase relations is important for setting constraints on the conditions of their formation.
This study aimed to experimentally investigate the phase relations of Fe-Mg spinel-structured phases at conditions of the deep upper mantle and transition zone. Exploring the stability of new post-spinel phases and their characterization were also major goals of this study. Approaching a pyrolitic mantle composition by adding amounts of SiO2 in the system allowed constraints on the relevance of Fe-Mg post-spinel phases coexisting with mantle silicates to be made. ...
As part of two drilling campaigns of the International Continental Scientific Drilling Program (ICDP), several geophysical borehole measurements were carried out by the Leibniz Institute for Applied Geophysics (LIAG) in two lakes. The acquired data was used to answer stratigraphic and paleoclimatic research questions, including the establishment of robust age-depth models and the construction of continuous lithological profiles.
Lake Towuti is located on Sulawesi (Indonesia), within the "Indo-Pacific Warm Pool" (IPWP), a globally important region for atmospheric heat and moisture budgets. The lake exists for approximately one million years, but its exact age is uncertain. We present the first agedepth model for the approximately 100 m continuous sediment sequence from the central part of the lake. The basis for this model is the magnetic susceptibility measured in the borehole and a tephra layer with an age of about 797 ka at 72 m depth. Our age-depth model is inferred from cyclostratigraphic analysis of borehole data and covers a period from 903 ± 11 to 131 ± 67 ka. We suggest that orbital eccentricity and/or changes between global cold and warm periods are responsible for hydroclimatic changes in the IPWP, that these changes affect sedimentation processes in Lake Towuti, and that we can measure and observe this effect in the sediment properties today. Additionally, we created a continuous artificial lithological profile from a series of different borehole data using cluster analysis. This provides information from parts of the borehole where no sediment is available due to core loss.
Lake Ohrid is 1.36 million years old and is located on the Balkan Peninsula on the border between Albania and North Macedonia. The primary hole 'DEEP' in the central part of the lake has been the subject of several investigations, but information about sediments of the marginal locations 'Pestani' and 'Cerava' have not been published yet. In our study, we use natural gamma radiation (GR) measured in the borehole to generate an age-depth model for DEEP. This is performed using the correlation of GR to the global LR04 reference record of Lisiecki and Raymo (2005).
The age information is then transferred via prominent seismic marker horizons to the other two sites, Pestani and Cerava, where it provides the first age-control points for the construction of age-depth models from correlation of GR to LR04. The generated age-depth models are tested using cyclostratigraphic methods, but the limits of this approach are revealed. At DEEP, sedimentation rates (SR) from the cyclostratigraphic method and the correlative approach differ by 2.8 %, at Pestani this difference is 16.7 %, and at Cerava the quality of the data does not allow a reliable evaluation of SR using the cyclostratigraphic approach. We used cluster analysis to construct artificial lithological profiles at all three sites and integrated them into the respective age-depth models. This enables us to determine which sediment types were deposited at what time, and we recognize the change between warm and cold periods in the sediment properties at all three locations. The analyses in this study were all performed on borehole and seismic data and thus do not involve sediment core data. Especially at Pestani and Cerava, new insights into the sedimentological history of Lake Ohrid could be obtained.
In the last part we discuss the occurrence of the half-precession (HP) signal in the European region during the last one million years. The focus is on Lake Ohrid, but a range of other proxies, from the eastern Mediterranean, across the European continent, up to Greenland are analyzed in regards to HP. Applying filters, we focus on the frequency range with a period of 13-8.5 ka and only HP remains in the records. We use correlative methods to determine the clarity of the HP signal in proxies distributed across the European realm. Additionally, we determined the development of HP over time. The HP signal is clearest in the southeast and decreases toward the north. It is further more pronounced in interglacial periods and in the younger part (<621 ka) of most proxies. We suggest that there are mechanisms that transmit the HP signal from its origin near the equator to higher latitudes via different processes. In this context, for instance, the African monsoon, the Nile River and the Mediterranean outflow via the Strait of Gibraltar can be important factors.
Analysing the composition of ambient ultrafine particles (UFPs) is a challenging task due to the low mass and chemical complexity of small particles, yet it is a prerequisite for the identification of particle sources and the assessment of potential health risks. Here, we show the molecular characterization of UFPs, based on cascade impactor (Nano-MOUDI) samples that were collected at an air quality monitoring station near one of Europe's largest airports, in Frankfurt, Germany. At this station, particle-size-distribution measurements show an enhanced number concentration of particles smaller than 50 nm during airport operating hours. We sampled the lower UFP fraction (0.010–0.018, 0.018–0.032, 0.032–0.056 µm) when the air masses arrived from the airport. We developed an optimized filter extraction procedure using ultra-high-performance liquid chromatography (UHPLC) for compound separation and a heated electrospray ionization (HESI) source with an Orbitrap high-resolution mass spectrometer (HRMS) as a detector for organic compounds. A non-target screening detected ∼200 organic compounds in the UFP fraction with sample-to-blank ratios larger than 5. We identified the largest signals as homologous series of pentaerythritol esters (PEEs) and trimethylolpropane esters (TMPEs), which are base stocks of aircraft lubrication oils. We unambiguously attribute the majority of detected compounds to jet engine lubrication oils by matching retention times, high-resolution and accurate mass measurements, and comparing tandem mass spectrometry (MS2) fragmentation patterns between both ambient samples and commercially available jet oils. For each UFP stage, we created molecular fingerprints to visualize the complex chemical composition of the organic fraction and their average carbon oxidation state. These graphs underline the presence of the homologous series of PEEs and TMPEs and the appearance of jet oil additives (e.g. tricresyl phosphate, TCP). Targeted screening of TCP confirmed the absence of the harmful tri-ortho isomer, while we identified a thermal transformation product of TMPE-based lubrication oil (trimethylolpropane phosphate, TMP-P). Even though a quantitative determination of the identified compounds is limited, the presented method enables the qualitative detection of molecular markers for jet engine lubricants in UFPs and thus strongly improves the source apportionment of UFPs near airports.
Analysing the composition of ambient ultrafine particles (UFP) is a challenging task due to the low mass and chemical complexity of small particles, yet it is a prerequisite for the identification ofparticle sources and the assessment of potential health risks. Here, we show the molecular characterization of UFP, based on cascade impactor (Nano-MOUDI) 10samples that were collected at an air quality monitoring station nearby one of Europe`s largest airports in Frankfurt, Germany. At this station, particle-size-distribution measurements show enhanced number concentration of particles smaller than 50nm during airport operating hours. We sampled the lower UFP fraction (0.010-0.018 μm; 0.018-0.032 μm; 0.032-0.056 μm) when the air masses arrived from the airport. We developed an optimized filter extraction procedure, used ultra-high performance liquid chromatography (UHPLC) for compound separation, and a heated electrospray ionization (HESI) source with an 15Orbitrap high-resolution mass spectrometer (HRMS) as a detector for organic compounds. A non-target screening detected ~200 organic compounds in the UFP fraction with sample-to-blank ratios larger than five. We identified the largest signals as homologous series of pentaerythritol esters (PEE) and trimethylolpropane esters (TMPE), which are base stocks of aircraft lubrication oils. We unambiguously attribute the majority of detected compounds to jet engine lubrication oils by matching retention times, high-resolution/accurate mass (HR/AM) measurements, and comparing MS/MS fragmentation patterns between both ambient samples and commercially available jet oils. For each UFP stage, we created molecular fingerprints to visualize the complex chemical composition ofthe organic fraction and their average carbon oxidation state. These graphs underline the presence of the homologous series of PEE and TMPE, and the appearance of jet oil additives (e.g. tricresyl phosphate (TCP)). Targeted screening on TCP confirmed the absence of the harmful tri-orthoisomer, while we identified a thermal transformation product of TMPE-based lubrication oil (trimethylolpropane phosphate (TMP-P)). Even though a quantitative determination of the identified compounds is limited, the presented method enables the qualitative detection of molecular markers for jet engine lubricants in UFP and thus strongly improves the source apportionment of UFP near airports.
Modelling short-term variability in carbon and water exchange in a temperate Scots pine forest
(2015)
The vegetation–atmosphere carbon and water exchange at one particular site can strongly vary from year to year, and understanding this interannual variability in carbon and water exchange (IAVcw) is a critical factor in projecting future ecosystem changes. However, the mechanisms driving this IAVcw are not well understood. We used data on carbon and water fluxes from a multi-year eddy covariance study (1997–2009) in a Dutch Scots pine forest and forced a process-based ecosystem model (Lund–Potsdam–Jena General Ecosystem Simulator; LPJ-GUESS) with local data to, firstly, test whether the model can explain IAVcw and seasonal carbon and water exchange from direct environmental factors only. Initial model runs showed low correlations with estimated annual gross primary productivity (GPP) and annual actual evapotranspiration (AET), while monthly and daily fluxes showed high correlations. The model underestimated GPP and AET during winter and drought events. Secondly, we adapted the temperature inhibition function of photosynthesis to account for the observation that at this particular site, trees continue to assimilate at very low atmospheric temperatures (up to daily averages of −10 °C), resulting in a net carbon sink in winter. While we were able to improve daily and monthly simulations during winter by lowering the modelled minimum temperature threshold for photosynthesis, this did not increase explained IAVcw at the site. Thirdly, we implemented three alternative hypotheses concerning water uptake by plants in order to test which one best corresponds with the data. In particular, we analyse the effects during the 2003 heatwave. These simulations revealed a strong sensitivity of the modelled fluxes during dry and warm conditions, but no single formulation was consistently superior in reproducing the data for all timescales and the overall model–data match for IAVcw could not be improved. Most probably access to deep soil water leads to higher AET and GPP simulated during the heatwave of 2003. We conclude that photosynthesis at lower temperatures than assumed in most models can be important for winter carbon and water fluxes in pine forests. Furthermore, details of the model representations of water uptake, which are often overlooked, need further attention, and deep water access should be treated explicitly.
Modelling short-term variability in carbon and water exchange in a temperate Scots pine forest
(2015)
The vegetation–atmosphere carbon and water exchange at one particular site can strongly vary from year to year, and understanding this interannual variability in carbon and water exchange (IAVcw) is a critical factor in projecting future ecosystem changes. However, the mechanisms driving this IAVcw are not well understood. We used data on carbon and water fluxes from a multi-year eddy covariance study (1997–2009) in a Dutch Scots pine forest and forced a process-based ecosystem model (Lund–Potsdam–Jena General Ecosystem Simulator; LPJ-GUESS) with local data to, firstly, test whether the model can explain IAVcw and seasonal carbon and water exchange from direct environmental factors only. Initial model runs showed low correlations with estimated annual gross primary productivity (GPP) and annual actual evapotranspiration (AET), while monthly and daily fluxes showed high correlations. The model underestimated GPP and AET during winter and drought events. Secondly, we adapted the temperature inhibition function of photosynthesis to account for the observation that at this particular site, trees continue to assimilate at very low atmospheric temperatures (up to daily averages of −10 °C), resulting in a net carbon sink in winter. While we were able to improve daily and monthly simulations during winter by lowering the modelled minimum temperature threshold for photosynthesis, this did not increase explained IAVcw at the site. Thirdly, we implemented three alternative hypotheses concerning water uptake by plants in order to test which one best corresponds with the data. In particular, we analyse the effects during the 2003 heatwave. These simulations revealed a strong sensitivity of the modelled fluxes during dry and warm conditions, but no single formulation was consistently superior in reproducing the data for all timescales and the overall model–data match for IAVcw could not be improved. Most probably access to deep soil water leads to higher AET and GPP simulated during the heatwave of 2003. We conclude that photosynthesis at lower temperatures than assumed in most models can be important for winter carbon and water fluxes in pine forests. Furthermore, details of the model representations of water uptake, which are often overlooked, need further attention, and deep water access should be treated explicitly.
Chlorine monoxide (ClO) plays a key role in stratospheric ozone loss processes at midlatitudes. We present two balloonborne in situ measurements of ClO conducted in northern hemisphere midlatitudes during the period of the maximum of total inorganic chlorine loading in the atmosphere. Both ClO measurements were conducted on board the TRIPLE balloon payload, launched in November 1996 in Le´on, Spain, and in May 1999 in Aire sur l’Adour, France. For both flights a ClO daylight and night time vertical profile could be derived over an altitude range of approximately 15–31 km. ClO mixing ratios are compared to model simulations performed with the photochemical box model version of the Chemical Lagrangian Model of the Stratosphere (CLaMS). Simulations along 24-h backward trajectories were performed to study the diurnal variation of ClO in the midlatitude lower stratosphere. Model simulations for the flight launched in Aire sur l’Adour 1999 show a good agreement with the ClO measurements. For the flight launched in Le´on 1996, a similar good agreement is found, except at around ~ 650 K potential temperature (~26km altitude). However, a tendency is found that for solar zenith angles greater than 86°–87° the simulated ClO mixing ratios substantially overestimate measured ClO by approximately a factor of 2.5 or more for both flights. Therefore we conclude that no indication can be deduced from the presented ClO measurements that substantial uncertainties exist in midlatitude chlorine chemistry of the stratosphere. An exception is the situation at solar zenith angles greater than 86°–87° where model simulations substantial overestimate ClO observations.
Chlorine monoxide (ClO) plays a key role in stratospheric ozone loss processes at midlatitudes. We present two balloon-borne in situ measurements of ClO conducted in northern hemisphere midlatitudes during the period of the maximum of total inorganic chlorine loading in the atmosphere. Both ClO measurements were conducted on board the TRIPLE balloon payload, launched in November 1996 in León, Spain, and in May 1999 in Aire sur l'Adour, France. For both flights a ClO daylight and night-time vertical profile was derived over an altitude range of approximately 15-35 km. ClO mixing ratios are compared to model simulations performed with the photochemical box model version of the Chemical Lagrangian Model of the Stratosphere (CLaMS). Simulations along 24-hour backward trajectories were performed to study the diurnal variation of ClO in the midlatitude lower stratosphere. Model simulations for the flight launched in Aire sur l'Adour 1999 show an excellent agreement with the ClO measurements. For the flight launched in León 1996, an overall good agreement is found, whereas the flight is characterized by a more complex dynamical situation due to a possible mixture of vortex and non-vortex air. We note that for both flights at solar zenith angles greater than 86°-87° simulated ClO mixing ratios are higher than observed ClO mixing ratios. However, the present findings indicate that no substantial uncertainties exist in midlatitude chlorine chemistry of the stratosphere.
Central and western Europe were affected by a compressional tectonic event in the Late Cretaceous, caused by the convergence of Iberia and Europe. Basement uplifts, inverted graben structures, and newly formed marginal troughs are the main expressions of crustal shortening. Although the maximum activity occurred during a short period of time between 90 and 75 Ma, the exact timing of this event is still unclear. Dating of the start and end of Late Cretaceous basin inversion gives very different results depending on the method applied. On the basis of borehole data, facies, and thickness maps, the timing of basin reorganization was reconstructed for several basins in central Europe. The obtained data point to a synchronous start of basin inversion at 95 Ma (Cenomanian), 5 Myr earlier than commonly assumed. The end of the Late Cretaceous compressional event is difficult to pinpoint in central Europe, because regional uplift and salt migration disturb the signal of shifting marginal troughs. Late Campanian to Paleogene strata deposited unconformably on inverted structures indicate slowly declining uplift rates during the latest Cretaceous. The differentiation of separate Paleogene inversion phases in central Europe does not appear possible at present.
Central Europe was affected by a compressional tectonic event in the Late Cretaceous, caused by the convergence of Iberia and Europe. Basement uplifts, inverted graben structures and newly formed marginal troughs are the main expressions of crustal shortening. Although the maximum activity occurred in a short period between 90 and 75 Ma, the exact timing of this event is still unclear. Dating of start and end of basin inversion is very different depending on the applied method. On the basis of borehole data, facies and thickness maps, the timing of basin re-organisation was reconstructed for several basins in Central Europe. The obtained data point to a synchronous start of basin inversion already at 95 Ma (Cenomanian), 5 Million years earlier than commonly assumed. The end of the Late Cretaceous compressional event is more difficult to pinpoint, because regional uplift and salt migration disturb the signal of shifting marginal troughs. Unconformities of Late Campanian to Paleogene age on inverted structures indicate slowly declining uplift rates.
To study the effect of galactic cosmic rays on aerosols and clouds, the Cosmics Leaving OUtdoor Droplets (CLOUD) project was established. Experiments are carried out at a 26.1 m3 tank at CERN (Switzerland). In the experiments, the effect of ionizing radiation on H2SO4 particle formation and growth is investigated. To evaluate the experimental configuration, the experiment was simulated using a coupled multidimensional computational fluid dynamics (CFD) – particle model. In the model the coupled fields of gas/vapor species, temperature, flow velocity and particle properties were computed to investigate mixing state and mixing times of the CLOUD tank's contents. Simulation results show that a 1-fan configuration, as used in first experiments, may not be sufficient to ensure a homogeneously mixed chamber. To mix the tank properly, two fans and sufficiently high fan speeds are necessary. The 1/e response times for instantaneous changes of wall temperature and saturation ratio were found to be in the order of few minutes. Particle nucleation and growth was also simulated and particle number size distribution properties of the freshly nucleated particles (particle number, mean size, standard deviation of the assumed log-normal distribution) were found to be distributed over the tank's volume similar to the gas species.
Numerical simulation of flow, H₂SO₄ cycle and new particle formation in the CERN CLOUD chamber
(2011)
To study the effect of galactic cosmic rays on aerosols and clouds, the Cosmic Leaving OUtdoor Droplets (CLOUD) project was established. Experiments are carried out at a 26 m3 tank at CERN (Switzerland). In the experiments, the effect of ionising particle radiation on H2SO4 particle formation and growth is investigated. To evaluate the experimental configuration, the experiment was simulated using a coupled multidimensional CFD – particle model (CLOUD-FPM). In the model the coupled fields of gas/vapour species, temperature, flow velocity and particle properties were computed to investigate the tank's mixing state and mixing times. Simulation results show that the mixing state of the tank's contents largely depends on the characteristics of the mixing fans and a 1-fan configuration, as used in first experiments, may not be sufficient to ensure a homogeneously mixed chamber. To mix the tank properly, 2 fans are necessary. The 1/e response times for instantaneous changes of wall temperature and saturation ratio inside the chamber were found to be in the order of few minutes. Particle nucleation and growth was also simulated and particle number size distribution properties of the freshly nucleated particles (particle number, mean size, standard deviation of the assumed log-normal distribution) were found to be mixed over the tank's volume similar to the gas species.
This thesis deals with the analysis of “presolar” silicates and oxides by high resolution mass spectrometry and electron microscopy techniques. This “stardust” was identified by its extreme oxygen isotopic anomalies, which point to nucleosynthetic reactions in stellar interiors, in the carbonaceous chondrite Acfer 094. Isotopic, chemical and mineralogical studies on these stardust grains therefore allow the testing of astrophysical questions on Earth, which are otherwise only accessible by spectroscopy and theoretical models. The class of presolar silicates has been identified only six years ago in 2002, although it was known already from spectroscopic observations that silicates represent the most abundant type of dust in the galaxy. The development of the “NanoSIMS” was a crucial step in this respect, because this ion probe with its superior spatial resolution of only 50 nm allowed the detection of the typically 300 nm sized presolar silicates. A total of 142 presolar silicates and 20 presolar oxides were identified within Acfer 094, whose matrix therefore contains 163 ± 14 ppm presolar silicates and 26 ± 6 ppm presolar oxides. This is among the highest amounts reported so far for any primitive solar system material. The majority of detected stardust grains derive from asymptotic giant branch stars of 1 – 2.5 Msun and close-to-solar or slightly lower-than-solar metallicity. However, by measuring the Si isotopic compositions of some enigmatic grains, it could be shown that there is a sub-class of presolar silicates characterized by an extreme enrichment of 17O and a moderate enhancement of 30Si relative to solar, whose origins might be explained by formation in binary stellar systems. About 10% of all grains exhibit an enrichment in 18O and some of them also of 28Si relative to solar, which most likely point to an origin in type II supernova explosions. The Si isotopic measurements also allowed to quantify the effect of the s-process on the Si isotopes in low-mass asymptotic giant branch stars. The results agree well with theoretical predictions. The grains were furthermore characterized by SEM and the chemistries of about half of the grains were determined by Auger electron spectroscopy. The majority of grain morphologies are consistent with what is expected from condensation experiments. However, a lot of grains are altered by Fe-rich minerals, which are either of primary condensation or of secondary ISM or solar nebula origin. Furthermore, complex presolar grains consisting of refractory Al-rich grains attached to silicate material could be identified, which have been predicted by condensation theory and observational evidence. Nine presolar silicates were analyzed by combined NanoSIMS/TEM studies. The majority of grains are Mg-rich and amorphous, which is in contrast to astrophysical evidence, which mainly postulate crystalline Mg-rich and amorphous Fe-rich circumstellar condensates. However, the grains might have been rendered amorphous by secondary processes in the ISM or could have condensed under non-equilibrium, low-temperature conditions in the circumstellar outflow. The grains are more likely characterized by a variable, pyroxene-like chemistry, which could be a result of sputtering in the ISM, which preferentially removes Mg. The detected crystalline presolar silicates in this study and in other work are all olivines, whereas grains with a pyroxene stoichiometry are all amorphous except one. This supports astrophysical models which point to different formation pathways for these two types of grains and therefore different crystallinity. However, the relatively high Fe content of three detected presolar olivines in this study and in other work is in contrast to astrophysical evidence and theoretical considerations, which predict essentially Fe-free crystalline grains. It is therefore possible that the infrared spectra might also be compatible with less Mg-rich olivines. The only crystalline presolar silicate with a pyroxene-like stoichiometry is the unusual grain 1_07: although it is chemically enstatite, the electron diffraction pattern could only be indexed to silicate perovskite, which is stable above ~23 GPa. The discovery of a high-pressure phase of presolar origin shows that dust grains encountering interstellar shocks might not necessarily be completely destroyed. In astrophysical models it is in principle also possible that a fraction of larger grains might survive such a shock wave encounter as a high-pressure modification, which is supported by this discovery.
Motivated by the question of whether and how wave–wave interactions should be implemented into atmospheric gravity-wave parametrizations, the modulation of triadic gravity-wave interactions by a slowly varying and vertically sheared mean flow is considered for a non-rotating Boussinesq fluid with constant stratification. An analysis using a multiple-scale WKBJ (Wentzel–Kramers–Brillouin–Jeffreys) expansion identifies two distinct scaling regimes, a linear off-resonance regime, and a nonlinear near-resonance regime. Simplifying the near-resonance interaction equations allows for the construction of a parametrization for the triadic energy exchange which has been implemented into a one-dimensional WKBJ ray-tracing code. Theory and numerical implementation are validated for test cases where two wave trains generate a third wave train while spectrally passing through resonance. In various settings, of interacting vertical wavenumbers, mean-flow shear, and initial wave amplitudes, the WKBJ simulations are generally in good agreement with wave-resolving simulations. Both stronger mean-flow shear and smaller wave amplitudes suppress the energy exchange among a resonantly interacting triad. Experiments with mean-flow shear as strong as in the vicinity of atmospheric jets suggest that internal gravity-wave dynamics are dominated in such regions by wave modulation. However, triadic gravity-wave interactions are likely to be relevant in weakly sheared regions of the atmosphere.
During this study clumped isotope analysis of carbonates was established at the Goethe University of Frankfurt, Germany. Therefore, preparation protocols and analytical parameters were elaborated to obtain precise and accurate Δ47 data. Briefly, analyte CO2 was cleaned cryogenically using glass extraction lines to remove traces of water that enable re-equilibration of C–O bonds in the gases. Furthermore, analyte CO2 was passed through a gas chromatograph (GC) to clean it from contaminants that produce isobaric interferences with m/z 47. Initially, phosphoric acid digestions of carbonates was conducted at 25 °C in McCrea-type reaction vessels. Afterwards samples were reacted at 90 °C using a common acid bath. Mass spectrometric analyses were performed using a MAT 253 equipped with a dual inlet system. Δ47 values were directly projected to the absolute scale using CO2 gases equilibrated at distinct temperatures.
In cooperation with Stefano Bernasconi and his research group at ETH Zurich we studied the non-linearity that occurs for the measurement of m/z 47. This effect results from secondary electrons created by the m/z 44 beam. These electrons cause a negative background on the m/z 47 collector. A correction procedure was proposed that relies on the determination of the negative background on the m/z 47 Faraday cup. This approach might reduce time-consuming analyses of heated gases which were used so far to account for the observed non-linearity. However, the suggested correction of the negative background on the m/z 47 cup is only applicable if the slit width of the m/z 44 beam is significantly wider than that of the m/z 47 beam.
This thesis, furthermore, presents a comparison of the different phosphoric acid digestion techniques which are commonly used for carbonate clumped isotope analysis. For calcitic and aragonitic material digested at 25 °C in McCrea-type vessels we observed that the sample size has an effect on Δ47 data: higher mean Δ47 values and a larger scatter of data were received for samples <7 mg than for larger aliquots. For carbonate samples digested at 90 °C in a common acid bath no sample size effect was determined. We assume that secondary re-equilibration of CO2 with water preferentially occurs at 25 °C producing the observed differences. However, a sample size effect can be avoided if reaction temperature is increased to 90 °C.
In order to make carbonate Δ47 data obtained from acid digestions at 90 °C comparable to Δ47 data received from reactions at 25 °C the difference of the acid fractionation factores (Δ47*25-90) between both temperatures has to be known. For the determination of the Δ47*25-90 value we have considered Δ47 data made at 25 °C from samples >7 mg only. For calicte and aragonite we obtained differences in fractionation factores of 0.075‰ and 0.066‰, respectively. These Δ47*25-90 values are coincident with the theoretical prediction of 0.069‰ proposed for calcite (Guo et al., 2009).
Moreover, this dissertation comprises a calibration study of the clumped isotope thermometer based on various natural calcites that grew between 9 and 38 °C. The samples include a brachiopod shell, a bivalve shell, an eggshell of an ostrich and foraminifera tests which formed from distinct biomineralizing processes. Furthermore we included an authigenic carbonate crystallized from biological-induced precipitation. The following linear relationship between 1/T2 and Δ47 was determined (with Δ47 in ‰ and T in K):
Δ47 = 0.0327 (± 0.0026) x 106 / T2 + 0.3030 (± 0.0308) (R2 = 0.9915)
This equation differs from the pioneering Ghosh et al. (2006a) calibration. However, our regression line is statistically indistinguishable from that of Henkes et al. (2013) which is based on aragonitic mollusks and calcitic brachiopod shells. Both studies have in common that calibration data were, at first, directly referenced to the absolute scale. In addition, both datasets rely on similar digestion techniques. Furthermore, the two calibrations are conform with the theoretical prediction of Guo et al. (2009).
The calcite calibration of the clumped isotope paleothermometer received in this study was applied to Δ47 data measured for Silurian brachiopods shells from Gotland/Sweden. Prior to isotopic analysis the fossils were intensively investigated for their preservation state (CL, SEM, trace elements). The lowest T(Δ47) values of ca. 28 to 33 °C were estimated from ultrastructurally well-preserved regions of some shells. For these samples also the lowest δ18Ow values of Silurian seawater were determined. These estimates of ca. −1‰ confirm the assumption that the δ18O value of the Silurian ocean was buffered to (0 ± 1)‰.
Nevertheless, most studied shells were characterized by a patchwork of pristine and altered shell portions resulting in elevated T(Δ47) values which plot mostly between 40 and 60 °C. Our results indicate that the clumped isotopic composition of the shells were altered at low water-rock ratios, not affecting the δ18O values. Δ47 and δ18O data of associated diagenetic phases (sparitic and micritic phases of the inner fillings of the fossils) provide evidence that the sparitic cements grew during several diagenetic events which occurred at different temperatures in fluid-buffered systems. We, furthermore, conclude that the micritic phases lithified at a very early diagenetic stage with the δ18O values being most probably close to a Silurian seawater composition
Over recent decades, the global population has been rapidly increasing and human activities have altered terrestrial water fluxes to an unprecedented extent. The phenomenal growth of the human footprint has significantly modified hydrological processes in various ways (e.g. irrigation, artificial dams, and water diversion) and at various scales (from a watershed to the globe). During the early 1990s, awareness of the potential for increased water scarcity led to the first detailed global water resource assessments. Shortly thereafter, in order to analyse the human perturbation on terrestrial water resources, the first generation of large-scale hydrological models (LHMs) was produced. However, at this early stage few models considered the interaction between terrestrial water fluxes and human activities, including water use and reservoir regulation, and even fewer models distinguished water use from surface water and groundwater resources. Since the early 2000s, a growing number of LHMs have incorporated human impacts on the hydrological cycle, yet the representation of human activities in hydrological models remains challenging. In this paper we provide a synthesis of progress in the development and application of human impact modelling in LHMs. We highlight a number of key challenges and discuss possible improvements in order to better represent the human–water interface in hydrological models.
AirCore samplers have been increasingly used to capture vertical profiles of trace gases reaching from the ground up to about 30 km, in order to validate remote sens- ing instruments and to investigate transport processes in the stratosphere. When deployed to a weather balloon, accu- rately attributing the trace gas measurements to the sampling altitudes is nontrivial, especially in the stratosphere. In this paper we present the CO-spiking experiment, which can be deployed to any AirCore on any platform in order to evalu- ate different computational altitude attribution processes and to experimentally derive the vertical resolution of the profile by injecting small volumes of signal gas at predefined GPS altitudes during sampling. We performed two CO-spiking flights with an AirCore from the Goethe University Frankfurt (GUF) deployed to a weather balloon in Traînou, France, in June 2019. The altitude retrieval based on an instantaneous pressure equilibrium assumption slightly overestimates the sampling altitudes, especially at the top of the profiles. For these two flights our altitude attribution is accurate within 250 m below 20 km. Above 20 km the positive bias becomes larger and reaches up to 1.2 km at 27 km altitude. Differences in descent velocities are shown to have a major impact on the altitude attribution bias. We parameterize the time lag between the theoretically attributed altitude and the actual CO-spike release altitude for both flights together and use it to empirically correct our AirCore altitude retrieval. Regard- ing the corrected profiles, the altitude attribution is accurate within ±120 m throughout the profile. Further investigations are needed in order to test for the scope of validity of this correction parameter regarding different ambient conditions and maximum flight altitudes. We derive the vertical resolu- tion from the CO spikes of both flights and compare it to the modeled vertical resolution. The modeled vertical resolution is too optimistic compared to the experimentally derived res- olution throughout the profile, albeit agreeing within 220 m. All our findings derived from the two CO-spiking flights are strictly bound to the GUF AirCore dimensions. The newly introduced CO-spiking experiment can be used to test differ- ent combinations of AirCore configurations and platforms in future studies.
Atmospheric nanoaerosols have extensive effects on the Earth’s climate and human health. This cumulative work focuses on the development and characterization of instrumentation for measuring various parameters of atmospheric nanoaerosols, and its use to understand new particle formation from organic precursors. The principal research question is, how the chemical composition of nanoaerosol particles can be measured and how atmospheric chemistry influences aerosol processes, especially new particle formation and growth. Therefore, nanoaerosols are investigated under various aspects. More specifically, an instrument is developed to analyze nanoparticles, and field as well as chamber studies are conducted.
The main project is the instrument development of the Thermal Desorption Differential Mobility Analyzer (TD-DMA, project 1, Wagner et al. (2018)). This instrument analyzes the chemical composition of small aerosol particles. By characterization and testing in chamber experiments, it is proven to be suitable for the analysis of freshly nucleated particles.
The second project (Wagner et al. (2017)) applies a broad spectrum of aerosol measurement instruments for the characterization of aerosol particles produced by a skyscraper blasting. A comprehensive picture of the particle population emitted by the demolition is obtained.
Project 3 (K¨urten et al. (2016)) is also an ambient aerosol measurement, focusing of new particle formation in a rural area in central Germany, and the ability of a negative nitrate CI-APi-TOF to detect various substances in atmosphere. Project 4 (Heinritzi et al. (2016)) is a characterization of the negative nitrate CI-APi-TOF used in projects 1, 3, 5, 6, 7 and 8. The following projects focus on understanding new particle formation from atmospherically abundant organic precursors. Key instruments comprise the negative nitrate CI-APiTOF for gas-phase measurements of the nucleating species, and various sizing and counting instruments for quantifying the particle formation and growth. Project 5 (Kirkby et al. (2016)) shows that biogenic organic compounds formed from alpha-pinene can nucleate on their own without the influence of e.g. sulfuric acid. Project 6 (Tr¨ostl et al. (2016)) describe the subsequent growth of these particles. Project 7 (Stolzenburg et al. (2018)) covers the temperature dependence of this growth and in project 8 (Heinritzi et al. (2018)), the suppressing influence of isoprene on the new particle formation is assessed.
We present a study characterizing aerosol particles resulting from a skyscraper blasting. High mass concentrations with a maximum of 844.9 μg m-3 were present for a short time period of approximately 15 minutes. They result in a day mean of 32.6 μg m-3 compared to a 27.6 μg m-3 background not exceeding the 50 μg m-3 EU maximum permissive value. The increase in particle number concentration was less pronounced with a maximum concentration of 6.9 ⋅ 104 cm-3 compared to the local background value of 1.8 ⋅ 104 cm-3. The size-resolved number concentration shows a single mode of ultrafine particles at approximately 93 nm. The spatial distribution of deposited dust was investigated with Bergerhoff glass collection vessels, showing a decrease with distance. In the deposited dust samples the concentrations of twelve metals was determined, non of them exceeded the regional background concentrations significantly. The chemical composition of individual particles emitted by the demolition was studied by Scanning Electron Microscopy. They were mainly concrete and steel particles, with 60% calcium carbonates, 19% calcium sulfates, 19% silicates and 2% steel. In energy-dispersive X-Ray Spectroscopy, no fibers like asbestos were observed. Using a broad spectrum of instruments and methods, we obtain comprehensive characterization of the particles emitted by the demolition.
A new method for size-resolved chemical analysis of nucleation mode aerosol particles (size range from ∼10 to ∼30 nm) is presented. The Thermal Desorption Differential Mobility Analyzer (TD-DMA) uses an online, discontinuous principle. The particles are charged, a specific size is selected by differential mobility analysis and they are collected on a filament by electrostatic precipitation. Subsequently, the sampled mass is evaporated in a clean carrier gas and analyzed by a chemical ionization mass spectrometer. Gas-phase measurements are performed with the same mass spectrometer during the sampling of particles. The characterization shows reproducible results, with a particle size resolution of 1.19 and the transmission efficiency for 15 nm particles being slightly above 50 %. The signal from the evaporation of a test substance can be detected starting from 0.01 ng and shows a linear response in the mass spectrometer. Instrument operation in the range of pg m−3 is demonstrated by an example measurement of 15 nm particles produced by nucleation from dimethylamine, sulfuric acid and water.
A new method for size resolved chemical analysis of nucleation mode aerosol particles (size range from ~10 to ~30 nm) is presented. The Thermal Desorption Differential Mobility Analyzer (TD-DMA) uses an online, discontinuous principle. The particles are charged, a specific size is selected by differential mobility analysis and they are collected on a filament by electrostatic precipitation. Subsequently, the sampled mass is evaporated in a clean carrier gas and analyzed by a chemical ionization mass spectrometer. Gas phase measurements are performed with the same mass spectrometer during the sampling of particles. The characterization shows reproducible results, with a particle size resolution of 1.19 and the transmission efficiency for 15 nm particles being slightly above 50 %. The signal from the evaporation of a test substance can be detected starting from 0.01 ng and shows a linear response in the mass spectrometer. Instrument operation in the range of pg/m3 is demonstrated by an example measurement of 15 nm particles produced by nucleation from dimethylamine, sulfuric acid and water.
The formation of secondary particles in the atmosphere accounts for more than half of global cloud condensation nuclei. Experiments at the CERN CLOUD (Cosmics Leaving OUtdoor Droplets) chamber have underlined the importance of ions for new particle formation, but quantifying their effect in the atmosphere remains challenging. By using a novel instrument setup consisting of two nano-particle counters, one of them equipped with an ion filter, we were able to further investigate the ion-related mechanisms of new particle formation. In autumn 2015, we carried out experiments at CLOUD on four systems of different chemical compositions involving monoterpenes, sulfuric acid, nitrogen oxides, and ammonia. We measured the influence of ions on the nucleation rates under precisely controlled and atmospherically relevant conditions. Our results indicate that ions enhance the nucleation process when the charge is necessary to stabilize newly formed clusters, i.e. in conditions where neutral clusters are unstable. For charged clusters that were formed by ion-induced nucleation, we were able to measure, for the first time, their progressive neutralization due to recombination with oppositely charged ions. A large fraction of the clusters carried a charge at 1.2 nm diameter. However, depending on particle growth rates and ion concentrations, charged clusters were largely neutralized by ion–ion recombination before they grew to 2.2 nm. At this size, more than 90 % of particles were neutral. In other words, particles may originate from ion-induced nucleation, although they are neutral upon detection at diameters larger than 2.2 nm. Observations at Hyytiälä, Finland, showed lower ion concentrations and a lower contribution of ion-induced nucleation than measured at CLOUD under similar conditions. Although this can be partly explained by the observation that ion-induced fractions decrease towards lower ion concentrations, further investigations are needed to resolve the origin of the discrepancy.
Following votes in the Coniacian Working Group, the Cretaceous Subcommission and the International Commission on Stratigraphy, on May 1st, 2021, the International Union of Geological Sciences voted unanimously to ratify the Global Stratotype Section and Point (GSSP) proposal for the base of the Coniacian Stage of the Upper Cretaceous Series and Cretaceous System. The lower boundary of the Coniacian Stage is placed at the base of Bed 46 of the Salzgitter-Salder section in northern Germany. The boundary is defined by the first appearance of the inoceramid bivalve species Cremnoceramus deformis erectus (Meek) and complemented by the Navigation carbon isotope event. Additional data include the bivalve genus Didymotis, foraminifera, ammonite, nannofossil and organic-walled dinoflagellate cyst events. Three auxiliary sections (Słupia Nadbrzeżna, central Poland; Střeleč, Czech Republic; El Rosario, NE Mexico) supplement the details of the boundary record in various facies, and in differing geographic and biogeographic contexts.
Invasive plant species are increasingly altering species composition and the functioning of ecosystems from a local to a global scale. The grass species Pennisetum setaceum has recently raised concerns as an invader on different archipelagos worldwide. Among these affected archipelagos are the Canary Islands, which are a hotspot of endemism. Consequently, conservation managers and stakeholders are interested in the potential spreading of this species in the archipelago. We identify the current extent of the suitable habitat for P. setaceum on the island of La Palma to assess how it affects island ecosystems, protected areas (PAs), and endemic plant species richness. We recorded in situ occurrences of P. setaceum from 2010 to 2018 and compiled additional ones from databases at a 500 m × 500 m resolution. To assess the current suitable habitat and possible distribution patterns of P. setaceum on the island, we built an ensemble model. We projected habitat suitability for island ecosystems and PAs and identified risks for total as well as endemic plant species richness. The suitable habitat for P. setaceum is calculated to cover 34.7% of the surface of La Palma. In open ecosystems at low to mid elevations, where native ecosystems are already under pressure by land use and human activities, the spread of the invader will likely lead to additional threats to endemic plant species. Forest ecosystems (e.g., broadleaved evergreen and coniferous forests) are not likely to be affected by the spread of P. setaceum because of its heliophilous nature. Our projection of suitable habitat of P. setaceum within ecosystems and PAs on La Palma supports conservationists and policymakers in prioritizing management and control measures and acts as an example for the potential threat of this graminoid invader on other islands.
Increasing atmospheric CO2 stimulates photosynthesis which can increase net primary production (NPP), but at longer timescales may not necessarily increase plant biomass. Here we analyse the four decade-long CO2-enrichment experiments in woody ecosystems that measured total NPP and biomass. CO2 enrichment increased biomass increment by 1.05 ± 0.26 kg C m−2 over a full decade, a 29.1 ± 11.7% stimulation of biomass gain in these early-secondary-succession temperate ecosystems. This response is predictable by combining the CO2 response of NPP (0.16 ± 0.03 kg C m−2 y−1) and the CO2-independent, linear slope between biomass increment and cumulative NPP (0.55 ± 0.17). An ensemble of terrestrial ecosystem models fail to predict both terms correctly. Allocation to wood was a driver of across-site, and across-model, response variability and together with CO2-independence of biomass retention highlights the value of understanding drivers of wood allocation under ambient conditions to correctly interpret and predict CO2 responses.
Simulation of global temperature variations and signal detection studies using neural networks
(1998)
The concept of neural network models (NNM) is a statistical strategy which can be used if a superposition of any forcing mechanisms leads to any effects and if a sufficient related observational data base is available. In comparison to multiple regression analysis (MRA), the main advantages are that NNM is an appropriate tool also in the case of non-linear cause-effect relations and that interactions of the forcing mechanisms are allowed. In comparison to more sophisticated methods like general circulation models (GCM), the main advantage is that details of the physical background like feedbacks can be unknown. Neural networks learn from observations which reflect feedbacks implicitly. The disadvantage, of course, is that the physical background is neglected. In addition, the results prove to be sensitively dependent from the network architecture like the number of hidden neurons or the initialisation of learning parameters. We used a supervised backpropagation network (BPN) with three neuron layers, an unsupervised Kohonen network (KHN) and a combination of both called counterpropagation network (CPN). These concepts are tested in respect to their ability to simulate the observed global as well as hemispheric mean surface air temperature annual variations 1874 - 1993 if parameter time series of the following forcing mechanisms are incorporated : equivalent CO2 concentrations, tropospheric sulfate aerosol concentrations (both anthropogenic), volcanism, solar activity, and ENSO (all natural). It arises that in this way up to 83% of the observed temperature variance can be explained, significantly more than by MRA. The implication of the North Atlantic Oscillation does not improve these results. On a global average, the greenhouse gas (GHG) signal so far is assessed to be 0.9 - 1.3 K (warming), the sulfate signal 0.2 - 0.4 K (cooling), results which are in close similarity to the GCM findings published in the recent IPCC Report. The related signals of the natural forcing mechanisms considered cover amplitudes of 0.1 - 0.3 K. Our best NNM estimate of the GHG doubling signal amounts to 2.1K, equilibrium, or 1.7 K, transient, respectively.
The climate system can be regarded as a dynamic nonlinear system. Thus, traditional linear statistical methods fail to model the nonlinearities of such a system. These nonlinearities render it necessary to find alternative statistical techniques. Since artificial neural network models (NNM) represent such a nonlinear statistical method their use in analyzing the climate system has been studied for a couple of years now. Most authors use the standard Backpropagation Network (BPN) for their investigations, although this specific model architecture carries a certain risk of over-/underfitting. Here we use the so called Cauchy Machine (CM) with an implemented Fast Simulated Annealing schedule (FSA) (Szu, 1986) for the purpose of attributing and detecting anthropogenic climate change instead. Under certain conditions the CM-FSA guarantees to find the global minimum of a yet undefined cost function (Geman and Geman, 1986). In addition to potential anthropogenic influences on climate (greenhouse gases (GHG), sulphur dioxide (SO2)) natural influences on near surface air temperature (variations of solar activity, explosive volcanism and the El Nino = Southern Oscillation phenomenon) serve as model inputs. The simulations are carried out on different spatial scales: global and area weighted averages. In addition, a multiple linear regression analysis serves as a linear reference. It is shown that the adaptive nonlinear CM-FSA algorithm captures the dynamics of the climate system to a great extent. However, free parameters of this specific network architecture have to be optimized subjectively. The quality of the simulations obtained by the CM-FSA algorithm exceeds the results of a multiple linear regression model; the simulation quality on the global scale amounts up to 81% explained variance. Furthermore the combined anthropogenic effect corresponds to the observed increase in temperature Jones et al. (1994), updated by Jones (1999a), for the examined period 1856–1998 on all investigated scales. In accordance to recent findings of physical climate models, the CM-FSA succeeds with the detection of anthropogenic induced climate change on a high significance level. Thus, the CMFSA algorithm can be regarded as a suitable nonlinear statistical tool for modeling and diagnosing the climate system.
Attribution and detection of anthropogenic climate change using a backpropagation neural network
(2002)
The climate system can be regarded as a dynamic nonlinear system. Thus traditional linear statistical methods are not suited to describe the nonlinearities of this system which renders it necessary to find alternative statistical techniques to model those nonlinear properties. In addition to an earlier paper on this subject (WALTER et al., 1998), the problem of attribution and detection of the observed climate change is addressed here using a nonlinear Backpropagation Neural Network (BPN). In addition to potential anthropogenic influences on climate (CO2-equivalent concentrations, called greenhouse gases, GHG and SO2 emissions) natural influences on surface air temperature (variations of solar activity, volcanism and the El Niño/Southern Oscillation phenomenon) are integrated into the simulations as well. It is shown that the adaptive BPN algorithm captures the dynamics of the climate system, i.e. global and area weighted mean temperature anomalies, to a great extent. However, free parameters of this network architecture have to be optimized in a time consuming trial-and-error process. The simulation quality obtained by the BPN exceeds the results of those from a linear model by far; the simulation quality on the global scale amounts to 84% explained variance. Additionally the results of the nonlinear algorithm are plausible in a physical sense, i.e. amplitude and time structure. Nevertheless they cover a broad range, e.g. the GHG-signal on the global scale ranges from 0.37 K to 1.65 K warming for the time period 1856-1998. However the simulated amplitudes are situated within the discussed range (HOUGHTON et al., 2001). Additionally the combined anthropogenic effect corresponds to the observed increase in temperature for the examined time period. In addition to that, the BPN succeeds with the detection of anthropogenic induced climate change on a high significance level. Therefore the concept of neural networks can be regarded as a suitable nonlinear statistical tool for modeling and diagnosing the climate system.
Titanite is a potentially powerful U–Pb petrochronometer that may record metamorphism, metasomatism, and deformation. Titanite may also incorporate significant inherited Pb, which may lead to inaccurate and geologically ambiguous U–Pb dates if a proper correction is not or cannot be applied. Here, we present laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS)-derived titanite U–Pb dates and trace element concentrations for two banded calcsilicate gneisses from south-central Maine, USA (SSP18-1A and SSP18-1B). Single spot common Pb-corrected dates range from 400 to 280 Ma with ±12–20 Ma propagated 2SE. Titanite grains in sample SSP18-1B exhibit regular core-to-rim variations in texture, composition, and date. We identify four titanite populations: (1) 397 ± 5 Ma (95% CL) low Y + HREE cores and mottled grains, (2) 370 ± 7 Ma high Y + REE mantles and cores, (3) 342 ± 6 Ma cores with high Y + REE and no Eu anomaly, and (4) 295 ± 6 Ma LREE-depleted rims. We interpret the increase in titanite Y + HREE between ca. 397 and ca. 370 Ma to constrain the timing of diopside fracturing and recrystallization and amphibole breakdown. Apparent Zr-in-titanite temperatures (803 ± 36°C at 0.5 ± 0.2 GPa) and increased XDi suggest a thermal maximum at ca. 370 Ma. Population 3 domains dated to ca. 342 Ma exhibit no Eu anomaly and are observed only in compositional bands dominated by diopside (>80 vol%), suggesting limited equilibrium between titanite and plagioclase. Finally, low LREE and high U/Th in Population 4 titanite dates the formation of hydrous phases, such as allanite, during high XH2O fluid infiltration at ca. 295 Ma. In contrast to the well-defined date–composition–texture relationships observed for titanite from SSP18-1B, titanite grains from sample SSP18-1A exhibit complex zoning patterns and little correlation between texture, composition, and date. We hypothesize that the incorporation of variable amounts of radiogenic Pb from dissolved titanite into recrystallized domains resulted in mixed dates spanning 380–330 Ma. Although titanite may reliably record multiple phases of metamorphism, these data highlight the importance of considering U–Pb data along with chemical and textural data to screen for inherited radiogenic Pb.
Sulfur in the slab: a sulfur-isotopes and thermodynamic-modeling perspective from exhumed terranes
(2022)
Sulfur is a key element in the subduction zone-volcanic arc system; however, the mechanism(s) that recycle sulfur from the slab into the overlying volcanic arc are debated. Here we summarize recent advances in quantifying this component of the deep sulfur cycle. First, primary metamorphic or inherited sulfides in oceanic-type eclogites are only rarely observed as inclusions and are typically absent from the rock matrix. Additionally, sulfides are relatively common in rocks metasomatized at the slab-mantle interface by slab-derived fluids during exhumation. Combined, these two observations suggest that sulfur loss from subducted mafic crust is relatively efficient. Thermodynamic modeling in Perple_X using the Holland and Powell (2011) database combined with the Deep Earth Water model suggests that the efficiency and speciation of sulfur loss varies depending on the degree of seafloor alteration prior to subduction and the geothermal gradient of the slab. In relatively cold subduction zones, such as Honshu, slab-fluids derived from subducted mafic crust are predicted to exhibit elevated concentrations of HSO4-, SO42-, HSO3-, and CaSO4(aq), whereas hot subduction zones, such as Cascadia, are predicted to produce slab fluids enriched in HS- and H2S at lower pressures. The oxidation of sulfur expelled from subducted pyrite is balanced by the reduction of Fe3+ to Fe2+, consistent with the low Fe3+/SFe of exhumed eclogites relative to blueschists and altered oceanic crust. Where oxidized S-bearing fluids are produced, they are anticipated to interact with more reduced rocks at the slab-mantle interface and within the mantle wedge, resulting in sulfide precipitation and significant isotopic fractionation. The δ34S values of slab fluids are estimated to fall between -11 and +8 ‰. Rayleigh fractionation during progressive fluid-rock interaction results in fractionations of tens of per mil as oxidized species are depleted and sulfides are precipitated, resulting in δ34S values of sulfides that easily span the -21.7 to +13.9 ‰ range observed in metasomatic sulfides in exhumed high-pressure rocks. However, in subduction zones where reduced species prevail, the S isotopic signature of slab fluids is expected to reflect their source and will exhibit a narrower range in δ34S values. As a result, the δ34S values measured in arc magmas may not always be a reliable indicator of the contribution of different components of the slab, such as sediments vs. AOC. Additionally, the impact of S recycling on the oxygen fugacity of arc magmas is expected to vary both spatially and temporally throughout Earth history.
The oxidation state of sulfur in slab fluids is controversial, with both dominantly oxidized and reduced species proposed. Here we use in situ X-ray absorption spectroscopy analysis of sulfur-in-apatite to monitor changes in the oxidation state of sulfur during high-P metasomatism by slab fluids in the subduction channel. Our samples include a 73 cm continuous transect of reaction zones between a metagabbroic eclogite block and serpentinite matrix from a mélange zone on the island of Syros, Greece. The block core consists of garnet, omphacite, phengite, paragonite, epidote-clinozoisite, and rutile. In this region, apatite is only observed as elongate inclusions in omphacite cores. From the core outwards micas are increasingly replaced by epidote-clinozoisite, garnets are smaller and more frequent, pyrite + bornite is observed as inclusions in recrystallized omphacite, and apatite is increasingly abundant in the matrix and inclusions in garnet. A major transition at 48 cm separates an assemblage of Ca-Na amphibole, omphacite, chlorite, pyrite, and apatite from the inner garnet-bearing eclogite assemblages. Omphacite disappears from the assemblage at ~56 cm and amphibole compositions sharply transition to tremolite at 59 cm. Finally, the assemblage tremolite + talc + pyrite is observed after ~70 cm.Apatites in the eclogite assemblages exclusively display S6+ peaks in their absorption spectra. This includes apatite inclusions in omphacite in the least altered lithology, as well as matrix apatite and isolated apatite inclusions in garnet in the outermost metasomatized eclogite zone. In the intermediate pyrite-rich (~1-5 vol %) amphibole + omphacite + chlorite zone, apatite displays a strong S1- absorption peak in most grains, with rare analyses showing mixed S1- and S6+. Finally, apatite in the outermost tremolite-bearing assemblages only displays a S6+ peak. The pyrite-rich zone at 48 cm occurs at the initial interface between the serpentinite matrix and eclogite block, characterized by a dramatic decrease in Na content and Mg#. Our data suggest that reduction of S6+ in infiltrating fluids to S1- in pyrite became focused as Fe diffused across the steep Mg# gradient, resulting in pyrite precipitation. In contrast, S reduction in the Mg-rich tremolite-dominant portions of the transect was limited by a lack of Fe, resulting in low modes of pyrite and fluid buffered S6+ in apatite. Finally, S6+-bearing apatite is also observed in reaction zone lithologies from elsewhere on Syros, suggesting our observations are not isolated.Two important conclusions are drawn from these data and observations: (1) In the case of Syros, slab fluids at eclogite-facies conditions carried oxidized S6+, and (2) The interaction of these fluids with eclogites composed of ferrous-Fe silicates resulted in extensive sulfide precipitation.
Analyzing the impact of streamflow drought on hydroelectricity production: a global-scale study
(2021)
Electricity production by hydropower is negatively affected by drought. To understand and quantify risks of less than normal streamflow for hydroelectricity production (HP) at the global scale, we developed an HP model that simulates time series of monthly HP worldwide and thus enables analyzing the impact of drought on HP. The HP model is based on a new global hydropower database (GHD), containing 8,716 geo-localized plant records, and on monthly streamflow values computed by the global hydrological model WaterGAP with a spatial resolution of 0.5°. The GHD includes 44 attributes and covers 91.8% of the globally installed capacity. The HP model can reproduce HP trends, seasonality, and interannual variability that was caused by both (de)commissioning of hydropower plants and hydrological variability. It can also simulate streamflow drought and its impact on HP reasonably well. Global risk maps of HP reduction were generated for both 0.5° grid cells and countries, revealing that 67 out of the 134 countries with hydropower suffer, in 1 out of 10 years, from a reduction of more than 20% of mean annual HP and 18 countries from a reduction of more than 40%. The developed HP model enables advanced assessments of drought impacts on hydroelectricity at national to international levels.
Measurement of iodine species and sulfuric acid using bromide chemical ionization mass spectrometers
(2021)
Iodine species are important in the marine atmosphere for oxidation and new-particle formation. Understanding iodine chemistry and iodine new-particle formation requires high time resolution, high sensitivity, and simultaneous measurements of many iodine species. Here, we describe the application of a bromide chemical ionization mass spectrometer (Br-CIMS) to this task. During the iodine oxidation experiments in the Cosmics Leaving OUtdoor Droplets (CLOUD) chamber, we have measured gas-phase iodine species and sulfuric acid using two Br-CIMS, one coupled to a Multi-scheme chemical IONization inlet (Br-MION-CIMS) and the other to a Filter Inlet for Gasses and AEROsols inlet (Br-FIGAERO-CIMS). From offline calibrations and intercomparisons with other instruments, we have quantified the sensitivities of the Br-MION-CIMS to HOI, I2, and H2SO4 and obtained detection limits of 5.8 × 106, 3.8 × 105, and 2.0 × 105 molec. cm−3, respectively, for a 2 min integration time. From binding energy calculations, we estimate the detection limit for HIO3 to be 1.2 × 105 molec. cm−3, based on an assumption of maximum sensitivity. Detection limits in the Br-FIGAERO-CIMS are around 1 order of magnitude higher than those in the Br-MION-CIMS; for example, the detection limits for HOI and HIO3 are 3.3 × 107 and 5.1 × 106 molec. cm−3, respectively. Our comparisons of the performance of the MION inlet and the FIGAERO inlet show that bromide chemical ionization mass spectrometers using either atmospheric pressure or reduced pressure interfaces are well-matched to measuring iodine species and sulfuric acid in marine environments.
A list of authors and their affiliations appears at the end of the paper New-particle formation is a major contributor to urban smog, but how it occurs in cities is often puzzling. If the growth rates of urban particles are similar to those found in cleaner environments (1–10 nanometres per hour), then existing understanding suggests that new urban particles should be rapidly scavenged by the high concentration of pre-existing particles. Here we show, through experiments performed under atmospheric conditions in the CLOUD chamber at CERN, that below about +5 degrees Celsius, nitric acid and ammonia vapours can condense onto freshly nucleated particles as small as a few nanometres in diameter. Moreover, when it is cold enough (below −15 degrees Celsius), nitric acid and ammonia can nucleate directly through an acid–base stabilization mechanism to form ammonium nitrate particles. Given that these vapours are often one thousand times more abundant than sulfuric acid, the resulting particle growth rates can be extremely high, reaching well above 100 nanometres per hour. However, these high growth rates require the gas-particle ammonium nitrate system to be out of equilibrium in order to sustain gas-phase supersaturations. In view of the strong temperature dependence that we measure for the gas-phase supersaturations, we expect such transient conditions to occur in inhomogeneous urban settings, especially in wintertime, driven by vertical mixing and by strong local sources such as traffic. Even though rapid growth from nitric acid and ammonia condensation may last for only a few minutes, it is nonetheless fast enough to shepherd freshly nucleated particles through the smallest size range where they are most vulnerable to scavenging loss, thus greatly increasing their survival probability. We also expect nitric acid and ammonia nucleation and rapid growth to be important in the relatively clean and cold upper free troposphere, where ammonia can be convected from the continental boundary layer and nitric acid is abundant from electrical storms.
New particle formation in the upper free troposphere is a major global source of cloud condensation nuclei (CCN)1,2,3,4. However, the precursor vapours that drive the process are not well understood. With experiments performed under upper tropospheric conditions in the CERN CLOUD chamber, we show that nitric acid, sulfuric acid and ammonia form particles synergistically, at rates that are orders of magnitude faster than those from any two of the three components. The importance of this mechanism depends on the availability of ammonia, which was previously thought to be efficiently scavenged by cloud droplets during convection. However, surprisingly high concentrations of ammonia and ammonium nitrate have recently been observed in the upper troposphere over the Asian monsoon region5,6. Once particles have formed, co-condensation of ammonia and abundant nitric acid alone is sufficient to drive rapid growth to CCN sizes with only trace sulfate. Moreover, our measurements show that these CCN are also highly efficient ice nucleating particles—comparable to desert dust. Our model simulations confirm that ammonia is efficiently convected aloft during the Asian monsoon, driving rapid, multi-acid HNO3–H2SO4–NH3 nucleation in the upper troposphere and producing ice nucleating particles that spread across the mid-latitude Northern Hemisphere.
his study aims at a detailed characterization of an ultra-fine aerosol particle counting system for operation on board the Russian high altitude research aircraft M-55 "Geophysica" (maximum ceiling of 21 km). The COndensation PArticle counting Systems (COPAS) consists of an aerosol inlet and two dual-channel continuous flow Condensation Particle Counters (CPCs).
The aerosol inlet, adapted for COPAS measurements on board the M-55 "Geophysica", is described concerning aspiration, transmission, and transport losses. The counting efficiencies of the CPCs using the chlorofluorocarbon FC-43 as the working fluid are studied experimentally at two pressure conditions, 300 hPa and 70 hPa. Three COPAS channels are operated with different temperature differences between the saturator and the condenser block yielding smallest detectable particle sizes (dp50 – as 50% detection "cut off" diameters) of 6 nm, 11 nm, and 15 nm, respectively, at ambient pressure of 70 hPa. The fourth COPAS channel is operated with an aerosol heating line (250°C) for a determination of the non-volatile number of particles. The heating line is experimentally proven to volatilize pure H2SO4-H2O particles for a particle diameter (dp) range of 11 nm<dp<200 nm.
Additionally this study includes investigation to exclude auto-nucleation of the working fluid inside the CPCs. An instrumental inter-comparison (cross-correlation) has been performed for several measurement flights and mission flights in the Arctic and the Tropics are discussed. Finally, COPAS measurements are used for an aircraft plume crossing analysis.
A characterization of the ultra-fine aerosol particle counter COPAS (COndensation PArticle counting System) for operation on board the Russian high altitude research aircraft M-55 Geophysika is presented. The COPAS instrument consists of an aerosol inlet and two dual-channel continuous flow Condensation Particle Counters (CPCs) operated with the chlorofluorocarbon FC-43. It operates at pressures between 400 and 50 hPa for aerosol detection in the particle diameter (dp) range from 6 nm up to 1 micro m. The aerosol inlet, designed for the M-55, is characterized with respect to aspiration, transmission, and transport losses. The experimental characterization of counting efficiencies of three CPCs yields dp50 (50% detection particle diameter) of 6 nm, 11 nm, and 15 nm at temperature differences (DeltaT) between saturator and condenser of 17°C, 30°C, and 33°C, respectively. Non-volatile particles are quantified with a fourth CPC, with dp50=11 nm. It includes an aerosol heating line (250°C) to evaporate H2SO4-H2O particles of 11 nm<dp<200 nm at pressures between 70 and 300 hPa. An instrumental in-flight inter-comparison of the different COPAS CPCs yields correlation coefficients of 0.996 and 0.985. The particle emission index for the M-55 in the range of 1.4–8.4×10 16 kg -1 fuel burned has been estimated based on measurements of the Geophysika's own exhaust.
We present the prototype of a regional climate system model based on the COSMO-CLM regional climate model coupled with several model components, analyze the performance of the couplings and present a strategy to find an optimum configuration with respect to computational costs and time to solution.
The OASIS3-MCT coupler is used to couple COSMO-CLM with two land surface models (CLM and VEG3D), a regional ocean model for the Mediterranean Sea (NEMO-MED12), two ocean models for the North and Baltic Sea (NEMO-NORDIC and TRIMNP+CICE) and the atmospheric component of an earth system model (MPI-ESM). We present a unified OASIS3-MCT interface which handles all couplings in a similar way, minimizes the model source code modifications and describes the physics and numerics of the couplings. Furthermore, we discuss solutions for specific regional coupling problems like handling of different domains, multiple usage of MCT interpolation library and efficient exchange of 3D fields.
A series of real-case simulations over Europe has been conducted and the computational performance of the couplings has been analyzed. The usage of the LUCIA tool of the OASIS3-MCT coupler enabled separation of the direct costs of: coupling, load imbalance and additional computations. The resulting limits for time to solution and costs are shown and the potential of further improvement of the computational efficiency is summarized for each coupling.
It was found that the OASIS3-MCT coupler keeps the direct coupling costs of communication and horizontal interpolation small in comparison with the costs of the additional computations and load imbalance for all investigated couplings. For the first time this could be demonstrated for an exchange of approximately 450 2D fields per time step necessary for the atmosphere-atmosphere coupling between COSMO-CLM and MPI-ESM.
A procedure for finding an optimum configuration for each of the couplings was developed considering the time to solution and costs of the simulations. The optimum configurations are presented for sequential and concurrent coupling layouts. The procedure applied can be regarded as independent on the specific coupling layout and coupling details.
In this dissertation, different aspects of turbulent transport and thermally driven flows over complex terrain are investigated. Two publications concentrate on the vertical heat and moisture exchange in the convective boundary layer over mountainous terrain. To study this, Large-Eddy Simulation (LES) is used. Both turbulent and advective transport mechanisms are evaluated over the simple orography of a quasi-two-dimensional, periodic valley with prescribed surface fluxes. Here, terrain elevation varies along only one of the horizontal coordinate axes. Even a relatively shallow orography, possibly unresolved in existing numerical weather prediction models, modifies the domain-averaged moisture and temperature profiles. For the analysis, the flow is decomposed into a local turbulent part, a local mean circulation, and a large-scale part. An analysis of the turbulent kinetic energy and turbulent heat and moisture flux budgets shows that the thermal circulation significantly contributes to the vertical transport. It is found that thermal upslope winds are important for the moisture transport from the valley to the mountain tops. In total, moisture export out of the valley is mostly accomplished by the mean circulation. On the temperature distribution, which is horizontally relatively homogeneous, the thermal circulation has a weaker impact. If an upper-level wind is present, it interacts with the thermal circulation. This weakens the vertical transport of moisture and thus reduces its export out of the valley. The heat transport is less affected by the upper-level wind because of its weaker dependence on the thermal circulation. These findings were corroborated in a more realistic experiment simulating the full diurnal cycle using radiation forcing and an interactive land surface model.
Based on these results, coherent turbulent structures in the convective boundary layer over non-flat terrain are studied in further detail. A conditional sampling method based on the concentration of a decaying passive tracer is implemented in order to identify the boundary-layer plumes objectively. Conditional sampling allows to quantify the contribution of plume structures to the vertical transport of heat and moisture. In case of the idealized valley, vertical transport by coherent structures is the dominant contribution to the turbulent components of both heat and moisture flux. It is comparable in magnitude to the advective transport by the mean slope-wind circulation, although it is more important for heat than for moisture transport. A set of less idealized simulations considers the flow over three-dimensional terrain. In this case, conditional sampling is carried out by using a simple domain-decomposition approach. We demonstrate that thermal updrafts are generally more frequent on hill tops than over the surroundings, but they are less persistent on the windward sides when large-scale winds are present in the free atmosphere.
The tools for flow decomposition and budget analysis are also applied in another idealized case with a quasi-two-dimensional valley featuring the stable boundary layer. Here, the formation of a low stratus cloud is investigated. The main driver for the cloud formation is radiative cooling due to outgoing longwave radiation. Despite a purely horizontal flow, the advection terms in the prognostic equations for heat and moisture produce vertical mixing across the upper cloud edge leading to a loss of cloud water content. However, this behavior is not due to any kind of thermally-driven circulation. Instead, this spurious mixing is caused by the diffusive error of the advection scheme in regions where the sloping surfaces of the terrain-following vertical coordinate intersect the cloud top. It is shown that the intensity of the (spurious) numerical diffusion strongly depends on the horizontal resolution, the order of advection, and the choice of the scalar advection scheme. A LES with 4 m horizontal resolution serves as a reference. For horizontal resolutions of a few hundred meters, carried out with a model setup as it is used in Numerical Weather Prediction, a strong reduction of the simulated liquid-water path is observed. In order to keep the (spurious) numerical diffusion at coarser resolutions small, at least a fifth-order advection scheme should be used. In the present case, a WENO scalar advection scheme turns out to increase the numerical diffusion along a sharp cloud edge compared to an upwind scheme. Furthermore, the choice of the vertical coordinate has a strong impact on the simulated liquid-water path over orography. With a modified definition of the terrain-following sigma coordinate, it is possible to produce cloud water where the classical sigma coordinate does not allow any cloud formation.
Since the discovery of the ozone hole [Farman et al., 1985], the dynamics of the stratosphere and the transport of anthropogenic trace gases from the surface to the higher atmosphere has come into the focus of interest. In the tropics, air rises high into the stratosphere and is transported poleward. Trace gases like the CFCs1, photochemically stable in the troposphere, are thus transported into regions where they are photolyzed. The products of the photolysis reactions (eg. Cl, Br) interact in the catalytic ozone cycles and lead to enhanced ozone depletion. Regarding the transport of trace gases, the so-called lowermost stratosphere (LMS) is a very interesting region, where the troposphere and the stratosphere directly interact and air masses out of both regions are mixed. It is the lowest part of the stratosphere between the tropopause and the 380 K isentrope. Tropospheric air can enter this region directly via isentropic transport across the extra-tropical tropopause whereas stratospheric air descends across the 380 K isentrope via the mean meridional circulation. Stratosphere-troposphere exchange (STE) controls the chemical composition of the LMS as well as of the tropopause region and thus has an important effect on the radiative and chemical balance of these regions and the climate system. STE exhibits a strong seasonality [Holton et al., 1995]. While downwelling of stratospheric air across the 380 K isentrope is the dominant process in winter, troposphere-to-stratosphere transport (TST) gains importance in summer, when the downwelling from the stratosphere is only weak. Isentropic transport across the extra-tropical tropopause occurs in regions where the tropopause is strongly disturbed and is connected to tropopause folds, streamer events, frontal zones, polar and subtropical jets, warm conveyor belts and cut-off low systems. A short introduction into STE, the LMS region, and methods to study atmospheric transport is given in Chapter 1. One useful tool to analyse the motions of air and transport processes are longlived trace gases. Since the lifetimes of these tracers are longer than the time scale of the transport processes they are involved in, the distribution of tracers in the atmosphere is mostly determined by dynamics. In the context of this thesis, measurements of such long-lived tracers were performed and used to study transport into the LMS region in the northern hemisphere. During the Vintersol/EuPLEx and ENVISAT validation campaigns in winter 2003, long-lived tracers such as N2O, CH4, CFC-12, CFC-11, H-1211, H2, SF6 and CO2 were measured with the High Altitude Gas Analyser (HAGAR), a two channel in-situ gas chromatograph combined with a CO2 instrument, based on nondispersive infrared absorption. Combined with measurements taken during campaigns in Forli/Italy (ENVISAT validation) in July and October 2002, tracer data were gathered from the tropopause up to altitudes around 20 km during 25 flights on board the Russian high-altitude aircraft M55 Geophysica. Thus, a substantial set of high quality tracer data has been obtained covering the polar vortex region as well as the mid latitudes of the northern hemisphere. Chapter 2 gives an overview of the HAGAR instrument and necessary improvements of the instrumental set up (implementing a CH4 channel) that were performed in the context of this thesis, and review data processing, the measurement campaigns. In order to study transport into the LMS it is assumed that air basically enters the LMS via three different pathways: a) quasi-isentropic transport from the troposphere, b) downward advection from the middle stratosphere through the 380 K surface and c) in the polar vortex region subsidence of air from of the polar vortex. Fractions of air originating in each of these source regions are determined with a simple mass balance calculation by using observations of a subset of the above species with distinct lifetimes (N2O, CH4, CFC-11, H-1211, H2 and O3) yielding complementary constraints on transport from each region. Details of the mass balance calculation and the results are presented in Chapter 3. During the mid-latitude measurement campaigns in Forlí the passing of a cut-off low system associated with an elongated streamer over Europe was observed. The impacts of this event on the trace gas mixing ratios in the LMS are examined in Chapter 4. Finally, a summary is given in Chapter 5.