Refine
Year of publication
Document Type
- Doctoral Thesis (139) (remove)
Has Fulltext
- yes (139)
Is part of the Bibliography
- no (139)
Keywords
- Stratosphäre (4)
- Eiskeime (3)
- Meteorologie (3)
- Monsun (3)
- Schwerewellen (3)
- Anisotropie (2)
- BTEX (2)
- Belize (2)
- Boden (2)
- Caribbean (2)
Institute
- Geowissenschaften (139) (remove)
During this study clumped isotope analysis of carbonates was established at the Goethe University of Frankfurt, Germany. Therefore, preparation protocols and analytical parameters were elaborated to obtain precise and accurate Δ47 data. Briefly, analyte CO2 was cleaned cryogenically using glass extraction lines to remove traces of water that enable re-equilibration of C–O bonds in the gases. Furthermore, analyte CO2 was passed through a gas chromatograph (GC) to clean it from contaminants that produce isobaric interferences with m/z 47. Initially, phosphoric acid digestions of carbonates was conducted at 25 °C in McCrea-type reaction vessels. Afterwards samples were reacted at 90 °C using a common acid bath. Mass spectrometric analyses were performed using a MAT 253 equipped with a dual inlet system. Δ47 values were directly projected to the absolute scale using CO2 gases equilibrated at distinct temperatures.
In cooperation with Stefano Bernasconi and his research group at ETH Zurich we studied the non-linearity that occurs for the measurement of m/z 47. This effect results from secondary electrons created by the m/z 44 beam. These electrons cause a negative background on the m/z 47 collector. A correction procedure was proposed that relies on the determination of the negative background on the m/z 47 Faraday cup. This approach might reduce time-consuming analyses of heated gases which were used so far to account for the observed non-linearity. However, the suggested correction of the negative background on the m/z 47 cup is only applicable if the slit width of the m/z 44 beam is significantly wider than that of the m/z 47 beam.
This thesis, furthermore, presents a comparison of the different phosphoric acid digestion techniques which are commonly used for carbonate clumped isotope analysis. For calcitic and aragonitic material digested at 25 °C in McCrea-type vessels we observed that the sample size has an effect on Δ47 data: higher mean Δ47 values and a larger scatter of data were received for samples <7 mg than for larger aliquots. For carbonate samples digested at 90 °C in a common acid bath no sample size effect was determined. We assume that secondary re-equilibration of CO2 with water preferentially occurs at 25 °C producing the observed differences. However, a sample size effect can be avoided if reaction temperature is increased to 90 °C.
In order to make carbonate Δ47 data obtained from acid digestions at 90 °C comparable to Δ47 data received from reactions at 25 °C the difference of the acid fractionation factores (Δ47*25-90) between both temperatures has to be known. For the determination of the Δ47*25-90 value we have considered Δ47 data made at 25 °C from samples >7 mg only. For calicte and aragonite we obtained differences in fractionation factores of 0.075‰ and 0.066‰, respectively. These Δ47*25-90 values are coincident with the theoretical prediction of 0.069‰ proposed for calcite (Guo et al., 2009).
Moreover, this dissertation comprises a calibration study of the clumped isotope thermometer based on various natural calcites that grew between 9 and 38 °C. The samples include a brachiopod shell, a bivalve shell, an eggshell of an ostrich and foraminifera tests which formed from distinct biomineralizing processes. Furthermore we included an authigenic carbonate crystallized from biological-induced precipitation. The following linear relationship between 1/T2 and Δ47 was determined (with Δ47 in ‰ and T in K):
Δ47 = 0.0327 (± 0.0026) x 106 / T2 + 0.3030 (± 0.0308) (R2 = 0.9915)
This equation differs from the pioneering Ghosh et al. (2006a) calibration. However, our regression line is statistically indistinguishable from that of Henkes et al. (2013) which is based on aragonitic mollusks and calcitic brachiopod shells. Both studies have in common that calibration data were, at first, directly referenced to the absolute scale. In addition, both datasets rely on similar digestion techniques. Furthermore, the two calibrations are conform with the theoretical prediction of Guo et al. (2009).
The calcite calibration of the clumped isotope paleothermometer received in this study was applied to Δ47 data measured for Silurian brachiopods shells from Gotland/Sweden. Prior to isotopic analysis the fossils were intensively investigated for their preservation state (CL, SEM, trace elements). The lowest T(Δ47) values of ca. 28 to 33 °C were estimated from ultrastructurally well-preserved regions of some shells. For these samples also the lowest δ18Ow values of Silurian seawater were determined. These estimates of ca. −1‰ confirm the assumption that the δ18O value of the Silurian ocean was buffered to (0 ± 1)‰.
Nevertheless, most studied shells were characterized by a patchwork of pristine and altered shell portions resulting in elevated T(Δ47) values which plot mostly between 40 and 60 °C. Our results indicate that the clumped isotopic composition of the shells were altered at low water-rock ratios, not affecting the δ18O values. Δ47 and δ18O data of associated diagenetic phases (sparitic and micritic phases of the inner fillings of the fossils) provide evidence that the sparitic cements grew during several diagenetic events which occurred at different temperatures in fluid-buffered systems. We, furthermore, conclude that the micritic phases lithified at a very early diagenetic stage with the δ18O values being most probably close to a Silurian seawater composition
Terrestrial climate and ecosystem evolution during ‘Greenhouse Earth’ phases of the early Paleogene remain incompletely known. Particularly, paleobotanical records from high southern latitudes are giving only limited insights into the Paleocene and early Eocene vegetation of the region. Hence, data from continuous well-calibrated sequences are required to make progress with the reconstruction of terrestrial climate and ecosystem dynamics from the southern latitudes during the early Paleogene.
In order to elucidate the terrestrial conditions from the high southern latitudes during the early Paleogene, terrestrial palynology was applied in the present study to two well-dated deep-marine sediment cores located at the Australo-Antarctic region: (i) IODP Site U1356 (Wilkes Land margin, East Antarctica) and (ii) ODP Site 1172 (East Tasman Plateau, southwest Pacific Ocean). The studied sequence from IODP Site U1356 comprises mid-shelfal sediments from the early to middle Eocene (53.9 – 46 million years ago [Ma]). For the ODP Site 1172, the studied succession is characterized by sediments deposited in shallow marine environments of the middle Paleocene to the early Eocene (60.7 – 54.2 Ma).
Based on the obtained pollen and spores (sporomorphs) results from the studied sequences of Site U1356 and Site 1172, this study aims to: (1) decipher the terrestrial climate conditions along the Australo-Antarctic region from the middle Paleocene to the middle Eocene; (2) evaluate the structure, diversity and compositional patterns of forests that throve in the Australo-Antarctic region during the early Paleogene; (3) understand the response of forests from the high southern latitudes to the climate dynamics from the early Paleogene; (4) establish a connection between the generated terrestrial palynomorph data and published Sea Surface Temperatures (SSTs) from the same cores.
To decipher the terrestrial climatic conditions on the Australo-Antarctic region, this study relies on the nearest living relative (NLR) concept that assumes that fossil taxa have similar climate requirements as their modern counterparts. This approach was applied to the sporomorph results of Site U1356 and Site 1172, following mainly the bioclimatic analysis. With regard to the structure and diversity patterns of the vegetation from the same region, the present study presents combined qualitative (i.e., reconstruction of the vegetation based mainly on the habitats of the known living relatives) and quantitative (i.e., application of ordination techniques, rarefaction and diversity indices) analyses of the fossil sporomorphs results.
The overall results from the paleoclimatic and vegetation reconstruction approaches applied in the present study, indicate that temperate and paratropical forests during the early Paleogene throve under different climatic conditions on the Wilkes Land margin and on Tasmania, at paleolatitudes of ∼70°S and ∼65°S, respectively.
Specifically, the sporomorph results from Site U1356, suggest that a highly diverse forest similar to present-day forests from New Caledonia was thriving on Antarctica during the early Eocene (53.9 – 51.9 Ma). These forests were characterized by the presence of termophilous taxa that are restricted today to tropical and subtropical settings, notably Bombacoideae, Strasburgeria, Beauprea, Spathiphyllum, Anacolosa and Lygodium. In combination with MBT/CBT paleotemperature results, they provide strong evidence for near-tropical warmth at least in the coastal lowlands along the Wilkes Land margin. The coeval presence of frost tolerant taxa such as Nothofagus, Araucariaceae and Podocarpaceae during the early Eocene on the same record suggests that paratropical forests were thriving along the Wilkes Land margin. Due to the presence of this kind of vegetation, it is possible to suggest that forests in this region were subject to a climatic gradient related to differences in elevation and/or the proximity to the coastline.
By the middle Eocene, the paratropical forests that characterized the vegetation of the early Eocene on the Wilkes Land margin were replaced by low diversity temperate forests dominated by Nothofagus, and similar to present-day cool-temperate forests from New Zealand. The dominance of these forests and the absence of thermophilous elements together with the lower temperatures suggested by the MBT/CBT and the sporomorph-based temperatures indicate consistently cooler conditions during this time interval.
With regard to the sporomorph results of Site 1172, this study suggests that three vegetation types were thriving on Tasmania from the middle Paleocene to the early Eocene under different climatic conditions. During the middle to late Paleocene, warm-temperate forests dominated by Podocarpaceae and Araucariaceae were the prevailing vegetation on Tasmania. The dominance of these forests was interrupted by the transient predominance of cool-temperate forests dominated by Nothofagus and Araucariaceae across the middle/late Paleocene transition interval (~59.5 to ~59.0 Ma). This cool-temperate forest was characterized by a lack of frost-sensitive elements (i.e., palms and cycads) indicating cooler conditions with harsher winters on Tasmania during this time interval. By the early Eocene, and linked with the Paleocene Eocene Thermal Maximum (PETM), Paleocene temperate forests dominated by gymnosperms were replaced by paratropical rainforests with the remarkable presence of the tropical mangrove palm Nypa during the PETM and the earliest Eocene. The overall results from Site U1356 and Site 1172, provide a new assessment of the terrestrial climatic conditions in the Australo-Antarctic region for validating climate models and understanding the response of high-latitude terrestrial ecosystems to the climate dynamics of the early Paleogene on southern latitudes.
The climatic conditions in the higher latitudes during the early Paleogene were further unravelled by comparing the obtained terrestrial and marine results. The integration of the obtained sporomorph data with previously published TEX86-based SSTs from Site 1172 documents that the vegetation dynamics were closely linked with the temperature evolution from the Australo-Antarctic region. Moreover, the comparison of TEX86-based SSTs and sporomorph-based climatic estimations from Site 1172 suggests a warm-season bias of both calibrations of TEX86 (i.e., TEX86Hand TEX86H), when this proxy is applied to high southern latitudes records of the early Paleogene.
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.
Die Seismizität des nördlichen Oberrheingrabens (ORG) ist aufgrund seines Potentials für die geothermische Nutzung und der damit möglicherweise verbundenen seismischen Risiken von allgemeinem Interesse. Detaillierte Kenntnisse der natürlichen Seismizität erlauben Rückschlüsse auf aktive Störungszonen und Spannungsverhältnisse im Untergrund. Sie liefert außerdem wichtige Hintergrundinformationen für die Abschätzung einer möglichen induzierten Seismizität. Untersuchungen zur Charakterisierung der natürlichen Seismizität, des Spannungsfeldes und der seismischen Gefährdung des nördlichen ORG sind Hauptbestandteil dieser Arbeit, die innerhalb des BMU/BMWi-Projektes SiMoN (Seismisches Monitoring im Zusammenhang mit der geothermischen Nutzung des Nördlichen Oberrheingrabens) entstanden ist. Aufzeichnungen eines Netzwerkes aus 13 seismischen Stationen dienen als Datengrundlage zur Charakterisierung der Seismizität innerhalb eines etwa 50 x 60 km2 großen Areals im dichtbesiedelten Rhein-Main Gebiet. Untersuchungen der Rauschbedingungen zur Bewertung der Eignung der Stationsorte für das Aufzeichnen der natürlichen Seismizität lieferten bei den Stationen auf felsigem Untergrund sehr gute spektrale Eigenschaften, während alle Stationen im Sediment des ORG deutlich höhere Rauschanteile aufzeigten. Anhand systematischer Messungen in flachen Bohrlöchern konnten laterale und vertikale Variationen des seismischen Rauschens beschrieben werden und dadurch eine Verbesserung der Detektionsschwelle beobachtet werden.
Es werden die Ergebnisse des seismischen Monitorings für den Zeitraum November 2010 bis Dezember 2014 dargestellt. Die Detektionsschwelle für das Netzwerk liegt bei einer Lokalmagnitude von etwa 0,5, die Vollständigkeitsmagnitude beträgt Mc = 1,2. Seit Beginn der Datenaufzeichnung konnten 243 Erdbeben im unmittelbaren Bereich des Stations-netzwerkes mit Magnituden im Bereich zwischen ML = -0,5 und ML = 4,2 lokalisiert werden. Die Epizentren liegen hauptsächlich entlang der östlichen Grabenschulter und im Graben; entlang der westlichen Grabenschulter ist die seismische Aktivität deutlich geringer. Eine weitere aktive Region konnte entlang der südlichen Ausläufer des Taunus im Nordwesten des Untersuchungsgebietes identifiziert werden. Die Seismizität erstreckt sich bis in eine Tiefe von 24 km mit einem Maximum der hypozentralen Tiefenverteilung im Bereich von 12-18 km. Im Graben ist die Seismizität dabei auf die tiefere Kruste im Bereich von 9-24 km beschränkt. Das Fehlen von seismischer Aktivität in der oberen Kruste bis ca. 9 km Tiefe im Graben könnte auf eine aseismische Deformation in diesem Tiefenbereich hindeuten. Seit Mai 2014 konnte südöstlich von Darmstadt bei der Ortschaft Ober-Ramstadt zum ersten Mal seit fast 150 Jahren eine Schwarmbebenaktivität im Bereich des nördlichen ORG registriert werden. Die Hypozentren sind in zwei Cluster unterteilt, die räumlich voneinander getrennt sind und unterschiedliche Aktivitätsraten aufweisen. Die Herdtiefen liegen im Bereich von 1-8 km.
Zusätzlich zu den Daten des SiMoN Netzwerkes wurden Aufzeichnungen der regionalen Erdbebendienste in Herdflächenanalysen für insgesamt 58 Erdbeben einbezogen. Die Herdflächenlösungen weisen überwiegend Blattverschiebungen (Strike-slip-Mechanismen) auf. Auf- und Abschiebungen spielen nur eine untergeordnete Rolle. Die berechneten Herdmechanismen bestätigen, dass sich das Spannungsfeld des nördlichen ORG transtensional verhält, im Vergleich zu früheren Studien konnte jedoch eine deutlich ausgeprägte Blattverschiebungskomponente identifiziert werden. Zur Bestimmung der Hauptspannungsachsen wurde eine Inversion der Herdflächenlösungen durchgeführt und die Richtung der maximalen horizontalen Spannung, welche hauptsächlich in N135°E orientiert ist, bestimmt.
Aufbauend auf den neu gewonnen Erkenntnissen zur natürlichen Seismizität und zum Spannungsfeld des nördlichen ORG wurde eine probabilistische seismische Gefährdungsanalyse durchgeführt. Um Unsicherheiten in den seismischen Quellregion-modellen zu berücksichtigen, wurden sechs unterschiedliche Modelle entwickelt. Für jede Quellregion wurden spezifische Parameter bestimmt. Ihre Unsicherheiten werden in einem logischen Baum behandelt. Auf der Grundlage eines neu zusammengestellten Momentmagnituden-basierten Erdbebenkatalogs wurden die Magnitudenhäufigkeits-parameter bestimmt. Unter Berücksichtigung des tektonischen Regimes in jeder Quelle wurden unterschiedliche Dämpfungsrelationen der Bodenbeschleunigung verwendet. Zur Quantifizierung der maximal zu erwartenden Magnitude in jeder Quelle wurden Wahrscheinlichkeitsdichtefunktionen berechnet. Die Resultate der Gefährdungsanalyse werden in Form von Karten der Spektralbodenbeschleunigungen und Spitzenboden-beschleunigungen für Wiederkehrperioden von 475 und 2475 Jahren und Antwortbeschleunigungsspektren dargestellt. Im Vergleich zu früheren Studien konnte eine erhöhte seismische Gefährdung für den nördlichen ORG festgestellt werden.
This study describes the Holocene sedimentary lagoonal deposition history, including event sedimentation and benthic foraminiferal analyzes, from about 10 kyrs BP until today. This is the first study describing the sedimentation of a Maldivian atoll lagoon in such detail. Thirty-nine sediment cores have been recovered from the deep Rasdhoo Atoll lagoon of the Maldives (4°N/73°W). Seventeen sediment cores were opened, described, and 296 sediment samples have been collected and analyzed. Different methods have been used to evaluate the coarse- and fine-grained carbonate components and a total of fifty-eight samples have been dated radiometrically by Beta Analytic Inc., Miami, Florida. In general, the Rasdhoo Atoll lagoon sediments can be divided into (1) a Late Pleistocene soil, (2) an early Holocene peat layer composed of mangrove deposits which mark the beginning inundation of the atoll lagoon by the rising Holocene sea-level at 10,320 ± 100 yrs BP, and (3) carbonate sediments starting to fill up the lagoon 7850 ± 140 yrs BP until today. The transition from peat to carbonate is characterized by a considerable hiatus. Six different carbonate sediment facies are classified by statistical analyses, listed in decreasing abundance:
(1) mollusk-coral-algal floatstone to rudstone (30%)
(2) mollusk-coral-red algae rudstone (23%)
(3) mollusk-coral-algal wackestone to floatstone (23%)
(4) mollusk-coral wackestone (13%)
(5) mollusk-coral mudstone to wackestone (9%)
(6) mollusk mudstone (2%)
Based on grain-sizes in combination with coral identification, the facies represent both lagoonal background sedimentation (mostly fine-grained sediments (matrix >50%)) and event sedimentation (coarse-grained sediment layers composing reefal components).
Six coarser grained layers in muddy background sediments of the Rasdhoo Atoll lagoon were interpreted as Holocene tsunami events, based on the increase of allochthonous skeletal material with shallow-water reef affinity such as fragments of shallow-water coral species, coralline red algae, and reef-dwelling foraminifera in these layers, as well as AMS dating:
• Event 1: 420 - 890 yrs BP (655 yrs BP)
• Event 2: 890 - 1560 yrs BP (1225 yrs BP)
• Event 3: 2040 - 2340 yrs BP (2190 yrs BP)
• Event 4: 2420 - 3380 yrs BP (2900 yrs BP)
• Event 5: 3890 - 4330 yrs BP (4110 yrs BP)
• Event 6: 5480 - 5760 yrs BP (5620 yrs BP)
Five of the six layers may be correlated to previously published tsunami events at adjacent coastal research sites. The mid-late Holocene atoll lagoon archive is incomplete though based on the assumption that major earthquakes at the Indonesian subduction zone generated more than six major tsunamis during the past 6.5 kyrs.
According to Gischler (2006), the sediments of the Rasdhoo Atoll lagoon can be divided into two areas: (1) a central to marginal deep lagoon with a lateral west-to-east gradient of sediment facies distribution, visible in sections <4 kyrs BP with sedimentary facies of mudstone to wackestone in the western part (e.g., cores 16, 18, and 34) and coarse-grained coral and algal-rich sediments in the eastern part of the lagoon (e.g., cores 30 and 31). (2) A northern enclosed and shallow area between the sand apron and the sand spit accumulating “sandy” sediments of wackestone facies (cores 2, 19, 25, and 26).
Comparing the sediment accumulation data of the lagoon with two reconstructed local sea-level curves, three different sequence-stratigraphical systems tracts are visible: (1) a lowstand systems tract (LST) >10 kyrs BP. Pleistocene brownish soil superposing subaerially exposed Pleistocene reef limestone. (2) A transgressive systems tract (TST) 10-6.5 kyrs BP. A peat layer marks the beginning of the inundation, and the carbonate sedimentation starts with very low sedimentation rates of 0.02 m/kyr. (3) A highstand systems tract (HST) 6.5-0 kyrs BP, further divided into three stages (6.5-3, 3-1, 1-0 kyrs BP). The sea-level rise slowed down, sedimentation rates are increasing continuously up to a maximum of 1.4 m/kyr, the sand spit developed some 4 kyrs BP, the lagoonal circulation got restricted, and the lateral west-to-east gradient of grain-size accumulation started. From 1-0 kyrs BP the sedimentation rates slowed down to modern mean sedimentation rates of 0.6 m/kyr.
Two cores, one core from the center of the lagoon (core 16) and one core from the northern margin of the lagoon (core 19), have been analyzed on diversity and assemblages of benthic foraminifera in high-resolution. The transitions of Ammonia spp. to a more even and diverse fauna marks a significant environmental change at 7.0 kyrs BP in core 16 (onset of a stable environment in the deep lagoon after the sea-level rise slowed down at HST stage 1) and at 4.0 kyrs BP in core 19. A continuing environmental change after 1.4 kyrs BP in core 16 caused the fauna to become more even, a recovery of diversity and a permanent decline of foraminiferal accumulation rate. The changes in the faunas at 4.0 kyrs BP and at 1.4 kyrs BP could be explained with the sand spit formation in the northwestern and western lagoon. The sand spit has apparently acted as an obstacle in lagoonal circulation and might have caused unstable environmental conditions due to a more rapid circulation at the shallow marine site of core 19 and a slowdown of bottom water circulation in the main lagoon (core 16) leading to higher residence times and to lower oxygen and higher nutrient concentrations.
Derivation and characterization of a new filter for nonlinear high-dimensional data assimilation
(2015)
Data assimilation (DA) combines model forecasts with real-world observations to achieve an optimal estimate of the state of a dynamical system. The quality of predictions in nonlinear and chaotic systems such as atmospheric or oceanic circulation is strongly sensitive to the initial conditions. Therefore, beyond the consistent reconstruction of past states, a primary relevance of advanced DA methods concerns the proper model initialization. The ensemble Kalman filter (EnKF) and its deterministic variants, mostly square root filters such as the ensemble transform Kalman filter (ETKF), represent a popular alternative to variational DA schemes. They are applied in a wide range of research and operations. Their forecast step employs an ensemble integration that fully respects the nonlinear nature of the analyzed system. In the analysis step, they implicitly assume the prior state and observation errors to be Gaussian. Consequently, in nonlinear systems, the mean and covariance of the analysis ensemble are biased and these filters remain suboptimal. In contrast, the fully nonlinear, non-Gaussian particle filter (PF) relies on Bayes' theorem without further assumptions, which guarantees an exact asymptotic behavior. However, it is exposed to weight collapse, particularly in higher-dimensional settings, known as the curse of dimensionality.
This work presents a new method to obtain an analysis ensemble with mean and covariance that exactly match the corresponding Bayesian estimates. This is achieved by a deterministic matrix square root transformation of the forecast ensemble, and subsequently a suitable random rotation that significantly contributes to filter stability while preserving the required second-order statistics. The forecast step remains as in the ETKF. The algorithm, which is fairly easy to implement and computationally efficient, is referred to as the nonlinear ensemble transform filter (NETF). The limitation with respect to fully-nonlinear filtering is that the NETF only considers the mean and covariance of the Bayesian analysis density, neglecting higher-order moments.
The properties and performance of the proposed algorithm are investigated via a set of experiments. The results indicate that such a filter formulation can increase the analysis quality, even for relatively small ensemble sizes, compared to other ensemble filters in nonlinear, non-Gaussian scenarios. They also confirm that localization enhances the applicability of this PF-inspired scheme in larger-dimensional systems. Finally, the novel filter is coupled to a large-scale ocean general circulation model with a realistic observation scenario. The NETF remains stable with a small ensemble size and shows a consistent behavior. Additionally, its analyses exhibit low estimation errors, as revealed by a comparison with a free ensemble integration and the ETKF. The results confirm that, in principle, the filter can be applied successfully and as simple as the ETKF in high-dimensional problems. No further modifications are needed, even though the algorithm is only based on the particle weights. Thus, it is able to overcome the curse of dimensionality, even in deterministic systems. This proves that the NETF constitutes a promising and user-friendly method for nonlinear high-dimensional DA.
In the past sixty years, excessive water consumption and dam construction have significantly influenced natural flow regimes and surface freshwater ecosystems throughout China, and thus resulted in serious environmental problems. In order to balance the competing water demands between human and environment and provide knowledge on sustainable water management, assessments on anthropogenic flow alterations and their impacts on aquatic and riparian ecosystems in China are needed.
In this study, the first evaluation on quantitative relationships between anthropogenic flow alterations and ecological responses in eleven river basins and watersheds in China was performed based on the data that could be obtained from published case studies. Quantitative relationships between changes in average annual discharge, seasonal low flow and seasonal high flow and changes in ecological indicators (fish diversity, fish catch and vegetation cover, etc.) were analyzed. The results showed that changes in riparian vegetation cover as well as changes in fish diversity and fish catch were strongly correlated with the changes in flow magnitude (r = 0.77, 0.66), especially with changes in average annual river discharge. In addition, more than half of the variations in vegetation cover could be explained by changes in average annual river discharge (r² = 0.63) and roughly 50 % changes in fish catch in arid and semi-arid region and 60% changes of fish catch in humid region could be related to alterations in average annual river discharge (r² = 0.53, 0.58).
In a supplementary analysis of this study, the first estimation on quantitative relationships between decreases in native fish species richness and anthropogenic flow alterations in 34 river basins and sub-basins in China was conducted. Linear relationships between losses of native fish species and five ecologically relevant flow indicators were analyzed by single and multiple regression models. For the single regression analysis, significant linear relationships were detected for the indicators of long-term average annual discharge (ILTA) and statistical low flow Q90 (IQ90). For the multiple regressions, no indicator other than ILTA has significant relationships with changes in number of fish species mainly due to collinearity. Two conclusions emerged from the analysis: 1) losses of fish species were positively correlated with changes in ILTA in China and 2) indicator of ILTA was dominant over other flow indicators included in this research for the given dataset. These results provide a guideline for the sustainable water resources management in rivers with high risk of fish extinction in China.
Bayesian Networks are computer-based environmental models that are frequently used to support decision-making under uncertainty. Under data scarce conditions, Bayesian Networks can be developed, parameterized, and run based on expert knowledge only. However, the efficiency of expert-based Bayesian Network modeling is limited by the difficulty in deriving model inputs in the time available during expert workshops. This thesis therefore aimed at developing a simple and robust method for deriving conditional probability tables from expert estimates in a time-efficient way. The design and application of this new elicitation and conversion method is demonstrated using a case study in Xinjiang, Northwest China. The key characteristics of this method are its time-efficiency and the approach to use different conversion tables based on varying levels of confidence. Although the method has its limitations, e.g. it can only be applied for variables with one conditioning variable; it provides the opportunity to support the parameterization of Bayesian Networks which would otherwise remain half-finished due to time constraints. In addition, a case study in the Murray-Darling Basin, Australia, is used to compare Bayesian Network types and software to improve the presentation clarity of large Bayesian Networks. Both case studies aimed at gaining insights on how to improve the applicability of Bayesian Networks to support environmental management.
The Late Cretaceous is known to be mostly affected by warm periods interrupted temporarily by a number of cooling events. The reconstruction of the paleoclimatic conditions during a period of high concentration of CO2 in the atmosphere is of great importance for the creation of future climate models. We applied the recently developed method reconstructing the SST from the TEX86 (TetraEther indeX of tetraethers consisting of 86 carbon atoms).
The sample material used for the present study was obtained from the tropical Late Cretaceous southern Tethys upwelling system (Negev/Israel), lasting from the Late Santonian to the Early Maastrichtian (~ 85 to 68 Ma). On the core samples from the Shefela basin, representing the outer belt of the upwelling system and the outcrop profile from the open mine Mishor Rotem (Efe Syncline), representing the inner belt, various bulk geochemical and biomarker studies were performed in this thesis.
Derived from TEX86 data, a significant long-term SST cooling trend from 36.0 to 29.3 °C is recognized during the Late Santonian and the Early Campanian in the southern Tethys margin. This is consistent with the opening and deepening of the Equatorial Atlantic Gateway (EAG) and the intrusion of cooler deep water from the southern Atlantic Ocean influencing the global SSTs and also the Tethys Ocean. Furthermore, the cooler near shore SST usually found in modern upwelling systems could be verified in case of the ancient upwelling system investigated in the present study. The calculated mean SST in the inner belt (27.7 °C) represented in the Efe Syncline was 1.5 °C cooler in comparison to the more seaward located outer belt (Shefela basin).
Moreover, geochemical and biomarker analyses were used to identify both the accumulation of high amounts of phosphate in the PM and good preservation of organic matter (OM) in the lower part of the OSM section. Total organic carbon (TOC) contents are highly variable over the whole profile reaching from 0.6 % in the MM, to 24.5 % in the OSM. Total iron (TFe) varies from 0.1 % in the PM to 3.3 % in the OSM and total sulfur (TS) varies between 0.1 % in the MM and 3.4 % in the OSM. Different correlations of TS, TOC and TFe were used to identify the conditions during the deposition of the different facies types. Natural sulfurization was found to play a key role in the preservation of the OM particularly in the lower part of the OSM. Samples from the OSM and the PM were deposited under dysoxic to anoxic conditions and iron limitation lasted during the deposition of the OSM and the PM, which effected the incorporation of sulfur into OM.
Phosphorus is highly accumulated in the sediments of the PM with a mean proportion of 11.5 % total phosphorus (TP), which is drastically reduced to a mean value of 0.9 % in the OSM and the MM. From the correlation of the bulk geochemical parameters TOC/TOCOR ratio and TP a major contribution of sulfate reducing bacteria to the phosphate deposition is concluded. This interrelation has previously been investigated in recent coastal upwelling systems off Peru, Chile, California and Namibia. This was further supported by the analysis of branched and monounsaturated fatty acids indicating the occurrence of sulfate reducing and sulfide oxidizing bacteria during the deposition.
According to the results from the analysis of n-alkanes and C27- to C29-steranes up to 95 % of the OM was of marine origin.
Organic sulfur compounds (OSC) were a major compound class in the aromatic hydrocarbon fraction and n-Alkyl and isoprenoid thiophenes were the most abundant, with highest amounts found for 2-methyl-5-tridecyl-thiophene (28 µg/g TOC). The relatively high abundance of ββ-C35 hopanoid thiophenes and epithiosteranes is equivalent to an incorporation of sulfur during the early stages of diagenesis.
Moreover, the geochemical parameters δ13Corg, δ15Norg, C/N and the pristane/phytane (Pr/Ph) ratio, were studied for reconstruction of seafloor and water column depositional environments. The high C/N ratio along with relatively low values of δ15Norg (4 ‰ to 6 ‰) and δ13Corg (-29 ‰ to -28 ‰) are consistent with a significant preferential loss of nitrogen-rich organic compounds during diagenesis. Oxygen-depleted conditions lasted during the deposition of the PM and the bottom of the OSM, reflected by the low Pr/Ph ratio of 0.11–0.7. In the upper part of the OSM and the MM the conditions changed from anoxic to dysoxic or oxic conditions. This environmental trend is consistent with co-occurring foraminiferal assemblages in the studied succession and implies that the benthic species in the Negev sequence were adapted to persistent minimum oxygen conditions by performing complete denitrification as recently found in many modern benthic foraminifera.
Furthermore, the anammox process could have influenced the nitrogen composition of the sediments. In this anaerobically process nitrite and ammonia are converted to molecular nitrogen.
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.
Die Anwendung der WKB-Theorie zur Simulation der schwach nichtlinearen Dynamik von Schwerewellen
(2014)
Es ist schon seit Längerem bekannt, dass Schwerewellen die Zirkulation der mittleren Atmosphäre beeinflussen. Sie werden fast ausschließlich in der Troposphäre durch Prozesse wie Gebirgsüberströmung, Konvektionen, Frontogenese etc. erzeugt. Sie propagieren von ihrem Entstehungsort in der Troposphäre in die höheren Schichten der Atmosphäre und transportieren dabei ihre Energie und ihren Impuls. Unter der Voraussetzung, dass die Energie von Schwerewellen erhalten bleibt und die Dichte der Atmosphäre mit der Höhe exponentiell abnimmt, wächst die Amplitude der Schwerewellen so stark an, dass sie brechen und ihren Impuls in Stratosphäre und Mesosphäre deponieren. Als Folge davon beeinflussen Schwerewellen die großräumige Zirkulation der Atmosphäre und sind damit ein wichtiges Bindeglied, welches die Troposphäre mit anderen Atmosphärenbereichen verbindet. Folglich ist es wichtig, dass die Klima- und Wettermodelle in der Lage sind, die Schwerewellendynamik zu beschreiben. Bedauerlicherweise können diese Modelle nicht das komplette Schwerewellenspektrum auflösen. Somit müssen Schwerewellen in den Modellen parametrisiert werden. Viele Parametrisierungsschemen basieren auf Wentzel-Kramer-Brillouin(WKB)-Theorie. Die WKB-Gleichungen, die sogenannten Strahlengleichungen, beschreiben die räumliche und zeitliche Variation der Welleneigenschaften wie Wellenzahl, Wellenamplitude und Wellenfrequenz entlang der Charakteristiken, welche durch die lokale Gruppengeschwindigkeit vorgegeben sind. Die numerische Modelle, die auf den Strahlengleichungen basieren, werden als Strahlenmodelle bezeichnet. In Strahlenmodellen werden Schwerewellen durch Wellenteilchen dargestellt. Zur Zeit verwenden die Strahlenmodelle stationäre Strahlengleichungen, da die Wechselwirkung eines zeitabhängigen Schwerewellenfeldes mit einem zeit- und ortsabhängigen Hintergrund zu Problemen in Strahlenmodellen führen kann. Die Strahlengleichungen basieren auf der Annahme, dass sich nie zwei Wellenteilchen mit den unterschiedlichen Welleneigenschaften an einer Position befinden können. Wenn an einer Position zwei Wellenteilchen mit den unterschiedlichen Wellenzahlen befinden, entsteht sogenannte Kaustik: ein Punkt im Raum, an dem sich mehrere Charakteristiken kreuzen. Wenn eine Kaustik entsteht, kann die Wellenamplitude nicht mehr bestimmt werden. Ziel der vorliegenden Arbeit ist es mithilfe der WKB-Theorie die Ausbreitung von Schwerewellenpaketen in einer raum- und zeitabhängigen Hintergrundströmung zu beschreiben und ein numerisches Modell zu entwickeln, welches die Schwerewellen parametrisieren und ihre Wechselwirkung mit der raum- und zeitabhängigen Hintergrundströmung beschreiben kann.
Einfachheitshalber wird in dieser Arbeit nur die Wechselwirkung zwischen horizontal periodischen, vertikal lokalisierten Schwerewellenpaketen und der raum- und zeitabhängigen Hintergrundströmung betrachtet.
Interpretation magnetotellurischer Messungen auf Island und 3D-Modellierungen des Island-Plumes
(2014)
Die Lage und der geologische Aufbau Islands sind eine einmalige geologische und geophysikalische Besonderheit - nur hier liegt der Mittelatlantische Rücken deutlich oberhalb des Meeresspiegels. Eine Theorie ist, dass unter Island ein sogenannter Hotspot existiert - ein aufsteigender Mantelplume (oder eine andere Temperatur- und Schmelzanomalie). Das bedeutet, dass die Tektonik, die Petrologie und die (geo)physikalischen Parameter von der Wechselwirkung der Meeresbodenspreizung mit dem aufsteigenden Mantelplume bestimmt werden. Die magnetotellurischen Untersuchungen auf Island und die magnetotellurischen Modelle der vorliegenden Arbeit lassen sich nach Zieltiefen unterteilen. Im ersten Teil steht vor allem die Untersuchung der Leitfähigkeitsstruktur des oberen Erdmantels im Vordergrund, wobei der Schwerpunkt auf Modellstudien liegt. Hier wird die Theorie des aufsteigenden Mantelplumes aufgegriffen und hinsichtlich der Möglichkeiten der Magnetotellurik (MT) diesen zu detektieren, untersucht. In dieser Arbeit werden deshalb verschiedene synthetische 3D-Modelle untersucht. Diese sollen klären, unter welchen Bedingungen, d.h. bei welcher Tiefenlage, Ausdehnung und Leitfähigkeit, ein Mantelplume in der Magnetotellurik messbare Signale liefert. Da Temperatur und Schmelzanteil hierbei eine entscheidende Rolle spielen, werden auch Ergebnisse von geodynamischen Modellierungen verwendet. Die verschiedenen Modellierungen machen deutlich, dass der Plume in Modellen mit maximal 1% Schmelzanteil nur sehr schlecht mit der MT aufgelöst werden kann. Ein weiteres Problem ist die relativ geringe laterale Ausdehnung des Bereichs hoher Leitfähigkeit. Insbesondere bei Berücksichtigung des Krustenleiters liegen die Effekte des Plumes im Bereich der Messfehler bei realen Daten. Erst bei 3% Schmelzanteil im Plume ergeben sich Charakteristiken in den Sondierungskurven, die eindeutig dem Plume zugeordnet werden können. Im zweiten Teil liegt das Augenmerk auf der Leitfähigkeitsstruktur der isländischen Kruste, für deren Untersuchung ca. 200 Datensätze aus zahlreichen deutsch-isländischen Messkampagnen zur Verfügung stehen. In dieser Arbeit werden sie ganzheitlich ausgewertet und interpretiert. Hier steht im Speziellen der in verschiedenen früheren Untersuchungen gefundene gute Leiter in ca. 10 km Tiefe im Vordergrund. Seine Ausdehnung, Tiefenlage und Leitfähigkeit sowie laterale Änderungen werden genauer untersucht. Dazu werden aus den gemessenen Daten und den daraus ermittelten magnetotellurischen Impedanztensoren nicht nur die Sondierungskurven, sondern weitere Parameter wie Phasentensoren und Rotationsinvarianten abgeleitet und diskutiert. Um einen Überblick über die räumliche Verteilung der Untergrundstrukturen zu erhalten, werden die einzelnen Messpunkte charakterisiert und Stationen mit ähnlichem Verhalten in Gruppen zusammengefasst. 1D- und 2D- Inversionsmodelle liefern einen direkten Hinweis auf die Leitfähigkeitsverteilung. Diese zeigt ausgeprägte laterale Änderungen. Aus den 2D-Modellen ergeben sich Hinweise, die darauf schließen lassen, dass diese Zone niedrigen Widerstands nicht durchgängig oder geschlossen sein muss, um die gemessenen Daten zu erklären. Der gute Leiter scheint weniger mit der Aufwölbung durch den vermuteten Mantelplume oder eine andere tiefer liegende Schmelzanomalie korreliert zu sein. Die lateralen und vertikalen Variationen lassen auf einen Zusammenhang mit den Vulkansystemen auf Island schließen. Darauf deuten auch die Induktionsvektoren hin, die lateral betrachtet stark variieren. Ein weiteres Ergebnis ist das Fehlen dieses Krustenleiters in einem mindestens 20km breiten Streifen entlang der Südküste.
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.
The objective of the present doctoral thesis was to investigate the occurrence, distribution, and behaviour of six hydrophilic ethers: ethyl tert-butyl ether (ETBE), 1,4-dioxane, ethylene glycol dimethyl ether (monoglyme), diethylene glycol dimethyl ether (diglyme), triethylene glycol dimethyl ether (triglyme), and tetraethylene glycol dimethyl ether (tetraglyme) in surface-, waste-, ground- and drinking water samples. Solid phase extraction and gas chromatography/mass spectrometry were used to analyze the six hydrophilic ethers. Altogether more than 150 surface water samples, almost 100 of each groundwater and wastewater samples, and 10 raw and drinking water samples were analyzed during the research project.
Initially, the method was validated in order to simultaneously determine the analytes of interest in various aquatic environments. A solid phase extraction method that uses coconut charcoal (Resprep® activated coconut charcoal, Restek) or carbon molecular sieve material (SupelcleanTM Envi-CarbTM Plus, Supelco) for analyte absorption were found suitable for determination of ETBE, 1,4-dioxane, and glymes in surface-, drinking-, ground- and wastewater samples. Precision and accuracy of both methods was demonstrated for all analytes of interest. The recovery of target compounds from the ultrapure water spiked at 1.0 µg L−1 was between 86.8 % and 98.2 %, with relative standard deviation below 6 %. The samples spiked at 10.0 µg L−1 gave slightly higher recovery of 90.6 % to 112.2 % with a relative standard deviation below 3.4 % for each analyte. Detection and quantification limits in ultrapure water and surface waters were furthermore established. The limit of quantitation (LOQ) in ultrapure water ranged between 0.024 µg L−1 to 0.057 µg L−1 using Restek cartridges, and 0.030 µg L−1 to 0.069 µg L−1 using Supelco cartridges. In the surface water samples the calculated LOQ was 0.032 µg L−1 to 0.067µg L−1 using coconut charcoal material and 0.032 µg L−1 to 0.052 µg L−1 using the carbon molecular sieve material. Moreover, stability of the unpreserved and preserved water samples as well as the extracts was determined. Preservation of samples with sodium bisulfate (at 1 gram per Liter) resulted in much better stability of the ethers in water samples. Subsequently, 27 samples obtained from seven surface water bodies in Germany (Rivers Rhine, Lippe, Main, Oder, Rur, Schwarzbach and Wesel-Datteln Canal) were analyzed for the six hydrophilic ethers. ETBE was present in only two surface waters (Rhine River and Wesel-Datteln Canal) with concentrations close to the LOQ (up to 0.065 µg L−1). 1,4-Dioxane was detected in all of the water samples at concentrations reaching 1.93 µg L–1. Monoglyme was identified only in the Main and Rhine Rivers at the maximum concentration of 0.114 µg L–1 and 0.427 µg L–1, respectively. Very high concentrations (up to 1.73 µg L−1) of diglyme, triglyme, and tetraglyme were detected in the samples from the Oder River. These glymes were also detected in the Rhine River; however the concentrations did not exceed 0.200 µg L–1. Furthermore, tetraglyme was detected in the Main River at an average concentration of 0.409 µg L–1 (n = 6) and in one sample from the Rur River at 0.192 µg L–1.
Four sampling campaigns were conducted at the Oderbruch polder between October 2009 and May 2012, in order to study the behavior of the hydrophilic ethers and organophosphates during riverbank filtration and in the anoxic aquifer. Moreover the suitability of these target compounds was assessed for their use as groundwater organic tracers. At the time of each sampling campaign, concentrations of triglyme and tetraglyme in the Oder River were between 20–185 ng L–1 (n = 4) and 273¬–1576 ng L–1 (n = 4). Monoglyme, diglyme, and 1,4-dioxane were analyzed only during the two last sampling campaigns. At that time, the concentration of diglyme in Oder River was 65¬–94 ng L-1 (n = 2) and 1,4-dioxane 1610¬–3290 ng L–1 (n = 2). In the drainage ditch, following bank filtration, concentrations of ethers ranged between 1090 ng L–1 and 1467 ng L–1 for 1,4-dioxane, 23¬ng L–1 and 41 ng L–1 for diglyme, 37 ng L–1 and 149 ng L–1 for triglyme, and 496 ng L–1 and 1403 ng L–1 for tetraglyme. In the anoxic aquifer, 1,4-dioxane showed the greatest persistence during the groundwater passage. At the distance of 1150 m from the river and an estimated groundwater age of 41.9 years, a concentration above 200 ng L−1 was detected. A positive correlation was found for the inorganic tracer chloride (Cl−) with 1,4-dioxane and tetraglyme. Similarities in the behavior of Cl− and the organic compound suggested that 1,4-dioxane and tetraglyme are controlled by the same hydraulic process and therefore can be used as additional tracers to study the dynamics of the groundwater system. These results show that high concentrations of ethers are present in the surface water and are not removed during bank filtration processes. Moreover, the hydrophilic ethers persist in the anoxic aquifer and little or no degradation is expected, supporting, their possible application as organic tracers.
A separate sampling project was conducted for 1,4-dioxane that focused primarily on its fate in the aquatic environment. This study provided missing information on the extent of water pollution with 1,4-dioxane is Germany. Numerous waste-, surface-, ground- and drinking water samples were collected in order to determine the persistence of 1,4-dioxane in the aquatic environment. The occurrence of 1,4-dioxane was determined in wastewater samples from four municipal sewage treatment plants (STP). The influent and effluent samples were collected during weekly campaigns. The average influent concentrations in all four plants ranged from 262 ± 32 ng L−1 to 834 ± 480 ng L−1, whereas the average effluents concentrations were between 267 ± 35 ng L−1 and 62,260 ± 36,000 ng L−1. The source of increased 1,4-dioxane concentrations in one of the effluents was identified to originate from impurities in the methanol used in the postanoxic denitrification process. Spatial and temporal distribution of 1,4-dioxane in the river Main, Rhine, and Oder was also examined. Concentrations reaching 2,200 ng L−1 in the Oder River, and 860 ng L−1 in both Main and Rhine River were detected. The average load during the sampling was estimated to be 6.5 kg d−1 in the Main, 34.1 kg d−1 in the Oder, and 134.5 kg d−1 in the Rhine River. In all of the sampled rivers, concentrations of 1,4-dioxane increased with distance from the mouth of the river and were found to negatively correlate with the discharge of the river. In order to determine if 1,4-dioxane can reach drinking water supplies, samples from a Rhine River bank filtration site and potable water from two drinking water production facilities were analyzed for the presence of 1,4-dioxane in the raw water and finished potable water. The raw water (following bank filtration) contained 650 ng L−1 to 670 ng L−1 of 1,4-dioxane, whereas the concentration in the finished drinking water fell only to 600 ng L−1 and 490 ng L−1, respectively.
During the final project, investigations of the source identification of high glyme concentrations in the Oder River were carried out. During four sampling campaigns between January, 2012 and April, 2013, 50 samples from the Oder River in the Oderbruch region and Poland were collected. During the first two samplings in the Oderbruch polder, glymes were detected at concentration reaching 0.07 µg L-1 (diglyme), 0.54 µg L−1 (triglyme) and 1.73 µg L−1 (tetraglyme) in the Oder River. The extensive sampling campaign of the Oder River (about 500 km) in Poland helped to identify the area of possible glyme entry into the river. During that sampling the maximum concentrations of triglyme and tetraglyme were 0.46 µg L−1 and 2.21 µg L−1, respectively. A closer investigation of the identified area of pollution, helped to determine the possible sources of glymes in the Oder River. Hence, the final sampling focused on the Kaczawa River, a left tributary of the Oder River and Czarna Woda, a left tributary of Kaczawa River. Moreover, samples from an industrial wastewater treatment plant were collected. Samples from Czarna Woda stream and Kaczawa River contained even higher concentrations of diglyme, triglyme, and tetraglyme, reaching 5.18 µg L−1, 12.87 µg L−1 and 80.81 µg L−1, respectively. Finally, three water samples from a wastewater treatment plant receiving influents from a copper smelter were analyzed. Diglyme, triglyme, and tetraglyme were present at an average concentration of 569 µg L−1, 4300 µg L−1, and 65900 µg L−1, respectively in the wastewater. Further research helped to identify the source of the glymes in the wastewater. The gas desulfurization process – Solinox implemented in the nearby copper smelter uses glymes as physical absorption medium for sulfur dioxide.
Results of this doctoral research provide important information about the occurrence, distribution, and behavior of hydrophilic ethers: 1,4-dioxane, monoglyme, diglyme, triglyme, and tetraglyme in the aquatic environment. A method capable of analyzing a wide range of ether compounds: from a volatile ETBE to a high molecular weight tetraglyme was validated. 1,4-Dioxane and tetraglyme were found to be applicable as organic tracers, since they are not easily attenuated during bank filtration and the anoxic groundwater passage. The extent of water pollution with 1,4-dioxane was shown in waste-, surface-, ground-, and drinking waters. One source of extremely high concentrations of 1,4-dioxane in a municipal sewage treatment plant applying postanoxic denitrification was identified, however more information is needed on the entry of 1,4-dioxane into surface waters. Moreover, 1,4-dioxane was present in drinking water samples from river bank filtration, which demonstrates its persistence in the aquatic environment and its low degradation potential during bank filtration and subsequent water treatment. Furthermore, this was the first study that focused primarily on identifying sources of glymes in surface waters. Glymes find a widespread use in industrial sectors, hence establishing their origin in the surface water is difficult (as with 1,4-dioxane). In this work, a gas desulphurization process was identified to be a dominating source of glyme pollution in the Oder River.
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.
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.
Forty two samples of the Late Eocene Kiliran oil shale, Central Sumatra Basin, Indonesia were collected from a 102 m long drill core. Palynofacies and geochemical analyses have been carried out to reconstruct the paleoenvironmental conditions and paleoecology during deposition of the oil shale. Amorphous organic matter (AOM) is very abundant (>76%). B. braunii palynomorph is present (3-16%) as the only autochtonous structured organic matter and generally more abundant in middle part of the profile. The stable carbon isotopic composition of organic matter (δ13C) varies from -27.0 to -30.5‰ and is generally more depleted in middle part of the profile. The ratio of total organic carbon to sulfur (TOC/S), used as salinity indicator, ranges from 2.5 to 15.8 and shows variations along the profile. Relatively less saline environments are observed in the middle part profile. Fungal remains are generally present only in middle part of the profile with distinct peak of abundances. The presence of fungal remains is regarded as an indication for a relatively warmer climate during deposition of middle part of the profile. The warmer climate is thought to influence the establishment of a thermocline, limiting the supply of recycled nutrients to the epilimnion. Consequently, the primary productivity in the Kiliran lake decreased during deposition of the middle part of the profile as indicated by the relatively depleted δ13C and the blooming of B. braunii. The chemocline was also shoaling during deposition of the middle part of the profile according to the higher abundance of isorenieratene derivatives of green sulfur bacteria origin. The warmer climate affected also to increase of water supply and thus less saline environments.
Tectonic subsidence is also thought to be a significant factor for the development of the Kiliran lake. The Zr/Rb ratio, an indicator for grain size, ranges from 0.4 to 1.3 and generally increases upwards along the profile. Three sudden decreases of the ratio are observed, indicating rapid change to finer grain size. These decreases are interpreted to indicate rapid deepening events of the lake due to mainly periodic subsidence. During deposition of lower part of the profile, the subsidence rates might have been relatively higher than sediment and water supply rates, resulting in a higher autochtonous fraction in the oil shale. During deposition of middle part of the profile, the sediment and water supply rates were relatively higher promoting distinct progradational sedimentation. Subsequently, the lake became more shallow and smaller during deposition of the upper part of the profile, leading to a relatively higher terrigenous input to the oil shale.
Norneohop-13(18)-ene and neohop-13(18)-ene derived from methanotrophic bacteria are the dominant hopanoid hydrocarbons. The sum of their concentrations varies from 40.6 to 360.0 μg/g TOC. The δ13C of these compounds are extremely depleted (-45.2 to -50.2‰). The occurrence of abundant bacteria including methanotrophic bacteria was responsible for the recycling of carbon below the chemocline of the lake. The effect of the recycling of carbon is observed by the presence of a concomitant depletion (about 7-9‰) in 13C of some specific biomarkers derived from organisms dwelling in the whole phototrophic zone.
4-Methylsterane and 4-methyldiasterene homologues occur in the oil shale as the predominant biomarkers. The sum of the concentrations of all homologues are about 40.3-1,009.2 μg/g TOC with generally higher values in uppermost and lower parts of the profile. Ca accounts as the predominant element in the oil shale, ranging from 5.0 to 16.7%. This element shows generally parallel variation with the 4-methylsterane homologues along the profile. This suggests that the 4-methylsteranes were derived from biological sources favoring more alkaline and more trophic environments. On the other hand, these compounds were less abundant in middle part of the profile which is consistent with less alkaline and less trophic environments promoting B. braunii to bloom.
The 4-methylsterane homologues are considered to originate from Dinoflagellates. Alternation between Dinoflagellates and B. braunii in Paleogene lake systems due to water chemistry changes are known from previous studies. Moreover, freshwater Dinoflagellates have been frequently reported to occur in the basin depocenters. In the present case, distinct alternation between B. braunii abundances and concentrations of 4-methylsterane homologues along the studied oil shale profile suggest that the 4-methylsterane homologues were derived from freshwater Dinoflagellates although dinosterane is not present in the sediment extracts. Water alkalinity and trophic level changes were most likely responsible for the alternation of Dinoflagellates and B. braunii blooming.
The mantle xenoliths collected by kimberlites indicate that the subcratonic mantle underneath the Archean crust is mostly a residue of high degrees of partial melting which was subsequently reenriched. The majority of the xenoliths show cryptic metasomatism and only few modal metasomatism.
Much effort has been put into deciphering different kinds of enrichment processes within the mantle. Here, we take the approach to look into the inventory of subcalcic garnets which stem from cpx-free harzburgites and dunites. These subcalcic garnets, commonly with sinusoidal REE patterns, carry the major budget of the trace elements of their host rock. Thus, they are promising objects to study both depletion and enrichment. Most importantly, the analysis of a single grain subcalcic garnetwill provide almost all important information of the bulk rock. Our aim is to gain detailed information mainly on metasomatism on a craton wide scale by combining major, trace elements and Lu-Hf and Sm-Nd isotopic signatures from subcalcic garnets. Eventually, we will summarize the metasomatic agent(s) and processes and possibly the timing of the enrichment within the lithospheric mantle underneath the Kaapvaal craton.