Soil application of Meliaceae derivatives: effect on carbon and nitrogen dynamics in the soil-plant system

The effect of soil incorporation of 7 Meliaceae derivatives (6 commercial neem cakes and leaves of Melia azedarach L.) on C and N dynamics and on nutrient availability to micropropagated GF677 rootstock was investigated. In a first laboratory incubation experiment the derivatives showed different N mineralization dynamics, generally well predicted by their C:N ratio and only partly by their initial N concentration. All derivatives increased microbial biomass C, thus representing a source of C for the soil microbial population. Soil addition of all neem cakes (8 g kg-1) and melia leaves (16 g kg-1) had a positive effect on plant growth and increased root N uptake and leaf green colour of micropropagated plants of GF677. In addition, the neem cakes characterized by higher nutrient concentration increased P and K concentration in shoot and leaves 68 days after the amendment. In another experiment, soil incorporation of 15N labeled melia leaves (16 g kg-1) had no effect on the total amount of plant N, however the percentage of melia derived-N of treated plants ranged between 0.8% and 34% during the experiment. At the end of the growing season, about 7% of N added as melia leaves was recovered in plant, while 70% of it was still present in soil. Real C mineralization and the priming effect induced by the addition of the derivatives were quantified by a natural 13C abundance method. The real C mineralization of the derivatives ranged between 22% and 40% of added-C. All the derivatives studied induced a positive priming effect and, 144 days after the amendment, the amount of C primed corresponded to 26% of added-C, for all the derivatives. Despite this substantial priming effect, the C balance of the soil, 144 days after the amendment, always resulted positive.

The Dynamics of Passive Torsional Fatigue Test Rigs: Innovative Applications of Universal Joints

The dynamics of a passive back-to-back test rig have been characterised, leading to a multi-coordinate approach for the analysis of arbitrary test configurations. Universal joints have been introduced into a typical pre-loaded back-to-back system in order to produce an oscillating torsional moment in a test specimen. Two different arrangements have been investigated using a frequency-based sub-structuring approach: the receptance method. A numerical model has been developed in accordance with this theory, allowing interconnection of systems with two-coordinates and closed multi-loop schemes. The model calculates the receptance functions and modal and deflected shapes of a general system. Closed form expressions of the following individual elements have been developed: a servomotor, damped continuous shaft and a universal joint. Numerical results for specific cases have been compared with published data in literature and experimental measurements undertaken in the present work. Due to the complexity of the universal joint and its oscillating dynamic effects, a more detailed analysis of this component has been developed. Two models have been presented. The first represents the joint as two inertias connected by a massless cross-piece. The second, derived by the dynamic analysis of a spherical four-link mechanism, considers the contribution of the floating element and its gyroscopic effects. An investigation into non-linear behaviour has led to a time domain model that utilises the Runge-Kutta fourth order method for resolution of the dynamic equations. It has been demonstrated that the torsional receptances of a universal joint, derived using the simple model, result in representation of the joint as an equivalent variable inertia. In order to verify the model, a test rig has been built and experimental validation undertaken. The variable inertia of a universal joint has lead to a novel application of the component as a passive device for the balancing of inertia variations in slider-crank mechanisms.

Characterization of structure and dynamics of submicrometric systems via fluorescence correlation spectroscopy

The presented thesis revolves around the study of thermally-responsive PNIPAAm-based hydrogels in water/based environments, as studied by Fluorescence Correlation Spectroscopy (FCS).rnThe goal of the project was the engineering of PNIPAAm gels into biosensors. Specifically, a gamma of such gels were both investigated concerning their dynamics and structure at the nanometer scale, and their performance in retaining bound bodies upon thermal collapse (which PNIPAAm undergoes upon heating above 32 ºC).rnFCS’s requirements, as a technique, match the limitations imposed by the system. Namely, the need to intimately probe a system in a solvent, which was also fragile and easy to alter. FCS, on the other hand, both requires a fluid environment to work, and is based on the observation of diffusion of fluorescents at nanomolar concentrations. FCS was applied to probe the hydrogels on the nanometer size with minimal invasivity.rnVariables in the gels were addressed in the project including crosslinking degree; structural changes during thermal collapse; behavior in different buffers; the possibility of decreasing the degree of inhomogeneity; behavior of differently sized probes; and the effectiveness of antibody functionalization upon thermal collapse.rnThe evidenced results included the heightening of structural inhomogeneities during thermal collapse and under different buffer conditions; the use of annealing to decrease the inhomogeneity degree; the use of differently sized probes to address different length scale of the gel; and the successful functionalization before and after collapse.rnThe thesis also addresses two side projects, also carried forward via FCS. One, diffusion in inverse opals, produced a predictive simulation model for diffusion of bodies in confined systems as dependent on the bodies’ size versus the characteristic sizes of the system. The other was the observation of interaction of bodies of opposite charge in a water solution, resulting in a phenomenological theory and an evaluation method for both the average residence time of the different bodies together, and their attachment likelihood.

Development of a field method to assess root parasites of woody plants - ecological analyses of the grapevine-phylloxera interaction

The development and the growth of plants is strongly affected by the interactions between roots, rootrnassociated organisms and rhizosphere communities. Methods to assess such interactions are hardly torndevelop particularly in perennial and woody plants, due to their complex root system structure and theirrntemporal change in physiology patterns. In this respect, grape root systems are not investigated veryrnwell. The aim of the present work was the development of a method to assess and predict interactionsrnat the root system of rootstocks (Vitis berlandieri x Vitis riparia) in field. To achieve this aim, grapernphylloxera (Daktulosphaira vitifoliae Fitch, Hemiptera, Aphidoidea) was used as a graperoot parasitizingrnmodel.rnTo develop the methodical approach, a longt-term trial (2006-2009) was arranged on a commercial usedrnvineyard in Geisenheim/Rheingau. All 2 to 8 weeks the top most 20 cm of soil under the foliage wallrnwere investigated and root material was extracted (n=8-10). To include temporal, spatial and cultivarrnspecific root system dynamics, the extracted root material was analyzed digitally on the morphologicalrnproperties. The grape phylloxera population was quantified and characterized visually on base of theirrnlarvalstages (oviparous, non oviparous and winged preliminary stages). Infection patches (nodosities)rnwere characterized visually as well, partly supported by digital root color analyses. Due to the knownrneffects of fungal endophytes on the vitality of grape phylloxera infested grapevines, fungal endophytesrnwere isolated from nodosity and root tissue and characterized (morphotypes) afterwards. Further abioticrnand biotic soil conditions of the vineyards were assessed. The temporal, spatial and cultivar specificrnsensitivity of single parameters were analyzed by omnibus tests (ANOVAs) and adjacent post-hoc tests.rnThe relations between different parameters were analyzed by multiple regression models.rnQuantitative parameters to assess the degeneration of nodosity, the development nodosity attachedrnroots and to differentiate between nodosities and other root swellings in field were developed. Significantrndifferences were shown between root dynamic including parameters and root dynamic ignoringrnparameters. Regarding the description of grape phylloxera population and root system dynamic, thernmethod showed a high temporal, spatial and cultivar specific sensitivity. Further, specific differencesrncould be shown in the frequency of endophyte morphotypes between root and nodosity tissue as wellrnas between cultivars. Degeneration of nodosities as well as nodosity occupation rates could be relatedrnto the calculated abundances of grape phylloxera population. Further ecological questions consideringrngrape root development (e.g. relation between moisture and root development) and grape phylloxerarnpopulation development (e.g. relation between temperature and population structure) could be answeredrnfor field conditions.rnGenerally, the presented work provides an approach to evaluate vitality of grape root systems. Thisrnapproach can be useful, considering the development of control strategies against soilborne pests inrnviticulture (e.g. grape phylloxera, Sorospheara viticola, Roesleria subterranea (Weinm.) Redhaed) as well as considering the evaluation of integrated management systems in viticulture.

A mesoscopic simulation method for electrolyte solutions

This thesis deals with the development of a novel simulation technique for macromolecules in electrolyte solutions, with the aim of a performance improvement over current molecular-dynamics based simulation methods. In solutions containing charged macromolecules and salt ions, it is the complex interplay of electrostatic interactions and hydrodynamics that determines the equilibrium and non-equilibrium behavior. However, the treatment of the solvent and dissolved ions makes up the major part of the computational effort. Thus an efficient modeling of both components is essential for the performance of a method. With the novel method we approach the solvent in a coarse-grained fashion and replace the explicit-ion description by a dynamic mean-field treatment. Hence we combine particle- and field-based descriptions in a hybrid method and thereby effectively solve the electrokinetic equations. The developed algorithm is tested extensively in terms of accuracy and performance, and suitable parameter sets are determined. As a first application we study charged polymer solutions (polyelectrolytes) in shear flow with focus on their viscoelastic properties. Here we also include semidilute solutions, which are computationally demanding. Secondly we study the electro-osmotic flow on superhydrophobic surfaces, where we perform a detailed comparison to theoretical predictions.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

The dynamics of ultra-thin polymer films in different confinements studied by resonance enhanced dynamic light scattering

Die Kontroverse über den Glasübergang im Nanometerbereich, z. B. die Glas¬über¬gangs-temperatur Tg von dünnen Polymerfilmen, ist nicht vollständig abgeschlossen. Das dynamische Verhalten auf der Nanoskala ist stark von den einschränkenden Bedingungen abhängig, die auf die Probe wirken. Dünne Polymerfilme sind ideale Systeme um die Dynamik von Polymerketten unter der Einwirkung von Randbedingungen zu untersuchen, wie ich sie in dieser Arbeit variiert habe, um Einblick in dieses Problem zu erhalten.rnrnResonanzverstärkte dynamische Lichtstreuung ist eine Methode, frei von z.B. Fluoreszenzmarkern, die genutzt werden kann um in dünnen Polymerfilmen dynamische Phänomene

Polymer nanoparticles : their functionalization for biomedical applications and dynamics near an electrode

In this work, a method for the functionalization of biocompatible, poly(lactic acid)-based nanoparticles with charged moieties or fluorescent labels is presented. Therefore, a miniemulsion solvent evaporation procedure is used in which prepolymerized poly(L-lactic acid) is used together with a previously synthesized copolymer of methacrylic acid or a polymerizable dye, respectively, and an oligo(lactic acid) macromonomer. Alternatively, the copolymerization has been carried out in one step with the miniemulsion solvent evaporation. The increased stability in salty solutions of the carboxyl-modified nanoparticles compared to nanoparticles consisting of poly(lactic acid) only has been shown in light scattering experiments. The properties of the nanoparticles that were prepared with the separately synthesized copolymer were almost identical to those in which the copolymerization and particle fabrication were carried out simultaneously. During the characterization of the fluorescently labeled nanoparticles, the focus was on the stable bonding between the fluorescent dye and the rest of the polymer chain to ensure that none of it is released from the particles, even after longer storage time or during lengthy experiments. In a fluorescence correlation spectroscopy experiment, it could be shown that even after two weeks, no dye has been released into the solvent. Besides biomedical research for which the above described, functionalized nanoparticles were optimized, nanoparticles also play a role in coating technology. One possibility to fabricate coatings is the electrophoretic deposition of particles. In this process, the mobility of nanoparticles near electrode interfaces plays a crucial role. In this thesis, the nanoparticle mobility has been investigated with resonance enhanced dynamic light scattering (REDLS). A new setup has been developed in which the evanescent electromagnetic eld of a surface plasmon that propagates along the gold-sample interface has been used as incident beam for the dynamic light scattering experiment. The gold layer that is necessary for the excitation of the plasmon doubles as an electrode. Due to the penetration depth of the surface plasmon into the sample layer that is limited to ca. 200 nm, insights on the voltage- and frequency dependent mobility of the nanoparticles near the electrode could be gained. Additionally, simultaneous measurements at four different scattering angles can be carried out with this setup, therefore the investigation of samples undergoing changes is feasible. The results were discussed in context with the mechanisms of electrophoretic deposition.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Circulation dynamics and its influence on European and Mediterranean January-April climate over the past half millennium: results and insights from instrumental data, documentary evidence and coupled climate models

We use long instrumental temperature series together with available field reconstructions of sea-level pressure (SLP) and three-dimensional climate model simulations to analyze relations between temperature anomalies and atmospheric circulation patterns over much of Europe and the Mediterranean for the late winter/early spring (January–April, JFMA) season. A Canonical Correlation Analysis (CCA) investigates interannual to interdecadal covariability between a new gridded SLP field reconstruction and seven long instrumental temperature series covering the past 250 years. We then present and discuss prominent atmospheric circulation patterns related to anomalous warm and cold JFMA conditions within different European areas spanning the period 1760–2007. Next, using a data assimilation technique, we link gridded SLP data with a climate model (EC-Bilt-Clio) for a better dynamical understanding of the relationship between large scale circulation and European climate. We thus present an alternative approach to reconstruct climate for the pre-instrumental period based on the assimilated model simulations. Furthermore, we present an independent method to extend the dynamic circulation analysis for anomalously cold European JFMA conditions back to the sixteenth century. To this end, we use documentary records that are spatially representative for the long instrumental records and derive, through modern analogs, large-scale SLP, surface temperature and precipitation fields. The skill of the analog method is tested in the virtual world of two three-dimensional climate simulations (ECHO-G and HadCM3). This endeavor offers new possibilities to both constrain climate model into a reconstruction mode (through the assimilation approach) and to better asses documentary data in a quantitative way.

Photodissociation dynamics of tert-butylnitrite following excitation to the S(1) and S(2) states. A study by velocity-map ion-imaging and 3D-REMPI spectroscopy

Excitation of tert-butylnitrite into the first and second UV absorption bands leads to efficient dissociation into the fragment radicals NO and tert-butoxy in their electronic ground states (2)Π and (2)E, respectively. Velocity distributions and angular anisotropies for the NO fragment in several hundred rotational and vibrational quantum states were obtained by velocity-map imaging and the recently developed 3D-REMPI method. Excitation into the well resolved vibronic progression bands (k = 0, 1, 2) of the NO stretch mode in the S(1) ← S(0) transition produces NO fragments mostly in the vibrational state with v = k, with smaller fractions in v = k - 1 and v = k - 2. It is concluded that dissociation occurs on the purely repulsive PES of S(1) without barrier. All velocity distributions from photolysis via the S(1)(nπ*) state are monomodal and show high negative anisotropy (β ≈ -1). The rotational distributions peak near j = 30.5 irrespective of the vibronic state S(1)(k) excited and the vibrational state v of the NO fragment. On average 46% of the excess energy is converted to kinetic energy, 23% and 31% remain as internal energy in the NO fragment and the t-BuO radical, respectively. Photolysis via excitation into the S(2) ← S(0) transition at 227 nm yields NO fragments with about equal populations in v = 0 and v = 1. The rotational distributions have a single maximum near j = 59.5. The velocity distributions are monomodal with positive anisotropy β ≈ 0.8. The average fractions of the excess energy distributed into translation, internal energy of NO, and internal energy of t-BuO are 39%, 23%, and 38%, respectively. In all cases ∼8500 cm(-1) of energy remain in the internal degrees of freedom of the t-BuO fragment. This is mostly assigned to rotational energy. An ab initio calculation of the dynamic reaction path shows that not only the NO fragment but also the t-BuO fragment gain large angular momentum during dissociation on the purely repulsive potential energy surface of S(2).

Towards dynamics in flood risk assessment

As a consequence of flood impacts, communities inhabiting mountain areas are increasingly affected by considerable damage to infrastructure and property. The design of effective flood risk mitigation strategies and their subsequent implementation is crucial for a sustainable development in mountain areas. The assessment of the dynamic evolution of flood risk is the pillar of any subsequent planning process that is targeted at a reduction of the expected adverse consequences of the hazard impact. Given these premises, firstly, a comprehensive method to derive flood hazard process scenarios for well-defined areas at risk is presented. Secondly, conceptualisations of a static and dynamic flood risk assessment are provided. These are based on formal schemes to compute the risk mitigation performance of devised mitigation strategies within the framework of economic cost-benefit analysis. In this context, techniques suitable to quantify the expected losses induced by the identified flood impacts are provided.

Open top chambers and infrared lamps : a comparison of heating efficacy and CO₂/CH₄ dynamics in a Lake Superior coastal peatland

Experimental warming provides a method to determine how an ecosystem will respond to increased temperatures. Northern peatland ecosystems, sensitive to changing climates, provide an excellent setting for experimental warming. Storing great quantities of carbon, northern peatlands play a critical role in regulating global temperatures. Two of the most common methods of experimental warming include open top chambers (OTCs) and infrared (IR) lamps. These warming systems have been used in many ecosystems throughout the world, yet their efficacy to create a warmer environment is variable and has not been widely studied. To date, there has not been a direct, experimentally controlled comparison of OTCs and IR lamps. As a result, a factorial study was implemented to compare the warming efficacy of OTCs and IR lamps and to examine the resulting carbon dioxide (CO2) and methane (CH4) flux rates in a Lake Superior peatland. IR lamps warmed the ecosystem on average by 1-2 #°C, with the majority of warming occurring during nighttime hours. OTC's did not provide any long-term warming above control plots, which is contrary to similar OTC studies at high latitudes. By investigating diurnal heating patterns and micrometeorological variables, we were able to conclude that OTCs were not achieving strong daytime heating peaks and were often cooler than control plots during nighttime hours. Temperate day-length, cloudy and humid conditions, and latent heat loss were factors that inhibited OTC warming. There were no changes in CO2 flux between warming treatments in lawn plots. Gross ecosystem production was significantly greater in IR lamp-hummock plots, while ecosystem respiration was not affected. CH4 flux was not significantly affected by warming treatment. Minimal daytime heating differences, high ambient temperatures, decay resistant substrate, as well as other factors suppressed significant gas flux responses from warming treatments.

Laser-induced transient grating analysis of dynamics of interaction between sensory rhodopsin II D75N and the HtrII transducer.

The interaction between sensory rhodopsin II (SRII) and its transducer HtrII was studied by the time-resolved laser-induced transient grating method using the D75N mutant of SRII, which exhibits minimal visible light absorption changes during its photocycle, but mediates normal phototaxis responses. Flash-induced transient absorption spectra of transducer-free D75N and D75N joined to 120 amino-acid residues of the N-terminal part of the SRII transducer protein HtrII (DeltaHtrII) showed only one spectrally distinct K-like intermediate in their photocycles, but the transient grating method resolved four intermediates (K(1)-K(4)) distinct in their volumes. D75N bound to HtrII exhibited one additional slower kinetic species, which persists after complete recovery of the initial state as assessed by absorption changes in the UV-visible region. The kinetics indicate a conformationally changed form of the transducer portion (designated Tr*), which persists after the photoreceptor returns to the unphotolyzed state. The largest conformational change in the DeltaHtrII portion was found to cause a DeltaHtrII-dependent increase in volume rising in 8 micros in the K(4) state and a drastic decrease in the diffusion coefficient (D) of K(4) relatively to those of the unphotolyzed state and Tr*. The magnitude of the decrease in D indicates a large structural change, presumably in the solvent-exposed HAMP domain of DeltaHtrII, where rearrangement of interacting molecules in the solvent would substantially change friction between the protein and the solvent.

Inverse modelling of the 14C bomb pulse in stalagmites to constrain the dynamics of soil carbon cycling at selected European cave sites

The decomposition of soil organic matter (SOM) is temperature dependent, but its response to a future warmer climate remains equivocal. Enhanced rates of decomposition of SOM under increased global temperatures might cause higher CO2 emissions to the atmosphere, and could therefore constitute a strong positive feedback. The magnitude of this feedback however remains poorly understood, primarily because of the difficulty in quantifying the temperature sensitivity of stored, recalcitrant carbon that comprises the bulk (>90%) of SOM in most soils. In this study we investigated the effects of climatic conditions on soil carbon dynamics using the attenuation of the 14C ‘bomb’ pulse as recorded in selected modern European speleothems. These new data were combined with published results to further examine soil carbon dynamics, and to explore the sensitivity of labile and recalcitrant organic matter decomposition to different climatic conditions. Temporal changes in 14C activity inferred from each speleothem was modelled using a three pool soil carbon inverse model (applying a Monte Carlo method) to constrain soil carbon turnover rates at each site. Speleothems from sites that are characterised by semi-arid conditions, sparse vegetation, thin soil cover and high mean annual air temperatures (MAATs), exhibit weak attenuation of atmospheric 14C ‘bomb’ peak (a low damping effect, D in the range: 55–77%) and low modelled mean respired carbon ages (MRCA), indicating that decomposition is dominated by young, recently fixed soil carbon. By contrast, humid and high MAAT sites that are characterised by a thick soil cover and dense, well developed vegetation, display the highest damping effect (D = c. 90%), and the highest MRCA values (in the range from 350 ± 126 years to 571 ± 128 years). This suggests that carbon incorporated into these stalagmites originates predominantly from decomposition of old, recalcitrant organic matter. SOM turnover rates cannot be ascribed to a single climate variable, e.g. (MAAT) but instead reflect a complex interplay of climate (e.g. MAAT and moisture budget) and vegetation development.