Refine
Year of publication
Document Type
- Doctoral Thesis (87) (remove)
Language
- English (87) (remove)
Has Fulltext
- yes (87)
Is part of the Bibliography
- no (87)
Keywords
- Boden (2)
- Caribbean (2)
- China (2)
- Roman (2)
- bank filtration (2)
- degradation (2)
- groundwater (2)
- 1,4-dioxane (1)
- 3D-Modellierung (1)
- Abbau (1)
Institute
- Geowissenschaften (87) (remove)
The biomarker record in two different lakes in central Europe, Lake Albano and Lake Constance, is used to reflect environmental changes and lake system response during the Late Glacial and Holocene. Extractable organic compounds in lake sediments, which can be assigned to their biological source (biomarkers) function as fingerprints of past aquatic or land plant organisms. Using gas chromatography coupled with mass spectrometry, 21 different biomarkers (predominantly steroids and triterpenoids) as well as a variety of n-alkanes, nalkanols, and n-alkanoic acids could be identified in the sediment records of Lake Albano and Lake Constance. In the Holocene sediments of Lake Albano, the distribution of biomarkers such as dinosterol (dinoflagellates), isoarborinol, and diplopterol (aquatic organisms) indicate three biomarker zones: The period between 0-3,800 years BP (zone 3) is characterized by high concentrations of these biomarkers and others such as tetrahymanol and diploptene. Conversely, zone 2 (3,800-6,500 years BP) shows very low concentrations of all autochthonous biomarkers. In zone 1 (6,500–11,480 years BP), dinosterol, isoarborinol, and diplopterol range on a relatively high level, whereas diploptene and tetrahymanol display comparatively low concentrations. The results suggest at least two distinct changes in the predominance of primary producers during the Holocene, which are related to changes in the lake system such as lake mixing and water column stratification. This interpretation is consistent with previous investigations of Lake Albano sediments including pigment and hydrogen index data (Ariztegui et al., 1996b; Guilizzoni et al., 2002). Allochthonous biomarkers such as long-chain n-alkanes, amyrenones and friedelin indicate a development from forest to a more open landscape from 6,000 and 5.000 years BP, respectively. After a period of high concentrations during the first half of the Holocene, all biomarkers derived from deciduous trees exhibit relatively low values until around 1,000 years BP. Again, this is consistent with results from previous pollen investigations (Ariztegui et al., 2000). The sediment core from Upper Lake Constance comprises the Late Glacial and Holocene. It was analysed for biomarkers and inorganic tracers in order to compare the biomarker results with other proxy data from the same core. Magnetic susceptibility (MS) was measured to get a high-resolution stratigraphic framework of the core and to obtain further information about changes of the proportions of allochthonous and autochthonous input. Enhanced concentrations and accumulation rates of dinosterol (biomarker for dinoflagellates) and biogenic calcite give evidence of increasing lake productivity at the beginning of the Holocene followed by a decrease in bioproductivity after around 7,000 years BP. Younger Dryas sediments are characterized by low amounts of both dinosterol and biogenic calcite indicating a low productivity. The comparison of the concentrations and accumulation rates of b-sitosterol and stigmastanol with parameters reflecting lake productivity suggests that both steroids in Lake Constance sediments are mainly derived from terrigenous sources. Biomarkers as well as concentrations and accumulation rates of allochthonous inorganic compounds such as titanium, magnesium and strontium indicate a slightly enhanced allochthonous input after 8,500 years BP. Significant increase of erosive matter input from enhanced soil erosion is not observed before 4,000 years BP. This can be attributed to the combined effects of precipitation increase as a result of climatic deterioration and anthropogenic deforestation which is consistent with observations from other lakes in Central Europe. The MS record of Lake Constance confirms these results by tracing the climatically induced shifts of more intense bioproduction (low MS caused by increased calcite deposition) during the ‘climatic optimum’. This is followed by increasing input of terrigenous sediment compounds during colder and wetter periods which lead to higher MS values in the lake sediments. The occurrence of tetrahymanol in Lake Constance sediments questions the unambiguous use of tetrahymanol as an indicator for water column stratification. Anaerobic organic macroaggregates within the oxygenated, photic zone of the water column have to be considered as a possible living space for anaerobic microorganisms containing tetrahymanol. The direct comparison of two very different lakes Albano and Constance with respect to biomarkers indicating climate or environmental change provides a contribution to the recent biomarker research for a better understanding of biomarkers in lacustrine sediments.
An eclogite barometer has profound importance in the study of upper mantle processes and potential application to diamond prospecting. Studies on the partitioning of Li between clinopyroxene (cpx) and garnet (grt) in natural samples have shown that this particular element is very sensitive to changes in pressure and could be calibrated as the barometer demanded for bimineralic eclogites. Experiments were performed from 4 to 13 GPa and 1100 to 1400°C in the CMAS (CaO, MgO, Al2O3, SiO2) system with Li added as Li3PO4 to quantify this pressure dependence into a barometer expressed in the following equation: P= (0.00255*T-lnKd)/0.2351 where P is in GPa, T is in °C and Kd is defined as the partition coefficient of Li (in ppm) between clinopyroxene and garnet. The experimental pressures are reproduced to ± 0.38 GPa (1σ) by this equation. This barometer is strictly applicable only to CMAS. Experiments at 1300°C, 8-12 GPa showed that Henry’s Law is fulfilled for Li partitioning between cpx and grt in the concentration range of approximately 0.01 – 1 wt% Li. Direct application of the equation to experiments in natural systems performed at 1300°C from 4 GPa to 13 GPa consistently overestimates pressures by approximately 2 GPa. Our previous experiments in the system CaO-MgO-Al2O3-SiO2 + Li3PO4 showed that the partitioning of Li between garnet and clinopyroxene is pressure dependent in eclogitic bulk compositions. This supports experimentally the hypothesis of Seitz et al. (2003), based on the analysis of Li in eclogitic xenoliths and inclusions in diamond, that the partitioning of this particular element between clinopyroxene and garnet is very sensitive to changes in pressure and could be calibrated as a barometer for bimineralic eclogites. In order to calibrate this pressure dependence into a barometer, experiments were performed in natural systems using starting materials sourced from a well preserved eclogitic xenolith from the Roberts Victor kimberlite pipe (South Africa) to extrapolate our findings in CMAS to natural systems. Sixteen multianvil experiments were performed from 4-13 GPa and 1100-1500°C. Our findings reinforced the general trend we observed in the CMAS system, that KdLi cpx-grt for Li decreases with increasing P, and that at P ≥ 12 GPa, garnet is able to incorporate more Li than clinopyroxene. Multiple linear regression was applied to our experimental results to create the barometer: P = (0.000963*T – ln KdLi cpx-grt + 1.581) / 0.252 Where P is pressure in GPa, T is temperature in °C and KdLi cpx-grt is defined as the partitioning coefficient of Li obtained by dividing the concentration of Li in cpx by the concentration of Li in garnet. This barometer reproduces the experimental conditions to ± 0.2 GPa. It is applicable to eclogitic xenoliths, to garnet pyroxenites and to peridotitic and eclogitic inclusions in diamond. Application of the barometer to diamond bearing xenoliths results in pressures in the diamond stability field. Clinopyroxene is easily corrupted in xenoliths and also preferentially takes in Li during short lived metasomatic processes. Care must be taken therefore to analyse primary, unaltered clinopyroxene. Our preliminary application to natural samples shows that the barometer can be applied beyond the experimental range to pressures down to 3 GPa. Seventeen eclogitic xenoliths were chosen from a sample set of greater than 200 for their fresh microscopic and macroscopic appearance and were analyzed for Li content in coexisting garnet (grt) and clinopyroxene (cpx). These samples can be subdivided into two groups on the basis of Mg in cpx (cpfu: cations per formula unit, based on 6 oxygens): Group 1 with Mg > 0.75, and Group 2 with Mg < 0.75. Group 1 xenoliths show lower Li contents in both grt and cpx compared to Group 2. The Li barom ter calibrated in Hanrahan et al. (2009b)/Chapter 3 was applied to these samples as well as available literature data to obtain pressures of provenance - Group 2 xenoliths often provide pressures that appear unrealistic for eclogitic xenoliths. In light of observed crystal chemical relations in the natural samples, a new fitting procedure was applied to the experimental data presented in Chapter 3. This new fit appears to be more realistic than the previous fit, although a strong relationship with Mg# remains present, suggesting that Li-barometry is, at present, only applicable to Mg-rich eclogites. Inclusions in diamond, with the exception of eclogitic inclusions of coexisting majorite and cpx, often yield pressures that are inconsistent with the pressures required for diamond formation. Although an interesting observation when comparing all of the data is that inclusions in diamond have significantly higher average Li concentrations compared to xenoliths, which suggests that Li is highly present in the fluids from which diamonds form in the mantle, an observation which was previously made for the deep mantle as a result of high Li in ferropericlase inclusions in diamond (Seitz et al. 2003).
Many natural minerals exist in the form of a solid solution. The systematic changes in structural and physical properties of oxide solid solutions are of geological importance and allow for wide applications. In order to understand the composition-structureproperty relations, substitutional solid solutions of CuxZn2−xTiO4, ZnxMg1−xTi2O5 and CuxMg1−xTi2O5 have been synthesised by mechanochemical activation assisted solid state synthesis. Self-propagating high-temperature synthesis has been employed to achieve the interstitial solid solutions of Ti5Si3Zx (Z refers to the element boron or oxygen).
The changes in the crystal structure and physical properties due to the formation of solid solutions are investigated by employing X-ray diffraction, neutron diffraction, Raman spectroscopy, low-temperature heat capacity, thermal expansion, scanning electron microscopy, UV-vis spectroscopy, plane-wave ultrasound spectroscopy and density functional theory calculations.
Enhanced aridification of Central Asia driven by the combined effects of orogenic surface uplift, Paratethys retreat, changes in atmospheric moisture transport and global cooling is one of the most prominent Cenozoic climate change events of the Northern Hemisphere. Deciphering regional long-term patterns of Central Asian hydrology is, therefore, a key element in understanding the role of Northern Hemisphere mid-latitude drying in the global hydrological system. Pedogenic carbonates record information of the paleoecosystem, the paleogeography, hydrology, tectonic and climatic conditions as well as the sedimentary regime during the time in which they formed. The calcrete-bearing paleosols in continental sedimentary basins yield the possibility for understanding these processes, mechanisms and controlling factors. This study characterizes long-term paleoenvironmental conditions between the late Oligocene and early Miocene in SE Kazakhstan based on stable isotopes, elemental geochemistry and laser ablation U-Pb geochronology from alluvial, fluvial and pedogenic deposits. The main topics addressed in this thesis are paleosol and calcrete formation processes in the light of geochemistry, tectonic and climatic influences on both, as well as the paleoenvironmental and hydrological conditions which led to these processes.
The sedimentological succession of the 14 km-long Kendyrlisai Valley in the Ili Basin, SE Kazakhstan in Central Asia, yields exceptional well-developed paleosols which provide an insight into the relationship between tectonics and sedimentation as well as soil formation processes. The valley accommodates more than 600 m of siliciclastic sediments deposited in a fluvial to alluvial environment in the lower part which grade into lacustrine to fluvial deposition with minor gypsum accumulation in the upper part. The yet undated sediments of Kendyrlisai Valley are compared with a biostratigraphically dated and well described succession in an adjacent location. The correlation of the two successions revealed a proximal-distal relation with the Kendyrlisai Valley succession represents the basin margin position with coarse alluvial-fluvial sediments. The calcrete nodules in Kendyrlisai Valley paleosols give the opportunity of U-Pb dating with LA-ICP-MS which uses the variation in both U and Pb to directly date the calcite. The U-Pb dating revealed a late Oligocene-early Miocene age for the investigated section. The analysis of facies and geochemistry of the paleosol profiles in Kendyrlisai Valley gives insight into calcretization stages, pedogenic processes and landscape stability. Kendyrlisai Valley paleosols show a low to moderate weathering intensity with the highest weathering intensity between 0.5 and 0.8 m depth in the paleosol profile. The comparison between acid leachable and non-acid leachable fraction indicate an incorporation of Ca and Ba in incipient calcrete calcite, whereas clay minerals adsorb Mg and Sr. During the evolution from early to more mature calcrete stages, i.e. calcrete nodules, Ba is lost by dissolution and subsequently adsorbed on clay minerals. The nodules consist almost exclusively of pure calcium carbonate with rare substitution of Ca by Mg indicating a constant supply of Ca by weathering and solution during calcrete formation. The occurrence of Mg-bearing clay minerals indicates high evaporative conditions with alkaline waters (pH >9).
Sampled calcrete nodule cross sections reveal more than one weathering and calcrete formation phase with a higher supply of Ca, Mg and Fe during early phases. The organic and inorganic carbon isotope composition of the calcrete nodules reflects C3 metabolism under occasionally moisture stress, resulting in higher δ13C values and lower respiration rates. This study also presents calculated atmospheric pCO2 values from the calcrete nodules with the equation from Cerling et al. (1999). The average calculated CO2 concentration for Kendyrlisai Valley paleosols is 313 ± 110 ppmv. The comparison with data from other studies during the late Oligocene–early Miocene reveal lower values for Kendyrlisai Valley paleosols, which may stem from an underestimation of the plant-respired CO2.
The knowledge of the variability within the paleosols and calcrete nodules enables a reliable paleoclimatic interpretation for the Kendyrlisai Valley and beyond the Ili Basin. Sedimentary facies and geochemical weathering indices suggest an increased surface and groundwater discharge fed by orographically-enhanced precipitation in the Tien Shan hinterland. In contrast, pedogenic stable isotope data and elevated rates of Mg fixation in clay minerals mirror enhanced rates of evaporation in the vadose zone due to protracted aridification. This study posits that pronounced surface uplift of the Tien Shan Mountains during the Oligocene–Miocene transition promoted regionally increased orographic precipitation and the development of fluvial discharge systems.
The comparison with other studies in adjacent basins creates a precipitation pattern for Central Asia during the late Oligocene to early Miocene. The westerlies supplied Central Asia with atmospheric moisture from the West, possibly from the Eastern Paratethys. The uplifting Tien Shan Mountain ranges captured this moisture on their luv-side, resulting in a pronounced orographic rainfall in the adjacent Ili and Issyk Kul Basins. The Tarim Basin and the Valley of Lakes experienced drier conditions due to the rain shadow effect on the lee-side of the Tien Shan and Altai Mountains. In course of this findings, the thesis highlights the crucial role the Tien Shan Mountains play and had been playing in former times as an orographic barrier for the distribution of atmospheric moisture across Central Asia.
The TTL is the transition layer between the tropical troposphere and stratosphere, and is the main region where tropospheric air enters the stratosphere. In this thesis different transport processes are studied by using in situ measurements of tracers. Long-lived tracers were measured with the High Altitude Gas Analyzer (HAGAR) on board the M55 Geophysica aircraft. The instrument was developed by the University of Frankfurt and measures the long-lived tracers CO2, N2O, CFC-12, CFC-11, H-1211, SF6, CH4 and H2 with two gas chromatographic channels and a CO2 sensor (LICOR). The measurements are supported by CO and O3 measurements of other instruments. Two campaigns were conducted to obtain measurements in the TTL: SCOUT-O3 (November/December 2005 in Darwin, Australia) and AMMA-SCOUT-O3 (August 2006 in Ouagadougou, Burkina Faso). After a general introduction of the thesis in chapters one and two, the third chapter describes the findings during this last campaign. Five local flights are analyzed to study the different transport processes that occur in the tropical tropopause layer above West-Africa: deep convection up to the level of main convective outflow, vertical mixing after overshooting of air in deep convection, horizontal inmixing from the extratropical lower stratosphere, and horizontal transport across the subtropical barrier. Main findings are that the TTL over West-Africa is mostly influenced by remote convection. The subtropical barrier is not a strong barrier but more a region of transition between the extratropical and the tropical stratosphere. Chapter 4 presents the results obtained during the SCOUT-O3 campaign. From the eight local flights the last four flights (051129, 051130a, 051130b, 051205) show enhanced values of ozone, CO and CO2 between 355 and 380 K potential temperature in comparison with the first four flights (051116, 051119, 051123, 051125). Horizontal inmixing from the extra-tropical stratosphere and influence of the local convective system Hector cannot explain the enhanced values of the two flights on 30 November Therefore, other possible explanations for these enhanced CO, CO2 and ozone levels are proposed. The first explanation is vertical mixing in the vicinity of the jet stream. However, the jet cannot explain the differences between the flights on 30 November and the flights on 29 November and 5 December. Another possible explanation is influence of polluted boundary layer air masses from the Indonesian region. Especially air sampled during the flights on November 30 crossed large parts of northern Indonesia between 8 and 10 days before the measurements. Convective uplift of biomass burning and other pollution plumes can transport CO and ozone precursors into the upper troposphere, where they can significantly enhance the ozone production. The last chapter deals with the vertical ascent rate in the TTL and uses measurements of both the SCOUT-O3 and AMMA-SCOUT-O3 campaign as well as data from previous aircraft campaigns (TROCCINOX and APE-THESEO). Time scales and residence times for mean vertical transport in the background TTL are estimated for different seasons and over different geographic regions using in situ observations of CO2 and long-lived tracers. The vertical transport time scales are constrained using the seasonal variation of CO2 in the tropical troposphere as a “tracer clock” for vertical ascent. Two methods are applied to calculate the residence time in the layer between 360 and 390 K potential temperature. The first method uses the slope of the CO2 index, the second method fits the CO2 index directly to the measurements assuming a constant ascent rate. The first method yields residence times for Australia,West Africa, and Brazil of the same order, 35-45 days to 380 K and 50 days to 390 K (where no value can be derived for Australia as the slope is changing approximately one month before the campaign). For APE-THESEO, the method does not yield reasonable results. The best estimates using the second method show moderate residence times between 360 and 390 K of 60±25 days SCOUT-O3 (NH autumn) and 43±8 days for AMMA/SCOUT-O3 (NH summer). These results agree well with the results calculated using the first method. For APE-THESEO and TROCCINOX the best fits yield shorter residence times of 23±7 and 40±10 days, respectively, both during winter. These results correspond well to the expectations based on the seasonal variation of the Brewer-Dobson circulation.
Physical soil properties feature high spatial variabilities which are known to affect geophysical measurements. However, these variations are not considered in most cases. The challenging task is to quantify the influence of soil heterogeneities on geophysical data. This question is analysed for DC resistivity and GPR measurements which are frequently used for near-surface explorations. To determine the pattern of electric soil properties in situ with the required high spatial resolution, geophysical measuring techniques are methodically enhanced. High-resolution dipole-dipole resistivity measurements are used to determine the electric conductivity distribution of the topsoil. Due to the small electrode separations, the actual electrode geometry has to be considered and an analytic expression for geometric factors is derived instead of assuming point electrodes. Two methods are used to determine soil permittivity with GPR:(i) the coefficient of reflection at the interface air-soil is measured with an air-launched horn antenna, (ii) the velocity of the groundwave is measured with a new setup using two receiver antennas enhancing the lateral resolution from in the best case 0.5 m for standard techniques to approximately 0.1 m with the new technique. With the optimised measuring techniques, the electric properties of sandy soils are determined in the field. Conductivity and permittivity show high spatial variability with correlation lengths of a few decimetres. Geostatistical simulation techniques are used to generate synthetic random media featuring the same statistical properties as in the field. FD calculations are carried out with this media to provide realistic synthetic data of resistivity and GPR measurements. Conductivity variations as determined in the field generate significant variations of simulated Schlumberger sounding curves resulting in uncertainties of the inverted models. Even in pedologically homogeneous sandy soil, moisture pattern and resulting permittivity variations cause strong GPR diffractions as demonstated by FD calculations. This influences the detectability of small objects such as e.g. landmines or of large reflectors as e.g. the groundwater table. Conductivity variations as typical for soils showed to have a minor effect on GPR measurements than variations of permittivity. In summary, geostatistical analysis and simulation provide a powerful tool to simulate geophysical measurements under field conditions including soil heterogeneity which can be used to quantify the uncertainty of field measurements by geologic noise.
The thesis is devoted to the study of the Antarctic polar vortex, mainly by analyzing data collected during APE-GAIA (1999) and ASHOE (1994) campaigns and recorded by the ADEOS satellite (1996-1997), and to improvement of the chromato-graphic processing schemes. A general introduction and overview of the campaigns and instruments relevant to the present work are given in Chapters 1 and 2. A relatively large part of the thesis (Chapters 3-5) is on improvement of the analysis of raw chromatographic data recorded during in-flight measurements of the trace gases. A Gaussian non-straight-base-line method, i.e. the Gaussian processing scheme (Chapter 3), is developed for better evaluation of the chromatographic peak size. Furthermore, a statistical cross-correlation method (Chapter 5) based on statistical behaviour of the whole chromatogram series fNchrg recorded, e.g., during a research flight or laboratory calibration, is developed and applied to measure the low-concentration trace gases. As demonstrated for HAGAR's chromatograms (HAGAR - High Altitude Gas Analyzer), the combination of the Gaussian fitting scheme for individual chromatograms and the statistical cross-correlation method for a series of subsequent chromatograms considerably improves and stabilizes quantitative analysis of in-flight chromatographic data. In this case, the detection accuracy of weak and noisy chromatographic signals can be improved by up to 40 %. A particular attention is paid to the in-flight two-standard calibration method. For this method, a special procedure, that allows to evaluate and effectively remove a weak background chromatographic signal associated with residual molecules in the carrier gas N2, is proposed and coded (Chapter 4). The developed approaches and methods are completely automized and, therefore, can be used for processing of in-flight chromatograms of recent and future field campaigns. The main part of the thesis (Chapters 6-8) deals with a two-dimensional quasi-Lagrangian coordinate system ... , based on a long-lived stratospheric trace gas i, and its systematic use for i = N2O in order to describe the structure of a well-developed Antarctic polar vortex, linearization and compactization of the tracer-tracer correlations in the polar vortex core (i.e. the stratospheric dynamics in this area), and the differential ozone losses in the Antarctic polar vortex area. In the coordinate system ... (...-method, Chapter 6), which refers to a well-developed polar vortex, the mixing ratio Âi is the vertical coordinate and ... = .... i is the reference profile in the vortex core) is the meridional coordinate. The quasi-Lagrangian coordinates ... are much more long-lived comparing with the standard quasi-isentropic coordinates, potential temperature ... and equivalent latitude ..e, do not require explicit reference to geographic space, and therefore well-suited for studying the dynamics of the Antarctic polar vortex and the relevant ozone loss processes. By using the introduced coordinate system ... to analyze the well-developed Antarctic vortex investigated in the APE-GAIA campaign, it is shown, in concurrence with the conclusion of A. M. Lee et al. (2001), that the Antarctic vortex area can be described in terms of the well-mixed and well-isolated vortex core, relatively wide vortex boundary region and adjoining surf zone. In this case, the reference profile ... i , which is compact in a well-developed and isolated polar vortex core [J. B. Greenblatt et al. (2002)], can be found by combining airborne (and/or balloon) data with high-altitude satellite measurements. A criterion, which uses the local in-situ measurements of Âi = Âi(£) and attributes the inner vortex edge to a rapid change (±-step) in the meridional pro¯le of the mixing ratio..., is developed in Chapter 6 to determine the (Antarctic) inner vortex edge. In turn, the outer vortex edge of a well-developed Antarctic vortex is proposed to attribute to the position of a local maximum of ...H2O in the polar vortex area. For a well-developed Antarctic vortex, the ...-parametrization of tracer-tracer correlations allows to distinguish the tracer-tracer inter-relationships in the vortex core, vortex boundary region and surf zone (Chapter 7). This is clearly illustrated by analyzing the tracer-tracer relationships Âi ¡ ÂN2O obtained from the in-situ data of the APE-GAIA campaign for i = CFCl3 (CFC-11), CF2Cl2 (CFC-12), CBrClF2 (H-1211) and SF6. The solitary anomalous points in the ...CFC11 ¡ ÂN2O correlation, observed in the Antarctic vortex core during the APE-GAIA and ASHOE campaigns, are interpreted in terms of small-scale localized differential descent. As detailed in Chapter 8, the quasi-Lagrangian coordinate system fÂN2O; ¢ÂN2Og is an effective tool for evaluation of the differential ozone losses in the polar vortex area. With this purpose, a two-parametric reference function ...O3 = F(...), which characterizes the unperturbed O3 distribution in the early winter polar vortex area, is introduced to separate and quantify in terms of the meridional coordinate ...2O the differential ozone losses in the vortex core and vortex boundary region. The method is applied to analyze the ozone depletion in the Antarctic stratosphere during the austral spring 1999 (APE-GAIA campaign). In Chapter 9, the main results of the thesis are summarized.
In this study, I investigate the crustal and upper mantle velocity structure beneath the Rwenzori Mountains in western Uganda. This mountain range is situated within the western branch of the East African Rift and reaches altitudes of more than 5000 m. I use four different approaches that belong to the travel-time tomography method. The first approach is based on the isotropic tomographic inversion of local data, which contain information about 2053 earthquakes recorded by a network of up to 35 stations covering an area of 140×90 km2. The LOTOS-09 algorithm described here is used to realize this approach. The second approach is based on the anisotropic tomographic inversion of the same local dataset. This method employs the tomographic code ANITA, developed with my participation, which provides 3D anisotropic P and isotropic S velocity distributions based on P and S travel-times from local seismicity. For the P anisotropic model, four parameters for each parameterization cell are determined. This represents an orthorhombic anisotropy with one vertically-oriented predefined direction. Three of the parameters describe slowness variations along three horizontal orientations with azimuths of 0°, 60° and 120°, and one is a perturbation along the vertical axis. The third approach is based on tomographic inversion of the teleseismic data, which contain information about the traveltimes of P-waves coming from 284 teleseismic events recorded by the seismic network stations. The TELELOTOS code, which is my own modification of the LOTOS-09 algorithm, is used in this approach. The TELELOTOS code is designed to iteratively invert the local and/or teleseismic datasets. Finally, I present the results of the new tomographic approach, which is based on the simultaneous inversion of the joint local and teleseismic data. The simultaneous use of these datasets for the tomographic inversion has several advantages. In this case, the velocity structure in the study area can be resolved as deep as in the teleseismic approach. At the same time, in the upper part of the study volume, the resolution of the obtained models is as good as in the local tomography. The TELELOTOS algorithm is used to perform the joint tomographic inversion. Special attention is paid in this work to synthetic testing. A number of different synthetic and real data tests are performed to estimate the resolution ability and robustness of the obtained models. In particular, synthetic tests have shown that the results of the anisotropic tomographic inversion of the local data have to be considered as unsatisfactory. For all approaches used in this study, I present synthetic models that reproduce the same pattern of anomalies as that obtained by inverting the real data. These models are used to interpret the results and estimate the real amplitudes of the obtained anomalies. The obtained models exhibit a relatively strong negative P anomaly (up to -10%) beneath the Rwenzori Mountains. Low velocities are found in the northeastern part of the array at shallower depths and are most likely related to sedimentary deposits, while higher velocities are found beneath the eastern rift shoulder and are thought to be related to old cratonic crust. The presence of low velocities in the northwestern part of the array may be caused by a magmatic intrusion beneath the Buranga hot springs. Relatively low velocities were observed within the lower crust and upper mantle in the western and southern parts of the study area (beneath the rift valley and the entire length of the Rwenzori range). The higher amplitude of the low-velocity anomaly in the south can be related to the thinner lithosphere in the southern part of the Albertine rift. In the center of the study area, a small negative anomaly is observed, with the intensity increasing with depth. This anomaly is presumably related to a fluids rising up from a plume branch in the deeper part of the mantle. According to the interpretation of the local earthquake distribution, the Rwenzori Mountains are located between two rift valleys with flanks marked by normal faults. The Rwenzori block is bounded by thrust faults that are probably due to compression.
Melt segregation inside the earth consists of two different processes: 1) Generation of partially molten rock and 2) separation of melt, produced from partially molten rock, from the solid residual matrix. This thesis focuses on the later process. The 2 phase flow dynamics combines the study of flow dynamics of melt and matrix. Several studies have given the background theoretical frameworks for the flow dynamics of melt inside the earth. [McKenzie, 1984] summarizes the studies of [Ahern and Turcotte, 1979; Frank, 1968; Sleep, 1975] and gives a complete set of governing equations for the 2-phase flow problem.
[Bercovici et al., 2001] gives a general formulation considering the univariate system of equations related to matrix and melt flow which includes the interfacial surface force. The assumption of melt having negligible viscosity compare to the matrix has been abandoned. Therefore, based on these formulations, we have constructed our numerical model and thereafter a fortran code PERCOL2D to get an insight of melt percolation process through porous media. Additionally, we have used the Helmhotz decomposition, which splits a smooth and rapidly decaying vector field into an irrotational vector field and an incompressible vector field [Srámek, 2007], for matrix and fluid viscosity in order to lower the number of linearly independent variables to minimize the computational complications. The melt residing at inter-granular areas of lithosphere, forms an interconnected network even at low porosity. Therefore, being less dense than the matrix, melt moves up through porous media due to its buoyancy. Compaction of matrix, which occurs to compensate the melt separation, is considered in this thesis, where the effective bulk and shear viscosity of matrix are function of melt fraction. We have effective bulk viscosity of matrix as inversely proportional to melt fraction. Porosity dependence of effective bulk and shear viscosity leads to stronger melt focusing in highly porous region like mid ocean ridges [Katz, 2008] since the ratio of bulk and shear viscosity is smaller (< 10) than the constant viscosity case for the porous waves having non dimensional amplitude 5% or higher. Moreover, it is observed in [Richard et al., 2012] that the solitary wave formed in porosity dependent viscous matrix settings are steeper than the one formed in the constant matrix viscosity setting.
Firstly some 1D numerical experiments with PERCOL2D have been carried out using fixed and periodic boundary conditions for zero source term (i.e. no melting or no freezing) and negligible surface tension.
3 series of model setups with different initial conditions have been carried out varying the width, non-dimensional amplitude and the background porosity value of the initial input of porous wave.
A mathematical derivation for 1D solitary wave solution for the two phase flow through porosity dependent compacting media, is obtained in this thesis which is different than the study of [Barcilon and Lovera, 1989; Barcilon and Richter, 1986; Scott and Stevenson, 1984; Spiegelman, 1993a,b] as the effective viscosity of matrix is constant there.
Although [Simpson and Spiegelman, 2011] gives the solitary wave solutions in 1D, 2D and 3D considering the porosity dependent effective viscosity of the matrix, but using the small background porosity approximation, they neglect the background porosity (i.e φ0) and therefore the effect of variation of compaction lengths, which causes variation in the shape and dynamics of the solitary wave. Therefore, the study [This thesis, Richard et al., 2012] can be used for more general purpose. Solitary waves in varying viscous medium, are steeper (cf fig.5.1) compared to the one in constant viscous medium and their speed decreases as an inverse function of the background porosity. Additionally, this analytical solution is used in our code PERCOL2D and also in FDCON for numerical benchmarking (1D) of PERCOL2D.
The role of melt grain contiguity is considered in the revised viscosity formulation [Schmeling et al., 2012] based on elastic moduli theory of a fluid filled poro-elastic medium. This formulation is used in this thesis to produce a comparative dispersion relationship between speed of the wave and the non dimensional amplitude of porous wave, based on both the viscosity formulations (fig. 6.20) where one can see that the model based on [Bercovici et al., 2001] formulation, converges to the same dispersion relationship obtained from [Simpson and Spiegelman, 2011]. Whereas, the dispersion relationship using [Schmeling et al., 2012] formulations, shows time-dependent decrease of phase velocity with increasing amplitude and it is not yet clear that whether these solutions converge to steady state porosity waves before the porosity becomes 1.