Refine
Year of publication
Document Type
- Doctoral Thesis (44) (remove)
Language
- English (34)
- German (9)
- Multiple languages (1)
Has Fulltext
- yes (44)
Is part of the Bibliography
- no (44) (remove)
Keywords
- Adaptation (1)
- Armenien (1)
- Biomineralization (1)
- COVID-19 lockdown (1)
- Carbonate (1)
- Carbonated apatite (1)
- Chemistry (1)
- Climate Change (1)
- Climate System, Climate Models, Information theory (1)
- Clumped isotopes (1)
Institute
- Geowissenschaften / Geographie (44) (remove)
Floodplains and other wetlands depend on seasonal river flooding and play an important role in the terrestrial water cycle. They influence evapotranspiration, water storage and river discharge dynamics, and they are the habitat of a large number of animals and plants. Thus, to assess the Earth’s system and its changes, a robust understanding of the dynamics of floodplain wetlands including inundated areas, water storages, and water flows is required.
This PhD thesis aims at improving the modeling of large floodplains and wetlands within the global-scale hydrological model WaterGAP, in order to better estimate water flows and water storage variations in different storage compartments. Within the scope of this thesis, I have developed a new approach to simulate dynamic floodplain inundation on a global-scale. This approach introduces an algorithm into WaterGAP, which has a spatial resolution of 0.5 degree (longitude and latitude) globally. The new approach uses subgrid-scale topography, based on high-resolution digital elevation models, to describe the floodplain elevation profile within each grid cell by applying a hypsographic curve. The approach comprises the modeling of a two-way river-floodplain interaction, the separate downstream water transport within the river and the floodplain – both with temporally and spatially different variable flow velocities – and the floodplain-groundwater interactions. The WaterGAP version that includes the floodplain algorithm, WaterGAP 2.2b_fpl, estimates floodplain and river water storage, inundated area and water table elevation, and also simulates backwater effects.
WaterGAP 2.2b_fpl was applied to model river discharge, river flow velocity, water storages, water heights and surface water extent on a global-scale. Model results were comprehensively validated against ground observations and remote sensing data. Overall, the modeled and observed data are in agreement. In comparison to the former version WaterGAP 2.2b, the model performance has improved significantly. The improvements are most remarkable in the Amazon River basin. However, the seasonal variation of surface water extent and total water storage anomalies are still too low in many regions on the globe when compared to observations. A detailed analysis of the simulated results suggests that in the Amazon River basin the introduction of backwater effects is important for realistically simulating water storages and surface water extent. Future efforts should focus on the simulation of water levels in order to better model the flow routing according to water slope. To further improve the model performance in specific regions, I recommend that the globally constant model parameters that affect inundation initiation, river-floodplain interaction, DEM correction for vegetation, and backwater amount at basin or subbasin-scale be adjusted.
Air-sea feedbacks between the Mediterranean Sea and the atmosphere on various temporal and spatial scales play a major role in the Mediterranean regional climate system and beyond. The Mediterranean Sea is a source of moisture due to excess evaporation and, on a long-term average, is associated with a warming of the lower atmosphere in contact with the sea surface due to heat loss at the air-sea interface. The complex air-sea interactions and feedbacks in the Mediterranean basin strongly modulate the sea surface fluxes and favor several cyclogenetic activities under certain meteorological conditions. Examples of such cyclonic activities are medicanes (Mediterranean hurricanes) and Vb-cyclones. Medicanes are mesoscale, marine, and warm-core Mediterranean cyclones that exhibit some similarities to tropical cyclones, while Vb-cyclones are extra-tropical cyclones, that propagate from the Western Mediterranean Sea and travel across the Eastern European Alps into the Central European region. Extremely strong winds and heavy precipitation associated with these cyclones can lead to severe destruction and flooding. Changes in the intensity and frequency of these cyclones are also projected under changing future climate conditions, where the Mediterranean region has been identified as a hotspot in terms of rising temperatures.
The development of high-resolution regional climate models (RCMs) has progressed our understanding of the processes characterizing the Mediterranean climate. However, large uncertainties still exist regarding the estimates of air-sea fluxes, which, in turn, affect the simulation of the Mediterranean climate. Several factors can be attributed to such discrepancies, such as data quality, temporal and spatial resolution, and the misrepresentation of physical processes. To overcome some of these inconsistencies and deficiencies of the existing climate simulations, a new high-resolution atmosphere-ocean regional coupled model (AORCM) has been developed to simulate the air-sea feedback mechanisms. This coupled model incorporates the coupling of RCM COSMO-CLM (CCLM) and the regional ocean model NEMO-MED12 for the Mediterranean Sea (MED) as well as NEMO-NORDIC for the North- and Baltic Sea (NORDIC). Several experiments were performed using both the coupled and uncoupled models to investigate the impact of air-sea interactions and feedbacks on sea surface heat fluxes, wind speed, and on the formation of Mediterranean cyclones (i.e., medicanes and Vb-cyclones). These experiments were performed using different horizontal atmospheric grid resolutions to analyze the effect of resolution on sea surface heat fluxes, wind speed, and the development of medicanes.
The results of the present study indicate that a finer atmospheric grid resolution ([is as appreciated as]9 vs. [is as appreciated as]50 km) improved the wind speed simulations (particularly near coastal areas) and subsequently improved the simulations of the turbulent heat fluxes. Both parameters were better simulated in the coupled simulations than in the uncoupled simulations, but coupling introduced a warm SST bias in winter. Radiation fluxes were slightly better represented in coarse-grid simulations than in fine-grid simulations. However, the higher-resolution coupled model could reproduce the observed net outgoing total surface heat flux over the Mediterranean Sea. In addition to that sub diurnal SST variations have a strong effect on sub-daily heat fluxes and wind speed but minor effects at longer timescales. Regarding the impact of atmospheric grid resolution ([is as appreciated as]50, 25, and [is as appreciated as]9 km) and ocean coupling on medicanes, it was detected that the coupled model with a finer atmospheric grid ([is as appreciated as]9 km) was able to not only reproduce most medicane events, but also improved the track length, warm core, and wind speed compared to the uncoupled model. The coupled model with the coarse-grid ([is as appreciated as]50 and [is as appreciated as]25 km) did not show any improvement in simulating medicanes compared to the uncoupled model. The spectral nudging technique, applied on the wind components above 850 hPa in the interior domain to keep large-scale circulation close to the driving data (i.e., ERAInterim reanalysis), improved the accuracy of the times and locations of generated medicanes, but no improvement was found in the track length and intensity.
Concerning the role of the Mediterranean Sea coupling on Vb cyclones, the investigation showed that atmosphere-ocean coupling had an overall positive impact, although with a strong case-by-case variation, on the trajectories and intensity of Vb-cyclones as a result of the variation in moisture source for each event. In general, all model configurations could replicate Vbcyclones, their trajectories, and associated precipitation fields. The average structure of the precipitation field was best represented in the coupled simulations. Coupling of the North- and Baltic Seas also showed an improvement in some of the simulated Vb-cyclones.
The atmosphere-ocean coupling showed an overall positive impact on the simulation of sea surface heat fluxes and Mediterranean cyclones (medicanes and Vb-cyclones). Moreover, the representation of sea surface heat fluxes, wind speed, and medicane features was more realistic when using a finer atmospheric grid resolution (less than 10 km). The present study suggests that the combination of a finer atmospheric grid resolution together with atmosphere-ocean coupling is advantageous in simulating the Mediterranean climate system.
One of the most important events in human history occurred during the Early Pleistocene: the dispersal of early hominins out of Africa and into Europe and Asia. In Western Europe, the earliest evidences of the genus Homo have been found in the Baza Basin, at the sites of Orce in the SE of the Iberian Peninsula. These sites contain fossils and lithic industry dated approximately as 1.4–1.3 Ma.While hominin remains and artifacts at Orce, as well as the accompanying fauna, have been extensively studied, the properties and evolution of the Early Pleistocene vegetation in the basin remain unknown. The general effect of climate change on the expansion of early hominins from Africa into Eurasia still remains unclear. It is not known if the Early Pleistocene climate changes and the development of glacials periods led to the extirpation of European communities, or if those communities were able to endure and persist through such adverse climatic periods. This open question highlights the need for climate and environmental analyses for the time before, during and after the first presence of Homo in Europe. This PhD thesis contributes to that need by the presentation of the first long pollen record of the Baza Basin, where the oldest hominin sites in Western Europe are found.
During my PhD, I was applying the clumped isotope technique to modern brachiopods and fossil belemnites, and I conducted methodological work. Carbonate clumped isotope thermometry is a tool to reconstruct carbonate precipitation temperatures. In contrast to oxygen isotope thermometry, i.e., the δ18O-thermometer, the carbonate clumped isotope thermometer does not require an estimate for the oxygen isotope composition of the seawater, as it considers the fractionation of isotopes exclusively amongst carbonate isotopologues. The ∆47 value of a carbonate expresses the abundance of the 13C–18O bond bearing carbonate isotopologue, within the carbonate, relative to its random distribution. In thermodynamic equilibrium, the ∆47 value of a given carbonate is solely a function of the carbonate precipitation temperature. However, kinetic isotope fractionations, i.e., vital effects, driven by diffusion, pH or incomplete oxygen isotope exchange between water and dissolved inorganic carbonate species can cause the carbonate to be precipitated with isotopic compositions that are offset from those predicted for thermodynamic equilibrium.
Brachiopods serve as important geochemical archives of past climate conditions. To investigate the nature and significance of kinetic controls on brachiopod shell δ18O and ∆47 values, in collaboration with the BASE-LiNE Earth ITN, I analysed the bulk and clumped isotope compositions of eighteen modern brachiopod shells, collected from different geographic locations and water depths that cover a substantial range of growth temperatures. Growth temperatures and seawater δ18O values for each brachiopod were independently determined. Most of the analysed brachiopods exhibit combined offsets from clumped and oxygen isotope equilibrium, and there is a significant negative correlation between the offset values. The observed correlation slope between offset ∆47 and offset δ18O point to the importance of kinetic effects associated with Knudsen diffusion and incomplete hydration and hydroxylation of CO2 (aq), occurring during biomineralisation. The correlations between the growth rates of the analysed brachiopods and both the offset ∆47 and the offset δ18O values provide further arguments for the presence of kinetic effects. In conclusion, the oxygen and clumped isotope composition of modern brachiopod shells are affected by growth rate-induced kinetic effects that hinder their use for palaeoceanography.
Diamant hat besondere physikalische und optische Eigenschaften sowie eine starke Resistenz gegenüber Strahlenschädigung. Diese Eigenschaften ermöglichen eine vielfältige Anwendung von Diamant in Wissenschaft und Technik, wie zum Beispiel als Sensormaterial in Strahlungsdetektoren.
Kubisches Zirconiumdioxid (ZrO2) wird aufgrund seiner mechanisch und optisch ähnlichen Eigenschaften unter anderem an Stelle von Diamant eingesetzt. Es ist ebenfalls ein geeignetes Material für viele technische Anwendungen und wird durch seine Strahlenresistenz in Strahlungsumgebungen verwendet. Da beide Materialien in diesem Anwendungsbereich hoher energetischer Strahlung ausgesetzt sind, sind Reaktionen auf die Bestrahlung wie etwa strukturelle Veränderungen oder die Änderungen von Materialeigenschaften von großem Interesse.
In der vorliegenden Arbeit wurde die Morphologie, Struktur und physikalischen Eigenschaften von Diamant und Yttriumoxid-stabilisiertem kubischem ZrO2 nach der Bestrahlung mit 14 MeV Au-Ionen und 1.6 GeV Au-Ionen untersucht. Die durch die Bestrahlung verursachten Veränderungen der Oberflächen und der bestrahlten Volumina wurden mit diversen komplementären analytischen Methoden charakterisiert, bewertet und für die verschiedenen Materialien und Ionenenergien verglichen.
Mittels Röntgenfluoreszenzmessungen wurde die Verteilung und Menge an implantiertem Au semi-quantitativ ermittelt. Die Oberflächen der Proben wurden mit optischer Mikroskopie, Rasterkraftmikroskopie, Rasterelektronenmikroskopie, Röntgenreflektometrie und Elektronenrückstreubeugung untersucht. Strukturelle Veränderungen wurden mit Raman-Spektroskopie analysiert. Der elektrische Widerstand, die Dichte, die Härte sowie das Ätzverhalten der bestrahlten Proben wurden ermittelt und geben Auskunft über die Änderung physikalischer Eigenschaften der Materialien.
Diamant und kubisches ZrO2 reagieren sehr unterschiedlich auf die Bestrahlung mit Au-Ionen gleicher Energien und Fluenzen. Die Diamantproben zeigen nach der Bestrahlung mit 14 MeV Au-Ionen deutliche Veränderungen und Schädigungen der Oberfläche sowie des bestrahlten Volumens. Es wird eine Änderung der Struktur, der Dichte, der Härte, des elektrischen Widerstands sowie des Ätzverhaltens der Proben beobachtet, was auf die Amorphisierung von Diamant zurückgeführt wird. Kubisches ZrO2 ist deutlich strahlungsresistenter gegenüber der Bestrahlung mit 14 MeV Au-Ionen. Es werden keine signifikanten strukturellen Änderungen im getesteten Fluenzbereich beobachtet.
Die mit 1.6 GeV Au-Ionen bestrahlten Diamanten zeigen nur geringe Schädigungen und keine deutliche Änderung der Struktur oder der physikalischen Eigenschaften. Die kubischen ZrO2 Proben sind als Folge der Bestrahlung mit 1.6 GeV Au-Ionen zerbrochen, was auf hohe interne Spannung durch Defektbildung zurückgeführt wird.
This thesis is focusing on the impact of Paratethys and Mediterranean water bodies over the Eurasian climate and the interplay between climate, tectonics and biosphere during the late Miocene. This target was the interval between 12.7 and 7.65 Ma for Paratethys, following the Eastern Paratethys restriction and isolation, and 7.2−6.5 Ma (the early Messinian) in Mediterranean, zooming on the effects of gateway restrictions over the eastern Mediterranean and the new born Aegean domain. In both cases restriction is overlapping with large scale climatic changes and tectonic reconfiguration, leading a sort of symbiotic relationship.
Paratethys was a giant epicontinental sea that covered a large part of Eurasia since Paleogene. Due to the Eurasia-Afro-Arabia collision and formation of the Alpine-Himalayan belt (Rögl, 1999; Popov et al., 2006), the Paratethys was divided during the late Miocene in smaller basins that in time were isolated of each other. The protracted isolation and intense continentalisation of paratethyan realm led to changes in humidity distribution, basin connectivity, sediment sources and salinity. These changes had in turn major consequences over water circulation, water availability, vegetation cover and biota. These changes are more intense after 11.6 Ma, when the Eastern Paratethys lost any sustained marine connection, evolving into an enclosed system with endemic fauna (Harzhauser and Piller, 2007).
Mediterranean Sea is a Mezozoic oceanic relic squeezed between Africa, Europe, Anatolia and Arabia, as Africa continued to subduct beneath the European plate. As opposed to Paratethys, it maintained the open connection with the ocean until Messinian, when the two Atlantic gateways (Betic and Rifian corridors) closed for a short time, isolating the basin. The cut off resulted in a dramatic drop down and onset of evaporitic precipitation in marginal basins, the event receiving the name of Messinian Salinity Crisis (5.97−5.55). The restriction affected all marine ecosystems, due to changes in salinity and stratification of water column.
The main objectives of this thesis were:
(1) build valid paleo-temperatures records for both basins based on biomarkers;
(2) reconstruct the hydrology for the late Miocene time interval;
(3) identify vegetation composition and changes;
(4) identify paleo-fires in the late Miocene sediment records;
(5) identify the biotic response to the overall climate and tectonic changes.
All the above objectives were attained with results published in specific journals (Chapters 5−7).
Based on Panagia section (Taman Peninsula, Russia) the longest Paratethys temperature record was completed (~5 Myr), covering the interval between 12.7 and 7.65 Ma. A comprehensive SST and MAT records was obtained, as well as soil pH and carbon (δ13C) and hydrogen (δ2H) stable isotopic compositions on n-alkanes and alkenones. The main findings are concentrated around three prolonged periods with severe droughts affecting the late Miocene circum-Paratethys region peaking at 9.65, 9.4 and 7.9 Ma, associated with a transition towards open land vegetation, intensification of fire activity and enhanced evaporation and aridity.
The time intervals with dryer conditions recorded in Panagia coincide with periods of mammal turnover and dispersal in Eurasia indicating that major environmental changes occurred in the circum-Paratethys region and Paratethys fragmentation had a great impact on the terrestrial ecosystems, when periods of prolonged droughts generated biotic crises and animal displacements across the Eurasian continent. The δ13CC29n-alkane values and charcoal morphologies from Panagia indicate an increased contribution of C4 plants adapted to drier conditions at 9.66 Ma. Similarly high δ13CC29n-alkane values continue until 9.4 Ma, when in Western Europe increased seasonality accelerated the demise of the evergreen subtropical woodlands and expansion of grasslands from Anatolia and Middle East to Europe.
As a result of basin fragmentation and climatic stress, the Eastern Paratethys sub-basins progressively lost their marine properties and turned into brackish-fresh water bodies fed primarily by riverine input. The shallower areas became in time emerged, obstructing connections and isolating the biota, inducing rapid adjusting or extinctions. Thus, the Paratethys harbored a highly endemic fauna (Rögl, 1999), such as dwarf whales, dolphins, seals (among mammals), as well as fish and other taxa (mollusks, ostracods, diatoms, foraminifera, algae, etc.).
Collectively the data structured and analyzed in chapter five support a model in which the Eastern Paratethys evolved as a largely (en)closed system, registering paleoenvironmental signals that are governed by interbasinal connectivity (or lack of it) and regional climate changes in the basin catchment. Acting as an important source of humidity for Western and Central Asia, the size and areal extent of the Paratethys water body is likely to have had a major impact on hydroclimate patterns in the Eurasian interior, with the cumulative fluctuations in both hydrology and surface temperature enhancing the aridity and seasonality, with different partition of moisture over the year. Our combined data suggests a decoupling of Paratethys from the global system as isolation advanced, dominated by regional tectonics and ultimately the Paratethys volume and areal extent reduction.
Atmospheric particles play an important role in the radiative balance of the Earth, as well as they affect human health and air quality. Hence, the chemical characterization constitutes a crucial task to determinate their properties, sources and fate. Particularly, the analysis of nanoparticles (d<100 nm) represents an analytical challenge, since these particles are abundant in number but have very little mass.
This accumulative thesis focuses on the chemical characterization of nanoparticles, performed in both laboratory and field studies. Here, I present four manuscripts, two of which are my main project as a lead author.
The first manuscript (Caudillo et al., 2021) focuses on the gas and the particle phase originated from biogenic precursor gases (α-pinene and isoprene). The experiments were performed in the CLOUD chamber at CERN to simulate pure biogenic new particle formation. Both gas and particle phases are measured with a nitrate CI-APi-TOF mass spectrometer, while the TD-DMA is coupled to it for particle-phase measurements, this setup allows a direct comparison as both measurements use the identical chemical ionization and detector. This study demonstrates the suitability of the TD-DMA for measuring newly formed nanoparticles and it confirms that isoprene suppresses new particle formation but contributes to the growth of newly formed particles.
The second manuscript (Caudillo et al., 2022) presents an intercomparison of four different techniques (including the TD-DMA) for measuring the chemical composition of SOA nanoparticles. The measurements were conducted in the CLOUD chamber. The intercomparison was done by contrasting the observed chemical composition, the calculated volatility, and the thermal desorption behavior (for the thermal desorption techniques). The methods generally agreed on the most important compounds that are found in the nanoparticles. However, they did see different parts of the organic spectrum. Potential explanations for these differences are suggested.
The third manuscript (Ungeheuer al., 2022) presents both laboratory and ambient measurements to investigate the ability of lubricant oil to form new particles. These new particles are an important source of ultrafine particles in the areas nearby large airports. The ambient measurements were performed downwind of Frankfurt International Airport, and it was found that the fraction of lubricant oil is largest in the smallest particles. In the laboratory, the main finding was that evaporated lubricant oil nucleates and forms new particles rapidly. The results suggest that nucleation of lubricant oil and subsequent particle growth can occur in the cooling exhaust plumes of aircraft-turbofans.
The fourth manuscript (Wang et al., 2022) is a new particle formation study in the CLOUD chamber at CERN. This study shows that nitric acid, sulfuric acid, and ammonia interact synergistically and rapidly form particles under upper free tropospheric conditions. These particles can grow by condensation (driven by the availability of ammonia) up to CCN sizes and INP particles. The ability of these particles to act as a CCN and INP was also investigated and it was found to be as efficient as for desert dust. This mechanism constitutes an important finding and it can account for previous observations of high concentrations of ammonia and ammonium nitrate over the Asia monsoon region.
First-principles modeling techniques offer the ability to simulate a wide range of systems under different physical conditions, such as temperature, pressure, and composition, without relying on empirical knowledge. Density functional theory (DFT), a quantum mechanical method, has become an exceptionally successful framework for materials science modeling. Employing DFT makes it possible to gain valuable insights into the fundamental state of a system, enabling the reliable determination of equilibrium crystal structures. Over time, DFT has become an essential tool that can be incorporated into various schemes for predicting the properties of a material related to its structure, insulating/metallic behavior, magnetism, and optics. DFT is regularly applied in numerous fields, spanning from fundamental subjects in condensed matter physics to the study of large-scale phenomena in geosciences. In the latter, the effectiveness of DFT stems from its ability to simulate the properties found on the Earth, other planets, and meteorites, which may pose challenges for their direct study or laboratory investigation.
In this thesis, a comprehensive examination of a family of monosulfides and a perovskite heterostructure was conducted. These materials are relevant for their potential applications in technology, energy harvesting, and in the case of monosulfides, their speculated abundance on the planet Mercury.
Firstly, a DFT approach was used to analyze two non-magnetic monosulfides, CaS and MgS. We determined their structural properties and then focused on the modeling of their reflectivity in the infrared region. The calculation of the reflectivity considered both harmonic and anharmonic contributions. In the harmonic limit, the non-analytic correction was employed to accurately determine the LO/TO splitting, which is necessary to delimit the retstrahlend band, that is, the maximum of the reflectivity. The anharmonic effects given by up to three-phonon and isotopic scatterings, which were included using perturbation theory, primarily smeared the reflectivity spectra edges in the high-wave region.
Secondly, four polymorphs of MnS were studied using a combination of first-principles methods to simulate their antiferromagnetic (AFM) and paramagnetic (PM) states. The integration of DFT+$U$ with special quasirandom structures (SQS) supercells, and occupation matrix control techniques was crucial for achieving convergence, structural optimization accuracy, and obtaining finite energy band gaps and local magnetic moments in the PM phases. The addition of the Hubbard $U$ correction was necessary to treat the highly-correlated Mn $d$-electrons. The success of our approach was clear based on our electronic structure predictions for the PM rock-salt B1-MnS polymorph. Experimentally this phase has been observed to be an insulator, but multiple \emph{ab initio} works resulted previously in metallic behavior. Our computations, on the other hand, predicted insulating and magnetic properties that compare well with available measurements. Additionally, the pressure-field stability of the four MnS polymorphs was studied. In the case of the PM phases, B1-MnS was identified to be the most stable up to about 21 GPa, then transforming into the B31-MnS polymorph. This finding was in close agreement with high-pressure experiments reporting a similar phase transformation. The optical properties of B1-, B4-, and B31-MnS were also simulated. The SQS technique was used to obtain soft-mode-free phonon band structures within the harmonic approximation. Then, the anharmonic effects were included, and the reflectivity was calculated for B1-MnS and B4-MnS. In both cases, a good agreement for the LO/TO splitting was achieved in comparison to experimental results.
Lastly, the oxygen-deficient heterostructure of LaAlO$_{3-\delta}$ /SrTiO$_{3-\delta}$ was investigated also employing DFT+$U$, with a particular emphasis on the potential impact of vacancy clustering at the interface. Six distinct configurations of pairs of vacancies were studied and their energies were compared to find the most stable one. The orbital reconstruction of Ti orbitals was also examined based on their location with respect to the vacancies and the local magnetic moments were calculated. The final results showed that linearly arranged vacancies located opposite to Ti ions give the most energetically stable configuration.
Melting inside earth is a common phenomenon and can be observed in many different regions where melt travels through the mantle and crust to eventually reach the surface where it crystallizes to build large volcanic provinces, whole stratigraphic layers of flood basalts, or even the oceanic crust. Often, melt reaching the surface is a good source of information. It can be used to achieve a better understanding about processes taking place in deeper regions inside the mantle and it is therefore essential to fundamentally understand melting and melt percolation processes. In order to achieve a deeper understanding, the aim of this thesis is to investigate processes that are connected to melting by using numerical models.
The physical model used is a so called two-phase flow model which describes the ability of melt to percolate through a viscously deforming, partially molten matrix. A famous feature of two-phase flow are solitary porosity waves, which are waves of locally higher porosity ascending through a partially molten background, keeping its shape constant, driven by decompaction and compaction of the solid matrix in front and behind the wave.
The viscosity law for shear- and volume viscosity was strongly simplified in most previous studies that modeled solitary waves. Often the porosity dependency is underestimated or its influence on the volume viscosity is even neglected, leading to too high viscosities. In this work more realistic laws are used that strongly decrease for small melt fractions. Those laws are incorporated into a 2D Finite Difference mantle convection code with two-phase flow to study the ascent of solitary porosity waves.
The model results show that an initial Gaussian-shaped wave rapidly evolves into a solitary wave with a certain amplitude, traveling upwards with constant velocity. Even though strongly weaker viscosities are used, the effect on dispersion curves and wave shape are only minor as long as the background porosity is rather small. The results are still in agreement to semi-analytical solutions which neglect shear stresses in the melt segregation equation. Higher background porosities and wave amplitudes lead to significant decrease in phase velocity and wave width, as the viscosity is strongly effected. However, the models show that solitary waves are still a possible mechanism for more realistic matrix viscosities.
While the ascending of porosity waves are mostly described by the movement of fluid melt, partially molten regions inside Earth trigger upwelling of both, solid and fluid phases, which can be called diapirism. While diapirs can have a wide variety of wavelengths, porosity waves are restricted to a few times the compaction length. The size of a melt perturbation in terms of compaction length therefore describes whether material is transported by diapirism or porosity waves. In this thesis we study the transition from diapiric rise to solitary porosity waves by systematically changing the size of a porosity perturbation from 1.8 to 120 times the compaction length. In case of a perturbation of the size of a few times the compaction length a single porosity wave will emerge, either with a positive or negative vertical matrix flux and if melt is not allowed to move relative to the matrix a diapir will emerge. In between these physical end members a regime can be observed where the partially molten perturbation will split up into numerous solitary waves, whose phase velocity is low compared to the Stokes velocity and the swarm of solitary waves will ascend jointly as a diapir, slowly elongating due to a higher amplitude main solitary wave.
Solitary waves will always emerge from a melt perturbation as long as two-phase flow is enabled, but the time for a solitary wave to emerge increases non-linearly with the perturbation radius in terms of compaction length. In nature, in many cases this time might be too long for solitary waves to emerge.
Another important feature when it comes to two-phase flow is the transport of trace elements in melt. Incompatible elements prefer to go into the melt, which eventually enriches the area where it crystallizes again. In order to model this redistribution, the code FDCON was extended to allow for fully consistent transport of elements in melt, including melting, freezing and re-equilibration with time. A 2D model, a simple representation of a volcanic back arc, is set up to investigate the behavior of trace elements. The influence of retention number and re-equilibration time is examined. Lava-lamp like convection can be observed in the lower part of the model, producing melt, that eventually leads to enrichment in trace elements in the upper high-viscous layer. The total enrichment in this layer approaches an asymptotic value and a 0D model is introduced to recreate this behavior.
As a cognitively-mediated response, autonomous adaptation at farm-gate levels constitutes reactionary actions by farmers against climate impacts. These actions are shaped by interacting factors such as household characteristics, livelihood scope and resources. It is driven by the goal of adapting cultivated farmlands to climate and for sustaining crop yields. Thus, interest in balancing adaptation goals with protection of vegetation conditions is less of a priority. Lack of research interest in understanding the gap between objectives of reactionary adaptation and protection of surface conditions (vegetation canopies) is a gap in research. In many studies, farm-gate level adaptation is described as a set of zero-feedback actions in response to climate impacts. This perception conceals the stress and impact-engendering attribute of reactionary adaptation. Inspired towards addressing this conceptual gap; this study investigates impact of farmers’ reactionary adaptation on vegetation cover in Keffi, Nasarawa, Nigeria. A twenty-year time-series NDVI and rainfall datasets are linearly regressed to examine the extent of NDVI-rainfall sensitivity. A weak linear relationship between NDVI and rainfall in Keffi for the period, 1999-2018 is observed. At a regression slope of 0.001, R squared, R2=0.129 (implying that only about 13% of the variability in NDVI in Keffi are explained by rainfall amount) and a bivariate regression coefficient, r=0.359; statistical evidence shows that rainfall amount are not significant predictors of NDVI in Keffi. In investigating the possible interference of non-rainfall factors on vegetation productivity (NDVI) in Keffi; a residual trend (RESTREND) analysis was carried out. Regression of residuals from NDVI-Rainfall linear regression produced a R=0.192 with a negative and downwards slope. The downward character of the RESTREND slope is suggestive of non-rainfall factors contained in the residuals. In validating the RESTREND analysis, a comparative analysis between observed and predicted NDVI derived from a reference NDVI value of 0.46 was carried out. The NDVI value of 0.46, is empirically assumed to be average NDVI value expected at a minimum rainfall amount of 850mm/year reported in tropical Savanna ecosystems. Using this empirical relationship, NDVI values were predicted for Keffi. Even at higher rainfall amounts≈1340mm/year, amounts were unable to produce corresponding higher NDVI values; rather a more plausible correlation between reference-derived predicted NDVI values and rainfall was obtained. A further analysis with predicted NDVI values, based on 1999 NDVI value in Keffi returned higher NDVI units than observed NDVI values. This strengthens the attribution of the possible interference of rainfall-NDVI sensitivity by non-rainfall factors like human activities on vegetation productivity. Surface soil analysis to exclude potential impacts of soil nutrients and moisture deficiency on vegetation productivity, showed that soil had insignificant effect on vegetation dynamics. Further inferential analysis, using the inter-annual NDVI and the reclassified bi-decadal NDVI maps showed that spatial vegetation distribution in Keffi were driven by farmers inter-annual rotational cultivation footprints than rainfall variability. With a three-class categorization, “gain, loss and significant loss”, the spatial distribution of vegetation in Keffi between (1999-2008) and (2009-2018) was assessed. Temporal condition (stressed and healthy) across the three classes supports the attribution of farmers’ reactionary adaptation and cultivation practices on the dynamic spatial vegetation distribution. Between 1999 -2018, an increase in areas with significant vegetation loss (42%), so with a decrease of -25% in areas with healthy vegetation was observed. The character of vegetation cover across the two decadal time slices, reflects landuse intensity and unsustainable farming practices. Preferences for modification of cultivation practices and changes in seed by farmers exerts positive feedbacks on vegetation cover. Higher statistical measures, 38.4% (yearly cropping) and 44% (shifting cultivation with less fallow periods) were observed in the chi-square analysis. These measures were higher than 2.0% relating to shifting cultivation with more fallow periods. While 11.6% farmers noted cultural practices as reasons for preferred cultivation methods, 48.4% farmers attributed climate as reason behind cultivation modification. This was higher than 24.4% who linked issues of tenure rights to cultivation practices. With preferences for yield- breaching strategies, the non-receding cultivation and shorter fallow practices in Keffi triggers feedback on vegetation dynamics. Evidence from this study shows that the NDVI-rainfall functional sensitivity in Keffi is plausibly dampened by effects of reactionary farm-gate level adaptation practices.