Geowissenschaften
Refine
Year of publication
- 2013 (46) (remove)
Document Type
- Article (36)
- Doctoral Thesis (7)
- Book (2)
- Part of Periodical (1)
Has Fulltext
- yes (46) (remove)
Is part of the Bibliography
- no (46)
Keywords
Barriereinseln und Außensände, geformt durch eine Kombination aus Wind, Wellen, Strömung und Küstenlängstransport gelten als morphologisch hoch aktive Küstenbereiche und variieren häufig in Ursprung, Genese und Entwicklung. Sie besitzen durch ihre dissipative Wirkungsweise eine bedeutende Schutzfunktion für rückwärtige Inseln, Halligen und Festlandbereiche und bilden vor der Küste Nordfrieslands die westliche Außengrenze des Wattenmeeres.
Ziel der nachfolgenden Studie ist es, anhand hochauflösender Georadarmessungen und sedimentologischer Daten aus Bohrungen die Landschaftsgeschichte an der Westküste Amrums sowie die Entstehung und interne sedimentäre Architektur des der Insel westlich angelagerten Kniepsandes zu untersuchen und die Prozesse, die zur Genese geführt haben, durch Datierungen zeitlich einzuordnen. Auf Grundlage der gewonnenen Daten wurden zwei stratigraphische Modelle entwickelt, welche die geologisch-geomorphologischen Prozesse und Sedimentationsbedingungen im Westküstenvorfeld erklären.
Ein Modell zeigt die Landschaftsentwicklung an der Westküste Amrums und beschreibt die Sedimentationsbedingungen, die im Vorfeld der Westküste herrschten, bevor der Kniepsand an die Insel heranwanderte. Auf der Landoberfläche des ertrinkenden saaleeiszeitlichen Geestkerns wurden zu Beginn der Flandrischen Transgression feinkörnige Sedimente eines Mischund Schlickwatts abgelagert. Es ist davon auszugehen, dass der damals noch weit vor der Küste Amrums liegende Kniepsand eine Barriere bildete und so an der heute hochenergetischen Westküste für strömungsberuhigte Sedimentationsbedingungen sorgte. Durch Erosion am Geestkern bildeten sich in unterschiedlichen Höhenpositionen fossile Kliffs, die dem damaligen Meeresspiegelstand entsprechen.
Ein weiteres Modell beschreibt den Andockmechanismus des Kniepsandes an die Insel Amrum. Durch die Anlagerung des ehemaligen Außensandes und den damit einhergehenden Sedimentinput wurden die Bedingungen für eine großflächige Dünenbildung geschaffen.
Comparing projections of future changes in runoff from hydrological and biome models in ISI-MIP
(2013)
Future changes in runoff can have important implications for water resources and flooding. In this study, runoff projections from ISI-MIP (Inter-sectoral Impact Model Intercomparison Project) simulations forced with HadGEM2-ES bias-corrected climate data under the Representative Concentration Pathway 8.5 have been analysed for differences between impact models. Projections of change from a baseline period (1981–2010) to the future (2070–2099) from 12 impacts models which contributed to the hydrological and biomes sectors of ISI-MIP were studied. The biome models differed from the hydrological models by the inclusion of CO2 impacts and most also included a dynamic vegetation distribution. The biome and hydrological models agreed on the sign of runoff change for most regions of the world. However, in West Africa, the hydrological models projected drying, and the biome models a moistening. The biome models tended to produce larger increases and smaller decreases in regionally averaged runoff than the hydrological models, although there is large inter-model spread. The timing of runoff change was similar, but there were differences in magnitude, particularly at peak runoff. The impact of vegetation distribution change was much smaller than the projected change over time, while elevated CO2 had an effect as large as the magnitude of change over time projected by some models in some regions. The effect of CO2 on runoff was not consistent across the models, with two models showing increases and two decreases. There was also more spread in projections from the runs with elevated CO2 than with constant CO2. The biome models which gave increased runoff from elevated CO2 were also those which differed most from the hydrological models. Spatially, regions with most difference between model types tended to be projected to have most effect from elevated CO2, and seasonal differences were also similar, so elevated CO2 can partly explain the differences between hydrological and biome model runoff change projections. Therefore, this shows that a range of impact models should be considered to give the full range of uncertainty in impacts studies.
The LPJ-GUESS dynamic vegetation model uniquely combines an individual- and patch-based representation of vegetation dynamics with ecosystem biogeochemical cycling from regional to global scales. We present an updated version that includes plant and soil N dynamics, analysing the implications of accounting for C-N interactions on predictions and performance of the model. Stand structural dynamics and allometric scaling of tree growth suggested by global databases of forest stand structure and development were well-reproduced by the model in comparison to an earlier multi-model study. Accounting for N cycle dynamics improved the goodness-of-fit for broadleaved forests. N limitation associated with low N mineralisation rates reduces productivity of cold-climate and dry-climate ecosystems relative to mesic temperate and tropical ecosystems. In a model experiment emulating free-air CO2 enrichment (FACE) treatment for forests globally, N-limitation associated with low N mineralisation rates of colder soils reduces CO2-enhancement of NPP for boreal forests, while some temperate and tropical forests exhibit increased NPP enhancement. Under a business-as-usual future climate and emissions scenario, ecosystem C storage globally was projected to increase by c. 10 %; additional N requirements to match this increasing ecosystem C were within the high N supply limit estimated on stoichiometric grounds in an earlier study. Our results highlight the importance of accounting for C-N interactions not only in studies of global terrestrial C cycling, but to understand underlying mechanisms on local scales and in different regional contexts.
Funded by the German Ministry for Education and Research (BMBF) a major research project called MiKlip (Mittelfristige Klimaprognose, Decadal Climate Prediction) was launched and global as well as regional predictive ensemble hindcasts have been generated. The aim of the project is to demonstrate for past climate change whether predictive models have the capability of predicting climate on time scales of decades. This includes the development of a decadal forecast system, on the one hand to support decision making for economy, politics and society for decadal time spans. On the other hand, the scientific aspect is to explore the feasibility and prospects of global and regional forecasts on decadal time scales. The focus of this paper lies on the description of the regional hindcast ensemble for Europe generated by COSMO-CLM and on the assessment of the decadal variability and predictability against observations. To measure decadal variability we remove the long term bias as well as the long term linear trend from the data. Further, we applied low pass filters to the original data to separate the decadal climate signal from high frequency noise. The decadal variability and predictability assessment is applied to temperature and precipitation data for the summer and winter half-year averages/sums. The best results have been found for the prediction of decadal temperature anomalies, i.e. we have detected a distinct predictive skill and reasonable reliability. Hence it is possible to predict regional temperature variability on decadal timescales, However, the situation is less satisfactory for precipitation. Here we have found regions showing good predictability, but also regions without any predictive skill.
To reconstruct ocean circulation changes during specific periods of Earth history, benthic and planktic foraminifera were used as proxies in the different parts of this thesis. Both studied time periods, the Late Cretaceous and the early Pleistocene, are characterized by long-term climate cooling and major changes in ocean circulation. The first part of this thesis concentrated in the Late Cretaceous. During the Late Cretaceous long-term cooling phase, benthic foraminiferal δ18O values show a positive shift lasting about 1.5 Myr (71.5–70 Ma). This shift can be observed on a global scale and has become known as the Campanian-Maastrichtian Boundary Event (CMBE). It is proposed that this δ18O excursion is influenced either by changing intermediate- to deep-water circulation or by temporal build-up of Antarctic ice sheets. Benthic foraminiferal assemblage counts from a southern high-latitudinal site near Antarctica (ODP Site 690) are analyzed to test if the influence of the CMBE on the benthic species composition. One of the two discussed hypotheses for the causation of the δ18O transition is a change in intermediate- to deep-water circulation from low-latitude to high-latitude water masses. This change would result in cooler temperatures, higher oxygen concentration, and possibly lower organic-matter flux at the seafloor, causing a major benthic foraminiferal assemblage change. Another possible explanation of the δ18O transition of the CMBE is significant ice formation on Antarctica. However no major benthic foraminiferal assemblage change would be expected in this case. The benthic foraminiferal assemblage of Site 690 shows a separation of the studied succession into two parts with significantly different species composition. The older part (73.0–70.5 Ma) is dominated by species, which are typical for lower bottom water oxygen concentration and more common in low-latitude assemblages. Species dominating the younger part (70.0–68.0 Ma) are indicators for well-oxygenated bottom waters and more common in high-latitude assemblages. This change in the benthic foraminiferal assemblages is interpreted to represent a shift of low-latitude toward high-latitude dominated intermediateto deep-water sources. A change in oceanic circulation was therefore at least a major component of the CMBE. The Pacific Ocean contributed significantly to the climatic development during the Late Cretaceous cooling period. The contribution of ocean circulation changes in the Pacific Ocean to the Late Cretaceous climatic development in general and the CMBE and Mid-Maastrichtian Event (MME) in particular, however, is poorly understood. Previously measured high resolution planktic and benthic stable isotope data and a neodymium (Nd) isotope record from the Pacific ODP Site 1210 (Shatsky Rise, tropical Pacific Ocean) for the Campanian to Maastrichtian (69.5 to 72.5 Ma) are used to reconstruct changes in surface- and bottom water temperatures as well as changes in the source region of deep- to intermediate waters [see Appendix 4; Jung et al. 2013]. The results of the benthic foraminiferal δ18O and Nd isotope records in combination with Nd isotope records from other studies indicate changes in the intensity of intermediate- to deep ocean circulation in the tropical Pacific across the Campanian-Maastrichtian interval [see Appendix 4; Jung et al. 2013]. During the early Maastrichtian (72.5 to 69.5 Ma), a three-million-year-long period of cooler conditions and a simultaneous change towards less radiogenic Nd isotope signatures is interpreted to represent a period of increased admixture and northward flow of deep waters from the Southern Ocean (Southern Component Water, SCW). This change was probably caused by an intensified formation of deep waters in the Southern Ocean. This was reduced again during the MME (69.5 to 68.5 Ma). This early Maastrichtian cold interval is similar to the CMBEδ13C fall and succeeding δ13C rise towards the MME and is therefore also interpreted to represent tectonically forced, long-term changes in the global carbon cycle and thus a tectonic forcing of the early Maastrichtian climate cooling. Overall, the Campanian-Maastrichtian Nd and stable isotope records of Shatsky Rise indicate changes in ocean circulation that are paralleled by global warming and cooling periods. The fluctuating strength of SCW contribution in the tropical Pacific points towards an increased respectively weakened ocean circulation, which is probably related to the strength of deep-water formation in the Southern Ocean [see Appendix 4; Jung et al. 2013]. For this study, the analysis of benthic foraminiferal assemblages of Site 1210 is carried out for the same time interval (69.5 to 72.5 Ma) as Nd and stable isotopes to evaluate the influence of intermediate- to deep ocean circulation changes on the benthic foraminiferal community. The possible reaction of benthic foraminiferal assemblages is compared to the results of stable isotope and neodymium isotopes. The observed changes in species abundances only partly reflect the circulation changes reconstructed with Nd and stable oxygen istopes. For example, Stensioina spp., Aragonia spp. and Lenticulina spp., cold-water preferring species, start to be increasingly abundant at the beginning of enhanced influence of SCW. However, their abundance pattern does not follow the varying strength of the cold SCW influence at Shatsky Rise. Other species prefer lesser oxygen concentrations and warmer bottom water, e.g. Paralabamina spp. and Globorotalites spp. Paralabamina spp. has its highest relativ abundance at the beginning of the studied succession, where the influence of SCW is small. However, this taxa occurs throughout the record, even though the influence of SCW increases. Globorotalites spp. is even most abundance after the CMBE, where bottom waters are till cold and influenced by SCW. This leads to the conclusion that the varying strength of SCW in the tropical Pacific at Shatsky Rise through the studied interval is not facilitating a significant faunal turnover as has been observed at the South Atlantic Site 690 (Chapter 3). These results of the benthic foraminiferal assemblage analysis suggest a rather minor influence of the SCW on the major environmental factors that are generally influencing benthic foraminiferal communities (e.g., oxygen concentration, organic matter flux to the sea floor, bottom-water temperature). The second major part of this thesis focused on the late Pliocene-earliest Pleistocene. The late Pliocene is characterized by a long-term global cooling trend resulting in a major increase of Arctic ice sheets from around 3 Ma onwards, culminating in the Plio-Pleistocene intensification of the Northern Hemisphere glaciation. At around 2.7 Ma, large amplitude glacial-interglacial excursions (~1‰ δ18O in benthic foraminiferal calcite) in benthic oxygen isotopes can be observed. Marine isotope stage (MIS) 100 at around 2.55 Ma is the first glacial, when widespread ice rafted debris has been found in sediments in the North Atlantic Ocean. To gain a deeper understanding of the climatic evolution of the latest Pliocene-early Pleistocene, it is necessary to improve the reconstructions of North Atlantic paleohydrography, as the North Atlantic provides a key region for global climate. The consequences of the intensification of Northern Hemisphere on the early Pleistocene North Atlantic thermocline stratification and intermediate waters are still poorly understood. However, surface hydrography, the history of the thermocline and development of North Atlantic intermediate waters are well-studied for the Last Glacial Maximum (LGM). These well-known mechanisms responsible for the LGM in comparison with the present-day interglacial North Atlantic are used as an analogue for te early Pleistocene glacialinterglacials cycles. In this study, suborbitally resolved stable oxygen and carbon isotope and Mg/Ca records are measured from a deep-dwelling planktic foraminifera (Globorotaliacrassaformis) from Integrated Ocean Drilling Program Site U1313 (North Atlantic, 41°N) covering marine oxygen isotope stages MIS 103 to 95 (early Pleistocene, 2.6 to 2.4 Ma). The results are interpreted to represent a change in intermediate-water masses on glacialinterglacial timescales. During glacials geochemical records in G. crassaformis (~500–1000 m) bear the imprint of Glacial North Atlantic Intermediate Water (GNAIW), while during interglacials this species reflects the signature of the influence of Mediterranean Outflow Water (MOW) in combination with the subtropical gyre. The comparison of this data with the published records from G. ruber from the same samples facilitates the reconstruction of glacial-interglacial stratification changes of the upper water column at Site U1313. The results show that larger gradients of temperature, salinity and δ13C prevailed during glacials, suggesting a stronger stratification of the upper water column. This can be seen to indicate glacial-interglacial changes in ntermediate water masses in the North Atlantic similar to those reconstructed for the latest Pleistocene. As an additional proxy, the clumped isotope paleothermometer is applied for the Late Cretaceous study as well as for the early Pleistocene. This proxy is commonly assumed to be independent of other factors than temperature. Clumped isotopes are measured for the Late Cretaceous Site 690 on the planktic foraminiferal species Archaeoglobigerina australis and compared to already existing stable oxygen isotopes of this species. This is assumed to enable the reconstruction of paleotemperature independent of ice volume and therefore contribute to the long-lasting discussion whether there was a temporal ice build-up on Antarctic during the Campanian-Maastrichtian cooling period. For the early Pleistocene, the planktic foraminiferal species G. crassaformis is used from Site U1313 from MIS 99 (interglacial) and MIS 98 (glacial). This provides the opportunity to separate ice volume, salinity and temperature effects on the measured δ18O record of G. crassaformis. The results of the clumped isotope measurements reveal comparatively large standard errors. For the Late Cretaceous the standard error of the clumped isotope measurements proved too large to allow any conclusions on the temperature component on the δ18O record of A. australis. For the early Pleistocene, the temperature difference is also too small to be reconstructed with the standard error of the clumped isotope measurements in this study. Measuring many replicates of one sample would minimize the standard error considerably. However, the amount necessary to measure replicates cannot be gained for either time period, as almost all foraminifera were picked from the respective samples. It is concluded that the respective questions may be solved with a different method of clumped isotope analysis requiring less sample material. This method is, for example, available at the ETH Zurich.
Mantle convection is the process by which heat from the Earth’s core is transferred upwards to the surface and it is accepted to explain the dynamics of the Earth’s interior. On geological time-scales, mantle material flows like a viscous fluid as a consequence of the buoyancy forces arising from thermal expansion. Indeed, mantel convection provides a framework which links together the major disciplines, such as seismology, mineral physics, geochemistry tectonic and geology. The numerical model has been applied to understand the dynamic, structure and evaluation of the Earth, and other terrestrial planets and the investigations continue to explore, different aspects of the mantle convection.
In fact, to model this phenomenon, two complementary approaches are possible. On the one hand, one can solve self-consistently the equations of thermal convection, including parameters and employing physical relationships derived from mineral physics. Our understanding of mantle convection depends ultimately upon the success of such fully self-consistent dynamic models in explaining observable features of the flow. Although, these models presently unable to predict the actual convection pattern of the Earth, they are extremely useful to investigate general characteristics of given physical systems. On the other hand, to permit comparison with specific observables associated with the flow, one can consider a more restricted problem. Instead of focusing on the time evolution of mantle flow, if we know a priori the temperature - and hence presumably the density - anomalies that drive the convection, we can try to build a snapshot of the present-day flow pattern, consistent with those anomalies, that can successfully predict the observables. As matter of fact, the aim of this study is to investigate both approaches in comparison with the main geophysical constraints on mantle structure. These constraints include the geoid anomalies, the dynamic surface and core-mantle boundary topography and tectonic plate motions.
The most appropriate mathematical basis functions for describing a bounded and continuous function on a spherical surface are spherical harmonics. We may therefore expand the geodynamic observables in terms of spherical harmonics. We have investigated two methods of the global spherical harmonic analysis by specific attention to the dynamic geoid computation of the geodynamic models. The first method is the quadrature method in which the loss of the orthogonality of the Legendre functions in transition from continues to discrete case is the major drawback to the method. Particularly, we showed that in the absence of the tesseral harmonics, quadrature formulation leads to obtain inaccurate results. The second method is the least-squares which can be considered as the best linear unbiased estimator that provides the exact results. We showed that even with a low resolution grid data it is possible to reconstruct the data and achieve an accurate result by using this method, which is extremely remarkable in three-dimensional global convection studies. However, special care has to be taken since there is some source of errors that might influence the efficiency of this method.
In general, to better understanding of the properties of the mantle, it is useful to assess observable characteristics of plumes in the mantle, including geoid, topography and heat flow anomalies. However, only few studies exist on geoid and topography for axi-symmetric convection and their models were restricted to isoviscous (or stratified) mantle and low Rayleigh numbers. We studied fully coupled depth and temperature dependent Arrhenius type of viscosity in axi-symmetric spherical shell geometry in order to investigate the shape of geoid anomalies and dynamic topography above a plume. Indeed, the topography and geoid anomalies produced from plumes are sensitive to rheology of the mantle and rheology of the plume; both have effects on shape and amplitude of the geoid anomalies. As results we are able to define different classes of plumes by their geoid signals.
Mainly depth-dependent viscosity models show a geoid with negative sign above the plume which can turn to the positive sign by decrease the viscosity contrast. This can be considered as a transition between the strongly depth dependent and the constant viscosity case. Our results basically support the idea by Morgan [1965] and McKenzie [1977]. They have shown the magnitude and even the sign of the total gravity anomaly depend on the spatial variation in effective viscosity. In addition, Hager [1984] has concluded that the total gravity field is depend on the radial distribution of effective viscosity, and a small change in viscosity contrast leads to varying sign of the response function.
In the case of temperature-dependent viscosity, the formation of an immobile lithosphere is a natural result, and the flow as well as the total geoid becomes strongly time dependent. When we increase the activation energy, all geoids associated with the first arriving plumes look like bell shaped whereas for typical plumes, after reaching a statistical steady state, bell-shaped geoids with decreasing amplitude as well as linear flank shaped geoids are observed. It is surprising that in spite of large differences in lateral and depth varying viscosities, the shapes of the geoid anomalies remained rather similar. We also identified different behaviors in the combined model with temperature-and pressure-dependent viscosity. In fact, in spite of the strongly different rheology, the geoid anomalies in all cases were surprisingly similar. Furthermore, we proposed a scaling law for the geoid which makes our results directly applicable to other planets. Moreover, we can apply the results of our calculation to find relations between different rheology and sub-lid temperature, since we know that the mantle temperature can change significantly with variation in pressure-temperature dependent viscosity. It is also possible to define a range of stagnant lid thickness related to the amplitude of the geoid which can be reasonable for study of the lid thickness in Venus or Mars.
Nevertheless, in these series of models, we simplified a number of complexities within the Earth. One of the most important of such simplification is the Boussinesq approximation. This approximation is valid if the temperature scale height (i.e. the depth over which temperature increases by a factor of “ ” due to adiabatic compression) is much greater than the convection depth. However, a temperature scale height in the Earth’s mantle is at best only slightly greater than the mantle depth. Hence, the Boussinesq approximation could mask some very important stratification and compressibility effects that influence both the spatial and temporal structure of the convection. Therefore, in more advance models we considered compressibility in our mantle convection models, assuming that density vary both radially and laterally, being determined as a function of pressure and temperature through an appropriate equation of the state. Moreover, thermodynamic properties assumed to be a function of depth.
We examined the details of the structure of the spherical axi-symmetric Anelastic Liquid Approximation model (ALA) with special attention to the Arrhenius rheology, and compare it to the cases of compressible convection without depth dependent thermodynamical properties, and to cases of the extended Boussinesq approximation. At the same time, the effects of the interaction between temperature and pressure-dependent viscosity and thermodynamic parameters in the compressible mantle convection on the geoid and topography have been studied. We showed that assuming compressible convection with depth-dependent thermodynamic properties strongly influence the geoid undulations. Using compressible convection with constant thermodynamic properties is physically inconsistent and may lead to spurious results for the geoid and convection pattern. Indeed, by a systematic study of different approaches of compressibility in the spherical shell convection for different Arrhenius viscosity laws we proved that only in the unrealistic case of zero activation energy the different compressibility modes result in comparable convection and geoid patterns. In all other rheological cases, large differences have been obtained, that stressing the important role of consistent compressible thermodynamic properties for mantle convection.
In addition, we examine the impact of compressibility as well as different rheologies on the power law relation that connects the Nusselt number to the Rayleigh number. We have discovered that the power law index of the relationship is controlled by the rheology, independent of which approximation is used. Instead, the bound of this relation is controlled by a combination of different approximation and rheology.
Next, instead of focusing on the time evolution of mantle flow, we have carried out three-dimensional spherical shell models of mantle circulation to investigate the effects of joint radial and lateral viscosity variations on the Earth’s non-hydrostatic geoid, surface and core-mantle boundary topographies. These models include realistic lateral viscosity variations (LVV) in the lithosphere, upper mantle and lower mantle in combination with different stratified viscosity structures. We have demonstrated that the contradictory results concerning the effects of LVV can be clarified by the most straight-forward problem in geoid modeling; namely, rather poorly known stratified viscosity structure. We explored three classes of dynamic geoid models due to lateral viscosity variations. In the first class, the LVV strongly improved the fit to the observed geoid. Indeed, when the viscosity contrast between lower and upper mantles is not large enough to produce a good fit to geoid the LVVs are able to perform this action by adjusting amplitudes, so that it becomes comparable with observation. In the second class, inducing the LVV moderately improved the fit. Actually, when the geoid induced by a stratified viscosity structure already has a good correlation with observation, then the LVV causes its amplitude to further improve. In the last class, if the viscosity contrast between upper and lower mantle would be high enough, inducing LVV deteriorate the fit to the observed geoid.. Indeed, depending on the stratified viscosity, inducing the LVV may take place in one of these categories.
We also quantified the effects of LVV in the mantle and lithosphere individually. We found that the presence of LVV in the mantle (upper and lower) improves the fit to the observed geoid regardless of stratified viscosity. While LVV in the lithosphere is a crucial parameter, and dependent of the stratified viscosity, may increase or decrease the geoid fit. In fact, when the lower mantle considers being viscous enough, it would support the negative buoyancy of subducting slabs. Thus, it transmits some of the stress back to the top boundary and causes a weak coupling between slab and surface. Therefore, by including the low viscous plate boundaries in this model, the slabs and overriding plates decouples and the fit to the observed geoid degrades. In contrast, when the lower mantle viscosity is not sufficiently stiff, the presence of the low viscous plate boundaries assists to weaken the strong mechanical coupling between slab and surface. Hence, a better fit achieved.
In this study we show how size-resolved measurements of aerosol particles and cloud condensation nuclei (CCN) can be used to characterize the supersaturation of water vapor in a cloud. The method was developed and applied for the investigation of a cloud event during the ACRIDICON-Zugspitze campaign (17 September to 4 October 2012) at the high-alpine research station Schneefernerhaus (German Alps, 2650 m a.s.l.). Number size distributions of total and interstitial aerosol particles were measured with a scanning mobility particle sizer (SMPS), and size-resolved CCN efficiency spectra were recorded with a CCN counter system operated at different supersaturation levels.
During the evolution of a cloud, aerosol particles are exposed to different supersaturation levels. We outline and compare different estimates for the lower and upper bounds (Slow, Shigh) and the average value (Savg) of peak supersaturation encountered by the particles in the cloud. For the investigated cloud event, we derived Slow ≈ 0.19–0.25%, Shigh ≈ 0.90–1.64% and Savg ≈ 0.38–0.84%. Estimates of Slow, Shigh and Savg based on aerosol size distribution data require specific knowledge or assumptions of aerosol hygroscopicity, which are not required for the derivation of Slow and Savg from the size-resolved CCN efficiency spectra.
Das Thema dieser Arbeit war die Untersuchung der natürlichen Variationen von den zwei primordialen Uranisotopen (238U und 235U) mit einem Schwerpunkt auf Proben, die (1) die kontinentale Kruste und ihre Verwitterungsprodukte (d.h. Granite, Shales und Flusswasser) repräsentieren, (2) Produkte der hydrothermalen Alteration vom mittelozeanischen Rücken widerspiegeln (d.h. alterierte Basalte, Karbonatgänge und hydrothermales Wasser) und (3) aus abgegrenzten euxinischen Becken (d.h. Proben aus der Wassersäule und den dazugehörigen Sedimenten) stammen. Das allgemeine Ziel war das Verständnis, unter welchen Bedingungen und Mechanismen eine Fraktionierung der zwei häufigsten Uranisotope (238U und 235U) in der Natur erfolgt, zu verbessern.
Die untersuchten Haupt- und Nebenflüsse unterscheiden sich sowohl in Ihrer Urankonzentration (c(U)) als auch in Ihrer Uranisotopenzusammensetzung (δ238U), wobei die Nebenflüsse eine geringere Urankonzentration (0.87 nmol/kg bis 3.08 nmol/kg) und eine schwerere Uranisotopenzusammensetzung aufweisen (-0.29 ‰ bis +0.01 ‰ im δ238U) im Vergleich zu den Hauptflüssen (c(U) = 5.19 nmol/kg bis 11.69 nmol/kg und d238U = -0.31 ‰ bis +0.13 ‰) aufweisen. Die untersuchten Gesteinsproben fallen alle in einen recht schmalen Bereich von δ238U, zwischen -0.45 ‰ und -0.21 ‰, mit einem Durchschnittswert von -0.30 ‰ ± 0.04 ‰ (doppelte Standardabweichung). Deren Uranisotopenvariationen sind unabhängig von der Urankonzentration (11.8 µg/g bis 1.3 µg/g), dem Alter (3.80 Ga bis 328 Ma), der Probenlokalität und Grad der Differenzierung. Basierend auf den Ergebnissen der Hauptflüsse, die die Uranhauptquelle für den Ozean darstellen, schlagen wir für zukünftige Berechnungen in der Massenbilanz des Urans einen neuen Wert als beste Abschätzung für die Quelle des Urans im Ozean vor, δ238U = -0.23 ‰.
Die Produkte der hydrothermalen Alteration, alterierte Basalte und Kalziumkarbonatgänge, zeigten etwas stärkere Isotopenvariationen (δ238U zwischen -0.63 ‰ und +0.27 ‰) als erwartet und die hydrothermalen Fluide wiesen eine etwas leichtere Uranisotopenzusammensetzung als Meerwasser ((-0.43 ± 0.25) ‰ vs. (-0.37 ± 0.03) ‰) auf. Diese Ergebnisse sind in Übereinstimmung mit einem Modell, dass annimmt, dass die beobachtete Isotopenfraktionierung hauptsächlich ein Ergebnis von Redoxprozessen ist, z.B. die partielle Reduktion von löslichem UVI aus dem Meerwasser während der hydrothermalen Alteration, was zu einer Anreicherung der schweren Uranisotope in der reduzierten Uranspezies (UIV) führt und 2) das bevorzugte Entfernen von UIV aus den hydrothermalen Fluid und der Einbau in die alterierte ozeanische Kruste. Durch diesen Prozess wird das hydrothermale Fluid an schweren Uranisotopen verarmt und somit würden auch die alterierten Basalte und Karbonate ein niedriges δ238U aufweisen, wenn sie mit dem isotopisch leichten hydrothermalen Fluid in Kontakt gekommen sind.
Die Untersuchung von Wasser- und Sedimentproben aus der Ostsee und dem anoxischen Kyllaren Fjord (Norwegen) auf deren Uran- und Mo-Isotopenzusammensetzung zeigte, dass die Uranisotopenzusammensetzung der Sedimente abhängt von (1) dem Ausmaß des Uranaustrags aus der Wassersäule (in einer ähnlichen Art und Weise wie bei den Molybdänisotopen) und (2) der Sedimentationsrate, d.h. der Fraktion von authigenem- relativ zum dedritischen Uran in den Sedimenten. Aufgrund der hohen Sedimentationsrate zeigen die Sedimente aus dem Kyllaren Fjord nur eine moderate authigene Urananreicherung und eine leichtere Uranisotopenzusammensetzung als Sedimente aus dem Schwarzen Meer. In den anoxischen Becken der Ostsee erfolgt dagegen eine starke Mo- und schwache U-Isotopenfraktionierung zwischen Wasser und Sediment. Durch die regelmäßigen auftretenden Spülereignisse mit sauerstoffreichem Wasser wurden vermutlich die ursprünglichen anoxischen Mo- und U-Isotopensignaturen der Sedimente verändert. Demzufolge müssen die Sedimente durchgehend anoxischen Bedingungen ausgesetzt sein, um eine Mo- und U-Isotopensignatur von den Redoxbedingungen während der Ablagerungen zu speichern.
Der Vergleich zwischen Molybdän- und Uranisotopen in der Ostsee und dem anoxischen Kyllaren Fjord zeigte, dass sich Uran- und Molybdänisotope in stark euxinischen Wassersäulen (c(H2S) > 11 µmol/L) entgegengesetzt verhalten. Dementsprechend ergänzen sich die beiden Isotopensysteme und können genutzt werden, um die Ablagerungsbedingungen in abgeschlossenen Becken und die Redoxentwicklung des Paläoozeans zu untersuchen.
Gridded maps of meteorological variables are needed for the evaluation of weather and climate models and for climate change monitoring. In order to produce them, values at locations where no observing stations are available need to be estimated from point-wise observations. For the interpolation of meteorological observations deterministic and stochastic methods are often combined. Deterministic methods can account for ancillary information such as elevation, continentality or satellite observations. Stochastic methods such as kriging reproduce observed values at the station locations and also account for spatial variability. In the first two studies of this thesis, a flexible interpolation method for the gridding of locally observed daily extreme temperatures is developed that also provides an optimal estimate of the interpolation ncertainty. In the third study, an observational dataset is created using this interpolation method and then applied to evaluate a climate simulation for Africa.
In the first study, the Regression-Kriging-Kriging (RKK) method is tested for the interpolation of daily minimum and maximum temperatures (Tmin and Tmax) in different regions in Europe. RKK accounts for elevation, continentality index and zonal mean temperature and is applicable in regions of differing station density and climate. The accuracy of RKK is compared to Inverse Distance Weighting, a common deterministic interpolation method, and to Ordinary Kriging, a common stochastic interpolation method. The first step in RKK is to use regression kriging, in which multiple linear regression accounts for topographical effects on the temperature field and kriging minimizes the regression error, to interpolate climatological means. In the second step daily deviations from the monthly climatology are interpolated using simple kriging. Owing to the large climatological differences across the investigation area the interpolation is performed in homogeneous subregions defined according to the Köppen-Geiger climate classification. Cross validation demonstrates the superiority of RKK over the simpler algorithms in terms of accuracy and preservation of spatial variability. The interpolation performance however strongly varies across Europe, being considerably higher over Central Europe (highest station density) than over Greenland (few stations along the coast line). This illustrates the strong impact of the station density on the accuracy of the interpolation result. Satellites provide comprehensive observations of climate variables such as land surface temperature (LST) and cloud cover (CC). However, LST is associated with high uncertainty (standard error ~ 1-2°C), preventing its direct application in meteorology and climatology. The second study investigates the usefulness of LST and CC as predictors for the gridding of daily Tmin and Tmax. The RKK algorithm is compared with similar interpolation methods that apply LST and CC in addition to the predictors used with the RKK algorithm. The investigation is conducted in two regions, Central Europe and the Iberian Peninsula, which differ strongly in average cloud cover (Central Europe is approximately 30% cloud free and the Iberian Peninsula approximately 60 % cloud free). RKKLST (in which monthly mean LST is used as an additional predictor) yields for Central Europe no clear improvement over RKK, yet it reduces the interpolation error over the Iberian Peninsula. This finding can be explained by the higher percentage of cloud free pixels over that region in summer which enables a more robust determination of monthly mean LST. Adding a regression step for daily anomalies (using the predictor CC) yields the RKRK method and improves the preservation of spatial variability over the Iberian Peninsula. Moreover, a successive reduction of the station number (from 140 to 10 stations) reveals an increasing superiority of RKKLST and RKRK over RKK in both regions.
The application of a gridded observational dataset for climate monitoring or climate model validation requires knowledge of the uncertainties associated with the dataset. The estimation of the interpolation uncertainty, here the inter quartile range is the used uncertainty measure, is therefore an important issue within the frame of this thesis. By means of cross validation it is shown that the largest uncertainties occur in regions of low station density (e.g. Greenland), in mountainous regions and along coastlines (in these regions model evaluation results should be interpreted carefully). The magnitude of the interpolation error mainly depends on the station density, while the complexity of terrain has substantially less influence. On average over all regions and investigation days the target precision of the uncertainty estimate is reached. However, on local scales and for single days it can be clearly over- or underestimated. The application of satellite-derived predictors (LST and CC) yields no noteworthy improvement of the uncertainty estimate.
In the last study two regional climate simulations for Africa using the ERA-Interim driven COSMO-CLM (CCLM) model at two different horizontal resolutions (0.22° and 0.44°) are validated. It is assessed whether observed patterns and statistical properties of daily Tmin and Tmax are correctly represented in the model. The ERA-Interim reanalysis and a specially created observational dataset are used as reference. The observational dataset is generated by applying the RKRK algorithm (developed within the second study). The investigations show an occasionally large bias in Tmin and Tmax. The hemispheric summers are generally too warm and the temporal variability in temperature is too high, particularly over extra tropical Africa. The diurnal temperature range is overestimated by about 2°C in the northern subtropics but underestimated by about 2°C over large parts of the African tropics. CCLM reproduces the observed frequency distribution of daily Tmin and Tmax in all African climate regions, and the extreme values in the lower percentiles (5, 10, 20%) for Tmin are well simulated. The higher percentiles (80, 90, 95%) for Tmax are however overestimated by 2-5°C. For both Tmin and Tmax the 0.22° simulation is on average 0.5°C warmer than the 0.44° simulation. Additionally, the higher percentiles are about 1°C warmer for both Tmin and Tmax in the higher resolution run, while the lower percentiles in both runs match very well. Although the temperature pattern is represented in more detail along the coastlines and in topographically complex regions, the higher resolution simulation yields no qualitative improvement.
To summarize, the choice of the appropriate algorithm mainly depends on the interpolation conditions. In cases where the station density is high across the target region and the predictor space is adequately covered by observing stations, the computationally less demanding RK algorithm should be preferred. In regions where the station density is low the more robust RKRK algorithm should be the first choice. Due to the strong physical relation of both CC and LST to Tmin and Tmax the missing information is at least partially compensated for. The estimation of the interpolation uncertainty could be improved by applying a normal score transformation to the data prior to a kriging step. This is because the kriging assumption that the increments of the variable of interest are second order stationary can be approximately met by a normal score transformation.
This study presents a method for adjusting long-term climate data records (CDRs) for the integrated use with near-real-time data using the example of surface incoming solar irradiance (SIS). Recently, a 23-year long (1983–2005) continuous SIS CDR has been generated based on the visible channel (0.45–1 μm) of the MVIRI radiometers onboard the geostationary Meteosat First Generation Platform. The CDR is available from the EUMETSAT Satellite Application Facility on Climate Monitoring (CM SAF). Here, it is assessed whether a homogeneous extension of the SIS CDR to the present is possible with operationally generated surface radiation data provided by CM SAF using the SEVIRI and GERB instruments onboard the Meteosat Second Generation satellites. Three extended CM SAF SIS CDR versions consisting of MVIRI-derived SIS (1983–2005) and three different SIS products derived from the SEVIRI and GERB instruments onboard the MSG satellites (2006 onwards) were tested. A procedure to detect shift inhomogeneities in the extended data record (1983–present) was applied that combines the Standard Normal Homogeneity Test (SNHT) and a penalized maximal T-test with visual inspection. Shift detection was done by comparing the SIS time series with the ground stations mean, in accordance with statistical significance. Several stations of the Baseline Surface Radiation Network (BSRN) and about 50 stations of the Global Energy Balance Archive (GEBA) over Europe were used as the ground-based reference. The analysis indicates several breaks in the data record between 1987 and 1994 probably due to artefacts in the raw data and instrument failures. After 2005 the MVIRI radiometer was replaced by the narrow-band SEVIRI and the broadband GERB radiometers and a new retrieval algorithm was applied. This induces significant challenges for the homogenisation across the satellite generations. Homogenisation is performed by applying a mean-shift correction depending on the shift size of any segment between two break points to the last segment (2006–present). Corrections are applied to the most significant breaks that can be related to satellite changes. This study focuses on the European region, but the methods can be generalized to other regions. To account for seasonal dependence of the mean-shifts the correction was performed independently for each calendar month. In comparison to the ground-based reference the homogenised data record shows an improvement over the original data record in terms of anomaly correlation and bias. In general the method can also be applied for the adjustment of satellite datasets addressing other variables to bridge the gap between CDRs and near-real-time data.
When studying new particle formation, the uncertainty in determining the "true" nucleation rate is considerably reduced when using condensation particle counters (CPCs) capable of measuring concentrations of aerosol particles at sizes close to or even at the critical cluster size (1–2 nm). Recently, CPCs able to reliably detect particles below 2 nm in size and even close to 1 nm became available. Using these instruments, the corrections needed for calculating nucleation rates are substantially reduced compared to scaling the observed formation rate to the nucleation rate at the critical cluster size. However, this improved instrumentation requires a careful characterization of their cut-off size and the shape of the detection efficiency curve because relatively small shifts in the cut-off size can translate into larger relative errors when measuring particles close to the cut-off size.
Here we describe the development of two continuous-flow CPCs using diethylene glycol (DEG) as the working fluid. The design is based on two TSI 3776 counters. Several sets of measurements to characterize their performance at different temperature settings were carried out. Furthermore, two mixing-type particle size magnifiers (PSM) A09 from Airmodus were characterized in parallel. One PSM was operated at the highest mixing ratio (1 L min−1 saturator flow), and the other was operated in a scanning mode, where the mixing ratios are changed periodically, resulting in a range of cut-off sizes. The mixing ratios are determined by varying the saturator flow, where the aerosol flow stays constant at 2.5 L min−1. Different test aerosols were generated using a nano-differential mobility analyser (nano-DMA) or a high-resolution DMA, to obtain detection efficiency curves for all four CPCs. One calibration setup included a high-resolution mass spectrometer (APi-TOF) for the determination of the chemical composition of the generated clusters. The lowest cut-off sizes were achieved with negatively charged ammonium sulfate clusters, resulting in cut-offs of 1.4 nm for the laminar flow CPCs and 1.2 and 1.1 nm for the PSMs. A comparison of one of the laminar-flow CPCs and one of the PSMs measuring ambient and laboratory air showed good agreement between the instruments.
In this study we report the set-up of a novel twin chamber technique that uses the comparative method and establishes an appropriate connection of atmospheric and laboratory methods to broaden the tools for investigations. It is designed to study the impact of certain parameters and gases on ambient processes, such as particle formation online, and can be applied in a large variety of conditions. The characterisation of both chambers proved that both chambers operate identically, with a residence time xT (COMPASS1) = 26.5 ± 0.3 min and xT (COMPASS2) = 26.6 ± 0.4 min, at a typical flow rate of 15 L min−1 and a gas leak rate of (1.6 ± 0.8) × 10−5 s−1. Particle loss rates were found to be larger (due to the particles' stickiness to the chamber walls), with an extrapolated maximum of 1.8 × 10−3 s−1 at 1 nm, i.e. a hundredfold of the gas leak rate. This latter value is associated with sticky non-volatile gaseous compounds, too. Comparison measurement showed no significant differences. Therefore operation under atmospheric conditions is trustworthy. To indicate the applicability and the benefit of the system, a set of experiments was conducted under different conditions, i.e. urban and remote, enhanced ozone and terpenes as well as reduced sunlight. In order to do so, an ozone lamp was applied to enhance ozone in one of two chambers; the measurement chamber was protected from radiation by a first-aid cover and volatile organic compounds (VOCs) were added using a small additional flow and a temperature-controlled oven. During the elevated ozone period, ambient particle number and volume increased substantially at urban and remote conditions, but by a different intensity. Protection of solar radiation displayed a clear negative effect on particle number, while terpene addition did cause a distinct daily pattern. E.g. adding β pinene particle number concentration rose by 13% maximum at noontime, while no significant effect was observable during darkness. Therefore, the system is a useful tool for investigating local precursors and the details of ambient particle formation at surface locations as well as potential future feedback processes.
The anthropogenic influence on climate and environment has increased strongly since industrialization about 150 yr ago. The consequences for the atmosphere became more and more apparent and nowadays affect our life quality on Earth progressively. Because of that it is very important to understand the atmospheric processes, on which these effects are based on, in detail. In this study we report the set-up of a novel twin chamber technique that uses the comparative method and establishes an appropriate connection of atmospheric and laboratory methods to broaden the tools for investigations. It is designed to study the impact of certain parameters and gases on ambient processes such as particle formation online and can be applied in a large variety of conditions. The characterisation of both chambers proved that both chambers operate identically with a residence time (xT (COMPASS 1) = 26.5 ± 0.3 min and xT (COMPASS 2) = 26.6 ± 0.4 min) at a typical flow rate of 15 L min−1 and a deposition rate (1.6 ± 0.8) × 10−5 s−1. Comparison measurement showed no significant differences. Therefore operation under atmospheric conditions is trustworthy. To indicate the applicability and the benefit of the system a set of experiments was conducted at different conditions, i.e. urban and remote, enhancing ozone and terpenes as well as reducing sunlight. In the ozone enhanced ambient particle number and volume increased substantially at urban and remote conditions in a different strength. Solar radiation displayed a clear positive effect on particle number as well as terpene addition did at remote conditions. Therefore the system is a useful tool to investigate local precursors, the details of ambient particle formation at surface locations as well as future feedback processes.
Planktonic Foraminifera are important marine calcifiers, and the ongoing change in the oceanic carbon system makes it essential to understand the influence of environmental factors on the biomineralisation of their shells. The amount of calcite deposited by planktonic Foraminifera during calcification has been hypothesized to reflect a range of environmental factors. However, it has never been assessed whether their calcification only passively responds to the conditions of the ambient seawater or whether it reflects changes in resource allocation due to physiological stress. To disentangle these two end-member scenarios, an experiment is required where the two processes are separated. A natural analogue to such an experiment occurred during the deposition of the Mediterranean sapropels, where large changes in surface water composition and stratification at the onset of the sapropel deposition were decoupled from local extinctions of planktonic Foraminifera species. We take advantage of this natural experiment and investigate the reaction of calcification intensity, expressed as size-normalized weight (SNW), of four species of planktonic Foraminifera to changing conditions during the onset of Sapropel S5 (126–121 ka) in a sediment core from the Levantine Basin. We observe a significant relationship between SNW and surface water properties, as reflected by stable isotopes in the calcite of Foraminifera shells, but we failed to observe any reaction of calcification intensity on ecological stress during times of decreasing abundance culminating in local extinction. The reaction of calcification intensity to surface water perturbation at the onset of the sapropel was observed only in surface dwelling species, but all species calcified more strongly prior to the sapropel deposition and less strongly within the sapropel than at comparable conditions during the present day. These results indicate that the high-salinity environment of the glacial Mediterranean Sea prior to sapropel deposition induced a more intense calcification, whereas the freshwater injection to the surface waters associated with sapropel deposition inhibited calcification. The results are robust to changes in carbonate preservation and collectively imply that changes in normalized shell weight in planktonic Foraminifera should reflect mainly abiotic forcing.
Planktonic Foraminifera are important marine calcifiers, and the ongoing change in the oceanic carbon system makes it essential to understand the influence of environmental factors on the biomineralization of their shells. The amount of calcite deposited by planktonic Foraminifera during calcification has been hypothesized to reflect a range of environmental factors. However, it has never been assessed whether their calcification only passively responds to the conditions of the ambient seawater or whether it reflects changes in resource allocation due to physiological stress. To disentangle these two end-member scenarios, an experiment is required where the two processes are separated. A natural analogue to such an experiment occurred during the deposition of the Mediterranean sapropels, where large changes in surface water composition and stratification at the onset of the sapropel deposition were decoupled from local extinctions of planktonic Foraminifera species. We took advantage of this natural experiment and investigated the reaction of calcification intensity, expressed as mean area density (MAD), of four species of planktonic Foraminifera to changing conditions during the onset of Sapropel S5 (126–121 ka) in a sediment core from the Levantine Basin. We observed a significant relationship between MAD and surface water properties, as reflected by stable isotopes in the calcite of Foraminifera shells, but we failed to observe any reaction of calcification intensity on ecological stress during times of decreasing abundance culminating in local extinction. The reaction of calcification intensity to surface water perturbation at the onset of the sapropel was observed only in surface-dwelling species, but all species calcified more strongly prior to the sapropel deposition and less strongly within the sapropel than at similar conditions during the present-day. These results indicate that the high-salinity environment of the glacial Mediterranean Sea prior to sapropel deposition induced a~more intense calcification, whereas the freshwater injection to the surface waters associated with sapropel deposition inhibited calcification. The results are robust to changes in carbonate preservation and collectively imply that changes in normalized shell weight in planktonic Foraminifera should reflect mainly abiotic forcing.
In old and heavily weathered soils, the availability of P might be so small that the primary production of plants is limited. However, plants have evolved several mechanisms to actively take up P from the soil or mine it to overcome this limitation. These mechanisms involve the active uptake of P mediated by mycorrhiza, biotic de-occlusion through root clusters, and the biotic enhancement of weathering through root exudation. The objective of this paper is to investigate how and where these processes contribute to alleviate P limitation on primary productivity. To do so, we propose a process-based model accounting for the major processes of the carbon, water, and P cycle including chemical weathering at the global scale. We use simulation experiments to assess the relative importance of the different uptake mechanisms to alleviate P limitation on biomass production. Implementing P limitation on biomass synthesis allows the assessment of the efficiencies of biomass production across different ecosystems.
We find that active P-uptake is an essential mechanism for sustaining P availability on long time scales, whereas biotic de-occlusion might serve as a buffer on time scales shorter than 10 000 yr. Although active P uptake is essential for reducing P losses by leaching, humid lowland soils reach P limitation after around 100 000 yr of soil evolution. Given the generalized modeling framework, our model results compare reasonably with observed or independently estimated patterns and ranges of P concentrations in soils and vegetation. Furthermore, our simulations suggest that P limitation might be an important driver of biomass production efficiency (the fraction of the gross primary productivity used for biomass growth), and that vegetation on older soils becomes P-limited leading to a smaller biomass production efficiency.
With this study, we provide a theoretical basis for investigating the responses of terrestrial ecosystems to P availability linking geological and ecological time scales under different environmental settings.
A satellite-based climate record of monthly mean surface solar irradiance (SIS) is investigated with regard to possible inhomogeneities in time. The data record is provided by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Climate Monitoring (CM SAF) for the period of 1983 to 2005, covering a disk area between ±70° in latitude and longitude. The Standard Normal Homogeneity Test (SNHT) and two other homogeneity tests are applied with and without the use of reference SIS data (from the Baseline Surface Radiation Network (BSRN) and from the ECMWF (European Centre for Medium-Range Weather Forecasts) ERA -Interim reanalysis. The focus is on the detection of break-like inhomogeneities, which may occur due to satellite or SIS retrieval algorithm changes. In comparison with the few suitable BSRN SIS observation series with limited extension in time (no data before 1992), the CM SAF SIS time series do not show significant inhomogeneities, even though slight discrepancies in the surface measurements appear. The investigation of the full CM SAF SIS domain reveal inhomogeneities related to most of the documented satellite and retrieval changes, but only for relatively small domain fractions (especially in mountainous desert-like areas in Africa). In these regions the retrieval algorithm is not capable of adjusting for the changes of the satellite instruments. For other areas, e.g., Europe, no such breaks in the time series are found. We conclude that the CM SAF SIS data record has to be further assessed and regionally homogenized before climate trend investigations can be conducted.
Emissions of halogenated very short-lived substances (VSLS) are poorly constrained. However, their inclusion in global models is required to simulate a realistic inorganic bromine (Bry) loading in both the troposphere, where bromine chemistry perturbs global oxidizing capacity, and in the stratosphere, where it is a major sink for ozone (O3). We have performed simulations using a 3-D chemical transport model (CTM) including three top-down and a single bottom-up derived emission inventory of the major brominated VSLS bromoform (CHBr3) and dibromomethane (CH2Br2). We perform the first concerted evaluation of these inventories, comparing both the magnitude and spatial distribution of emissions. For a quantitative evaluation of each inventory, model output is compared with independent long-term observations at National Oceanic and Atmospheric Administration (NOAA) ground-based stations and with aircraft observations made during the NSF HIAPER Pole-to-Pole Observations (HIPPO) project. For CHBr3, the mean absolute deviation between model and surface observation ranges from 0.22 (38%) to 0.78 (115%) parts per trillion (ppt) in the tropics, depending on emission inventory. For CH2Br2, the range is 0.17 (24%) to 1.25 (167%) ppt. We also use aircraft observations made during the 2011 "Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere" (SHIVA) campaign, in the tropical West Pacific. Here, the performance of the various inventories also varies significantly, but overall the CTM is able to reproduce observed CHBr3 well in the free troposphere using an inventory based on observed sea-to-air fluxes. Finally, we identify the range of uncertainty associated with these VSLS emission inventories on stratospheric bromine loading due to VSLS (BryVSLS). Our simulations show BryVSLS ranges from ~ 4.0 to 8.0 ppt depending on the inventory. We report an optimised estimate at the lower end of this range (~ 4 ppt) based on combining the CHBr3 and CH2Br2 inventories which give best agreement with the compilation of observations in the tropics.
Emissions of halogenated very short-lived substances (VSLS) are poorly constrained. However, their inclusion in global models is required to simulate a realistic inorganic bromine (Bry) loading in both the troposphere, where bromine chemistry perturbs global oxidising capacity, and in the stratosphere, where it is a major sink for ozone (O3). We have performed simulations using a 3-D chemical transport model (CTM) including three top-down and a single bottom-up derived emission inventory of the major brominated VSLS bromoform (CHBr3) and dibromomethane (CH2Br2). We perform the first concerted evaluation of these inventories, comparing both the magnitude and spatial distribution of emissions. For a quantitative evaluation of each inventory, model output is compared with independent long-term observations at National Oceanic and Atmospheric Administration (NOAA) ground-based stations and with aircraft observations made during the NSF (National Science Foundation) HIAPER Pole-to-Pole Observations (HIPPO) project. For CHBr3, the mean absolute deviation between model and surface observation ranges from 0.22 (38%) to 0.78 (115%) parts per trillion (ppt) in the tropics, depending on emission inventory. For CH2Br2, the range is 0.17 (24%) to 1.25 (167%) ppt. We also use aircraft observations made during the 2011 Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere (SHIVA) campaign, in the tropical western Pacific. Here, the performance of the various inventories also varies significantly, but overall the CTM is able to reproduce observed CHBr3 well in the free troposphere using an inventory based on observed sea-to-air fluxes. Finally, we identify the range of uncertainty associated with these VSLS emission inventories on stratospheric bromine loading due to VSLS (BryVSLS). Our simulations show BryVSLS ranges from ~4.0 to 8.0 ppt depending on the inventory. We report an optimised estimate at the lower end of this range (~4 ppt) based on combining the CHBr3 and CH2Br2 inventories which give best agreement with the compilation of observations in the tropics.
The SAWNUC (Sulphuric Acid Water NUCleation) microphysical aerosol nucleation model is used to study the effect of reactor walls on the interpretation of nucleation experiments with respect to nucleation theory. This work shows that loss processes, such as wall losses, influence the interpretation of nucleation experiments, especially at low growth rates and short lifetimes of freshly nucleated particles. In these cases the power dependency of the formation rates, determined at a certain particle size, with respect to H2SO4 does not correspond to the approximate number of H2SO4 molecules in the critical cluster as expected by the first nucleation theorem. Observed ∂log(J)/∂log([H2SO4]) therefore can vary widely for identical nucleation conditions but different sink terms.
We investigate the contribution of oceanic methyl iodide (CH3I) to the stratospheric iodine budget. Based on CH3I measurements from three tropical ship campaigns and the Lagrangian transport model FLEXPART, we provide a detailed analysis of CH3I transport from the ocean surface to the cold point in the upper tropical tropopause layer (TTL). While average oceanic emissions differ by less than 50% from campaign to campaign, the measurements show much stronger variations within each campaign. A positive correlation between the oceanic CH3I emissions and the efficiency of CH3I troposphere–stratosphere transport has been identified for some cruise sections. The mechanism of strong horizontal surface winds triggering large emissions on the one hand and being associated with tropical convective systems, such as developing typhoons, on the other hand, could explain the identified correlations. As a result of the simultaneous occurrence of large CH3I emissions and strong vertical uplift, localized maximum mixing ratios of 0.6 ppt CH3I at the cold point have been determined for observed peak emissions during the SHIVA (Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere)-Sonne research vessel campaign in the coastal western Pacific. The other two campaigns give considerably smaller maxima of 0.1 ppt CH3I in the open western Pacific and 0.03 ppt in the coastal eastern Atlantic. In order to assess the representativeness of the large local mixing ratios, we use climatological emission scenarios to derive global upper air estimates of CH3I abundances. The model results are compared with available upper air measurements, including data from the recent ATTREX and HIPPO2 aircraft campaigns. In the eastern Pacific region, the location of the available measurement campaigns in the upper TTL, the comparisons give a good agreement, indicating that around 0.01 to 0.02 ppt of CH3I enter the stratosphere. However, other tropical regions that are subject to stronger convective activity show larger CH3I entrainment, e.g., 0.08 ppt in the western Pacific. Overall our model results give a tropical contribution of 0.04 ppt CH3I to the stratospheric iodine budget. The strong variations in the geographical distribution of CH3I entrainment suggest that currently available upper air measurements are not representative of global estimates and further campaigns will be necessary in order to better understand the CH3I contribution to stratospheric iodine.
As part of the CLACE-6 campaign we performed size-resolved CCN measurements for a~supersaturation range of S = 0.079 % to 0.66% at the high-alpine research station Jungfraujoch, Switzerland, in March~2007. The derived effective hygroscopicity parameter κ describing the influence of particle composition on CCN activity was on average 0.23–0.30 for Aitken (50–100 nm) and 0.32–0.43 for accumulation mode particles (100–200 nm). The campaign average value of κ = 0.3 is similar to the average value of κ for other continental locations. When air masses came from southeasterly directions crossing the Po Valley in Italy, particles were much more hygroscopic (κ ≈ 0.42) due to large sulfate mass fractions. The κ values obtained at S = 0.079 % exhibited a good negative correlation with the organic mass fractions derived from PM1 aerosol mass spectrometer (AMS) measurements. Applying a simple mixing rule the organic and inorganic mass fractions observed by the AMS could be used to reproduce the temporal fluctuations of the hygroscopicity of accumulation mode particles quite well.
We show how during a cloud event the aerosol particles were activated as cloud droplets and then removed from the air by precipitation leaving behind only a small amount of accumulation mode particles consisting mainly of weakly CCN-active particles, most likely externally mixed unprocessed soot particles.
During the campaign we had the opportunity to directly compare two DMT CCN counters for a certain time. The total CCN concentration (NCCN,tot) obtained by the two instruments at equal supersaturations agreed well for both possible operating modes: detecting NCCN,tot directly by sampling the polydisperse aerosol with the CCNC, or indirectly by combining size-resolved measurements of the activated fraction with parallel measurements of the particle size distribution (e.g., by SMPS). However, some supersaturation setpoints differed between the two CCNCs by as much as 20% after applying the instrument calibrations, which resulted in differences of the corresponding NCCN,tot of up to 50%. This emphasizes that it is extremely important to carefully calibrate the supersaturation of the instrument, especially at low S.
The link between atmospheric radicals and newly formed particles at a spruce forest site in Germany
(2013)
It has been claimed for more than a century that atmospheric new particle formation is primarily influenced by the presence of sulphuric acid. However, the activation process of sulphuric acid related clusters into detectable particles is still an unresolved topic. In this study we focus on the PARADE campaign measurements conducted during August/September 2011 at Mt. Kleiner Feldberg in central Germany. During this campaign a set of radicals, organic and inorganic compounds and oxidants and aerosol properties were measured or calculated. We compared a range of organic and inorganic nucleation theories, evaluating their ability to simulate measured particle formation rates at 3 nm in diameter (J3) for a variety of different conditions. Nucleation mechanisms involving only sulphuric acid tentatively captured the observed noon-time daily maximum in J3, but displayed an increasing difference to J3 measurements during the rest of the diurnal cycle. Including large organic radicals, i.e. organic peroxy radicals (RO2) deriving from monoterpenes and their oxidation products in the nucleation mechanism improved the correlation between observed and simulated J3. This supports a recently proposed empirical relationship for new particle formation that has been used in global models. However, the best match between theory and measurements for the site of interest was found for an activation process based on large organic peroxy radicals and stabilized Criegee intermediates (sCI). This novel laboratory derived algorithm simulated the daily pattern and intensity of J3 observed in the ambient data. In this algorithm organic derived radicals are involved in activation and growth and link the formation rate of smallest aerosol particles with OH during daytime and NO3 during nighttime. Because of the RO2s lifetime is controlled by HO2 and NO we conclude that peroxy radicals and NO seem to play an important role for ambient radical chemistry not only with respect to oxidation capacity but also for the activation process of new particle formation. This is supposed to have significant impact of atmospheric radical species on aerosol chemistry and should to be taken into account when studying the impact of new particles in climate feedback cycles.
This thesis aims to analyse in a first step the physical and chemical properties of soil profiles along pedomorphological transects in different land used conditions (protected, partly protected as well as cultivated and pastured areas) in North West Benin and in South East Burkina Faso. The information about soils, which are carried out in consideration of the pedogenesis processes like weathering types, saprolitisation, formation of laterite crusts and denudation within the planation surfaces are therefore correlated in a second step with the structure and dynamic of woody plant around individual soil profiles. The relationship soil properties and woody plant is investigated in order to assess the reciprocal influence between the diversity of woody plants and soil characteristics within a small scale study and under different land use conditions.
A common vertical and lateral differentiation of physical and chemical properties regardless of the partly protected, protected and cultivated status of the sites can be noticed. Thus, in the cultivated site of Kikideni and in the partly protected zone of Natiabouani (South East Burkina Faso) sandy loam and sandy clay loam soil surfaces are widespread because of the occurrence of similar erosion processes like sheet wash, rill and gully erosion while in the central part of the Pendjari National Park loamy soil textures are prevailing. In fact, the steepness of the relief and the length of the slopes in the Pendjari Park seem to limit the development of some erosion forms as gully. Furthermore, the classification of soils reflects the variation of pedological processes along the transects and thus the occurrence of different soil types. The status of the sites may play an insignificant role in the differentiation of soil properties within the scale of small pedomorphological transects. A direct comparison of the vegetation type in the land use respectively partly protected and in the total protected sites (National Park of Pendjari) reveals a transition from the shrub savanna to the tree savanna. In conclusion it is important to insist on the fact that the variations of soil parameters within small slopes and the different sites are more conditioned by varying erosion processes and drainage conditions than the status protected or land use sites while the composition and diversity of plants is influenced by the status of the sites, the prevailing management tools, the pedogenetic conditions as well as the presence of wild animals like elephants. The ordination diagram shows that the organic matter is better correlated to the subgroup representing principally the sites of the hunting zone of the Pendjari Park and might be an explaining factor to the distribution of these sample sites groups. CEC ratios in the partly protected site of Natiabouani represent the highest measured in all sites. Nevertheless, statistical analysis of the CCA (canonical correspondence analysis) indicates generally a low correlation. This tendency is consolidated by the Monte Carlo test (p=0.14) which is a good indicator of species and environmental conditions. The detailed analysis of soil properties and the vegetation dynamic as well as their relationship within small pedomorphological transects represent an important pedological and botanical data collection involving different compartments. This thesis contributes to the better understanding of the savanna landscapes of West Africa and may provide essential scientific background for each development project directed towards interdisciplinary and integrative researches.
Nucleation of aerosol particles from trace atmospheric vapours is thought to provide up to half of global cloud condensation nuclei. Aerosols can cause a net cooling of climate by scattering sunlight and by leading to smaller but more numerous cloud droplets, which makes clouds brighter and extends their lifetimes. Atmospheric aerosols derived from human activities are thought to have compensated for a large fraction of the warming caused by greenhouse gases. However, despite its importance for climate, atmospheric nucleation is poorly understood. Recently, it has been shown that sulphuric acid and ammonia cannot explain particle formation rates observed in the lower atmosphere. It is thought that amines may enhance nucleation, but until now there has been no direct evidence for amine ternary nucleation under atmospheric conditions. Here we use the CLOUD (Cosmics Leaving Outdoor Droplets) chamber at CERN and find that dimethylamine above three parts per trillion by volume can enhance particle formation rates more than 1,000-fold compared with ammonia, sufficient to account for the particle formation rates observed in the atmosphere. Molecular analysis of the clusters reveals that the faster nucleation is explained by a base-stabilization mechanism involving acid–amine pairs, which strongly decrease evaporation. The ion-induced contribution is generally small, reflecting the high stability of sulphuric acid–dimethylamine clusters and indicating that galactic cosmic rays exert only a small influence on their formation, except at low overall formation rates. Our experimental measurements are well reproduced by a dynamical model based on quantum chemical calculations of binding energies of molecular clusters, without any fitted parameters. These results show that, in regions of the atmosphere near amine sources, both amines and sulphur dioxide should be considered when assessing the impact of anthropogenic activities on particle formation.
In dieser Arbeit wird die Richtungsabhängigkeit seismischer Geschwindigkeiten im Erdmantel unterhalb Deutschlands und angrenzender Gebiete durch die Analyse der teleseismischen Kernphase SKS auf Doppelbrechung untersucht (Scherwellen-Splitting). Die Anisotropie wird durch die Splittingparameter Φ und δt beschrieben und erlaubt Rückschlüsse auf geodynamische Prozesse.
Untersucht werden Aufzeichnungen des Deutschen Seismologischen Regionalnetzes (GRSN) und assoziierter Stationen aus dem Zeitraum von 1993 bis 2009. Für drei Stationen des Gräfenberg-Arrays (GRF-Array) sind Wellenformen ab 1976 verfügbar, welche damit einen weltweit einmaligen Datensatz liefern.
Auf Grund des stetigen Ausbaus der seismologischen Netze und des langen Beobachtungszeitraumes können über 3.000 Seismogramme ausgewertet werden. Der Hauptteil dieser Arbeit besteht daher in der Entwicklung einer automatischen Methodik zur Analyse von SKS-Splitting: ADORE ("Automatische Bestimmung von DOppelbrechnungsparametern in REgionalseismischen Netzwerken"). Für regionale Netze wie das GRSN gewährleistet ADORE eine objektive Bestimmung der Splittingparameter. Zunächst wird das seismologische Netzwerk als seismisches Array aufgefasst, um durch eine Frequenz-Wellenzahl-Analyse den Einsatz der SKS-Phase ohne manuellen Eingriff zu bestimmen. Die Berechnung der Splittingparameter erfolgt durch eine Inversion nach der Methode der Minimierung des transversalen Energieanteils. Automatisch wird das optimale Fenster um den SKS-Einsatz positioniert, für jede Beben-Stations-Kombination werden dazu 3.600 Einzelinversionen durchgeführt.
Um diese Vielzahl von Auswertungen in akzeptabler Zeit zu berechnen, nutzt ADORE moderne Rechnerarchitekturen aus, verteilt die Berechnungen auf mehrere Computer im lokalen Netzwerk und erzielt damit eine Beschleunigung um einen Faktor 60.
Die Analyse des gesamten Datensatzes ergibt folgende Ergebnisse: An allen analysierten Stationen wurde ein Scherwellen-Splitting festgestellt, der Stationsuntergrund weist somit überall Anisotropie auf. Für 240 Erdbeben können insgesamt 494 Wertepaare mit höchster Qualität bestimmt werden.
Unter der Annahme einer homogenen ungeneigten anisotropen Schicht unterhalb der jeweiligen Station können die Einzelmessungen pro Station gemittelt werden. Damit sind Regionen mit ähnlichen Merkmalen gut zu identifizieren: Im Norden Deutschlands herrschen NW-SO-, in der Mitte W-O-Richtungen und im Süden SW-NO-Richtungen vor.
Die Verzögerungszeiten liegen im Bereich zwischen 1.0 (Station Taunus) und 2.2 Sekunden (Tannenbergsthal, TANN). Auf Grund des hohen Wertes sind die Ursachen für die hier beobachteten Zeiten dem Erdmantel und nicht der Kruste zuzuordnen. Die bevorzugte Ausrichtung von anisotropen Kristallen auf Grund von Fließprozessen von Mantelmaterial ist Quelle der beobachteten Anisotropie. Rezente Fließprozesse von Mantelmaterial sind vor allem an der Unterkante der Lithosphäre wahrscheinlich. Durch Gebirgsbildungsprozesse, vorhandene Gebirgswurzeln oder regionale Veränderungen in der Mächtigkeit der Lithosphäre entstehen Barrieren für viskoses Mantelmaterial.
Als tektonische Ursachen für die hier gemessenen Orientierungen ist im Norden die Tornquist-Teisseyre-Linie (TTZ), in der Mitte die Variszische Gebirgsbildung und im Süden Einflüsse des Alpenbogens anzusehen. Ausnahmen bilden die Stationen Clausthal-Zellerfeld (CLZ), Rügen und Black-Forest-Observatory (BFO). Während bei letzterer ein Einfluss der Spreizungszone des Oberrheingrabens zu vermuten ist, scheint die Intrusion des Brockengranits die Beobachtungen an CLZ zu prägen. Rügen liegt in einer Übergangszone zwischen Sorgenfrei-Tornquist-Zone und TTZ.
Durch die Vielzahl von vorhandenen Einzelmessungen lassen sich an manchen Stationen komplexe Modelle untersuchen. Dazu zählen neben Gradientmodellen auch die geneigte Schicht und Zwei-Schicht-Modelle. Für sechs Stationen kann ein Zwei-Schicht-Modell erstellt werden: BFO, Gräfenberg A1, Fürstenfeldbruck (FUR), Rüdersdorf (RUE), TANN und Unterbreitzbach (UBBA). Die Interpretation der Richtungen von oberer und unterer Schicht gelingt für einen Teil der genannten Stationen: An BFO liegt die Orientierung der unteren Schicht parallel zur Vorzugsrichtung der variszischen Gebirgsbildung, jene der obere Schicht antiparallel zur Spreizungsrichtung des Rheingrabens. Für die Station FUR ist eine Überlagerung mit der Streichrichtung des Alpenmassivs zu beobachten. An GRA1 wird die untere Schicht offenbar durch rezente oder eingefrorene Anisotropie des Böhmischen Massivs bzw. des Eger-Riftsystems beeinflusst. Eine vergleichbare Wirkung ist durch die TTZ an der Station RUE zu erkennen.
ADORE wurde weiterhin auf einen Datensatz des temporären RIFTLINK-Projektes angewandt.
Sulphuric acid, ammonia, amines, and oxidised organics play a crucial role in nanoparticle formation in the atmosphere. In this study, we investigate the composition of nucleated nanoparticles formed from these compounds in the CLOUD (Cosmics Leaving Outdoor Droplets) chamber experiments at CERN (Centre européen pour la recherche nucléaire). The investigation was carried out via analysis of the particle hygroscopicity, ethanol affinity, oxidation state, and ion composition. Hygroscopicity was studied by a hygroscopic tandem differential mobility analyser and a cloud condensation nuclei counter, ethanol affinity by an organic differential mobility analyser and particle oxidation level by a high-resolution time-of-flight aerosol mass spectrometer. The ion composition was studied by an atmospheric pressure interface time-of-flight mass spectrometer. The volume fraction of the organics in the particles during their growth from sizes of a few nanometers to tens of nanometers was derived from measured hygroscopicity assuming the Zdanovskii–Stokes–Robinson relationship, and compared to values gained from the spectrometers. The ZSR-relationship was also applied to obtain the measured ethanol affinities during the particle growth, which were used to derive the volume fractions of sulphuric acid and the other inorganics (e.g. ammonium salts). In the presence of sulphuric acid and ammonia, particles with a mobility diameter of 150 nm were chemically neutralised to ammonium sulphate. In the presence of oxidation products of pinanediol, the organic volume fraction of freshly nucleated particles increased from 0.4 to ~0.9, with an increase in diameter from 2 to 63 nm. Conversely, the sulphuric acid volume fraction decreased from 0.6 to 0.1 when the particle diameter increased from 2 to 50 nm. The results provide information on the composition of nucleated aerosol particles during their growth in the presence of various combinations of sulphuric acid, ammonia, dimethylamine and organic oxidation products.
The international research project RECONCILE has addressed central questions regarding polar ozone depletion, with the objective to quantify some of the most relevant yet still uncertain physical and chemical processes and thereby improve prognostic modelling capabilities to realistically predict the response of the ozone layer to climate change. This overview paper outlines the scope and the general approach of RECONCILE, and it provides a summary of observations and modelling in 2010 and 2011 that have generated an in many respects unprecedented dataset to study processes in the Arctic winter stratosphere. Principally, it summarises important outcomes of RECONCILE including (i) better constraints and enhanced consistency on the set of parameters governing catalytic ozone destruction cycles, (ii) a better understanding of the role of cold binary aerosols in heterogeneous chlorine activation, (iii) an improved scheme of polar stratospheric cloud (PSC) processes that includes heterogeneous nucleation of nitric acid trihydrate (NAT) and ice on non-volatile background aerosol leading to better model parameterisations with respect to denitrification, and (iv) long transient simulations with a chemistry-climate model (CCM) updated based on the results of RECONCILE that better reproduce past ozone trends in Antarctica and are deemed to produce more reliable predictions of future ozone trends. The process studies and the global simulations conducted in RECONCILE show that in the Arctic, ozone depletion uncertainties in the chemical and microphysical processes are now clearly smaller than the sensitivity to dynamic variability.
The International Halocarbons in Air Comparison Experiment (IHALACE) was conducted to document relationships between calibration scales among various laboratories that measure atmospheric greenhouse and ozone depleting gases. Six stainless steel cylinders containing natural and modified natural air samples were circulated among 19 laboratories. Results from this experiment reveal relatively good agreement among commonly used calibration scales for a number of trace gases present in the unpolluted atmosphere at pmol mol−1 (parts per trillion) levels, such as chlorofluorocarbons (CFCs), hydrochlorofluorocarbons (HCFCs), and hydrofluorocarbons (HFCs). Some scale relationships were found to be consistent with those derived from bi-lateral experiments or from analysis of atmospheric data, while others revealed discrepancies. The transfer of calibration scales among laboratories was found to be problematic in many cases, meaning that measurements tied to a common scale may not, in fact, be compatible. These results reveal substantial improvements in calibration over previous comparisons. However there is room for improvement in communication and coordination of calibration activities with respect to the measurement of halogenated and related trace gases.
Numerous reactive volatile organic compounds (VOCs) are emitted into the atmosphere by vegetation. Most biogenic VOCs are highly reactive towards the atmosphere's most important oxidant, the hydroxyl (OH) radical. One way to investigate the chemical interplay between biosphere and atmosphere is through the measurement of total OH reactivity, the total loss rate of OH radicals. This study presents the first determination of total OH reactivity emission rates (measurements via the comparative reactivity method) based on a branch cuvette enclosure system mounted on a Norway spruce (Picea abies) throughout spring, summer and autumn 2011. In parallel VOC emission rates were monitored by a second proton-transfer-reaction mass spectrometer (PTR-MS), and total ozone (O3) loss rates were obtained inside the cuvette. Total OH reactivity emission rates were in general temperature and light dependent, showing strong diel cycles with highest values during daytime. Monoterpene emissions contributed most, accounting for 56–69% of the measured total OH reactivity flux in spring and early summer. However, during late summer and autumn the monoterpene contribution decreased to 11–16%. At this time, a large missing fraction of the total OH reactivity emission rate (70–84%) was found when compared to the VOC budget measured by PTR-MS. Total OH reactivity and missing total OH reactivity emission rates reached maximum values in late summer corresponding to the period of highest temperature. Total O3 loss rates within the closed cuvette showed similar diel profiles and comparable seasonality to the total OH reactivity fluxes.
Total OH reactivity fluxes were also compared to emissions from needle storage pools predicted by a temperature-only-dependent algorithm. Deviations of total OH reactivity fluxes from the temperature-only-dependent emission algorithm were observed for occasions of mechanical and heat stress. While for mechanical stress, induced by strong wind, measured VOCs could explain total OH reactivity emissions, during heat stress they could not. The temperature-driven algorithm matched the diel variation of total OH reactivity emission rates much better in spring than in summer, indicating a different production and emission scheme for summer and early autumn. During these times, unmeasured and possibly unknown primary biogenic emissions contributed significantly to the observed total OH reactivity flux.
In this study, two different methods were applied to derive daily and monthly sunshine duration based on high-resolution satellite products provided by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Climate Monitoring using data from Meteosat Second Generation (MSG) SEVIRI (Spinning Enhanced Visible and Infrared Imager). The satellite products were either hourly cloud type or hourly surface incoming direct radiation. The satellite sunshine duration estimates were not found to be significantly different using the native 15-minute temporal resolution of SEVIRI. The satellite-based sunshine duration products give additional spatial information over the European continent compared with equivalent in situ-based products. An evaluation of the satellite sunshine duration by product intercomparison and against station measurements was carried out to determine their accuracy. The satellite data were found to be within ±1 h/day compared to high-quality Baseline Surface Radiation Network or surface synoptic observations (SYNOP) station measurements. The satellite-based products differ more over the oceans than over land, mainly because of the treatment of fractional clouds in the cloud type-based sunshine duration product. This paper presents the methods used to derive the satellite sunshine duration products and the performance of the different retrievals. The main benefits and disadvantages compared to station-based products are also discussed.
Reduction of greenhouse gas (GHG) emissions to minimize climate change requires very significant societal effort. To motivate this effort, it is important to clarify the benefits of avoided emissions. To this end, we analysed the impact of four emissions scenarios on future renewable groundwater resources, which range from 1600 GtCO2 during the 21st century (RCP2.6) to 7300 GtCO2 (RCP8.5). Climate modelling uncertainty was taken into account by applying the bias-corrected output of a small ensemble of five CMIP5 global climate models (GCM) as provided by the ISI-MIP effort to the global hydrological model WaterGAP. Despite significant climate model uncertainty, the benefits of avoided emissions with respect to renewable groundwater resources (i.e. groundwater recharge (GWR)) are obvious. The percentage of projected global population (SSP2 population scenario) suffering from a significant decrease of GWR of more than 10% by the 2080s as compared to 1971–2000 decreases from 38% (GCM range 27–50%) for RCP8.5 to 24% (11–39%) for RCP2.6. The population fraction that is spared from any significant GWR change would increase from 29% to 47% if emissions were restricted to RCP2.6. Increases of GWR are more likely to occur in areas with below average population density, while GWR decreases of more than 30% affect especially (semi)arid regions, across all GCMs. Considering change of renewable groundwater resources as a function of mean global temperature (GMT) rise, the land area that is affected by GWR decreases of more than 30% and 70% increases linearly with global warming from 0 to 3 ° C. For each degree of GMT rise, an additional 4% of the global land area (except Greenland and Antarctica) is affected by a GWR decrease of more than 30%, and an additional 1% is affected by a decrease of more than 70%.
Environmental change impacts on the C- and N-cycle of European forests: a model comparison study
(2013)
Forests are important components of the greenhouse gas balance of Europe. There is considerable uncertainty about how predicted changes to climate and nitrogen deposition will perturb the carbon and nitrogen cycles of European forests and thereby alter forest growth, carbon sequestration and N2O emission. The present study aimed to quantify the carbon and nitrogen balance, including the exchange of greenhouse gases, of European forests over the period 2010–2030, with a particular emphasis on the spatial variability of change. The analysis was carried out for two tree species: European beech and Scots pine. For this purpose, four different dynamic models were used: BASFOR, DailyDayCent, INTEGRATOR and Landscape-DNDC. These models span a range from semi-empirical to complex mechanistic. Comparison of these models allowed assessment of the extent to which model predictions depended on differences in model inputs and structure. We found a European average carbon sink of 0.160 ± 0.020 kgC m−2 yr−1 (pine) and 0.138 ± 0.062 kgC m−2 yr−1 (beech) and N2O source of 0.285 ± 0.125 kgN ha−1 yr−1 (pine) and 0.575 ± 0.105 kgN ha−1 yr−1 (beech). The European average greenhouse gas potential of the carbon sink was 18 (pine) and 8 (beech) times that of the N2O source. Carbon sequestration was larger in the trees than in the soil. Carbon sequestration and forest growth were largest in central Europe and lowest in northern Sweden and Finland, N. Poland and S. Spain. No single driver was found to dominate change across Europe. Forests were found to be most sensitive to change in environmental drivers where the drivers were limiting growth, where changes were particularly large or where changes acted in concert. The models disagreed as to which environmental changes were most significant for the geographical variation in forest growth and as to which tree species showed the largest rate of carbon sequestration. Pine and beech forests were found to have differing sensitivities to environmental change, in particular the response to changes in nitrogen and precipitation, with beech forest more vulnerable to drought. There was considerable uncertainty about the geographical location of N2O emissions. Two of the models BASFOR and LandscapeDNDC had largest emissions in central Europe where nitrogen deposition and soil nitrogen were largest, whereas the two other models identified different regions with large N2O emission. N2O emissions were found to be larger from beech than pine forests and were found to be particularly sensitive to forest growth.
When studying new particle formation, the uncertainty in determining the "true" nucleation rate is considerably reduced when using Condensation Particle Counters (CPCs) capable of measuring concentrations of aerosol particles at sizes close to or even at the critical cluster size (1–2 nm). Recently CPCs, able to reliably detect particles below 2 nm in size and even close to 1 nm became available. The corrections needed to calculate nucleation rates are substantially reduced compared to scaling the observed formation rate to the nucleation rate at the critical cluster size. However, this improved instrumentation requires a careful characterization of their cut-off size and the shape of the detection efficiency curve because relatively small shifts in the cut-off size can translate into larger relative errors when measuring particles close to the cut-off size.
Here we describe the development of two continuous flow CPCs using diethylene glycol (DEG) as the working fluid. The design is based on two TSI 3776 counters. Several sets of measurements to characterize their performance at different temperature settings were carried out. Furthermore two mixing-type Particle Size Magnifiers (PSM) A09 from Airmodus were characterized in parallel. One PSM was operated at the highest mixing ratio (1 L min−1 saturator flow), and the other was operated in a scanning mode, where the mixing ratios are changed periodically, resulting in a range of cut-off sizes. Different test aerosols were generated using a nano-Differential Mobility Analyzer (nano-DMA) or a high resolution DMA, to obtain detection efficiency curves for all four CPCs. One calibration setup included a high resolution mass spectrometer (APi-TOF) for the determination of the chemical composition of the generated clusters. The lowest cut-off sizes were achieved with negatively charged ammonium sulphate clusters, resulting in cut-offs of 1.4 nm for the laminar flow CPCs and 1.2 and 1.1 nm for the PSMs. A comparison of one of the laminar-flow CPCs and one of the PSMs measuring ambient and laboratory air showed good agreement between the instruments.
The SAWNUC microphysical aerosol nucleation model is used to study the effect of reactor walls on the interpretation of nucleation experiments with respect to nucleation theory. This work shows that loss processes, such as wall losses, influence the interpretation of nucleation experiments, especially at low growth rates and short lifetime of freshly nucleated particles. In these cases the power dependency of the formation rates, determined at a certain particle size, with respect to H2SO4 does not correspond to the approximate number of H2SO4 molecules in the critical cluster as expected by the first nucleation theorem. Observed ∂log(J)/∂log([H2SO4]) therefore can vary widely for identical nucleation conditions but different sink terms.
We investigate the contribution of oceanic methyl iodide (CH3I) to the stratospheric iodine budget. Based on CH3I measurements during three tropical ship campaigns and the Lagrangian transport model FLEXPART we provide a detailed analysis of CH3I transport from the ocean surface to the cold point in the upper tropical tropopause layer (TTL). While average oceanic emissions differ by less than 50% from campaign to campaign, the measurements show much stronger variations within each campaign. A positive correlation between the oceanic CH3I emissions and the efficiency of CH3I troposphere–stratosphere transport has been identified for some cruise sections. The mechanism of strong horizontal surface winds triggering large emissions on the one hand and being associated with tropical convective systems, such as developing typhoons, on the other hand, could explain the identified correlations. As a result of the simultaneous occurrence of large CH3I emissions and strong vertical uplift, localized maximum mixing ratios of 0.6 ppt CH3I at the cold point have been determined for observed peak emissions during the SHIVA-Sonne campaign in the coastal West Pacific. The other two campaigns give considerable smaller maxima of 0.1 ppt CH3I for the TransBrom campaign in the open West Pacific and 0.03 ppt for emissions from the coastal East Atlantic during the DRIVE campaign. In order to assess the representativeness of the large local mixing ratios we use climatological emission scenarios to derive global upper air estimates of CH3I abundances. The model results are compared to available upper air measurements including data from the recent ATTREX and HIPPO2 aircraft campaigns. In the East Pacific region, the location of the available measurement campaigns in the upper TTL, the comparisons give a good agreement indicating that around 0.01 to 0.02 ppt of CH3I enter the stratosphere. However, other tropical regions, which are subject to stronger convective activity show larger CH3I entrainment, e.g., 0.08 ppt in the West Pacific. The strong variations in the geographical distribution of CH3I entrainment suggest that currently available upper air measurements are not representative of global estimates and further campaigns will be necessary in order to better understand the CH3I contribution to stratospheric iodine.
Water footprints have been proposed as sustainability indicators, relating the consumption of goods like food to the amount of water necessary for their production and the impacts of that water use in the source regions. We have further developed the existing water footprint methodology by globally resolving virtual water flows and import and source regions at 5 arc minutes spatial resolution, and by assessing local impacts of export production. Applying this method to three exemplary cities, Berlin, Delhi and Lagos, we find major differences in amounts, composition, and origin of green and blue virtual water imports, due to differences in diets, trade integration and crop water productivities in the source regions. While almost all of Delhi's and Lagos' virtual water imports are of domestic origin, Berlin on average imports from more than 4000 km distance, in particular soy (livestock feed), coffee and cocoa. While 42% of Delhi's virtual water imports are blue water based, the fractions for Berlin and Lagos are 2% and 0.5%, respectively, roughly equal to local drinking water abstractions of these cities. Some of the external source regions of Berlin's virtual water imports appear to be critically water scarce and/or food insecure. However for deriving recommendations on sustainable consumption and trade, further analysis of context-specific costs and benefits associated with export production will be required.
Estimates of the recovery time of stratospheric ozone heavily rely on the exact knowledge of the processes that lead to the decomposition of the relevant halogenated source gases. Crucial parameters in this context are fractional release factors (FRFs) as well as stratospheric lifetimes and ozone depletion potentials (ODPs). We here present data from the analysis of air samples collected between 2009 and 2011 on board research aircraft flying in the mid- and high-latitude stratosphere and infer the above-mentioned parameters for ten major source gases: CFCl3 (CFC-11), CF2Cl2 (CFC-12), CF2ClCFCl2 (CFC-113), CCl4 (carbon tetrachloride), CH3CCl3 (methyl chloroform), CHF2Cl (HCFC-22), CH3CFCl2 (HCFC-141b), CH3CF2Cl (HCFC-142b), CF2ClBr (H-1211), and CF3Br (H-1301). The inferred correlations of their FRFs with mean ages of air reveal less decomposition as compared to previous studies for most compounds. When using the calculated set of FRFs to infer equivalent stratospheric chlorine, we find a reduction of more than 20% as compared to the values inferred in the most recent Scientific Assessment of Ozone Depletion by the World Meteorological Organisation (WMO, 2011). We also note that FRFs and their correlations with mean age are not generally time-independent as often assumed. The stratospheric lifetimes were calculated relative to that of CFC-11. Within our uncertainties the ratios between stratospheric lifetimes inferred here agree with the values in recent WMO reports except for CFC-11, CFC-12 and CH3CCl3. Finally, we calculate lower ODPs than recommended by WMO for six out of ten compounds, with changes most pronounced for the three HCFCs. Collectively these newly calculated values may have important implications for the severity and recovery time of stratospheric ozone loss.
Whereas fossil evidence indicates extensive treeless vegetation and diverse grazing megafauna in Europe and northern Asia during the last glacial, experiments combining vegetation models and climate models have to-date simulated widespread persistence of trees. Resolving this conflict is key to understanding both last glacial ecosystems and extinction of most of the mega-herbivores. Using a dynamic vegetation model (DVM) we explored the implications of the differing climatic conditions generated by a general circulation model (GCM) in “normal” and “hosing” experiments. Whilst the former approximate interstadial conditions, the latter, designed to mimic Heinrich Events, approximate stadial conditions. The “hosing” experiments gave simulated European vegetation much closer in composition to that inferred from fossil evidence than did the “normal” experiments. Given the short duration of interstadials, and the rate at which forest cover expanded during the late-glacial and early Holocene, our results demonstrate the importance of millennial variability in determining the character of last glacial ecosystems.
Since the 1980s, advances in wastewater treatment technology have led to considerably improved surface water quality in the urban areas of many high income countries. However, trace concentrations of organic wastewater-associated contaminants may still pose a key environmental hazard impairing the ecological quality of surface waters. To identify key impact factors, we analyzed the effects of a wide range of anthropogenic and environmental variables on the aquatic macroinvertebrate community. We assessed ecological water quality at 26 sampling sites in four urban German lowland river systems with a 0–100% load of state-of-the-art biological activated sludge treated wastewater. The chemical analysis suite comprised 12 organic contaminants (five phosphor organic flame retardants, two musk fragrances, bisphenol A, nonylphenol, octylphenol, diethyltoluamide, terbutryn), 16 polycyclic aromatic hydrocarbons, and 12 heavy metals. Non-metric multidimensional scaling identified organic contaminants that are mainly wastewater-associated (i.e., phosphor organic flame retardants, musk fragrances, and diethyltoluamide) as a major impact variable on macroinvertebrate species composition. The structural degradation of streams was also identified as a significant factor. Multiple linear regression models revealed a significant impact of organic contaminants on invertebrate populations, in particular on Ephemeroptera, Plecoptera, and Trichoptera species. Spearman rank correlation analyses confirmed wastewater-associated organic contaminants as the most significant variable negatively impacting the biodiversity of sensitive macroinvertebrate species. In addition to increased aquatic pollution with organic contaminants, a greater wastewater fraction was accompanied by a slight decrease in oxygen concentration and an increase in salinity. This study highlights the importance of reducing the wastewater-associated impact on surface waters. For aquatic ecosystems in urban areas this would lead to: (i) improvement of the ecological integrity, (ii) reduction of biodiversity loss, and (iii) faster achievement of objectives of legislative requirements, e.g., the European Water Framework Directive.
This research was conducted in the Rwenzori Region of the Western Branch, East African Rift System (EARS). The EARS is a tectonic structure extending over a length of more than 3000 km from the Afar Triple Junction, in Ethiopia, to Lake Malawi in the south. The Western Rift System is a roughly NE to ENE trending sector of the EARS, which runs along the western boundary of Uganda and the neighboring Democratic Republic of Congo (D.R.C). It stretches 2100 km from Nimule, NW on Uganda-Sudan border, extending to Lake Malawi in the SE of Africa. The unusual uplift of the Rwenzori Mountains within an extensional regime and the mechanisms associated with the high frequency of seismic activity in the region was hardly understood and therefore, had remained a subject of contention that needed to be critically addressed in detail. To my knowledge, this was probably the first study to be performed and documented in great depth within the domains of seismic noise variation, seismic anisotropy and b value analyses beneath the Rwenzori Region. After about six years of operation (2006-2012), the seismology group of the RIFTLINK Research Project (www.riftlink.org) acquired a vast amount of high-quality, digital data that were collected using a seismic network of well calibrated seismic equipment. The project was divided into two phases. Phase I, that operated between February 2006 - September 2007, consisted of thirty-two temporary seismic stations, which were selectively spread out in the Rwenzori Region on the Ugandan side, to detect and record extremely weak as well as strong naturally occurring earthquakes. The seismic equipment used included EDL and REFTEK digitizers, which were coupled with Güralp and MARK sensors respectively (REFTEKS: only short-period MARK sensors, EDLs: short-period MARK plus few broadband Güralp Sensors). Exactly 22375 earthquakes were recorded. The data were processed using the SEISAN software package. About 14413 earthquakes were carefully localized using the velocity model of Bram (1975) that implements a Vp=Vs ratio fixed at 1.74. Phase II, that extended between 2009-2012 consisted of thirty-two seismic stations, which were spread out around the Rwenzori Mountains, both on the Ugandan side and the neighboring D.R.C. Only Taurus digitizers that were coupled with Trillium sensors were used in the D.R.C. On the Ugandan side however, both EDL and Taurus digitizers, which were coupled with Trillium and Güralp sensors were used. ...
Projections of future changes in runoff can have important implications for water resources and flooding. In this study, runoff projections from ISI-MIP (Inter-sectoral Impact Model Intercomparison Project) simulations forced with HadGEM2-ES bias-corrected climate data under the Representative Concentration Pathway 8.5 have been analysed. Projections of change from the baseline period (1981–2010) to the future (2070–2099) from a number of different ecosystems and hydrological models were studied. The differences between projections from the two types of model were looked at globally and regionally. Typically, across different regions the ecosystem models tended to project larger increases and smaller decreases in runoff than the hydrological models. However, the differences varied both regionally and seasonally. Sensitivity experiments were also used to investigate the contributions of varying CO2 and allowing vegetation distribution to evolve on projected changes in runoff. In two out of four models which had data available from CO2 sensitivity experiments, allowing CO2 to vary was found to increase runoff more than keeping CO2 constant, while in two models runoff decreased. This suggests more uncertainty in runoff responses to elevated CO2 than previously considered. As CO2 effects on evapotranspiration via stomatal conductance and leaf-area index are more commonly included in ecosystems models than in hydrological models, this may partially explain some of the difference between model types. Keeping the vegetation distribution static in JULES runs had much less effect on runoff projections than varying CO2, but this may be more pronounced if looked at over a longer timescale as vegetation changes may take longer to reach a new state.