Refine
Year of publication
- 2022 (73) (remove)
Document Type
- Article (53)
- Doctoral Thesis (9)
- Postprint (7)
- Report (2)
- Habilitation Thesis (1)
- Master's Thesis (1)
Is part of the Bibliography
- yes (73)
Keywords
- climate change (6)
- Klimawandel (3)
- Lake Malawi (3)
- Shire River (3)
- convolutional neural network (3)
- drought intensity (3)
- drought projections (3)
- floods (3)
- meteorological drought (3)
- validation (3)
Institute
- Institut für Umweltwissenschaften und Geographie (73) (remove)
After a century of semi-restricted floodplain development, Southern Alberta, Canada, was struck by the devastating 2013 Flood. Aging infrastructure and limited property-level floodproofing likely contributed to the $4-6 billion (CAD) losses. Following this catastrophe, Alberta has seen a revival in flood management, largely focused on structural protections. However, concurrent with the recent structural work was a 100,000+ increase in Calgary's population in the 5 years following the flood, leading to further densification of high-hazard areas. This study implements the novel Stochastic Object-based Flood damage Dynamic Assessment (SOFDA) model framework to quantify the progression of the direct-damage flood risk in a mature urban neighborhood after the 2013 Flood. Five years of remote-sensing data, property assessment records, and inundation simulations following the flood are used to construct the model. Results show that in these 5 years, vulnerability trends (like densification) have increased flood risk by 4%; however, recent structural mitigation projects have reduced overall flood risk by 47% for this case study. These results demonstrate that the flood management revival in Southern Alberta has largely been successful at reducing flood risk; however, the gains are under threat from continued development and densification absent additional floodproofing regulations.
Phytoliths in particulate matter released by wind erosion on arable land in La Pampa, Argentina
(2022)
Silicon (Si) is considered a beneficial element in plant nutrition, but its importance on ecosystems goes far beyond that. Various forms of silicon are found in soils, of which the phytogenic pool plays a decisive role due to its good availability. This Si returns to the soil through the decomposition of plant residues, where they then participate in the further cycle as biogenic amorphous silica (bASi) or so-called phytoliths. These have a high affinity for water, so that the water holding capacity and water availability of soils can be increased even by small amounts of ASi. Agricultural land is a considerable global dust source, and dust samples from arable land have shown in cloud formation experiments a several times higher ice nucleation activity than pure mineral dust. Here, particle sizes in the particulate matter fractions (PM) are important, which can travel long distances and reach high altitudes in the atmosphere. Based on this, the research question was whether phytoliths could be detected in PM samples from wind erosion events, what are the main particle sizes of phytoliths and whether an initial quantification was possible.Measurements of PM concentrations were carried out at a wind erosion measuring field in the province La Pampa, Argentina. PM were sampled during five erosion events with Environmental Dust Monitors (EDM). After counting and classifying all particles with diameters between 0.3 and 32 mu m in the EDMs, they are collected on filters. The filters were analyzed by Scanning Electron Microscopy and Energy Dispersive X-Ray analysis (SEM-EDX) to investigate single or ensembles of particles regarding composition and possible origins.The analyses showed up to 8.3 per cent being phytoliths in the emitted dust and up to 25 per cent of organic origin. Particles of organic origin are mostly in the coarse dust fraction, whereas phytoliths are predominately transported in the finer dust fractions. Since phytoliths are both an important source of Si as a plant nutrient and are also involved in soil C fixation, their losses from arable land via dust emissions should be considered and its specific influence on atmospheric processes should be studied in detail in the future.
We present a chronology framework named LegacyAge 1.0 containing harmonized chronologies for 2831 pollen records (downloaded from the Neotoma Paleoecology Database and the supplementary Asian datasets) together with their age control points and metadata in machine-readable data formats.
All chronologies use the Bayesian framework implemented in Bacon version 2.5.3. Optimal parameter settings of priors (accumulation.shape, memory.strength, memory.mean, accumulation.rate, and thickness) were identified based on information in the original publication or iteratively after preliminary model inspection.
The most common control points for the chronologies are radiocarbon dates (86.1 %), calibrated by the latest calibration curves (IntCal20 and SHCal20 for the terrestrial radiocarbon dates in the Northern Hemisphere and Southern Hemisphere and Marine20 for marine materials).
The original publications were consulted when dealing with outliers and inconsistencies. Several major challenges when setting up the chronologies included the waterline issue (18.8% of records), reservoir effect (4.9 %), and sediment deposition discontinuity (4.4 %).
Finally, we numerically compare the LegacyAge 1.0 chronologies to those published in the original publications and show that the reliability of the chronologies of 95.4% of records could be improved according to our assessment.
Our chronology framework and revised chronologies provide the opportunity to make use of the ages and age uncertainties in synthesis studies of, for example, pollen-based vegetation and climate change.
The LegacyAge 1.0 dataset, including metadata, datings, harmonized chronologies, and R code used, is openaccess and available at PANGAEA (https://doi.org/10.1594/PANGAEA.933132; Li et al., 2021) and Zenodo (https://doi.org/10.5281/zenodo.5815192; Li et al., 2022), respectively.
Earthquake site responses or site effects are the modifications of surface geology to seismic waves. How well can we predict the site effects (average over many earthquakes) at individual sites so far? To address this question, we tested and compared the effectiveness of different estimation techniques in predicting the outcrop Fourier site responses separated using the general inversion technique (GIT) from recordings. Techniques being evaluated are (a) the empirical correction to the horizontal-to-vertical spectral ratio of earthquakes (c-HVSR), (b) one-dimensional ground response analysis (GRA), and (c) the square-root-impedance (SRI) method (also called the quarter-wavelength approach). Our results show that c-HVSR can capture significantly more site-specific features in site responses than both GRA and SRI in the aggregate, especially at relatively high frequencies. c-HVSR achieves a "good match" in spectral shape at similar to 80%-90% of 145 testing sites, whereas GRA and SRI fail at most sites. GRA and SRI results have a high level of parametric and/or modeling errors which can be constrained, to some extent, by collecting on-site recordings.
The 2020s are an essential decade for achieving the 2030 Agenda and its Sustainable Development Goals (SDGs). For this, SDG research needs to provide evidence that can be translated into concrete actions. However, studies use different SDG data, resulting in incomparable findings. Researchers primarily use SDG databases provided by the United Nations (UN), the World Bank Group (WBG), and the Bertelsmann Stiftung & Sustainable Development Solutions Network (BE-SDSN). We compile these databases into one unified SDG database and examine the effects of the data selection on our understanding of SDG interactions. Among the databases, we observed more different than similar SDG interactions. Differences in synergies and trade-offs mainly occur for SDGs that are environmentally oriented. Due to the increased data availability, the unified SDG database offers a more nuanced and reliable view of SDG interactions. Thus, the SDG data selection may lead to diverse findings, fostering actions that might neglect or exacerbate trade-offs.
The fluxes of water and solutes in the subsurface compartment of the Critical Zone are temporally dynamic and it is unclear how this impacts microbial mediated nutrient cycling in the spatially heterogeneous subsurface. To investigate this, we undertook numerical modeling, simulating the transport in a wide range of spatially heterogeneous domains, and the biogeochemical transformation of organic carbon and nitrogen compounds using a complex microbial community with four (4) distinct functional groups, in water saturated subsurface compartments. We performed a comprehensive uncertainty analysis accounting for varying residence times and spatial heterogeneity. While the aggregated removal of chemical species in the domains over the entire simulation period was approximately the same as that in steady state conditions, the sub-scale temporal variation of microbial biomass and chemical discharge from a domain depended strongly on the interplay of spatial heterogeneity and temporal dynamics of the forcing. We showed that the travel time and the Damkohler number (Da) can be used to predict the temporally varying chemical discharge from a spatially heterogeneous domain. In homogeneous domains, chemical discharge in temporally dynamic conditions could be double of that in the steady state conditions while microbial biomass varied up to 75% of that in steady state conditions. In heterogeneous domains, the interquartile range of uncertainty in chemical discharge in reaction dominated systems (log(10)Da > 0) was double of that in steady state conditions. However, high heterogeneous domains resulted in outliers where chemical discharge could be as high as 10-20 times of that in steady state conditions in high flow periods. And in transport dominated systems (log(10)Da < 0), the chemical discharge could be half of that in steady state conditions in unusually low flow conditions. In conclusion, ignoring spatio-temporal heterogeneities in a numerical modeling approach may exacerbate inaccurate estimation of nutrient export and microbial biomass. The results are relevant to long-term field monitoring studies, and for homogeneous soil column-scale experiments investigating the role of temporal dynamics on microbial redox dynamics.
Despite the amount of research focussed on the Alpine orogen, different hypotheses still exist regarding varying spatial seismicity distribution patterns throughout the region. Previous measurement-constrained regional 3D models of lithospheric density distribution and thermal field facilitate the generation of a data-based rheological model of the region.
In this study, we compute the long-term lithospheric strength and compare its spatial variation to observed seismicity patterns. We demonstrate how strength maxima within the crust (similar to 1 GPa) and upper mantle (> 2 GPa) occur at temperatures characteristic of the onset of crystal plasticity in those rocks (crust: 200-400 degrees C; mantle: similar to 600 degrees C), with almost all seismicity occurring in these regions. Correlation in the northern and southern forelands between crustal and lithospheric strengths and seismicity show different patterns of event distribution, reflecting their different tectonic settings. Seismicity in the plate boundary setting of the southern foreland corresponds to the integrated lithospheric strength, occurring mainly in the weaker domains surrounding the strong Adriatic plate. In the intraplate setting of the northern foreland, seismicity correlates to modelled crustal strength, and it mainly occurs in the weaker and warmer crust beneath the Upper Rhine Graben.
We, therefore, suggest that seismicity in the upper crust is linked to weak crustal domains, which are more prone to localise deformation promoting failure and, depending on the local properties of the fault, earthquakes at relatively lower levels of accumulated stress than their neighbouring stronger counterparts. Upper mantle seismicity at depths greater than modelled brittle conditions, can be either explained by embrittlement of the mantle due to grain-size sensitive deformation within domains of active or recent slab cooling, or by dissipative weakening mechanisms, such as thermal runaway from shear heating and/or dehydration reactions within an overly ductile mantle.
Results generated in this study are available for open access use to further discussions on the region.
Hundreds of basaltic plateau margins east of the Patagonian Cordillera are undermined by numerous giant slope failures. However, the overall extent of this widespread type of plateau collapse remains unknown and incompletely captured in local maps. To detect giant slope failures consistently throughout the region, we train two convolutional neural networks (CNNs), AlexNet and U-Net, with Sentinel-2 optical data and TanDEM-X topographic data on elevation, surface roughness, and curvature. We validated the performance of these CNNs with independent testing data and found that AlexNet performed better when learned on topographic data, and UNet when learned on optical data. AlexNet predicts a total landslide area of 12,000 km2 in a study area of 450,000 km2, and thus one of Earth's largest clusters of giant landslides. These are mostly lateral spreads and rotational failures in effusive rocks, particularly eroding the margins of basaltic plateaus; some giant landslides occurred along shores of former glacial lakes, but are least prevalent in Quaternary sedimentary rocks. Given the roughly comparable topographic, climatic, and seismic conditions in our study area, we infer that basalts topping weak sedimentary rocks may have elevated potential for large-scale slope failure. Judging from the many newly detected and previously unknown landslides, we conclude that CNNs can be a valuable tool to detect large-scale slope instability at the regional scale. However, visual inspection is still necessary to validate results and correctly outline individual landslide source and deposit areas.
Geostatistics as a subfield of statistics accounts for the spatial correlations encountered in many applications of, for example, earth sciences. Valuable information can be extracted from these correlations, also helping to address the often encountered burden of data scarcity. Despite the value of additional data, the use of geostatistics still falls short of its potential. This problem is often connected to the lack of user-friendly software hampering the use and application of geostatistics. We therefore present GSTools, a Python-based software suite for solving a wide range of geostatistical problems. We chose Python due to its unique balance between usability, flexibility, and efficiency and due to its adoption in the scientific community. GSTools provides methods for generating random fields; it can perform kriging, variogram estimation and much more. We demonstrate its abilities by virtue of a series of example applications detailing their use.
The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.
Wind erosion of agricultural soils affects their stock of essential elements for plants, like phosphorus (P). It is known that the composition of the eroded sediments varies with height, according to the size and density of the transported substances. Aim of this study was to analyze the concentration and enrichment ratios of P forms in sediments transported by the wind. A wind-tunnel study was performed on a sandy-and a sandy loam soil in order to measure P forms concentrations in the saltating sediments. P concentrations were also measured in the particulate matter (PM) of each soil, gained with the Easy Dust Generator. In both soils, inorganic-(Pi) and organic P (Po) were preferentially transported in PM, with enrichment ratios of 1.8 and 5.5, respectively. Nevertheless, a Pi/Po of 0.9 indicated that the accumulation of the minor Po in PM was more pronounced than Pi. This agrees with P-rich light and easily erodible organic compounds, almost exclusively accumulated in PM, and in relatively heavy and less erodible minerals, like apatites, in lower height sediments. Labile P (Pl) was preferentially transported in saltating sediments of both soils. This was attributed to the selective Bray & Kurtz I's extraction of the abundant inorganic P forms of these sediments. Total P (Pt) copied the transport trends of Pi, the major form. According to the transporting trends, Pi and Po would be re-sedimented at longer distances from the source than Pl. Outcomes become useful for modeling the influence of wind erosion on P cycling.
Boreal forests of Siberia play a relevant role in the global carbon cycle. However, global warming threatens the existence of summergreen larch-dominated ecosystems, likely enabling a transition to evergreen tree taxa with deeper active layers. Complex permafrost-vegetation interactions make it uncertain whether these ecosystems could develop into a carbon source rather than continuing atmospheric carbon sequestration under global warming. Consequently, shedding light on the role of current and future active layer dynamics and the feedbacks with the apparent tree species is crucial to predict boreal forest transition dynamics and thus for aboveground forest biomass and carbon stock developments. Hence, we established a coupled model version amalgamating a one-dimensional permafrost multilayer forest land-surface model (CryoGrid) with LAVESI, an individual-based and spatially explicit forest model for larch species (Larix Mill.), extended for this study by including other relevant Siberian forest species and explicit terrain. <br /> Following parameterization, we ran simulations with the coupled version to the near future to 2030 with a mild climate-warming scenario. We focus on three regions covering a gradient of summergreen forests in the east at Spasskaya Pad, mixed summergreen-evergreen forests close to Nyurba, and the warmest area at Lake Khamra in the southeast of Yakutia, Russia. Coupled simulations were run with the newly implemented boreal forest species and compared to runs allowing only one species at a time, as well as to simulations using just LAVESI. Results reveal that the coupled version corrects for overestimation of active layer thickness (ALT) and soil moisture, and large differences in established forests are simulated. We conclude that the coupled version can simulate the complex environment of eastern Siberia by reproducing vegetation patterns, making it an excellent tool to disentangle processes driving boreal forest dynamics.
The Ca Mau peninsula (CMP) is a key economic region in southern Vietnam. In recent decades, the high demand for water has increased the exploitation of groundwater, thus lowering the groundwater level and leading to risks of degradation, depletion, and land subsidence, as well as salinity intrusion in the groundwater of the whole Mekong Delta region. By using a finite element groundwater model with boundary expansion to the sea, we updated the latest data on hydrogeological profiles, groundwater levels, and exploitation. The basic model setup covers seven aquifers and seven aquitards. It is determined that the inflow along the coastline to the mainland is 39% of the total inflow. The exploitation of the study area in 2019 was 567,364 m(3)/day. The most exploited aquifers are the upper-middle Pleistocene (qp(2-3)) and the middle Pliocene (n(2)(2)), accounting for 63.7% and 24.6%, respectively; the least exploited aquifers are the upper Pleistocene and the upper Miocene, accounting for 0.35% and 0.02%, respectively. In the deeper aquifers, qp(2-3) and n(2)(2), the change in storage is negative due to the high exploitation rate, leading to a decline in the reserves of these aquifers. These groundwater model results are the calculations of groundwater reserves from the coast to the mainland in the entire system of aquifers in the CMP. This makes groundwater decision managers, stakeholders, and others more efficient in sustainable water resources planning in the CMP and Mekong Delta (MKD).
Channel steepness index, k(s), is a metric derived from the stream power model that, under certain conditions, scales with relative rock uplift rate. Channel steepness index is a property of rivers, which can be relatively easily extracted from digital elevation models (DEMs). As DEM data sets are widely available for Earth and are becoming more readily available for other planetary bodies, channel steepness index represents a powerful tool for interpreting tectonic processes. However, multiple approaches to calculate channel steepness index exist. From this several important questions arise; does choice of approach change the values of channel steepness index, can values be so different that choice of approach can influence the findings of a study, and are certain approaches better than others? With the aid of a synthetic river profile and a case study from the Sierra Nevada, California, we show that values of channel steepness index vary over orders of magnitude according to the methodology used in the calculation. We explore the limitations, advantages and disadvantages of the key approaches to calculating channel steepness index, and find that choosing an appropriate approach relies on the context of a study. Given these observations, it is important that authors acknowledge the methodology used to calculate channel steepness index, to ensure that results can be contextualised and reproduced.
Study region:
Ca Mau Province (CMP), Mekong Delta (MD), Vietnam.
Study focus:
Groundwater from deep aquifers is the most reliable source of freshwater in the MD but extensive overexploitation in the last decades led to the drop of hydraulic heads and negative environmental impacts. Therefore, a comprehensive groundwater investigation was conducted to evaluate its composition in the context of Quaternary marine transgression and regression cycles, geochemical processes as well as groundwater extraction.
New hydrological insights for the region:
The abundance of groundwater of Na-HCO3 type and distinct ion ratios, such as Na+/Cl-, indicate extensive freshwater intrusion in an initially saline hydrogeological system, with decreasing intensity from upper Pleistocene to deeper Miocene aquifers, most likely during the last marine regression phase 60-12 ka BP. Deviations from the conservative mixing line between the two endmembers seawater and freshwater are attributed to ion-exchange processes on mineral surfaces, making ion ratios in combination with a customized water type analysis a useful tool to distinguish between salinization and freshening processes. Elevated salinity in some areas is attributed to HCO3- generation by organic matter decomposition in marine sediments rather than to seawater intrusion. Nevertheless, a few randomly distributed locations show strong evidence of recent salinization in an early stage, which may be caused by the downwards migration of saline Holocene groundwater through natural and anthropogenic pathways into deep aquifers.
Controls on the deformation pattern (shortening mode and tectonic style) of orogenic forelands during lithospheric shortening remain poorly understood. Here, we use high-resolution 2D thermomechanical models to demonstrate that orogenic crustal thickness and foreland lithospheric thickness significantly control the shortening mode in the foreland. Pure-shear shortening occurs when the orogenic crust is not thicker than the foreland crust or thick, but the foreland lithosphere is thin (<70-80 km, as in the Puna foreland case). Conversely, simple-shear shortening, characterized by foreland underthrusting beneath the orogen, arises when the orogenic crust is much thicker. This thickened crust results in high gravitational potential energy in the orogen, which triggers the migration of deformation to the foreland under further shortening. Our models present fully thick-skinned, fully thin-skinned, and intermediate tectonic styles in the foreland. The first tectonics forms in a pure-shear shortening mode whereas the others require a simple-shear mode and the presence of thick (>similar to 4 km) sediments that are mechanically weak (friction coefficient <similar to 0.05) or weakened rapidly during deformation. The formation of fully thin-skinned tectonics in thick and weak foreland sediments, as in the Subandean Ranges, requires the strength of the orogenic upper lithosphere to be less than one-third as strong as that of the foreland upper lithosphere. Our models successfully reproduce foreland deformation patterns in the Central and Southern Andes and the Laramide province.
Building and changing a microbiome at will and maintaining it over hundreds of generations has so far proven challenging. Despite best efforts, complex microbiomes appear to be susceptible to large stochastic fluctuations. Current capabilities to assemble and control stable complex microbiomes are limited. Here, we propose a looped mass transfer design that stabilizes microbiomes over long periods of time. Five local microbiomes were continuously grown in parallel for over 114 generations and connected by a loop to a regional pool. Mass transfer rates were altered and microbiome dynamics were monitored using quantitative high-throughput flow cytometry and taxonomic sequencing of whole communities and sorted subcommunities. Increased mass transfer rates reduced local and temporal variation in microbiome assembly, did not affect functions, and overcame stochasticity, with all microbiomes exhibiting high constancy and increasing resistance. Mass transfer synchronized the structures of the five local microbiomes and nestedness of certain cell types was eminent. Mass transfer increased cell number and thus decreased net growth rates mu'. Subsets of cells that did not show net growth mu'SCx were rescued by the regional pool R and thus remained part of the microbiome. The loop in mass transfer ensured the survival of cells that would otherwise go extinct, even if they did not grow in all local microbiomes or grew more slowly than the actual dilution rate D would allow. The rescue effect, known from metacommunity theory, was the main stabilizing mechanism leading to synchrony and survival of subcommunities, despite differences in cell physiological properties, including growth rates.
Provisioning a sufficient stable source of food requires sound knowledge about current and upcoming threats to agricultural production. To that end machine learning approaches were used to identify the prevailing climatic and soil hydrological drivers of spatial and temporal yield variability of four crops, comprising 40 years yield data each from 351 counties in Germany. Effects of progress in agricultural management and breeding were subtracted from the data prior the machine learning modelling by fitting smooth non-linear trends to the 95th percentiles of observed yield data. An extensive feature selection approach was followed then to identify the most relevant predictors out of a large set of candidate predictors, comprising various soil and meteorological data. Particular emphasis was placed on studying the uniqueness of identified key predictors. Random Forest and Support Vector Machine models yielded similar although not identical results, capturing between 50% and 70% of the spatial and temporal variance of silage maize, winter barley, winter rapeseed and winter wheat yield. Equally good performance could be achieved with different sets of predictors. Thus identification of the most reliable models could not be based on the outcome of the model study only but required expert's judgement. Relationships between drivers and response often exhibited optimum curves, especially for summer air temperature and precipitation. In contrast, soil moisture clearly proved less relevant compared to meteorological drivers. In view of the expected climate change both excess precipitation and the excess heat effect deserve more attention in breeding as well as in crop modelling.
Climate change and increasing water demand in urban environments necessitate planning water utility companies' finances. Traditionally, methods to estimate the direct water utility business interruption costs (WUBIC) caused by droughts have not been clearly established. We propose a multi-driver assessment method. We project the water yield using a hydrological model driven by regional climate models under radiative forcing scenarios. We project water demand under stationary and non-stationary conditions to estimate drought severity and duration, which are linked with pricing policies recently adopted by the Sao Paulo Water Utility Company. The results showed water insecurity. The non-stationary trend imposed larger differences in the drought resilience financial gap, suggesting that the uncertainties of WUBIC derived from demand and climate models are greater than those associated with radiative forcing scenarios. As populations increase, proactively controlling demand is recommended to avoid or minimize reactive policy changes during future drought events, repeating recent financial impacts.
The detection of auto-fluorescence in phytogenic, hydrated amorphous silica depositions (phytoliths) has been found to be a promising approach to verify if phytoliths were burnt or not, especially in archaeological contexts. However, it is unknown so far at what temperature and how auto-fluorescence is induced in phytoliths. We used fluorescence microscopy, scanning electron microscope-energy dispersive X-ray spectroscopy (SEM-EDX), and Fourier transform infrared spectroscopy to analyze auto-fluorescence in modern phytoliths extracted from plant samples or in intact leaves of winter wheat. Leaves and extracted phytoliths were heated at different temperatures up to 600 degrees C. The aims of our experiments were i) to find out what temperature is needed to induce auto-fluorescence in phytoliths, ii) to detect temperature-dependent changes in the molecular structure of phytoliths related to auto-fluorescence, and iii) to derive a mechanistic understanding of auto-fluorescence in phytoliths. We found organic compounds associated with phytoliths to cause auto-fluorescence in phytoliths treated at temperatures below approx. 400 degrees C. In phytoliths treated at higher temperatures, i.e., 450 and 600 degrees C, phytolith auto-fluorescence was mainly caused by molecular changes of phytolith silica. Based on our results we propose that auto-fluorescence in phytoliths is caused by clusterization-triggered emissions, which are caused by overlapping electron clouds forming non-conventional chromophores. In phytoliths heated at temperatures above about 400 degrees C dihydroxylation and the formation of siloxanes result in oxygen clusters that serve as non-conventional chromophores in fluorescence events. Furthermore, SEM-EDX analyses revealed that extractable phytoliths were dominated by lumen phytoliths (62%) compared to cell wall phytoliths (38%). Our findings might be not only relevant in archaeological phytolith-based examinations, but also for studies on the temperature-dependent release of silicon from phytoliths and the potential of long-term carbon sequestration in phytoliths.