Doctoral Thesis
Refine
Document Type
- Doctoral Thesis (8) (remove)
Language
- English (8) (remove)
Has Fulltext
- yes (8)
Is part of the Bibliography
- no (8)
Keywords
- Modellierung (8) (remove)
In terms of climate change and climate change mitigation, the quantitative knowledge of global carbon pools is important information. On the one hand, knowledge on the amount of carbon cycling among – and stored in – global pools (i.e. Atmosphere, Biosphere, Cryosphere, Hydrosphere, and Lithosphere) may improve the reliability of models predicting atmospheric CO2 concentrations in terms of fossil fuel combustion. On the other hand, the carbon sequestration potential of specific ecosystems allows for estimating their feasibility regarding carbon trade mechanisms such as the Clean Development Mechanism or the Reducing Emissions from Deforestation and Degradation Program (REDD+). However, up to date, the majority of terrestrial carbon assessments have focused on forests and peatlands, leaving a data gap open regarding the remaining ecosystems. This data gap is likely to be explained by the relatively high carbon densities and/or productivities of forests and peatlands. Nevertheless, to get a precise as possible global picture, information on carbon pools and sequestration of other ecosystems is needed. Although desert ecosystems generally express low carbon densities, they may absolutely store a remarkable amount of carbon due to their large areal extent. In this context, Central Asian Deserts (in particular within the Turanian Deserts, i.e. Karakum, Kysylkum, Muyunkum) likely inhibit comparably high carbon pools as they express a sparse vegetation cover due to an exceptionally high annual precipitation if compared to the World’s deserts. In this dissertation, three important woody plant species – Populus euphratica and Haloxylon aphyllum and Haloxylon persicum – of Central Asian Deserts were investigated for their carbon pools and carbon sequestration potential. These species were chosen as they I) locally express high carbon densities, II) are dominant species, III) have a rather large spatial distribution, and IV) have experienced a strong degradation throughout the 20th century. Thus, they likely show a remarkable potential for carbon re-sequestration through restoration and thus for an application of carbon trade mechanisms (CHAPTER I). P. euphratica was investigated in the nature reserve Kabakly at the Amu Darya, Turkmenistan and in Iminqak at the Tarim He, Xinjiang, China. The assessment of Haloxylon species was restricted to the Turanian deserts west of the Tain Shan. To achieve a first scientific basis for large scale estimates, different methodologies, ranging from allometric formulas, over dendrochronology to remote sensing were combined (CHAPTERS II-V). In CHAPTER II allometric formulas were successfully developed for Haloxylon aphyllum and Haloxylon persicum and applied to six study sites distributed over the Turanian Deserts to represent the allometric variability of Haloxylon species in Central Asia. CHAPTER III derives another allometric formula (only based on canopy area) for H. aphyllum and combines it with a remote sensing analysis from the nature reserve Repetek. Thereby, a first large scale estimate covering the Northeastern Karakum Desert of carbon pools related to mono specific H. aphyllum stands is achieved. CHAPTER IV describes the wood structure of Populus euphratica forests in the nature reserve Kabakly (Turkmenistan) and in Iminqak (Xinjiang, China). In CHAPTER V a dendrochronological approach derives models for predicting the Net Primary Productivity (NPP) and the age of P. euphratica in the nature reserve Kabakly. Thereby, a first feasibility assessment regarding remote sensing analyses and the upscaling of the obtained NPP results is carried out. First estimates based on these local studies (CHAPTER VI), reveal carbon densities ranging from 0.1 – 26.3 t C ha 1 for the three investigated species. Highest maximum and median carbon densities were found for P. euphratica, but Haloxylon aphyllum expressed remarkable maximum carbon densities (13.1 t C ha-1), too. The total carbon pools were estimated at 6480 kt C for P. euphratica, 520 kt C for H. aphyllum stands and 6900 kt C for Haloxylon persicum shrubland. Accounting for the extent of degraded areas, the total re-sequestration potentials of the respective species were estimated at 4320 kt C, 1620 kt C and 21900 kt C, this highlighting the remarkable absolute re-sequestration potential of H. persicum shrubland despite its low average carbon densities. In the end, the main results were put into a broader context (CHAPTER VI), discussing the general feasibility of reforestations both in ecological terms as well as in terms of carbon trade mechanisms. A short example highlights the strong connection between the feasibility of reforestations and the global carbon market. Finally, open research questions are brought forth revealing the yet large research potential of Central Asian Desert ecosystems in general and in terms of carbon sequestration.
In this work, we theoretically investigate both aspects of charge-transferring atom-surface collisions: local-moment-type correlations and emission of secondary electrons from surfaces. Ideally, one chooses an approach that keeps as many electronic and lattice degrees of freedom at an ab-initio level as possible. In practice, however, this sophistication is hard to maintain. In this work, we do not aim to perform a description from first principles which could utilize density functional theory or quantum-chemical techniques. Instead, we keep only the most important degrees of freedom of the scattering process and use effective models for them. These are basically the Anderson-impurity model leading to time-dependent Anderson-Newns Hamiltonians and Gadzuk’s semiempirical approach to describe the projectile-target interaction from classical image shifts. In direct comparison with the description from first principles, the semiempirical approach offers a flexible basis for the modeling of a great variety of projectile-target combinations. The addition of further effective models to increase the general quality of the results is possible since the approach is very modular. The clear physical interpretation of each effective model, as well as the requirement for only a few and generally available parameters are further advantages of this approach. Rewritten in terms of Coleman’s pseudo-particle operators, the model is then numerically analyzed. This is done within a non-crossing approximation for the hybridization self-energies which are utilized by contour-ordered Green functions for each relevant electronic state of the projectile.
High resolution palaeo-ecological analysis of an Arctic ice-wedge polygon mire (Kytalyk, NE Siberia)
(2020)
Ice-wedge polygon mires are typical features of the Artic and therefore especially affected by climate change. They show, caused by soil-ice action, an amazing regular polygonal structure in meter dimension of higher and lower elevated dry and wet parts, and to this microtopography adapted vegetation. Polygon mires play, analogous to other mires, an important role in carbon sequestration, water balance, wildlife habitat and archive value with local to global significance. By storing enormous amounts of the global soil carbon polygon mires are crucial for our climate. Despite this relevance by covering large areas, polygon mires are comparatively poorly scientifically investigated and understood. It is still difficult to make forecasts on how polygon mires will develop under a changing climate in the Arctic, especially because internal factors and self-organisation complicate the understanding of their functioning. Therefore the investigation of modern and past polygon mires is necessary. This dissertation presents high resolution palaeo-ecological studies of a Northeast Siberian model polygon: ice-wedge polygon Lhc11 located in the Indigirka Lowlands at the scientific station Kytalyk. During field work in July 2011 the study site, covering an area of 26 × 21 m was divided into 546 plots, in which vegetation composition and microtopographical elevation characteristics were assessed and surface samples were collected. For palaeoecological analysis a 105.5 cm long peat section was excavated from the same site. Cluster analysis revealed five plant communities, which are clearly separated with respect to ground surface height, frost surface height and coverages of open water and vegetation, confirming the pattern already identified in other studies of Arctic ice-wedge polygons. The correct recognition of these patterns is crucial in palaeoecological studies in order to reconstruct landscape elements and their dynamics. This recognition requires insight in the short-distance relationships between surface elevation/wetness, vegetation and pollen deposition. The applied pollen-vegetation reference study shows that in general modern pollen deposition in polygon Lhc11 corresponds well with actual vegetation, allowing accurate reconstruction of local site conditions from fossil palynomorph sequences, including the reconstruction of the dynamics of closely spaced microtopographical elements. We conducted an evaluation of common palaeo proxies to compare their wetness reconstruction potential. The analysed proxies macrofossils, pollen, testate amoebae, geochemistry and sediment properties show similar wetness trends. Macrofossils provided the most detailed wetness reconstruction, spanning several wetness classes from very dry to wet, because they could be identified to genus or species level. However, as the proxies sometimes show contradictory results, a multi-proxy approach is preferable over a single proxy interpretation as it allows the reconstruction of environmental development in a broader palaeoecological context. For a better understanding of polygon dynamics and former greenhouse gas fluxes, more detailed and better quantified palaeo-microtopographical information is required. Therefore we developed a new transfer approach for modelling past Ground Surface Heights (GSH) in polygon mires from plant fossils. Based on the composition of modern vegetation we constructed two sets of potential fossil types (plant macrofossils and pollen), an extensive and a more restricted one. We applied Canonical Correspondence Analysis to model the relationships between potential fossil types and measured GSH. Both models show a strong relationship between modelled and measured GSH values and a high accuracy in prediction. Finally, we used the models to predict GSH values for Holocene peat samples. We found a fair correspondence with expert-based multi-proxy reconstruction of wetness conditions, even though only a minor part of the encountered fossils were represented in the GSH models, illustrating the robustness of the approach. The method can thus be used to reconstruct palaeoenvironmental conditions in a more objective way and can serve as a template for further palaeoecological studies. The 4000 years lasting history of the Lhc11 polygon site started with the establishment of a low-centre polygon in a drained thermokarst lake basin. Polygon Lhc11 formed part of a low-centre polygon for about 2000 years, experiencing enormous environmental influences discernible by incidence of silt, charred detritus, change of fossils composition and strongly declined peat accumulation rates and finally developed into a mature and degradation stage, into a low-high-centre polygon, currently characterized by high elevation differences. In the context of less studied but large-scale polygon mire occurrence, the high-resolution analysed ice-wedge polygon Lhc11 delivers insights into state and dynamics of a representative Siberian polygon site, in terms of modern and past vegetation and elevation characteristics. Furthermore the present study provides facilities for palaeoecological polygon studies including a new quantitative elevation modelling approach and provides valuable datasets for future research, e.g. greenhouse gas emissions and therefore contributes to a better understanding of these climate relevant ecosystems.
The importance of ion propulsion devices as an option for in-space propulsion of space
crafts and satellites continues to grow. They are more efficient than conventional chemi-
cal thrusters, which rely on burning their propellant, by ionizing the propellant gas in a
discharge channel and emitting the heavy ions at very high velocities. The ion emission
region of a thruster is called the plume and extends several meters axially and radially
downstream from the exit of a thruster. This region is particularly important for the effi-
ciency of a thruster, because it determines energy and angular distribution of the emitted
ions. It also determines the interaction with the carrier space craft by defining the electric
potential shape and the fluxes and energies of the emitted high energy ions, which are the
key parameters for sputter erosion of satellite components such as solar panels. Developing
new ion thrusters is expensive because of the high number of prototypes and testing cycles
required. Numerical modeling can help to reduce the costs in thruster development, but
the vastly differing length and time scales of the system, particularly the large differences of
scales between the discharge chamber and the plume, make a simulation challenging. Often
both regions are considered to be decoupled and are treated with different models to make
their simulation technically feasible. The coupling between channel and plume plasmas and
its influence on each other is disregarded, because there is no interaction between the two
regions. Therefore, this thesis investigates the physical effects which arise from this cou-
pling as well as models suitable for an integrated simulation of the whole coupled problem
of channel and plume plasmas. For this purpose the High Efficiency Multistage Plasma
Thruster (HEMP-T) ion thruster is considered.
For the discharge channel plasma, a fully kinetic model is required and the Particle-in-Cell
(PIC) method is applied. The PIC method requires very high spatial and temporal resolu-
tions which makes it computationally costly. As a result, only the discharge channel and the
near-field plume close to the channel exit can be simulated. In the channel, the results show
that electrons are magnetized and follow the magnetic field lines. The orientation of the
magnetic field there is mostly parallel to the symmetry axis and the channel walls which re-
sults in a high parallel electron transport and leads to a flat electric potential and a reduced
plasma-wall sheath. Only at the magnetic cusps, which are characteristic of HEMP-Ts the
electrons are guided towards the wall, with ions following due to quasineutrality, where a
classical plasma-wall sheath develops. The ion-wall contact is thus limited to the cusp re-
gion. The small radial drop of the potential towards the wall gives rather low energies of
ions impinging at the wall and minimizes erosion in the HEMP-T.
In the near-field plume, which extends from the thruster exit plane to some centimeters
downstream, the ion emission characteristics is defined. The ratio of radial and axial elec-
tric field components in this region determines the ion emission angle which should be
minimized for maximum thruster efficiency. The plasma discharge in the channel produces
high plasma densities and the subsequent drop from plasma to vacuum potential occurs
further downstream for higher densities. This increases the ratio of radial and axial electric
field components because the plasma expands radially outside of the confinement from the
dielectric discharge channel walls. The potential structure in the near-field plume impacts
also the supply of electrons for the channel discharge because the electrons enter the channel
from the plume. An effect which arises from this coupling is the breathing mode oscilla-
tion. It is an oscillation which is observed in all plasma quantities and is located near the
thruster exit. The oscillation frequency measured in the simulation is in good agreement
with a predator-prey estimate which validates this ansatz. However, the electron tempera-
ture, assumed constant in the predator-prey model, correlates inversely with the oscillation,
i.e. it is minimal at the current maximum and vice versa, which contributes to the observed
oscillations. Because of the oscillation of the plasma number density, the potential drop also
oscillates in the exit region and thus the ratio of radial to axial electric field components,
which results in the oscillation of the mean ion emission angle.
Regarding suitable models for a combined simulation of channel and plume plasmas, the
PIC model for channel and near-field plume is explicitly coupled to a hybrid fluid-PIC
model for the plume. The latter treats the electrons as a fluid, hence increasing the effective
spatial and temporal resolutions which can be applied in the plume simulations at the cost
of reduced accuracy of the electron model. Plasma densities decrease by two orders of
magnitude two meters downstream from the channel exit. The explicitly coupled kinetic
and hybrid PIC models are well suited for the computation of a HEMP-T and its plume
expansion, but they disregard the coupling of channel and plume plasmas for which other
methods are necessary. For this purpose a new approach is presented with a proof-of-
principle validation. The limited spatial resolution in the plume can be overcome with the
mesh-coarsening method, which increases the resolution in regions of low plasma density
without numerical artifacts. Sub-cycling for the electrons in the plume can then be used
to increase the temporal resolution in the plume. The combination of both methods, called
the sub-cycling mesh-coarsening (SMC) algorithm in the scope of this work, promises high
savings in computational cost which can make a combined simulation of plume and channel
plasmas feasible.
Kinetic modeling and infrared spectroscopy of charge carriers across the plasma-wall interface
(2022)
In this thesis, charge transport at the plasma-wall interface is investigated theoretically, on a semiclassical, microscopic level. Based on the Boltzmann and Poisson equations a set of equations is derived and numerically solved to model charge carriers both within a semiconducting wall and a gaseous plasma in front of it. While the plasma is considered collision-free, within the solid, phonon collisions, as well as recombination processes between conduction band electrons and valence band holes are considered. This results, for the first time, in a self-consistent modeling of both the gaseous electron-ion plasma and the electron-hole plasma in the solid on the same footing. Utilizing specific approximations for different physical scenarios, numerical solutions are presented both for the floating and the electronically contacted (biased) interface. In the latter case, the current voltage characteristic is calculated and shown to heavily depend on the charge kinetics within the wall.
Furthermore, we present optical methods to measure the wall charge noninvasively. These utilize the influence of the deposited surplus charges on the optical reflection coefficient of the surface. By calculating the optical response of these charges, we show that the magnitude of the surface charge can be inferred from the change in the reflectivity of the surface caused by the presence of the plasma. While nonlocal effects are considered, it is shown analytically and numerically that these can be neglected at the scales of the considered physical systems.
Gram-negative bacteria secrete lipopolysaccharides (LPS), leading to a host immune
response of proinflammatory cytokine secretion. Those proinflammatory cytokines are
TNF-α and IFN-γ, which induce the production of indoleamine 2,3-dioxygenase (IDO). IDO production is increased during severe sepsis, and septic shock. High IDO
levels are associated with increased mortality. This enzyme catalyzes the degradation of tryptophan (TRP) to kynurenine (KYN) along the kynurenine pathway (KP).
KYN is further degraded to kynurenic acid (KYNA). Increased IDO levels accompany
with increased levels of KYNA, which is associated with immunoparalysis.
Due to its central role, the KP is a potential target of therapeutic intervention.
The degradation of TRP to KYN by IDO was intervened by 1-Methyltryptophan (1-
MT), which is assumed to inhibit IDO. By administering 1-MT, the survival of
1-MT-treated mice suffering from sepsis increased compared to mice not treated with
1-MT. The levels of downstream metabolites such as KYN and KYNA were
expected to be decreased. Surprisingly, in healthy mice and pigs, an increase in KYNA
after 1-MT administration was reported. Those unexpected metabolite alterations after 1-MT administration, and the mode of action, were not the focus of recent
research. Hence, there is no explanation for KYNA increase, while KYN did not change.
This thesis aims to postulate a possible degradation pathway of 1-MT along the KP
with the help of ordinary differential equation (ODE) systems.
Moreover, the developed ODE models were used to determine the ability of 1-MT to
inhibit IDO in vivo. Therefore, a multiplicity of ODE models were developed, including
a model of the KP, an extension by lipopolysaccharide (LPS) administration, and 1-MT
administration.
Moreover, seven ODE models were developed, all considering possible degradation pathways of 1-MT. The most likely degradation pathway was combined with the ODE model
of LPS administration, including the inhibitory effects of 1-MT.
Those models consist of several dependent equations describing the dynamics of the KP.
For each component of the KP, one equation describes the alterations over time. Equations for TRP, KYN, KYNA, and quinolinic acid (QUIN) were developed.
Moreover, the alterations of serotonin (SER) were also included. All together belong
to the TRP metabolism. They include the degradation of TRP to SER and to KYN,
which is further degraded to KYNA and QUIN. Every degradation is catalyzed by an enzyme. Therefore, Michaelis-Menten (MM) equations were used employing the substrate
constant Km and the maximal degradation velocity Vmax. To reduce the complexity of
parameter calculation, Km values of the different enzymes were fixed to literature values.
The remaining parameters of the equations were determined so that the trajectories of
the calculated metabolite levels correspond to data. The parameters of different models were determined. To propose a degradation pathway of 1-MT leading to increased
KYNA levels, seven models were developed and compared. The most likely model was
extended to test whether the inhibitory effects of 1-MT on IDO can be determined.
Three different approaches determined the ODE model parameters of the different hypothesis of 1-MT degradation. In the first approach, ODE model parameters were fixed
to values fitted to an independent data set. In the second approach, parameters were
fitted to a subset of the data set, which was used for simulations of the different hypotheses. The third approach calculated ODE model parameters 100 times without
fixed parameters. The parameter set ending up in trajectories of the TRP metabolites,
which have the smallest distance to the data, was assumed to be the most likely. The
ODE model parameters were fitted to data measured in pigs. Two different
experimental models delivered data used in this thesis. The first experimental model
activates IDO by LPS administration in pigs. The second one combines the IDO
activation by LPS with the administration of 1-MT in pigs.
The most likely hypothesis, according to approach 1 was the degradation of 1-MT to
KYNA and TRP. For the second data set the most likely one was the direct degradation of 1-MT to KYNA. With approach 2 the most likely degradation pathways were
the combination of all degradation pathways and the degradation of 1-MT to TRP and
TRP to KYNA. With approach 3 the most likely way of KYNA increase was given by
the direct degradation of 1-MT to KYNA. In summary, the three approaches revealed
hypothesis 2, the direct degradation of 1-MT to KYNA most frequently. A cell-free
assay validated this result. This experiment combined 1-MT or TRP with or without
the enzyme kynurenine aminotransferase (KAT). KAT was already shown to degrade
TRP directly to KYNA. The levels of TRP, KYN and KYNA were measured. The
highest KYNA levels were yielded with an assay adding KAT to 1-MT, corresponding
to hypothesis 2. The models describing the inhibitory effects of 1-MT revealed that
the model without inhibitory effects of 1-MT on IDO was more likely for all three approaches.
The correctness of hypothesis 2 has to be confirmed by further in vitro experiments. It
also has to be investigated which reactions promote the degradation of 1-MT to KYNA.
The missing inhibitory properties of 1-MT on IDO, determined by the in silico ODE
models, align with previous research. It was shown that the saturation of 1-MT was too
low, e.g. in pigs, to inhibit IDO efficiently.
In this study, the first possible degradation pathway of 1-MT along the KP is proposed.
The reliability of the results depends on the quality of the experimental data, and the
season, when data were measured. Moreover, the results vary between the different
approaches of parameter fitting. Different approaches of parameter fitting have to be
included in the analysis to get more evidence for the correctness of the results.
In this doctoral thesis, algorithms are presented that are designed for the investigation in the mesopause region between the upper Mesosphere and Lower Thermosphere (MLT). The photochemical models are proposed and applied to represent the oxygen airglow and the oxygen photochemistry in the MLT. Atomic oxygen, O, in the ground state, O(3P), is of special interest because it is a reactive trace gas actively contributing to the Earth’s airglow. The retrievals of O(3P) concentrations, [O(3P)], are based on the nightglow time series of the green line emission measured remotely as in the first part of this thesis and the individual profiles of multiple nightglow emissions of O and molecular oxygen (O2) measured in situ as in the second part of this thesis. To process the complete spectral time series measured by using the satellite-borne instrument SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY), an intricate set of algorithms is developed and applied with the regularized total least squares minimization approach to estimate a set of the optimal regularization parameters and to retrieve a corresponding set of vertical Volume Emission Rate (VER) profiles. Furthermore, these algorithms take emissions of another origin and the Earth's shape into account. Considering not identified states of O2, the established photochemical models are adjusted resulting in two model modifications. Both model modifications are employed to retrieve the [O(3P)] time series on the basis of the VER time series in the MLT. The model input parameters vary in the atmosphere that motivated to propose these two model modifications and to employ available sources of the input parameters. One semi-empirical model, one general circulation model and the satellite-borne instrument SABER (Sounding of the Atmosphere using Broadband Emission Radiometry) are employed as sources of the reference [O(3P)] and input parameters time series. The SABER instrument employed as a source of the input parameters is preferred according to the comparison of the retrieved and reference [O(3P)] time series. Studying the impact of the 11-year solar cycle on O(3P) in the MLT, an algorithm is developed and applied with the Levenberg-Marquardt algorithm to estimate the optimal fit parameters step-wise. The results of the O(3P) sensitivity analysis obtained with respect to the solar activity forcing at the 11 year and 27 day time scales and the lunar gravitational forcing agree with the reference model simulations. The hypothesis regarding vertical shifts between different of Meinel bands at least partly caused by the hydroxyl radical (OH*) quenching with O(3P) is confirmed experimentally. Based on the conclusion drawn in the first part of this thesis that the data sets’ self-consistency is high as for the averaged SABER and SCIAMACHY measurements, a comprehensive set of available data with a higher level of the data sets’ self-consistency is employed in the second part of this thesis. Multiple airglow emissions measured in situ during four campaigns are employed to propose the Multiple Airglow Chemistry (MAC) model. Processed emissions are the Herzberg I, Chamberlain, Atmospheric and Infrared Atmospheric band emissions of O2 and the green line emission of O. Considering all widely known and additionally complemented reactions, the MAC model is proposed to represent the oxygen airglow and the oxygen photochemistry in the MLT. The presented MAC model is based on the hypothesis of Slanger et al. (2004) stating that higher excited states of O2 are coupled with each other through vibronic de-excitation caused by collisions among molecules of this group of O2 states in the MLT. This hypothesis is modified excluding the singlet Herzberg state of O2 from the group of O2 states considered by Slanger et al. (2004). The MAC calculations are carried out sequentially starting with higher excited O2 states to provide the retrieved output concentrations of these O2 states as the input concentrations to the next calculation steps. The final step is only based on concentrations of all species, whereas each of the earlier steps is based on a corresponding VER profile besides of the input concentrations. The oxygen photochemistry in the MLT is represented by all species considered at the final step that makes it possible to adopt the MAC reactions in a general circulation model. Four modifications of the MAC model, i.e. including or excluding the triplet Herzberg states of O2 and including or excluding ozone and odd hydrogen (hydrogen, OH* and hydroperoxy radical), lead to negligible differences in the retrieved [O(3P)] profiles. Based on the MAC calculations verified and validated on the basis of the four rocket campaigns, one of the effective modifications of the MAC model (excluding the triplet Herzberg states of O2, ozone and odd hydrogen) is further reduced to the most effective modification. This implies that for the [O(3P)] retrieval only the O2 Atmospheric band emission, temperature and concentrations of molecular nitrogen (N2) and O2 are sufficient to apply. Calculations carried out by using the most effective modification of the MAC model are verified and validated on the basis of self-consistent in situ measurements obtained simultaneously. The MAC model enables identifying precursors of (1) the three lowest O2 valence states and (2) the second excited O state responsible for (1) the Atmospheric and Infrared Atmospheric band emissions of O2 and (2) the green line emission of O, respectively. Particularly, the singlet Herzberg state of O2 is identified as the major precursor of the second excited O state resulting in the green line emission. In focus of potential further research is an extension of the MAC model with vibrationally excited states of O2 and ionized species.
Underground hard coal mining operations irreversibly disrupt the pre-existing mechanical equilibrium of the geological media. The employment of high-recovery methods modifies the stress field of the sedimentary sequence, generating movement and faulting of the rock layers above and below mined seams. These new fracture zones do affect the original conditions of the hydrogeological system by modifying flow pathways and increasing the permeability of the rock sequence. Moreover, the surface area of rock exposed to air and water is increased, conditioning the water-rock interaction. Despite this rather clear conceptualization, flow and reactive transport processes in fractured overburdens are rarely modeled simultaneously. Discrete setups that consider fractures and porous matrix require extensive characterization of both media, which is impractical for regional case studies. As a result, most post-mining models explicitly ignore fracture structures by employing the equivalent porous approach or even both media with lumped parameter models. However, omitting either medium represents a delicate simplification, considering that mining-related fractures control the rate and direction of water flow within moderately permeable but relatively highly porous rock sequences.
In this dissertation, the specific contribution of fractured and matrix continua to the transient discharge and water quality of a post-mining coal zone is quantified and evaluated. For this purpose, dual and multiple interacting continua models are employed to simulate fluid flow and reactive mass transport in fractured and variable water-saturated rock sequences. The effectiveness of the models is evaluated by simulating the origin, generation and transport of acid mine drainage (i.e., water with elevated concentrations of hydrogen, iron, sulfate and chloride) within the shallow overburden of the Ibbenbüren Westfield. Compared to other coal districts in Germany, this area is strongly delimited by the local geology and topography, resulting in a well-defined hydrogeological system to test the models. Petrographic and chemical analyses performed on core samples from the area show the strong influence of mining-derived fractures on the water-rock interaction within the Carboniferous sequence. The presence of oxidized pyrite along with amorphous iron hydroxide phases in weathering fronts on both sides of the fractures demonstrates the exchange of solutes and gases between the fractured and the porous matrix media.
Based on the previous evidence, the TOUGHREACT software is employed to characterize flow and reactive transport processes in the Westfield. However, each of the two processes is simulated at separate stages to have more control in the adjustment of sensitive parameters for which little information is available. For the flow component, a dual continuum model, with Richard’s equations is used to characterize the unsaturated water flow in both fractured and matrix media. Under this approach, the model adequately reproduces the bimodal flow behavior of the discharges measured in the mine drainage for the years 2008 and 2017. Simulation results show how the fractured continuum generates intense discharge events during the winter months while the rock matrix controls smooth discharge limbs in summer, when water is slowly released back to the fractures. With the flow component calibrated, the second part of the study incorporates the geochemical processes into the model based on actual data from the rock samples. Their simulation requires extending the two-continuum setup to a multiple continua model with five nested block strings: one for the fractures and four for the rock matrix. This further subdivision prevents under-representations of kinetic reactions with short equilibrium length scales and numerical instabilities due to lack of chemical and flow gradients. As a result, the new multiple continua model provides good agreement with respect to long- and short-term concentrations and discharge trends measured in the mine drainage. The flow of oxygen and meteoric water through the fractured continuum leads to a high and steady release of hydrogen, iron and sulfate ions derived from pyrite oxidation in the matrix continua closest to the fractures. Moreover, high chloride concentrations result from the mixing and gradual release of relatively immobile solutes in the matrix as they interact with percolating water in the fracture. Both findings are equally congruent with the reactive pyrite oxidation and iron hydroxide precipitation fronts identified in the fractured core samples.
In the end, the multiple continua models, the simulation procedure and the results of the benchmark and sensitivity analysis scenarios developed for the Westfield pave the way for the application of the approach in other mining zones. The first candidate emerges in the Ibbenbüren Eastfield, where a coupled elemental-isotopic approach included in this thesis has confirmed that water-conducting fracture zones are primary elements for solute generation and transport in the first 300 meters of the overburden. In the latter case, calibration and verification of the models can be complemented with measurements of δ34S in sulfates and δ18O, δ2H, and Tritium in water.