550 Geowissenschaften
Refine
Year of publication
- 2021 (91)
- 2016 (78)
- 2020 (56)
- 2019 (48)
- 2022 (46)
- 2011 (44)
- 2015 (44)
- 2010 (42)
- 2013 (42)
- 2014 (40)
- 2009 (39)
- 2017 (38)
- 2008 (35)
- 2012 (31)
- 2018 (31)
- 2005 (20)
- 2023 (18)
- 2007 (15)
- 2006 (11)
- 2003 (10)
- 2002 (9)
- 2024 (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 (672)
- Doctoral Thesis (116)
- Conference Proceeding (14)
- Working Paper (11)
- Book (6)
- Part of a Book (4)
- Master's Thesis (3)
- Preprint (2)
- Review (2)
- Bachelor Thesis (1)
Language
- English (832) (remove)
Has Fulltext
- yes (832)
Is part of the Bibliography
- no (832)
Keywords
- climate change (11)
- Climate change (7)
- COSMO-CLM (6)
- precipitation (6)
- Atmospheric chemistry (5)
- Palaeoclimate (5)
- Biogeochemistry (4)
- Geochemistry (4)
- Palaeoceanography (4)
- Bewässerung (3)
Institute
- Geowissenschaften (616)
- Geowissenschaften / Geographie (126)
- Biodiversität und Klima Forschungszentrum (BiK-F) (61)
- Senckenbergische Naturforschende Gesellschaft (45)
- Geographie (22)
- Biowissenschaften (21)
- Institut für Ökologie, Evolution und Diversität (6)
- Institut für sozial-ökologische Forschung (ISOE) (5)
- Frankfurt Institute for Advanced Studies (FIAS) (4)
- Extern (3)
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.
Strong seasonal variability of hygric and thermal soil conditions are a defining environmental feature in northern Australia. However, how such changes affect the soil–atmosphere exchange of nitrous oxide (N2O), nitric oxide (NO) and dinitrogen (N2) is still not well explored. By incubating intact soil cores from four sites (three savanna, one pasture) under controlled soil temperatures (ST) and soil moisture (SM) we investigated the release of the trace gas fluxes of N2O, NO and carbon dioxide (CO2). Furthermore, the release of N2 due to denitrification was measured using the helium gas flow soil core technique. Under dry pre-incubation conditions NO and N2O emissions were very low (<7.0 ± 5.0 μg NO-N m−2 h−1; <0.0 ± 1.4 μg N2O-N m−2 h−1) or in the case of N2O, even a net soil uptake was observed. Substantial NO (max: 306.5 μg N m−2 h−1) and relatively small N2O pulse emissions (max: 5.8 ± 5.0 μg N m−2 h−1) were recorded following soil wetting, but these pulses were short lived, lasting only up to 3 days. The total atmospheric loss of nitrogen was generally dominated by N2 emissions (82.4–99.3% of total N lost), although NO emissions contributed almost 43.2% to the total atmospheric nitrogen loss at 50% SM and 30 °C ST incubation settings (the contribution of N2 at these soil conditions was only 53.2%). N2O emissions were systematically higher for 3 of 12 sample locations, which indicates substantial spatial variability at site level, but on average soils acted as weak N2O sources or even sinks. By using a conservative upscale approach we estimate total annual emissions from savanna soils to average 0.12 kg N ha−1 yr−1 (N2O), 0.68 kg N ha−1 yr−1 (NO) and 6.65 kg N ha−1 yr−1 (N2). The analysis of long-term SM and ST records makes it clear that extreme soil saturation that can lead to high N2O and N2 emissions only occurs a few days per year and thus has little impact on the annual total. The potential contribution of nitrogen released due to pulse events compared to the total annual emissions was found to be of importance for NO emissions (contribution to total: 5–22%), but not for N2O emissions. Our results indicate that the total gaseous release of nitrogen from these soils is low and clearly dominated by loss in the form of inert nitrogen. Effects of seasonally varying soil temperature and moisture were detected, but were found to be low due to the small amounts of available nitrogen in the soils (total nitrogen <0.1%).
Strong seasonal variability of hygric and thermal soil conditions are a defining environmental feature in Northern Australia. However, how such changes affect the soil–atmosphere exchange of nitrous oxide (N2O), nitric oxide (NO) and dinitrogen (N2) is still 5 not well explored. By incubating intact soil cores from four sites (3 savanna, 1 pasture) under controlled soil temperatures (ST) and soil moisture (SM) we investigated the release of the trace gas fluxes of N2O, NO and carbon dioxide (CO2). Furthermore, the release of N2 due to denitrification was measured using the helium gas flow soil core technique. Under dry pre-incubation conditions NO and N2O emission were very low (< 7.0± 5.0 μgNO-Nm−2 h−1; < 0.0± 1.4 μgN2O-Nm−2 h−1) or in case of N2O, even a net soil uptake was observed. Substantial NO (max: 306.5 μgNm−2 h−1) and relatively small N2O pulse emissions (max: 5.8±5.0 μgNm−2 h−1) were recorded following soil wetting, but these pulses were short-lived, lasting only up to 3 days. The total atmospheric loss of nitrogen was dominated by N2 emissions (82.4–99.3% of total N lost), although NO emissions contributed almost 43.2% at 50% SM and 30 °C ST. N2O emissions were systematically higher for 3 of 12 sample locations, which indicates substantial spatial variability at site level, but on average soils acted as weak N2O sources or even sinks. Emissions were controlled by SM and ST for N2O and CO2, ST and pH for NO, and SM and pH for N2.