Refine
Year of publication
Document Type
- Article (2598)
- Doctoral Thesis (407)
- Postprint (138)
- Other (74)
- Review (51)
- Preprint (17)
- Monograph/Edited Volume (13)
- Conference Proceeding (12)
- Habilitation Thesis (11)
- Master's Thesis (6)
Language
- English (3329) (remove)
Keywords
- climate change (51)
- Holocene (44)
- erosion (27)
- Himalaya (26)
- permafrost (26)
- Climate change (23)
- remote sensing (23)
- Tibetan Plateau (22)
- Earthquake source observations (21)
- Pollen (20)
Institute
- Institut für Geowissenschaften (3329) (remove)
Profundal lake sediment cores are often interpreted in line with diverse and detailed sedimentological processes to infer paleoenvironmental conditions. The effects of frozen lake surfaces on terrigenous sediment deposition and how climate changes on the Tibetan Plateau are reflected in these lakes, however, is seldom discussed. A lake sediment core from Hala Lake (590 km(2)), northeastern Tibetan Plateau spanning the time interval from the Last Glacial Maximum to the present was investigated using high-resolution grain-size composition of lacustrine deposits. Seismic analysis along a north-south profile across the lake was used to infer the sedimentary setting within the lake basin. Periods of freezing and melting processes on the lake surface were identified by MODIS (MOD10A1) satellite data. End-member modeling of the grain size distribution allowed the discrimination between lacustrine, eolian and fluvial sediments. The dominant clay sedimentation (slack water type) during the global Last Glacial Maximum (LGM) reflects ice interceptions in long cold periods, in contrast to abundant eolian input during abrupt cold events. Therefore, fluvial and slack water sedimentation processes can indicate changes in the local paleoclimate during periods of the lake being frozen, when eolian input was minor. Inferred warm (i.e., similar to 22.7 and 19.5 cal. ka BP) and cold (i.e., similar to 11-9 and 3-1.5 cal. ka BP) spells have significant environmental impacts, not only in the regional realm, but they are also coherent with global-scale climate events. The eolian input generally follows the trend of the mid-latitude westerly wind dynamics in winter, contributing medium-sized sand to the lake center, deposited within the ice cover during icing and melting phases. Enhanced input was dominant during the Younger Dryas, Heinrich Event 1 and at around 8.2 ka, equivalent to the well-known events of the North Atlantic realm. (C) 2015 Elsevier Ltd. All rights reserved.
To find out the future of nowadays reef ecosystem turnover under the environmental stresses such as global warming and ocean acidification, analogue studies from the geologic past are needed. As a critical time of reef ecosystem innovation, the Permian-Triassic transition witnessed the most severe demise of Phanerozoic reef builders, and the establishment of modern style symbiotic relationships within the reef-building organisms. Being the initial stage of this transition, the Middle Permian (Capitanian) mass extinction coursed a reef eclipse in the early Late Permian, which lead to a gap of understanding in the post-extinction Wuchiapingian reef ecosystem, shortly before the radiation of Changhsingian reefs. Here, this thesis presents detailed biostratigraphic, sedimentological, and palaeoecological studies of the Wuchiapingian reef recovery following the Middle Permian (Capitanian) mass extinction, on the only recorded Wuchiapingian reef setting, outcropping in South China at the Tieqiao section.
Conodont biostratigraphic zonations were revised from the Early Permian Artinskian to the Late Permian Wuchiapingian in the Tieqiao section. Twenty main and seven subordinate conodont zones are determined at Tieqiao section including two conodont zone below and above the Tieqiao reef complex. The age of Tieqiao reef was constrained as early to middle Wuchiapingian.
After constraining the reef age, detailed two-dimensional outcrop mapping combined with lithofacies study were carried out on the Wuchiapingian Tieqiao Section to investigate the reef growth pattern stratigraphically as well as the lateral changes of reef geometry on the outcrop scale. Semi-quantitative studies of the reef-building organisms were used to find out their evolution pattern within the reef recovery. Six reef growth cycles were determined within six transgressive-regressive cycles in the Tieqiao section. The reefs developed within the upper part of each regressive phase and were dominated by different biotas. The timing of initial reef recovery after the Middle Permian (Capitanian) mass extinction was updated to the Clarkina leveni conodont zone, which is earlier than previous understanding. Metazoans such as sponges were not the major components of the Wuchiapingian reefs until the 5th and 6th cycles. So, the recovery of metazoan reef ecosystem after the Middle Permian (Capitanian) mass extinction was obviously delayed. In addition, although the importance of metazoan reef builders such as sponges did increase following the recovery process, encrusting organisms such as Archaeolithoporella and Tubiphytes, combined with microbial carbonate precipitation, still played significant roles to the reef building process and reef recovery after the mass extinction.
Based on the results from outcrop mapping and sedimentological studies, quantitative composition analysis of the Tieqiao reef complex were applied on selected thin sections to further investigate the functioning of reef building components and the reef evolution after the Middle Permian (Capitanian) mass extinction. Data sets of skeletal grains and whole rock components were analyzed. The results show eleven biocommunity clusters/eight rock composition clusters dominated by different skeletal grains/rock components. Sponges, Archaeolithoporella and Tubiphytes were the most ecologically important components within the Wuchiapingian Tieqiao reef, while the clotted micrites and syndepositional cements are the additional important rock components for reef cores. The sponges were important within the whole reef recovery. Tubiphytes were broadly distributed in different environments and played a key-role in the initial reef communities. Archaeolithoporella concentrated in the shallower part of reef cycles (i.e., the upper part of reef core) and was functionally significant for the enlargement of reef volume.
In general, the reef recovery after the Middle Permian (Capitanian) mass extinction has some similarities with the reef recovery following the end-Permian mass extinction. It shows a delayed recovery of metazoan reefs and a stepwise recovery pattern that was controlled by both ecological and environmental factors. The importance of encrusting organisms and microbial carbonates are also similar to most of the other post-extinction reef ecosystems. These findings can be instructive to extend our understanding of the reef ecosystem evolution under environmental perturbation or stresses.
Following the Middle Permian (Capitanian) mass extinction there was a global ‘reef eclipse’, and this event had an important role in the Paleozoic-Mesozoic transition of reef ecosystems. Furthermore, the recovery pattern of reef ecosystems in the Wuchiapingian of South China, before the radiation of Changhsingian reefs, is poorly understood. Here, we present a detailed sedimentological account of the Tieqiao section, South China, which records the only known Wuchiapingian reef setting from South China. Six reef growing phases were identified within six transgressive-regressive cycles. The cycles represent changes of deposition in a shallow basin to a subtidal outer platform setting, and the reefal build-ups are recorded in the shallowest part of the cycles above wave base in the euphotic zone. Our results show that the initial reef recovery started from the shallowing up part of the 1st cycle, within the Clarkina leveni conodont zone, which is two conodont zones earlier than previously recognized. In addition, even though metazoans, such as sponges, do become important in the development of the reef bodies, they are not a major component until later in the Wuchiapingian in the 5th and 6th transgressive-regressive cycles. This suggests a delayed recovery of metazoan reef ecosystems following the Middle Permian extinction. Furthermore, even though sponges do become abundant within the reefs, it is the presence and growth of the encrusters Archaeolithoporella and Tubiphytes and abundance of microbial micrites that play an important role in stabilizing the reef structures that form topographic highs.
Molecular Paleoclimate Reconstructions over the Last 9 ka from a Peat Sequence in South China
(2016)
To achieve a better understanding of Holocene climate change in the monsoon regions of China, we investigated the molecular distributions and carbon and hydrogen isotope compositions delta C-13 and delta D values) of long-chain n-alkanes in a peat core from the Shiwangutian SWGT) peatland, south China over the last 9 ka. By comparisons with other climate records, we found that the delta C-13 values of the long-chain n-alkanes can be a proxy for humidity, while the dD values of the long-chain n-alkanes primarily recorded the moisture source dD signal during 9-1.8 ka BP and responded to the dry climate during 1.8-0.3 ka BP. Together with the average chain length ACL) and the carbon preference index CPI) data, the climate evolution over last 9 ka in the SWGT peatland can be divided into three stages. During the first stage 9-5 ka BP), the delta C-13 values were depleted and CPI and Paq values were low, while ACL values were high. They reveal a period of warm and wet climate, which is regarded as the Holocene optimum. The second stage 5-1.8 ka BP) witnessed a shift to relatively cool and dry climate, as indicated by the more positive delta C-13 values and lower ACL values. During the third stage 1.8-0.3 ka BP), the delta C-13, delta D, CPI and Paq values showed marked increase and ACL values varied greatly, implying an abrupt change to cold and dry conditions. This climate pattern corresponds to the broad decline in Asian monsoon intensity through the latter part of the Holocene. Our results do not support a later Holocene optimum in south China as suggested by previous studies.
Knowledge of the possible impacts of early diagenesis on lipid biomarkers in geologic settings is important for robust applications of lipid proxies for paleoclimate reconstructions. In this study, molecular distributions and carbon isotopic compositions of lipids were compared in two particle-size fractions (<0.3 mm and >0.3 mm) of twelve surface peat samples collected from Dajiuhu peatland, central China. The average chain length (ACL) values of long-chain n-alkanes, n-fatty alcohols, n-fatty acids and n-alkan-2-ones show no significant differences between the finer and coarser fractions. In contrast, the carbon preference index values of long-chain n-alkanes, n-fatty alcohols and n-alkan-2-ones have relatively smaller values in the finer fractions than in the coarser ones. Stanols were also more abundant in the finer fractions. In addition, the delta C-13 values of odd-numbered n-alkanes (C-23-C-33) were generally less negative in the finer fractions. Our results indicate that (1) the finer fractions probably experienced stronger degradation than the coarser fractions; (2) the less negative delta(CC)-C-13 values of odd-numbered n-alkanes (C-23-C-33) in the finer fractions were largely a result of greater heterotrophic reworking during degradation; (3) ACL values of long-chain n-alkyl lipids (n-alkanes, n-fatty alcohols and n-fatty acids, n-alkan-2-ones) appear to be reliable proxies to trace lipid sources and their associated paleoenvironmental signals in peat deposits.
Pollen records from large lakes have been used for quantitative palaeoclimate reconstruction, but the influences that lake size (as a result of species-specific variations in pollen dispersal patterns that smaller pollen grains are more easily transported to lake centre) and taphonomy have on these climatic signals have not previously been systematically investigated. We introduce the concept of pollen source area to pollen-based climate calibration using the north-eastern Tibetan Plateau as our study area. We present a pollen data set collected from large lakes in the arid to semi-arid region of central Asia. The influences that lake size and the inferred pollen source areas have on pollen compositions have been investigated through comparisons with pollen assemblages in neighbouring lakes of various sizes. Modern pollen samples collected from different parts of Lake Donggi Cona (in the north-eastern part of the Tibetan Plateau) reveal variations in pollen assemblages within this large lake, which are interpreted in terms of the species-specific dispersal and depositional patterns for different types of pollen, and in terms of fluvial input components. We have estimated the pollen source area for each lake individually and used this information to infer modern climate data with which to then develop a modern calibration data set, using both the multivariate regression tree (MRT) and weighted-averaging partial least squares (WA-PLS) approaches. Fossil pollen data from Lake Donggi Cona have been used to reconstruct the climate history of the north-eastern part of the Tibetan Plateau since the Last Glacial Maximum (LGM). The meanannual precipitation was quantitatively reconstructed using WA-PLS: extremely dry conditions are found to have dominated the LGM, with annual precipitation of around 100 mm, which is only 32% of present-day precipitation. A gradually increasing trend in moisture conditions during the Late Glacial is terminated by an abrupt reversion to a dry phase that lasts for about 1000 yr and coincides with "Heinrich event 1" in the North Atlantic region. Subsequent periods corresponding to the Bolling/Allerod interstadial, with annual precipitation (P-ann) of about 350 mm, and the Younger Dryas event (about 270 mm P-ann) are followed by moist conditions in the early Holocene, with annual precipitation of up to 400 mm. A drier trend after 9 cal. ka BP is followed by a second wet phase in the middle Holocene, lasting until 4.5 cal. ka BP. Relatively steady conditions with only slight fluctuations then dominate the late Holocene, resulting in the present climatic conditions. The climate changes since the LGM have been primarily driven by deglaciation and fluctuations in the intensity of the Asian summer monsoon that resulted from changes in the Northern Hemisphere summer solar insolation, as well as from changes in the North Atlantic climate through variations in the circulation patterns and intensity of the westerlies.
The past climate in central Asia, and especially on the Tibetan Plateau (TP), is of great importance for an understanding of global climate processes and for predicting the future climate. As a major influence on the climate in this region, the Asian Summer Monsoon (ASM) and its evolutionary history are of vital importance for accurate predictions. However, neither the evolutionary pattern of the summer monsoon nor the driving mechanisms behind it are yet clearly understood. For this research, I first synthesized previously published Late Glacial to Holocene climatic records from monsoonal central Asia in order to extract the general climate signals and the associated summer monsoon intensities. New climate and vegetation sequences were then established using improved quantitative methods, focusing on fossil pollen records recovered from Tibetan lakes and also incorporating new modern datasets. The pollen-vegetation and vegetation-climate relationships on the TP were also evaluated in order to achieve a better understanding of fossil pollen records. The synthesis of previously published moisture-related palaeoclimate records in monsoonal central Asia revealed generally different temporal patterns for the two monsoonal subsystems, i.e. the Indian Summer Monsoon (ISM) and East Asian Summer Monsoon (EASM). The ISM appears to have experienced maximum wet conditions during the early Holocene, while many records from the area affected by the EASM indicate relatively dry conditions at that time, particularly in north-central China where the maximum moisture levels occurred during the middle Holocene. A detailed consideration of possible driving factors affecting the summer monsoon, including summer solar insolation and sea surface temperatures, revealed that the ISM was primarily driven by variations in northern hemisphere solar insolation, and that the EASM may have been constrained by the ISM resulting in asynchronous patterns of evolution for these two subsystems. This hypothesis is further supported by modern monsoon indices estimated using the NCEP/NCAR Reanalysis data from the last 50 years, which indicate a significant negative correlation between the two summer monsoon subsystems. By analogy with the early Holocene, intensification of the ISM during coming decades could lead to increased aridification elsewhere as a result of the asynchronous nature of the monsoon subsystems, as can already be observed in the meteorological data from the last 15 years. A quantitative climate reconstruction using fossil pollen records was achieved through analysis of sediment core recovered from Lake Donggi Cona (in the north-eastern part of the TP) which has been dated back to the Last Glacial Maximum (LGM). A new data-set of modern pollen collected from large lakes in arid to semi-arid regions of central Asia is also presented herein. The concept of "pollen source area" was introduced to modern climate calibration based on pollen from large lakes, and was applied to the fossil pollen sequence from Lake Donggi Cona. Extremely dry conditions were found to have dominated the LGM, and a subsequent gradually increasing trend in moisture during the Late Glacial period was terminated by an abrupt reversion to a dry phase that lasted for about 1000 years and coincided with the first Heinrich Event of the northern Atlantic region. Subsequent periods corresponding to the warm Bølling-Allerød period and the Younger Dryas cold event were followed by moist conditions during the early Holocene, with annual precipitation of up to about 400 mm. A slightly drier trend after 9 cal ka BP was then followed by a second wet phase during the middle Holocene that lasted until 4.5 cal ka BP. Relatively steady conditions with only slight fluctuations then dominated the late Holocene, resulting in the present climatic conditions. In order to investigate the relationship between vegetation and climate, temporal variations in the possible driving factors for vegetation change on the northern TP were examined using a high resolution late Holocene pollen record from Lake Kusai. Moving-window Redundancy Analyses (RDAs) were used to evaluate the correlations between pollen assemblages and individual sedimentary proxies. These analyses have revealed frequent fluctuations in the relative abundances of alpine steppe and alpine desert components, and in particular a decrease in the total vegetation cover at around 1500 cal a BP. The climate was found to have had an important influence on vegetation changes when conditions were relatively wet and stable. However, after the 1500 cal a BP threshold in vegetation cover was crossed the vegetation appears to have been affected more by extreme events such as dust storms or fluvial erosion than by the general climatic trends. In addition, pollen spectra over the last 600 years have been revealed by Procrustes analysis to be significantly different from those recovered from older samples, which is attributed to an increased human impact that resulted in unprecedented changes to the composition of the vegetation. Theoretical models that have been developed and widely applied to the European area (i.e. the Extended R-Value (ERV) model and the Regional Estimates of Vegetation Abundance from Large Sites (REVEALS) model) have been applied to the high alpine TP ecosystems in order to investigate the pollen-vegetation relationships, as well as for quantitative reconstructions of vegetation abundance. The modern pollen–vegetation relationships for four common pollen species on the TP have been investigated using Poaceae as the reference taxa. The ERV Submodel 2 yielded relatively high PPEs for the steppe and desert taxa (Artemisia Chenopodiaceae), and low PPEs for the Cyperaceae that are characteristic of the alpine Kobresia meadows. The plant abundances on the central and north-eastern TP were quantified by applying these PPEs to four post-Late Glacial fossil pollen sequences. The reconstructed vegetation assemblages for the four pollen sequences always yielded smaller compositional species turnovers than suggested by the pollen spectra, indicating that the strength of the previously-reported vegetation changes may therefore have been overestimated. In summary, the key findings of this thesis are that (a) the two ASM subsystems show asynchronous patterns during both the Holocene and modern time periods, (b) fossil pollen records from large lakes reflect regional signals for which the pollen source areas need to be taken into account, (c) climate is not always the main driver for vegetation change, and (d) previously reported vegetation changes on the TP may have been overestimated because they ignored inter-species variations in pollen productivity.
The concept of a Global Monsoon (GM) has been proposed based on modern precipitation observations, but its application over a wide range of temporal scales is still under debate. Here, we present a synthesis of 268 continental paleo-moisture records collected from monsoonal systems in the Eastern Hemisphere, including the East Asian Monsoon (EAsM), the Indian Monsoon (IM), the East African Monsoon (EAfM), and the Australian Monsoon (AuM) covering the last 18,000 years. The overall pattern of late Glacial to Holocene moisture change is consistent with those inferred from ice cores and marine records. With respect to the last 10,000 years (10 ka), i.e. a period that has high spatial coverage, a Fuzzy c-Means clustering analysis of the moisture index records together with "Xie-Beni" index reveals four clusters of our data set. The paleoclimatic meaning of each cluster is interpreted considering the temporal evolution and spatial distribution patterns. The major trend in the tropical AuM, EAfM, and IM regions is a gradual decrease in moisture conditions since the early Holocene. Moisture changes in the EAsM regions show maximum index values between 8 and 6 ka. However, records located in nearby subtropical areas, i.e. in regions not influenced by the intertropical convergence zone, show an opposite trend compared to the tropical monsoon regions (AuM, EAfM and IM), i.e. a gradual increase. Analyses of modern meteorological data reveal the same spatial patterns as in the paleoclimate records such that, in times of overall monsoon strengthening, lower precipitation rates are observed in the nearby subtropical areas. We explain this pattern as the effect of a strong monsoon circulation suppressing air uplift in nearby subtropical areas, and hence hindering precipitation. By analogy to the modern system, this would mean that during the early Holocene strong monsoon period, the intensified ascending airflows within the monsoon domains led to relatively weaker ascending or even descending airflows in the adjacent subtropical regions, resulting in a precipitation deficit compared to the late Holocene. Our conceptual model therefore integrates regionally contrasting moisture changes into the Global Monsoon hypothesis. (C) 2017 Elsevier Ltd. All rights reserved.
Previous studies based on fossil pollen data have reported significant changes in vegetation on the alpine Tibetan Plateau during the Holocene. However, since the relative proportions of fossil pollen taxa are largely influenced by individual pollen productivities and the dispersal characteristics, such inferences on vegetation have the potential to be considerably biased. We therefore examined the modern pollen-vegetation relationships for four common pollen species on the Tibetan Plateau, using Extended R-value (ERV) models. Assuming an average radius of 100 m for the sampled lakes, we estimated the relevant source area of pollen (RSAP) to be 2200 m (which represents the distance from the lake). Using Poaceae as the reference taxa (Pollen Productivity Estimate, PPE = 1), ERV Submodel 2 derived relative high PPEs for the steppe and desert taxa: 2.079 +/- 0.432 for Artemisia and 5.379 +/- 1.077 for Chenopodiaceae. Low PPEs were estimated for the Cyperaceae (1.036 +/- 0.012). whose plants are characteristic of the alpine Kobresia meadows. Applying these PPEs to four fossil pollen sequences since the Late Glacial, the plant abundances on the central and north-eastern Tibetan Plateau were quantified using the "Regional Estimates of Vegetation Abundance from Large Sites" (REVEALS) model. The proportions of Artemisia and Chenopodiaceae were greatly reduced compared to their original pollen percentages in the reconstructed vegetation, owing to their high productivities and their dispersal characteristics, while Cyperaceae showed a relative increase in the vegetation reconstruction. The reconstructed vegetation assemblages of the four pollen sequence sites always yielded smaller compositional species turnovers than suggested by the pollen spectra, as revealed by Detrended Canonical Correspondence Analyses (DCCA) of the Holocene sections. The strength of the previously reported vegetation changes may therefore have been overestimated, which indicates the importance of taking into account pollen-vegetation relationships when discussing the potential drivers (such as climate, land use, atmospheric CO(2) concentrations) and implications (such as for land surface-climate feedbacks, carbon storage, and biodiversity) of vegetation change.
Lacustrine sediments have been widely used to investigate past climatic and environmental changes on millennial to seasonal time scales. Sedimentary archives of lakes in mountainous regions may also record non-climatic events such as earthquakes. We argue herein that a set of 64 annual laminae couplets reconciles a stratigraphically inconsistent accelerator mass spectrometry (AMS) C-14 chronology in a similar to 4-m-long sediment core from Lake Mengda, in the north-eastern Tibetan Plateau. The laminations suggest the lake was formed by a large landslide, triggered by the 1927 Gulang earthquake (M = 8.0). The lake sediment sequence can be separated into three units based on lithologic, sedimentary, and isotopic characteristics. Starting from the bottom of the sequence, these are: (1) unweathered, coarse, sandy valley-floor deposits or landslide debris that pre-date the lake, (2) landslide-induced, fine-grained soil or reworked landslide debris with a high organic content, and (3) lacustrine sediments with low organic content and laminations. These annual laminations provide a high-resolution record of anthropogenic and environmental changes during the twentieth century, recording enhanced sediment input associated with two phases of construction activities. The high mean sedimentation rates of up to 4.8 mm year(-1) underscore the potential for reconstructing such distinct sediment pulses in remote, forested, and seemingly undisturbed mountain catchments.
Temporally changing drivers for late-Holocene vegetation changes on the northern Tibetan Plateau
(2012)
Fossil pollen records have been widely used as indicators of past changes in vegetation and variations in climate. The driving mechanisms behind these vegetation changes have, however, remained unclear. In order to evaluate vegetation changes that have occurred in the northern part of the Tibetan Plateau and the possible drivers behind these changes, we have applied a moving-window Redundancy Analysis (RDA) to high resolution (10-15 years) pollen and sedimentary data from Lake Kusai covering the last 3770 years. Our analyses reveal frequent fluctuations in the relative abundances of alpine steppe and alpine desert components. The sedimentary proxies (including total organic carbon content, total inorganic carbon content, and "end-member" indices from grain-size analyses) that explain statistically some of the changes in the pollen assemblage vary significantly with time, most probably reflecting multiple underlying driving processes. Climate appears to have had an important influence on vegetation changes when conditions were relatively wet and stable. However, a gradual decrease in vegetation cover was identified after 1500 cal a BP, after which the vegetation appears to have been affected more by extreme events such as dust-storms or fluvial erosion than by general climatic trends. Furthermore, pollen spectra over the last 600 years are shown by Procrustes analysis to be statistically different from those recovered from older samples, which we attribute to increased human impact that resulted in unprecedented changes to the vegetation composition. Overall, changes in vegetation and climate on the northern part of the Tibetan Plateau appear to have roughly followed the evolution of the Asian Summer Monsoon. After taking into account the highly significant millennial (1512 years) periodicity revealed by time-series analysis, the regional vegetation and climate changes also show variations that appear to match variations in the mid-latitude westerlies.
Inorganic minerals form a major component of lacustrine sediments and have the potential to reveal detailed information on previous climatic and hydrological conditions. The ability to extract such information however, has been restricted by a limited understanding of the relationships between minerals and the environment. In an attempt to fill in this gap in our knowledge, 146 surface sediment samples have been investigated from 146 lakes on the Tibetan Plateau. The mineral compositions derived from these samples by X-Ray Diffraction (XRD) were used to examine the relationships between mineral compositions and the environmental variables determined for each site. Statistical techniques including Multivariate regression trees (MRT) and Redundancy Analysis (RDA), based on the mineral spectra and environmental variables, reveal that the electrical conductivity (EC) and Mg/Ca ratios of lake water are the most important controls on the composition of endogenic minerals. No endogenic minerals precipitate under hyper-fresh water conditions (EC lower than 0.13 mS/cm), with calcite commonly forming in water with EC values above 0.13 mS/cm. Between EC values of 0.13 and 26 mS/cm the mineral composition of lake sediments can be explained in terms of variations in the Mg/Ca ratio: calcite dominates at Mg/Ca ratios of less than 33, whereas aragonite commonly forms when the ratio is greater than 33. Where EC values are between 26 and 39 mS/cm, monohydrocalcite precipitates together with calcite and aragonite; above 39 mS/cm, gypsum and halite commonly form. Information on the local geological strata indicates that allogenic (detrital) mineral compositions are primarily influenced by the bedrock compositions within the catchment area. By applying these relationships to the late glacial and Holocene mineral record from Chaka Salt Lake, five lake stages have been identified and their associated EC conditions inferred. The lake evolved from a freshwater lake during the late glacial (before 11.4 cal. ka BP) represented by the lowest EC values (< 0.13 mS/cm), to a saline lake with EC values slightly higher than 39 mS/cm during the early and mid Holocene (ca. 11.4-5.3 cal. ka BP), and finally to a salt lake (after 5.3 cal. ka BP). These results illustrate the utility of our mineral-environmental model for the quantitative reconstruction of past environmental conditions from lake sediment records.
We hypothesized that at the very beginning of terrestrial ecosystem development, airborne testate amoebae play a pivotal role in facilitating organismic colonization and related soil processes. We, therefore, analyzed size and quantity of airborne testate amoebae and immigration and colonization success of airborne testate amoebae on a new land surface (experimental site "Chicken Creek", artificial post-mining water catchment). Within an altogether 91-day exposure of 70 adhesive traps, 12 species of testate amoebae were identified to be of airborne origin. Phryganella acropodia (51% of all individuals found, diameter about 35-45 mu m) and Centropyxis sphagnicola (23% of all individuals found, longest axis about 55-68 mu m), occurred most frequently in the adhesive traps. We extrapolated an aerial amoeba deposition of 61 individuals d(-1) m(-2) (living and dead individuals combined). Although it would be necessary to have a longer sequence (some additional years), our analysis of the "target substrate" of aerial immigration (catchment site) may point to a shift from a stochastic (variable) beginning of community assembly to a more deterministic (stable) course. This shift was assigned to an age of seven years of initial soil development. Although experienced specialists are necessary to conduct these time-consuming studies, the presented data suggest that terrestrial amoebae are suitable indicators for initial ecosystem development and utilization.
Testate amoebae with self-secreted siliceous shell platelets ("idiosomes") play an important role in terrestrial silicon (Si) cycles. In this context, Si-dependent culture growth dynamics of idiosomic testate amoebae are of interest. Clonal cultures of idiosomic testate amoebae were analyzed under three different Si concentrations: low (50 mu mol L-1), moderate/site-specific (150 mu mol L-1) and high Si supply (500 mu mol L-1). Food (Saccharomyces cerevisiae) was provided in surplus. (i) Shell size of four different clones of idiosomic testate amoebae either decreased (Trinema galeata, Euglypha filifera cf.), increased (E. rotunda cf.), or did not change (E. rotunda) under the lowest Si concentration (50 mu mol Si L-1). (ii) Culture growth of idiosomic Euglypha rotunda was dependent on Si concentration. The more Si available in the culture medium, the earlier the entry into exponential growth phase. (iii) Culture growth of idiosomic Euglypha rotunda was dependent on origin of inoculum. Amoebae previously cultured under a moderate Si concentration revealed highest sustainability in consecutive cultures. Amoebae derived from cultures with high Si concentrations showed rapid culture growth which finished early in consecutive cultures. (iv) Si (diluted in the culture medium) was absorbed by amoebae and fixed in the amoeba shells resulting in decreased Si concentrations. (C) 2016 Elsevier GmbH. All rights reserved.
Process life cycle assessment (PLCA) is widely used to quantify environmental flows associated with the manufacturing of products and other processes. As PLCA always depends on defining a system boundary, its application involves truncation errors. Different methods of estimating truncation errors are proposed in the literature; most of these are based on artificially constructed system complete counterfactuals. In this article, we review the literature on truncation errors and their estimates and systematically explore factors that influence truncation error estimates. We classify estimation approaches, together with underlying factors influencing estimation results according to where in the estimation procedure they occur. By contrasting different PLCA truncation/error modeling frameworks using the same underlying input-output (I-O) data set and varying cut-off criteria, we show that modeling choices can significantly influence estimates for PLCA truncation errors. In addition, we find that differences in I-O and process inventory databases, such as missing service sector activities, can significantly affect estimates of PLCA truncation errors. Our results expose the challenges related to explicit statements on the magnitude of PLCA truncation errors. They also indicate that increasing the strictness of cut-off criteria in PLCA has only limited influence on the resulting truncation errors. We conclude that applying an additional I-O life cycle assessment or a path exchange hybrid life cycle assessment to identify where significant contributions are located in upstream layers could significantly reduce PLCA truncation errors.
We present the first high-resolution trace element (Mg/Ca, Sr/Ca, Ba/Ca) record from a stalagmite in southwestern Romania covering the last 3.6 ka, which provides the potential for quantitative climate reconstruction. Precise age control is based on three independent dating methods, in particular for the last 250 yr, where chemical lamina counting is combined with the identification of the 20th century radiocarbon bomb peak and Th-230/U dating. Long-term cave monitoring and model simulations of drip water and speleothem elemental variability indicate that precipitation-related processes are the main drivers of speleothem Mg/Ca ratios. Calibration against instrumental climate data shows a significant anti-correlation of speleothem Mg/Ca ratios with autumn/winter (October to March) precipitation (r = -0.61, p < 0.01), which is statistically robust when considering age uncertainties and auto-correlation. This relationship is used to develop a quantitative reconstruction of autumn/winter precipitation. During the late Holocene, our data suggest a heterogeneous pattern of past regional winter hydroclimate in the Carpathian/Balkan realm, along with intermittent weakening of the dominant influence of North Atlantic forcing. In agreement with other regional paleo-hydrological reconstructions, the observed variability reveals periodically occurring strong NW-SE hydro-climate gradients. We hypothesize, that this pattern is caused by shifts of the eastern edge of the area of influence of the NAO across central eastern Europe due to the confluence of North Atlantic forcing, and other climatic features such as the East Atlantic/Western Russia (EAWR) pattern. (C) 2018 Elsevier B.V. All rights reserved.
In this study we investigated the effect of different land use options (wildlife versus livestock) on species richness of plants and reptiles in the protected Kgalagadi Transfrontier Park (KTP) versus adjacent non-protected farmland within the same savanna habitat type (Aoub dune veld) in the southern Kalahari, South Africa. Our results show that both plant and reptile species richness as well as plant cover and reptile abundance was significantly higher in the protected KTP than in the non-protected farmland. The higher proportion of shrub but lower proportions of perennial grass cover, herb cover, and herb species richness in the farmland can be explained by higher stocking rates and the differences in feeding behaviour between native wild ungulates (e.g. Antidorcas marsupialis, Oryx gazella) and livestock (mainly sheep). The reptile's prey availability and microhabitats (perennial grass tussocks and rodent burrows) for thermoregulation and protection against predators were significantly lower in the farmland. To conclude, our results clearly show that long term effects of different land use options (wildlife in protected KTP versus extensive livestock production in the non-protected farmland) even within the same habitat type have led to significant changes in vegetation composition, availability of microhabitat structures and in the reptile community.
Although the effects of grazing-induced savannah degradation on animal diversity are well documented, knowledge of how they affect space use or responding behaviour remains poor. In this study, we analysed space use of the spotted sand lizard (Pedioplanis l. lineoocellata) in degraded versus nondegraded habitats of southern Kalahari savannah habitats. Lizards were radio tracked, daily movement distances recorded and home range sizes calculated. In degraded Kalahari savannah habitats where plant diversity and perennial grass cover are low but shrub cover high, P. lineoocellata moves larger distances (40.88 +/- 6.42 m versus 27.43 +/- 5.08 m) and occupies larger home ranges (646.64 +/- 244.84 m(2) versus 209.15 +/- 109.84 m(2)) than in nondegraded habitats (high plant diversity, high perennial grass cover and low shrub cover). We assume that this increase in daily movement distances and home range sizes is a behavioural plasticity to limited food resources in degraded savannah habitats. Although P. lineoocellata is able to adjust to resource-poor savannah habitats, the increase in the lizard's movement activities is likely to result in a higher predation risk. This is supported by the lower availability of protective vegetation i.e. perennial grass cover. Hence, we conclude that despite behavioural plasticity of P. lineoocellata, overgrazing has a severe negative impact on the space use of P. lineoocellata.
Volcano seismology
(2001)
Ambient vibration techniques are promising methods for assessing the subsurface structure, in particular the shear-wave velocity profile (V-s). They are based on the dispersion property of surface waves in layered media. Therefore, the penetration depth is intrinsically linked to the energy content of the sources. For ambient vibrations, the spectral content extends in general to lower frequency when compared to classical artificial sources. Among available methods for processing recorded signals, we focus here on the spatial autocorrelation method. For stationary wavefields, the spatial autocorrelation is mathematically related to the frequency-dependent wave velocity c(omega). This allows the determination of the dispersion curve of traveling surface waves, which, in turn, is linked to the V-s profile. Here, we propose a direct inversion scheme for the observed autocorrelation curves to retrieve, in a single step, the V-s profile. The powerful neighborhood algorithm is used to efficiently search for all solutions in an n- dimensional parameter space. This approach has the advantage of taking into account the existing uncertainty over the measured curves, thus generating all V-s profiles that fit the data within their experimental errors. A preprocessing tool is also developed to estimate the validity of the autocorrelation curves and to reject parts of them if necessary before starting the inversion itself. We present two synthetic cases to test the potential of the method: one with ideal autocorrelation curves and another with autocorrelation curves computed from simulated ambient vibrations. The latter case is more realistic and makes it possible to figure out the problems that may be encountered in real experiments. The V-s profiles are correctly retrieved up to the depth of the first major velocity contrast unless low-velocity zones are accepted. We demonstrate that accepting low-velocity zones in the parameterization has a dramatic influence on the result of the inversion, with a considerable increase in the nonuniqueness of the problem. Finally, a real data set is processed with the same method
The variation of Rayleigh ellipticity versus frequency is gaining popularity in site characterization. It becomes a necessary observable to complement dispersion curves when inverting shear wave velocity profiles. Various methods have been proposed so far to extract polarization from ambient vibrations recorded on a single three-component station or with an array of three-component sensors. If only absolute values were recovered 10 yr ago, new array-based techniques were recently proposed with enhanced efficiencies providing also the ellipticity sign. With array processing, higher-order modes are often detected even in the ellipticity domain. We suggest to explore the properties of a high-resolution beamforming where radial and vertical components are explicitly included. If N is the number of three-component sensors, 2N x 2N cross-spectral density matrices are calculated for all presumed directions of propagation. They are built with N radial and N vertical channels. As a first approach, steering vectors are designed to fit with Rayleigh wave properties: the phase shift between radial and vertical components is either -Pi/2 or Pi/2. We show that neglecting the ellipticity tilt due to attenuation has only minor effects on the results. Additionally, we prove analytically that it is possible to retrieve the ellipticity value from the usual maximization of the high-resolution beam power. The method is tested on synthetic data sets and on experimental data. Both are reference sites already analysed by several authors. A detailed comparison with previous results on these cases is provided.
The specific chemical composition of monazite in shear zones is controlled by the syndeformation dissolution-precipitation reactions of the rock-forming minerals. This relation can be used for dating deformation, even when microfabric characteristics like shape preferred orientation or intracrystalline deformation of monazite itself are missing. Monazite contemporaneously formed in and around the shear zones may have different compositions. These depend on the local chemical context rather than reflecting successive crystallization episodes of monazite.
This is demonstrated in polymetamorphic, mylonitic high-pressure (HP) garnet-kyanite granulites of the Alpine Sidironero Complex (Rhodope UHP terrain, Northern Greece). The studied mylonitic rocks escaped from regional migmatization at 40-36 Ma and from subsequent shearing through cooling until 36 Ma. In-situ laser-ablation split-stream inductively-coupled plasma mass spectrometry (LASS) analyses have been carried out on monazite from micro-scale shear zones, from pre-mylonitic microlithons as well as of monazite inclusions in relictic minerals complimented by U-Pb data on rutile and Rb-Sr data of biotite.
Two major metamorphic episodes, Mesozoic and Cenozoic, are constrained. Chemical compositions, isotopic characteristics and apparent ages systematically vary among monazite of four different microfabric domains (I-IV). Within three pre-mylonitic domains (inclusions in (I) pre-mylonitic kyanite and (II) garnet porphyroclasts, and (III) in pre-mylonitic microlithons) monazite yields ages of ca. 130-150 Ma for HP-granulite metamorphism, in line with previous geochronological results in the area. Patchy alteration of the pre-mylonitic monazite by intra-grain dissolution-precipitation processes variably increased negative Eu anomaly and reduced the HREE contents. The apparent age of this altered monazite is reduced.
Monazite in the syn-mylonitic shear bands (IV) differs in chemical composition from unaltered and altered monazite of the three pre-mylonitic domains by having a significantly more pronounced negative Eu anomaly, a flatter HREE pattern, and high Th content. These compositional characteristics are linked with syn-mylonitic formation of plagioclase and resorption of garnet in the shear bands under amphibolite fades conditions. The absence of pre-mylonitic monazite in the shear zones, in contrast to the other domains, suggests complete dissolution of old and formation of new monazite. This probably results from an increased alkalinity and reactivity of the fluid that again is controlled by syn-mylonitic interaction with feldspar and apatite in the shear zones. There, the deformation was accommodated by dissolution precipitation creep at ca. 690 +/- 50 degrees C and 6-7.5 kbar. Growth of monazite at 55 +/- 1 Ma dates this deformation, which precedes the regional migmatization of the Sidironero Complex, whereas rutile and biotite ages reflect these later stages. This new pressure-temperature-time constraint for a relictic deformation structure provides insight into the still missing parts of the overall metamorphic, deformation and exhumation processes of the UHP units in the Rhodope. (C) 2015 Elsevier B.V. All rights reserved.
The correlation of deformation fabrics and metamorphic reactions with geochronologic data of UHP metamorphic rocks demonstrate that the multistage subduction and exhumation history of the Central Dabie Complex requires rapid subduction and rapid initial exhumation. Moreover, these data show that volume diffusion is not the major resetting mechanism of radiogenic isotope systems. Thus, our age data do not simply reflect a thermal/cooling history. In the investigated section, the maximum age for UHP is given by the 244 +/- 3 Ma (2 sigma) U-Pb age of a pre-UHP titanite phenocryst that survived UHP metamorphism and subsequent tectonometamorphic events. A minimum age for UHP is set by the 238 +/- 1 Ma (2 sigma) U-238-Pb-206 mineral isochron age of titanite and cogenetic epidote. These minerals formed from local partial melts during ascent and their age suggests fast exhumation and emplacement in the middle crust. In the period of ca. 238-218 Ma, the UHP terrain records HT metamorphism, local partial melting, and extensive pervasive strain below the eclogite (jd+grt) stability field. Exhumation was polyphase with a first phase of fast exhumation, succeeded by episodes of HT metamorphism and concomitant deformation at deep/mid crustal level between 238 and 218 Ma. Slow exhumation related to the final emplacement of tectonic units along greenschist facies shear zones did not cease before ca. 209-204 Ma. The resetting and homogenization of radiogenic isotope systems were aided by dissolution precipitation creep, which was the dominant deformation mechanism in quartz-feldspar rocks, in combination with fluid influx. (c) 2005 Elsevier B.V. All rights reserved
A ground motion logic tree for seismic hazard analysis in the stable cratonic region of Europe
(2020)
Regions of low seismicity present a particular challenge for probabilistic seismic hazard analysis when identifying suitable ground motion models (GMMs) and quantifying their epistemic uncertainty. The 2020 European Seismic Hazard Model adopts a scaled backbone approach to characterise this uncertainty for shallow seismicity in Europe, incorporating region-to-region source and attenuation variability based on European strong motion data. This approach, however, may not be suited to stable cratonic region of northeastern Europe (encompassing Finland, Sweden and the Baltic countries), where exploration of various global geophysical datasets reveals that its crustal properties are distinctly different from the rest of Europe, and are instead more closely represented by those of the Central and Eastern United States. Building upon the suite of models developed by the recent NGA East project, we construct a new scaled backbone ground motion model and calibrate its corresponding epistemic uncertainties. The resulting logic tree is shown to provide comparable hazard outcomes to the epistemic uncertainty modelling strategy adopted for the Eastern United States, despite the different approaches taken. Comparison with previous GMM selections for northeastern Europe, however, highlights key differences in short period accelerations resulting from new assumptions regarding the characteristics of the reference rock and its influence on site amplification.
Probabilistic assessment of seismic hazard and risk over a geographical region presents the modeler with challenges in the characterization of the site amplification that are not present in site-specific assessment. Using site-to-site residuals from a ground motion model fit to observations from the Japanese KiK-net database, correlations between measured local amplifications and mappable proxies such as topographic slope and geology are explored. These are used subsequently to develop empirical models describing amplification as a direct function of slope, conditional upon geological period. These correlations also demonstrate the limitations of inferring 30-m shearwave velocity from slope and applying them directly into ground motion models. Instead, they illustrate the feasibility of deriving spectral acceleration amplification factors directly from sets of observed records, which are calibrated to parameters that can be mapped uniformly on a regional scale. The result is a geologically calibrated amplification model that can be incorporated into national and regional seismic hazard and risk assessment, ensuring that the corresponding total aleatory variability reflects the predictive capability of the mapped site proxy.
The chemical composition of eclogites, found as boulders in a Tertiary conglomerate from the Guajira Peninsula, Colombia suggests that these rocks are mainly metamorphosed basaltic andesites. They are depleted in LILE elements compared to MORB, have a negative Nb-anomaly and flat to enriched REE patterns, suggesting that their protoliths evolved in a subduction related tectonic setting. They show island-arc affinities and are similar to primitive island-arc rocks described in the Caribbean. The geochemical characteristics are comparable to low-grade greenschists from the nearby Etpana Terrane, which are interpreted as part of a Cretaceous intra-oceanic arc. These data support evidence that the eclogites and the Etpana terrane rocks formed from the same volcano-sedimentary sequence. Part of this sequence was accreted onto the margin and another was incorporated into the subduction channel and metamorphosed at eclogite facies conditions. Ar-40-Ar-39 ages of 79.2 +/- 1.1Ma and 82.2 +/- 2.5Ma determined on white micas, separated from two eclogite samples, are interpreted to be related to the cooling of the main metamorphic event. The formation of a common volcano-sedimentary protolith and subsequent metamorphism of these units record the ongoing Late Cretaceous continental subduction of the South American margin within the Caribbean intra-oceanic arc subduction zone. This gave way to an arc-continent collision between the Caribbean and the South American plates, where this sequence was exhumed after the Campanian.
To address one of the central questions of plate tectonics-How do large transform systems work and what are their typical features?-seismic investigations across the Dead Sea Transform (DST), the boundary between the African and Arabian plates in the Middle East, were conducted for the first time. A major component of these investigations was a combined reflection/ refraction survey across the territories of Palestine, Israel and Jordan. The main results of this study are: (1) The seismic basement is offset by 3-5 km under the DST, (2) The DST cuts through the entire crust, broadening in the lower crust, (3) Strong lower crustal reflectors are imaged only on one side of the DST, (4) The seismic velocity sections show a steady increase in the depth of the crust-mantle transition (Moho) from 26 km at the Mediterranean to 39 km under the Jordan highlands, with only a small but visible, asymmetric topography of the Moho under the DST. These observations can be linked to the left-lateral movement of 105 km of the two plates in the last 17 Myr, accompanied by strong deformation within a narrow zone cutting through the entire crust. Comparing the DST and the San Andreas Fault (SAF) system, a strong asymmetry in subhorizontal lower crustal reflectors and a deep reaching deformation zone both occur around the DST and the SAF. The fact that such lower crustal reflectors and deep deformation zones are observed in such different transform systems suggests that these structures are possibly fundamental features of large transform plate boundaries
Fault zones are the locations where motion of tectonic plates, often associated with earthquakes, is accommodated. Despite a rapid increase in the understanding of faults in the last decades, our knowledge of their geometry, petrophysical properties, and controlling processes remains incomplete. The central questions addressed here in our study of the Dead Sea Transform (DST) in the Middle East are as follows: (1) What are the structure and kinematics of a large fault zone? (2) What controls its structure and kinematics? (3) How does the DST compare to other plate boundary fault zones? The DST has accommodated a total of 105 km of left-lateral transform motion between the African and Arabian plates since early Miocene (similar to 20 Ma). The DST segment between the Dead Sea and the Red Sea, called the Arava/Araba Fault (AF), is studied here using a multidisciplinary and multiscale approach from the mu m to the plate tectonic scale. We observe that under the DST a narrow, subvertical zone cuts through crust and lithosphere. First, from west to east the crustal thickness increases smoothly from 26 to 39 km, and a subhorizontal lower crustal reflector is detected east of the AF. Second, several faults exist in the upper crust in a 40 km wide zone centered on the AF, but none have kilometer-size zones of decreased seismic velocities or zones of high electrical conductivities in the upper crust expected for large damage zones. Third, the AF is the main branch of the DST system, even though it has accommodated only a part (up to 60 km) of the overall 105 km of sinistral plate motion. Fourth, the AF acts as a barrier to fluids to a depth of 4 km, and the lithology changes abruptly across it. Fifth, in the top few hundred meters of the AF a locally transpressional regime is observed in a 100-300 m wide zone of deformed and displaced material, bordered by subparallel faults forming a positive flower structure. Other segments of the AF have a transtensional character with small pull-aparts along them. The damage zones of the individual faults are only 5-20 m wide at this depth range. Sixth, two areas on the AF show mesoscale to microscale faulting and veining in limestone sequences with faulting depths between 2 and 5 km. Seventh, fluids in the AF are carried downward into the fault zone. Only a minor fraction of fluids is derived from ascending hydrothermal fluids. However, we found that on the kilometer scale the AF does not act as an important fluid conduit. Most of these findings are corroborated using thermomechanical modeling where shear deformation in the upper crust is localized in one or two major faults; at larger depth, shear deformation occurs in a 20-40 km wide zone with a mechanically weak decoupling zone extending subvertically through the entire lithosphere.
The motion of tectonic plates is accommodated at fault zones. One of the unanswered questions about fault zones relates to the role they play in controlling shallow and local hydrology. This study focuses on the Arava/Araba Fault (AF) zone, the southern portion of the Dead Sea Transform (DST) in the Middle East. We combine seismic and electromagnetic methods (EM) to image the geometry and map the petro-physical properties and water occurrence in the top 100 m of this active fault. For three profiles, P-velocity and resistivity images were derived independently. Using a neural network cluster analysis three classes with similar P-velocity and resistivities could then be determined from these images. These classes correspond to spatial domains of specific material and wetness. The first class occurs primarily east of the fault consisting of 'wet' sand (dunes) and brecciated sediments, whereas the second class composed of similar material located west of the fault is 'dry'. The third class lies at depth below ca. 50 m and is composed of highly deformed and weathered Precambrian rocks that constitute the multi-branch fault zone of the AF at this location. The combination of two independent measurements like seismics and EM linked by a stringent mathematical approach has thus shown the potential to delineate the interplay of lithology and water near active faults.
Studies of seismic tomography have been highly successful at imaging the deep structure of subduction zones. In a study complementary to these tomographic studies, we use array seismology and reflected waves to image a stagnant slab in the mantle transition zone. Using P and S (SH) waves we find a steeply dipping reflector centred at ca. 400 km depth and ca. 550 km west of the present Mariana subduction zone (at 20N, 140E). The discovery of this anomaly in tomography and independently in array seismology (this paper) helps in understanding the evolution of the Mariana margin. The reflector/stagnant slab may be the remains of the hypothetical North New Guinea Plate, which was theorized to have subducted ca. 50 Ma.
A key question for the development of geothermal plants is the seismic detection and monitoring of fluid injections at several kilometers depth. The detection and monitoring limits are controlled by several parameters, for example, the strength of seismic sources, number of receivers, vertical stacking, and noise conditions. For a known reference reflector at 2.66 km depth at a geothermal site in northern Germany the results of a simple surface seismic experiment were therefore combined with numerical forward modeling for different injection scenarios at 3.8 km depth. The underlying idea is that changes of reflectivity from the injection at 3.8 km must be larger than the variance of the measurements to be observable. Assuming that the injection at 3.8 km depth would produce a subhorizontal disklike target with a fracture porosity of 2% or 5% (the critical porosity) the water injection volume has to be at least 443 and 115 m(3), respectively, to be detectable from the surface. If the injection on the other hand does not create subhorizontal but subvertical pathways or only reduces the seismic velocities via the increased pore pressure in the immediate vicinity of the bore hole, the injection is undetectable from the surface. The most promising approach is therefore to move sources and/or receivers closer to the target, that is, the use of borehole instrumentation
We remelted and analyzed crystallized silicate melt inclusions in quartz from a porphyritic albitezinnwaldite microgranite dike to determine the composition of highly evolved, shallowly intruded, Li- and F-rich granitic magma and to investigate the role of crystal fractionation and aqueous fluid exsolution in causing the extreme extent of magma differentiation. This dike is intimately associated with tin- and tungsten-mineralized granites of Zinnwald, Erzgebirge, Germany. Prior research on Zinnwald granite geochemistry was limited by the effects of strong and pervasive greisenization and alkali-feldspar metasomatism of the rocks. These melt inclusions, however, provide important new constraints on magmatic and mineralizing processes in Zinnwald magmas. The mildly peraluminous granitic melt inclusions are strongly depleted in CAFEMIC constituents (e.g., CaO, FeO, MgO, TiO2), highly enriched in lithophile trace elements, and highly but variably enriched in F and Cl. The melt inclusions contain up to several thousand ppm Cl and nearly 3 wt% F, on average; several inclusions contain more than 5 wt% F. The melt inclusions are geochemically similar to the corresponding whole-rock sample, except that the former contain much more F and less CaO, FeO, Zr, Nb, Sr, and Ba. The Sr and Ba abundances are very low implying the melt inclusions represent magma that was more evolved than that represented by the bulk rock. Relationships involving melt constituents reflect increasing lithophile-element and halogen abundances in residual melt with progressive magma differentiation. Modeling demonstrates that differentiation was dominated by crystal fractionation involving quartz and feldspar and significant quantities of topaz and F-rich zinnwaldite. The computed abundances of the latter phases greatly exceed their abundances in the rocks, suggesting that the residual melt was separated physically from phenocrysts during magma movement and evolution. Interactions of aqueous fluids with silicate melt were also critical to magma evolution. To better understand the role of halogen-charged, aqueous fluids in magmatic differentiation and in subsequent mineralization and metasomatism of the Zinnwald granites, Cl-partitioning experiments were conducted with a F-enriched silicate melt and aqueous fluids at 2,000 bar (200 MPa). The results of the experimentally determined partition coefficients for Cl and F, the compositions of fluid inclusions in quartz and other phenocrysts, and associated geochemical modeling point to an important role of magmatic-hydrothermal fluids in influencing magma geochemistry and evolution. The exsolution of halogen-charged fluids from the Li- and F- enriched Zinnwald granitic magma modified the Cl, alkali, and F contents of the residual melt, and may have also sequestered Li, Sri, and W from the melt. Many of these fluids contained strongly elevated F concentrations that were equivalent to or greater than their Cl abundances. The exsolution of F-, Cl-, Li-, +/- W- and Sn-bearing hydrothermal fluids from Zinnwald granite magmas was important in effecting the greisenizing and alkali-feldspathizing metasomatism of the granites and the concomitant mineralization
The Yukon Coast in Canada is an ice-rich permafrost coast and highly sensitive to changing environmental conditions. Retrogressive thaw slumps are a common thermoerosion feature along this coast, and develop through the thawing of exposed ice-rich permafrost on slopes and removal of accumulating debris. They contribute large amounts of sediment, including organic carbon and nitrogen, to the nearshore zone.
The objective of this study was to 1) identify the climatic and geomorphological drivers of sediment-meltwater release, 2) quantify the amount of released meltwater, sediment, organic carbon and nitrogen, and 3) project the evolution of sediment-meltwater release of retrogressive thaw slumps in a changing future climate.
The analysis is based on data collected over 18 days in July 2013 and 18 days in August 2012. A cut-throat flume was set up in the main sediment-meltwater channel of the largest retrogressive thaw slump on Herschel Island. In addition, two weather stations, one on top of the undisturbed tundra and one on the slump floor, measured incoming solar radiation, air temperature, wind speed and precipitation. The discharge volume eroding from the ice-rich permafrost and retreating snowbanks was measured and compared to the meteorological data collected in real time with a resolution of one minute.
The results show that the release of sediment-meltwater from thawing of the ice-rich permafrost headwall is strongly related to snowmelt, incoming solar radiation and air temperature. Snowmelt led to seasonal differences, especially due to the additional contribution of water to the eroding sediment-meltwater from headwall ablation, lead to dilution of the sediment-meltwater composition. Incoming solar radiation and air temperature were the main drivers for diurnal and inter-diurnal fluctuations. In July (2013), the retrogressive thaw slump released about 25 000 m³ of sediment-meltwater, containing 225 kg dissolved organic carbon and 2050 t of sediment, which in turn included 33 t organic carbon, and 4 t total nitrogen. In August (2012), just 15 600 m³ of sediment-meltwater was released, since there was no additional contribution from snowmelt. However, even without the additional dilution, 281 kg dissolved organic carbon was released. The sediment concentration was twice as high as in July, with sediment contents of up to 457 g l-1 and 3058 t of sediment, including 53 t organic carbon and 5 t nitrogen, being released.
In addition, the data from the 36 days of observations from Slump D were upscaled to cover the main summer season of 1 July to 31 August (62 days) and to include all 229 active retrogressive thaw slumps along the Yukon Coast. In total, all retrogressive thaw slumps along the Yukon Coast contribute a minimum of 1.4 Mio. m³ sediment-meltwater each thawing season, containing a minimum of 172 000 t sediment with 3119 t organic carbon, 327 t nitrogen and 17 t dissolved organic carbon. Therefore, in addition to the coastal erosion input to the Beaufort Sea, retrogressive thaw slumps additionally release 3 % of sediment and 8 % of organic carbon into the ocean. Finally, the future evolution of retrogressive thaw slumps under a warming scenario with summer air temperatures increasing by 2-3 °C by 2081-2100, would lead to an increase of 109-114% in release of sediment-meltwater.
It can be concluded that retrogressive thaw slumps are sensitive to climatic conditions and under projected future Arctic warming will contribute larger amounts of thawed permafrost material (including organic carbon and nitrogen) into the environment.
Germany and the United Kingdom have domestic shale gas reserves which they may exploit in the future to complement their national energy strategies. However gas production releases volatile organic compounds (VOC) and nitrogen oxides (NOx), which through photochemical reaction form ground-level ozone, an air pollutant that can trigger adverse health effects e.g. on the respiratory system. This study explores the range of impacts of a potential shale gas industry in these two countries on local and regional ambient ozone. To this end, comprehensive emission scenarios are used as the basis for input to an online-coupled regional chemistry transport model (WRF-Chem). Here we simulate shale gas scenarios over summer (June, July, August) 2011, exploring the effects of varying VOC emissions, gas speciation, and concentration of NOx emissions over space and time, on ozone formation. An evaluation of the model setup is performed, which exhibited the model’s ability to predict surface meteorological and chemical variables well compared with observations, and consistent with other studies. When different shale gas scenarios were employed, the results show a peak increase in maximum daily 8-hour average ozone from 3.7 to 28.3 μg m–3. In addition, we find that shale gas emissions can force ozone exceedances at a considerable percentage of regulatory measurement stations locally (up to 21% in Germany and 35% in the United Kingdom) and in distant countries through long-range transport, and increase the cumulative health-related metric SOMO35 (maximum percent increase of ~28%) throughout the region. Findings indicate that VOC emissions are important for ozone enhancement, and to a lesser extent NOx, meaning that VOC regulation for a future European shale gas industry will be of especial importance to mitigate unfavorable health outcomes. Overall our findings demonstrate that shale gas production in Europe can worsen ozone air quality on both the local and regional scales.
The advantages of remote sensing using Unmanned Aerial Vehicles (UAVs) are a high spatial resolution of images, temporal flexibility and narrow-band spectral data from different wavelengths domains. This enables the detection of spatio-temporal dynamics of environmental variables, like plant-related carbon dynamics in agricultural landscapes. In this paper, we quantify spatial patterns of fresh phytomass and related carbon (C) export using imagery captured by a 12-band multispectral camera mounted on the fixed wing UAV Carolo P360. The study was performed in 2014 at the experimental area CarboZALF-D in NE Germany. From radiometrically corrected and calibrated images of lucerne (Medicago sativa), the performance of four commonly used vegetation indices (VIs) was tested using band combinations of six near-infrared bands. The highest correlation between ground-based measurements of fresh phytomass of lucerne and VIs was obtained for the Enhanced Vegetation Index (EVI) using near-infrared band b(899). The resulting map was transformed into dry phytomass and finally upscaled to total C export by harvest. The observed spatial variability at field- and plot-scale could be attributed to small-scale soil heterogeneity in part.
The advantages of remote sensing using Unmanned Aerial Vehicles (UAVs) are a high spatial resolution of images, temporal flexibility and narrow-band spectral data from different wavelengths domains. This enables the detection of spatio-temporal dynamics of environmental variables, like plant-related carbon dynamics in agricultural landscapes. In this paper, we quantify spatial patterns of fresh phytomass and related carbon (C) export using imagery captured by a 12-band multispectral camera mounted on the fixed wing UAV Carolo P360. The study was performed in 2014 at the experimental area CarboZALF-D in NE Germany. From radiometrically corrected and calibrated images of lucerne (Medicago sativa), the performance of four commonly used vegetation indices (VIs) was tested using band combinations of six near-infrared bands. The highest correlation between ground-based measurements of fresh phytomass of lucerne and VIs was obtained for the Enhanced Vegetation Index (EVI) using near-infrared band b(899). The resulting map was transformed into dry phytomass and finally upscaled to total C export by harvest. The observed spatial variability at field- and plot-scale could be attributed to small-scale soil heterogeneity in part.
The growing body of research on large-scale mass wasting events so far has only scarcely investigated the sedimentology of chaotic deposits from non-volcanic terrestrial landslides such that any overarching and systematic terminological framework remains elusive. Yet recent work has emphasized the need for better understanding the internal structure and composition of rockslide deposits as a means to characterise the mechanics during the final stages of runout and emplacement. We offer a comprehensive overview on the occurrence of rock fragmentation and frictional melt both at different geographic locations, and different sections within large (>10(6) m(3)) rockslide masses. We argue that exposures of pervasively fragmented and interlocked jigsaw-cracked rock masses; basal melange containing rip-up clasts and phantom blocks; micro-breccia; and thin bands of basal frictionite are indispensable clues for identifying deposits from giant rockslides that may remain morphologically inconspicuous otherwise. These sedimentary assemblages are diagnostic tools for distinguishing large rockslide debris from macro and microscopically similar glacial deposits, tectonic fault-zone breccias, and impact breccias, and thus help avoid palaeoclimatic and tectonic misinterpretations, let alone misestimates of the hazard from giant rockslides. Moreover, experimental results from Mossbauer spectroscopy of frictionite samples support visual interpretations of thin sections, and demonstrate that short-lived (<10 s) friction-induced partial melting at temperatures >1500 degrees C in the absence of water occurred at the base of several giant moving rockslides. This finding supports previous theories of dry excess runout accompanied by comminution of rock masses down to gm-scale, and indicates that catastrophic motion of large fragmenting rock masses does not require water as a potential lubricant.
A 3-D crustal shear wave velocity model and Moho map below the Semail Ophiolite, eastern Arabia
(2022)
The Semail Ophiolite in eastern Arabia is the largest and best-exposed slice of oceanic lithosphere on land. Detailed knowledge of the tectonic evolution of the shallow crust, in particular during and after ophiolite obduction in Late Cretaceous times is contrasted by few constraints on physical and compositional properties of the middle and lower continental crust below the obducted units. The role of inherited, pre-obduction crustal architecture remains therefore unaccounted for in our understanding of crustal evolution and the present-day geology. Based on seismological data acquired during a 27-month campaign in northern Oman, Ambient Seismic Noise Tomography and Receiver Function analysis provide for the first time a 3-D radially anisotropic shear wave velocity (V-S) model and a consistent Moho map below the iconic Semail Ophiolite. The model highlights deep crustal boundaries that segment the eastern Arabian basement in two distinct units. The previously undescribed Western Jabal Akhdar Zone separates Arabian crust with typical continental properties and a thickness of similar to 40-45 km in the northwest from a compositionally different terrane in the southeast that is interpreted as a terrane accreted during the Pan-African orogeny in Neoproterozoic times. East of the Ibra Zone, another deep crustal boundary, crustal thickness decreases to 30-35 km and very high lower crustal V-S suggest large-scale mafic intrusions into, and possible underplating of the Arabian continental crust that occurred most likely during Permian breakup of Pangea. Mafic reworking is sharply bounded by the (upper crustal) Semail Gap Fault Zone, northwest of which no such high velocities are found in the crust. Topography of the Oman Mountains is supported by a mild crustal root and Moho depth below the highest topography, the Jabal Akhdar Dome, is similar to 42 km. Radial anisotropy is robustly resolved in the upper crust and aids in discriminating dipping allochthonous units from autochthonous sedimentary rocks that are indistinguishable by isotropic V-S alone. Lateral thickness variations of the ophiolite highlight the Haylayn Ophiolite Massif on the northern flank of Jabal Akhdar Dome and the Hawasina Window as the deepest reaching unit. Ophiolite thickness is similar to 10 km in the southern and northern massifs, and <= 5 km elsewhere.
X-ray emission and X-ray Raman scattering spectroscopy are powerful tools to investigate the local electronic and atomic structure of high and low Z elements in situ. Notably, these methods can be applied for in situ spectroscopy at high pressure and high temperature using resistively or laser-heated diamond anvil cells in order to achieve thermodynamic conditions which appear in the Earth's interior. We present a setup for combined X-ray emission and X-ray Raman scattering studies at beamline P01 of PETRA III using a portable wavelength-dispersive von Hamos spectrometer together with the permanently installed multiple-analyzer Johann-type spectrometer. The capabilities of this setup are exemplified by investigating the iron spin crossover of siderite FeCO3 up to 49.3 GPa by measuring the Fe M2,3-edge and the Fe Kβ1,3 emission line simultaneously. With this setup, the Fe valence-to-core emission can be detected together with the Kβ1,3 emission line providing complementary information on the sample's electronic structure. By implementing a laser-heating device, we demonstrate the strength of using a von Hamos type spectrometer for spin state mapping at extreme conditions. Finally, we give different examples of low Z elements' absorption edges relevant for application in geoscience that are accessible with the Johann-type XRS spectrometer. With this setup new insights into the spin transition and compression mechanisms of Earth's mantle materials can be obtained of importance for comprehension of the macroscopic physical and chemical properties of the Earth's interior.
Deformation associated with plate convergence at subduction zones is accommodated by a complex system involving fault slip and viscoelastic flow. These processes have proven difficult to disentangle. The 2010 M-w 8.8 Maule earthquake occurred close to the Chilean coast within a dense network of continuously recording Global Positioning System stations, which provide a comprehensive history of surface strain. We use these data to assemble a detailed picture of a structurally controlled megathrust fault frictional patchwork and the three-dimensional rheological and time-dependent viscosity structure of the lower crust and upper mantle, all of which control the relative importance of afterslip and viscoelastic relaxation during postseismic deformation. These results enhance our understanding of subduction dynamics including the interplay of localized and distributed deformation during the subduction zone earthquake cycle.
High-resolution surface velocities and strain for Anatolia from Sentinel-1 InSAR and GNSS data
(2020)
Measurements of present-day surface deformation are essential for the assessment of long-term seismic hazard. The European Space Agency's Sentinel-1 satellites enable global, high-resolution observation of crustal motion from Interferometric Synthetic Aperture Radar (InSAR). We have developed automated InSAR processing systems that exploit the first similar to 5 years of Sentinel-1 data to measure surface motions for the similar to 800,000-km(2) Anatolian region. Our new 3-D velocity and strain rate fields illuminate deformation patterns dominated by westward motion of Anatolia relative to Eurasia, localized strain accumulation along the North and East Anatolian Faults, and rapid vertical signals associated with anthropogenic activities and to a lesser extent extension across the grabens of western Anatolia. We show that automatically processed Sentinel-1 InSAR data can characterize details of the velocity and strain rate fields with high resolution and accuracy over large regions. These results are important for assessing the relationship between strain accumulation and release in earthquakes. <br /> Plain Language Summary Satellite-based measurements of small rates of motion of the Earth's surface made at high spatial resolutions and over large areas are important for many geophysical applications including improving earthquake hazard models. We take advantage of recent advances in geodetic techniques in order to measure surface velocities and tectonic strain accumulation across the Anatolia region, including the highly seismogenic and often deadly North Anatolian Fault. We show that by combining Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data with Global Navigation Satellite System (GNSS) measurements we can enhance our view of surface deformation associated with active tectonics, the earthquake cycle, and anthropogenic processes.
To safeguard the sustainable use of ecosystems and their services, early detection of potentially damaging changes in functional capabilities is needed. To support a proper ecosystem management, the analysis of an ecosystem’s vulnerability provide information on its weaknesses as well as on its capacity to recover after suffering an impact. However, the application of the vulnerability concept to ecosystems is still an emerging topic. After providing background on the vulnerability concept, we summarize existing ecosystem vulnerability research on the basis of a systematic literature review with a special focus on ecosystem type, disciplinary background, and more detailed definition of the ecosystem vulnerability components. Using the Web of ScienceTM Core Collection, we overviewed the literature from 1991 onwards but used the 5 years from 2011 to 2015 for an in-depth analysis, including 129 articles. We found that ecosystem vulnerability analysis has been applied most notably in conservation biology, climate change research, and ecological risk assessments, pinpointing a limited spreading across the environmental sciences. It occurred primarily within marine and freshwater ecosystems. To avoid confusion, we recommend using the unambiguous term ecosystem vulnerability rather than ecological, environmental, population, or community vulnerability. Further, common ground has been identified, on which to define the ecosystem vulnerability components exposure, sensitivity, and adaptive capacity. We propose a framework for ecosystem assessments that coherently connects the concepts of vulnerability, resilience, and adaptability as different ecosystem responses. A short outlook on the possible operationalization of the concept by ecosystem vulnerabilty indices, and a conclusion section complete the review.
This review analyzes the potential role and long-term effects of field perennial polycultures (mixtures) in agricultural systems, with the aim of reducing the trade-offs between provisioning and regulating ecosystem services. First, crop rotations are identified as a suitable tool for the assessment of the long-term effects of perennial polycultures on ecosystem services, which are not visible at the single-crop level. Second, the ability of perennial polycultures to support ecosystem services when used in crop rotations is quantified through eight agricultural ecosystem services. Legume-grass mixtures and wildflower mixtures are used as examples of perennial polycultures, and compared with silage maize as a typical crop for biomass production. Perennial polycultures enhance soil fertility, soil protection, climate regulation, pollination, pest and weed control, and landscape aesthetics compared with maize. They also score lower for biomass production compared with maize, which confirms the trade-off between provisioning and regulating ecosystem services. However, the additional positive factors provided by perennial polycultures, such as reduced costs for mineral fertilizer, pesticides, and soil tillage, and a significant preceding crop effect that increases the yields of subsequent crops, should be taken into account. However, a full assessment of agricultural ecosystem services requires a more holistic analysis that is beyond the capabilities of current frameworks.
Isolated extreme habitats are ideally suited to investigate pivotal ecological processes such as niche use, local adaptation and dispersal. Extremophilic animals living in isolated habitats face the problem that dispersal is limited through the absence of suitable dispersal corridors, which in turn facilitates local adaptation. We used five rotifer isolates from extremely acidic mining lakes with a pH of below 3 as model organisms to test whether these isolates are acidotolerant or acidophilic, whether they survive and reproduce at their niche edges (here pH 2 and circum-neutral pH) and whether local adaptation has evolved. To evaluate potential dispersal limitation, we tested whether animals and their parthenogenetic eggs survive and remain reproductive or viable at unfavourable pH-conditions. All five isolates were acidophilic with a pH-optimum in the range of 4-6, which is well above the pH (< 3) of their lakes of origin. At unfavourable high pH, in four out of the five isolates parthenogenetic females produced a high number of non-viable eggs. Females and eggs produced at favourable pH (4) remained vital at an otherwise unfavourable pH of 7, indicating that for dispersal no acidic dispersal corridors are necessary. Common garden experiments revealed no clear evidence for local adaptation in any of the five isolates. Despite their acidophilic nature, all five isolates can potentially disperse via circum-neutral water bodies as long as their residence time is short, suggesting a broader dispersal niche than their realized niche. Local adaptation might have been hampered by the low population sizes of the rotifers in their isolated habitat and the short time span the mining lakes have existed.
Putting the plan into practice implementation of proposals for measures of local landscape plans
(2012)
The knowledge of the effectiveness of local landscape planning in Germany is in the main limited to particular cases and derives mostly from qualitative single case studies. This applies especially to the implementation of measures defined by landscape plans. To fill that gap, the paper focuses on the implementation of those measures. Furthermore, it discusses the factors and framework conditions which are crucial for this implementation. The potential factors and conditions of influence were derived from theory and compiled in 20 investigation hypotheses. In order to gain information on the execution of the measures, 28 randomly selected plans were first analysed, then interviews were carried out with administration representatives. It can be stated that landscape planning has positively influenced the development of nature and landscape in the investigated municipalities. A considerable number of measures had been implemented, although landscape planning as a supply-side instrument proposes generally a very large number of measures. Factors with a positive effect on the implementation of landscape planning measures are pointed out.
An innovative approach to improve SRTM DEM using multispectral imagery and artificial neural network
(2016)
Although the Shuttle Radar Topography Mission [SRTM) data are a publicly accessible Digital Elevation Model [DEM) provided at no cost, its accuracy especially at forested area is known to be limited with root mean square error (RMSE) of approx. 14 m in Singapore's forested area. Such inaccuracy is attributed to the 5.6 cm wavelength used by SRTM that does not penetrate vegetation well. This paper considers forested areas of central catchment of Singapore as a proof of concept of an approach to improve the SRTM data set. The approach makes full use of (1) the introduction of multispectral imagery (Landsat 8), of 30 m resolution, into SRTM data; (2) the Artificial Neural Network (ANN) to flex its known strengths in pattern recognition and; (3) a reference DEM of high accuracy (1 m) derived through the integration of stereo imaging of worldview-1 and extensive ground survey points. The study shows a series of significant improvements of the SRTM when assessed with the reference DEM of 2 different areas, with RMSE reduction of ∼68% (from 13.9 m to 4.4 m) and ∼52% (from 14.2 m to 6.7 m). In addition, the assessment of the resulting DEM also includes comparisons with simple denoising methodology (Low Pass Filter) and commercially available product called NEXTMap® World 30™.
One main challenge in constructing a reliable recurrence plot (RP) and, hence, its quantification [recurrence quantification analysis (RQA)] of a continuous dynamical system is the induced noise that is commonly found in observation time series. This induced noise is known to cause disrupted and deviated diagonal lines despite the known deterministic features and, hence, biases the diagonal line based RQA measures and can lead to misleading conclusions. Although discontinuous lines can be further connected by increasing the recurrence threshold, such an approach triggers thick lines in the plot. However, thick lines also influence the RQA measures by artificially increasing the number of diagonals and the length of vertical lines [e.g., Determinism (DET) and Laminarity (LAM) become artificially higher]. To take on this challenge, an extended RQA approach for accounting disrupted and deviated diagonal lines is proposed. The approach uses the concept of a sliding diagonal window with minimal window size that tolerates the mentioned deviated lines and also considers a specified minimal lag between points as connected. This is meant to derive a similar determinism indicator for noisy signal where conventional RQA fails to capture. Additionally, an extended local minima approach to construct RP is also proposed to further reduce artificial block structures and vertical lines that potentially increase the associated RQA like LAM. The methodology and applicability of the extended local minima approach and DET equivalent measure are presented and discussed, respectively.
As an effort to reduce parameter uncertainties in constructing recurrence plots, and in particular to avoid potential artefacts, this paper presents a technique to derive artefact-safe region of parameter sets. This technique exploits both deterministic (incl. chaos) and stochastic signal characteristics of recurrence quantification (i.e. diagonal structures). It is useful when the evaluated signal is known to be deterministic. This study focuses on the recurrence plot generated from the reconstructed phase space in order to represent many real application scenarios when not all variables to describe a system are available (data scarcity). The technique involves random shuffling of the original signal to destroy its original deterministic characteristics. Its purpose is to evaluate whether the determinism values of the original and the shuffled signal remain closely together, and therefore suggesting that the recurrence plot might comprise artefacts. The use of such determinism-sensitive region shall be accompanied by standard embedding optimization approaches, e.g. using indices like false nearest neighbor and mutual information, to result in a more reliable recurrence plot parameterization.
Aquifer thermal energy storage (ATES) as a complement to fluctuating renewable energy systems is a reliable technology to guarantee continuous energy supply for heating and air conditioning. We investigated a high-temperature (HT) mono-well system (c. 100 degrees C), where the well screens are separated vertically within the aquifer, as an alternative to conventional doublet ATES systems for an underground storage in northern Oman. We analysed the impact of thermal inference between injection and extraction well screens on the heat recovery factor (HRF) in order to define the optimal screento-screen distance for best possible systems efficiency. Two controlling interference parameters were considered: the vertical screen-to-screen distance and aquifer heterogeneities. The sensitivity study shows that with decreasing screen-to-screen distances, thermal interference increases storage performance. A turning point is reached if the screen distance is too close, causing either water breakthrough or negative thermal interference between the screens. Our simulations show that a combined heat plume with spherical geometry results in the highest heat recovery factors due to the lowest surface area to volume ratios. Thick aquifers for mono-well HT-ATES are thus not mandatory Our study shows that a HT-ATES mono-well system is a feasible storage design with high heat recovery factors for continuous cooling or heating purposes.
The cultivation of plants for use as energy resources is an agricultural and industrial sector with potentially synergistic benefits related to protecting the environment and generating income. Against the background of increasing land-use changes and new agricultural approaches to the production of energy crops, we present a method for identifying future-oriented crop rotations that supports both the economic and environmental components of decision-making strategies with respect to agriculture-related policy decisions (regional mission statements). The conflicting aspects of these objectives can be addressed with the analytic hierarchy process (AHP), a multi-attribute decision-making method that was integrated here. Three models are used to generate simulations of the defined objectives over a planning period of 30 years under the current climate scenario and provide input data for the multi-attribute assessment of several crop rotations. Based on the entire evaluation process, dimensionless global priority vectors are used to indicate how well the crop rotations meet the requirements of the defined mission statement. The method is tested in a municipality in NE Germany. (C) 2014 Elsevier Ltd. All rights reserved.
Individual great earthquakes are posited to release the elastic strain energy that has accumulated over centuries by the gradual movement of tectonic plates(1,2). However, knowledge of plate deformation during a complete seismic cycle-two successive great earthquakes and the intervening interseismic period-remains incomplete(3). A complete seismic cycle began in south-central Chile in 1835 with an earthquake of about magnitude 8.5 (refs 4,5) and ended in 2010 with a magnitude 8.8 earthquake(6). During the first earthquake, an uplift of Isla Santa Maria by 2.4 to 3m was documented(4,5). In the second earthquake, the island was uplifted(7) by 1.8 m. Here we use nautical surveys made in 1804, after the earthquake in 1835 and in 1886, together with modern echo sounder surveys and GPS measurements made immediately before and after the 2010 earthquake, to quantify vertical deformation through the complete seismic cycle. We find that in the period between the two earthquakes, Isla Santa Maria subsided by about 1.4 m. We simulate the patterns of vertical deformation with a finite-element model and find that they agree broadly with predictions from elastic rebound theory(2). However, comparison with geomorphic and geologic records of millennial coastline emergence(8,9) reveal that 10-20% of the vertical uplift could be permanent.
Groundwater is critical in supporting current and future reliable water supply throughout Africa. Although continental maps of groundwater storage and recharge have been developed, we currently lack a clear understanding on how the controls on groundwater recharge vary across the entire continent. Reviewing the existing literature, we synthesize information on reported groundwater recharge controls in Africa. We find that 15 out of 22 of these controls can be characterised using global datasets. We develop 11 descriptors of climatic, topographic, vegetation, soil and geologic properties using global datasets, to characterise groundwater recharge controls in Africa. These descriptors cluster Africa into 15 Recharge Landscape Units for which we expect recharge controls to be similar. Over 80% of the continents land area is organized by just nine of these units. We also find that aggregating the Units by similarity into four broader Recharge Landscapes (Desert, Dryland, Wet tropical and Wet tropical forest) provides a suitable level of landscape organisation to explain differences in ground-based long-term mean annual recharge and recharge ratio (annual recharge / annual precipitation) estimates. Furthermore, wetter Recharge Landscapes are more efficient in converting rainfall to recharge than drier Recharge Landscapes as well as having higher annual recharge rates. In Dryland Recharge Landscapes, we found that annual recharge rates largely varied according to mean annual precipitation, whereas recharge ratio estimates increase with increasing monthly variability in P-PET. However, we were unable to explain why ground based estimates of recharge signatures vary across other Recharge Landscapes, in which there are fewer ground based recharge estimates, using global datasets alone. Even in dryland regions, there is still considerable unexplained variability in the estimates of annual recharge and recharge ratio, stressing the limitations of global datasets for investigating ground-based information.
Convergence between the Indian and Asian plates has reshaped large parts of Asia, changing regional climate and biodiversity, yet geodynamic models fundamentally diverge on how convergence was accommodated since the India-Asia collision. Here we report palaeomagnetic data from the Burma Terrane, which is at the eastern edge of the collision zone and is famous for its Cretaceous amber biota, to better determine the evolution of the India-Asia collision. The Burma Terrane was part of a Trans-Tethyan island arc and stood at a near-equatorial southern latitude at similar to 95 Ma, suggesting island endemism for the Burmese amber biota. The Burma Terrane underwent significant clockwise rotation between similar to 80 and 50 Ma, causing its subduction margin to become hyper-oblique. Subsequently, it was translated northward on the Indian Plate by an exceptional distance of at least 2,000 km along a dextral strike-slip fault system in the east. Our reconstructions are only compatible with geodynamic models involving an initial collision of India with a near-equatorial Trans-Tethyan subduction system at similar to 60 Ma, followed by a later collision with the Asian margin.
Late Quaternary landscapes of unglaciated Beringia were largely shaped by ice-wedge polygon tundra. Ice Complex (IC) strata preserve such ancient polygon formations. Here we report on the Yukagir IC from Bol'shoy Lyakhovsky Island in northeastern Siberia and suggest that new radioisotope disequilibria (230Th/U) dates of the Yukagir IC peat confirm its formation during the Marine Oxygen Isotope Stage (MIS) 7a–c interglacial period. The preservation of the ice-rich Yukagir IC proves its resilience to last interglacial and late glacial–Holocene warming. This study compares the Yukagir IC to IC strata of MIS 5, MIS 3, and MIS 2 ages exposed on Bol'shoy Lyakhovsky Island. Besides high intrasedimental ice content and syngenetic ice wedges intersecting silts, sandy silts, the Yukagir IC is characterized by high organic matter (OM) accumulation and low OM decomposition of a distinctive Drepanocladus moss-peat. The Yukagir IC pollen data reveal grass-shrub-moss tundra indicating rather wet summer conditions similar to modern ones. The stable isotope composition of Yukagir IC wedge ice is similar to those of the MIS 5 and MIS 3 ICs pointing to similar atmospheric moisture generation and transport patterns in winter. IC data from glacial and interglacial periods provide insights into permafrost and climate dynamics since about 200 ka.
Geochemical processes such as mineral dissolution and precipitation alter the microstructure of rocks, and thereby affect their hydraulic and mechanical behaviour. Quantifying these property changes and considering them in reservoir simulations is essential for a sustainable utilisation of the geological subsurface. Due to the lack of alternatives, analytical methods and empirical relations are currently applied to estimate evolving hydraulic and mechanical rock properties associated with chemical reactions. However, the predictive capabilities of analytical approaches remain limited, since they assume idealised microstructures, and thus are not able to reflect property evolution for dynamic processes. Hence, aim of the present thesis is to improve the prediction of permeability and stiffness changes resulting from pore space alterations of reservoir sandstones.
A detailed representation of rock microstructure, including the morphology and connectivity of pores, is essential to accurately determine physical rock properties. For that purpose, three-dimensional pore-scale models of typical reservoir sandstones, obtained from highly resolved micro-computed tomography (micro-CT), are used to numerically calculate permeability and stiffness. In order to adequately depict characteristic distributions of secondary minerals, the virtual samples are systematically altered and resulting trends among the geometric, hydraulic, and mechanical rock properties are quantified. It is demonstrated that the geochemical reaction regime controls the location of mineral precipitation within the pore space, and thereby crucially affects the permeability evolution. This emphasises the requirement of determining distinctive porosity-permeability relationships
by means of digital pore-scale models. By contrast, a substantial impact of spatial alterations patterns on the stiffness evolution of reservoir sandstones are only observed in case of certain microstructures, such as highly porous granular rocks or sandstones comprising framework-supporting cementations. In order to construct synthetic granular samples a process-based approach is proposed including grain deposition and diagenetic cementation. It is demonstrated that the generated samples reliably represent the microstructural complexity of natural sandstones. Thereby, general limitations of imaging techniques can be overcome and various realisations of granular rocks can be flexibly produced. These can be further altered by virtual experiments, offering a fast and cost-effective way to examine the impact of precipitation, dissolution or fracturing on various petrophysical correlations.
The presented research work provides methodological principles to quantify trends in permeability and stiffness resulting from geochemical processes. The calculated physical property relations are directly linked to pore-scale alterations, and thus have a higher accuracy than commonly applied analytical approaches. This will considerably improve the predictive capabilities of reservoir models, and is further relevant to assess and reduce potential risks, such as productivity or injectivity losses as well as reservoir compaction or fault reactivation. Hence, the proposed method is of paramount importance for a wide range of natural and engineered subsurface applications, including geothermal energy systems, hydrocarbon reservoirs, CO2 and energy storage as well as hydrothermal deposit exploration.
Predicting macroscopic elastic rock properties requires detailed information on microstructure
(2017)
Predicting variations in macroscopic mechanical rock behaviour due to microstructural changes, driven by mineral precipitation and dissolution is necessary to couple chemo-mechanical processes in geological subsurface simulations. We apply 3D numerical homogenization models to estimate Young’s moduli for five synthetic microstructures, and successfully validate our results for comparable geometries with the analytical Mori-Tanaka approach. Further, we demonstrate that considering specific rock microstructures is of paramount importance, since calculated elastic properties may deviate by up to 230 % for the same mineral composition. Moreover, agreement between simulated and experimentally determined Young’s moduli is significantly improved, when detailed spatial information are employed.
Geochemical processes change the microstructure of rocks and thereby affect their physical behaviour at the macro scale. A micro-computer tomography (micro-CT) scan of a typical reservoir sandstone is used to numerically examine the impact of three spatial alteration patterns on pore morphology, permeability and elastic moduli by correlating precipitation with the local flow velocity magnitude. The results demonstrate that the location of mineral growth strongly affects the permeability decrease with variations by up to four orders in magnitude. Precipitation in regions of high flow velocities is characterised by a predominant clogging of pore throats and a drastic permeability reduction, which can be roughly described by the power law relation with an exponent of 20. A continuous alteration of the pore structure by uniform mineral growth reduces the permeability comparable to the power law with an exponent of four or the Kozeny-Carman relation. Preferential precipitation in regions of low flow velocities predominantly affects smaller throats and pores with a minor impact on the flow regime, where the permeability decrease is considerably below that calculated by the power law with an exponent of two. Despite their complete distinctive impact on hydraulics, the spatial precipitation patterns only slightly affect the increase in elastic rock properties with differences by up to 6.3% between the investigated scenarios. Hence, an adequate characterisation of the spatial precipitation pattern is crucial to quantify changes in hydraulic rock properties, whereas the present study shows that its impact on elastic rock parameters is limited. The calculated relations between porosity and permeability, as well as elastic moduli can be applied for upscaling micro-scale findings to reservoir-scale models to improve their predictive capabilities, what is of paramount importance for a sustainable utilisation of the geological subsurface.
Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations
(2018)
The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution.
Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations
(2018)
The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution.
Diagenetic trends of synthetic reservoir sandstone properties assessed by digital rock physics
(2021)
Quantifying interactions and dependencies among geometric, hydraulic and mechanical properties of reservoir sandstones is of particular importance for the exploration and utilisation of the geological subsurface and can be assessed by synthetic sandstones comprising the microstructural complexity of natural rocks. In the present study, three highly resolved samples of the Fontainebleau, Berea and Bentheim sandstones are generated by means of a process-based approach, which combines the gravity-driven deposition of irregularly shaped grains and their diagenetic cementation by three different schemes. The resulting evolution in porosity, permeability and rock stiffness is examined and compared to the respective micro-computer tomographic (micro-CT) scans. The grain contact-preferential scheme implies a progressive clogging of small throats and consequently produces considerably less connected and stiffer samples than the two other schemes. By contrast, uniform quartz overgrowth continuously alters the pore space and leads to the lowest elastic properties. The proposed stress-dependent cementation scheme combines both approaches of contact-cement and quartz overgrowth, resulting in granulometric, hydraulic and elastic properties equivalent to those of the respective micro-CT scans, where bulk moduli slightly deviate by 0.8%, 4.9% and 2.5% for the Fontainebleau, Berea and Bentheim sandstone, respectively. The synthetic samples can be further altered to examine the impact of mineral dissolution or precipitation as well as fracturing on various petrophysical correlations, which is of particular relevance for numerous aspects of a sustainable subsurface utilisation.
Whilst sophisticated multiphase fluid flow models are routinely employed to understand behaviour of oil and gas reservoirs, high-resolution data describing the three-dimensional (3D) distribution of rock characteristics is rarely available to populate models. We present a new approach to developing a quantitative understanding of the effect of individual controls on the distribution of petrophysical properties and their impact on fluid flow. This involves simulating flow through high-detail permeability architectures generated by forward modelling of the coupled depositional-diagenetic evolution of isolated platforms using CARB3D(+). This workflow is exemplified by an investigation of interactions between subsidence and climate, and their expression in spatial variations in reservoir quality in an isolated carbonate platform of similar size and subsidence history to the Triassic Latemar Platform.
Dissolutional lowering during subaerial exposure controls platform-top graininess via platform top hydrodynamics during the subsequent transgression. Dissolved carbonate is reprecipitated as cements by percolating meteoric waters. However, associated subsurface meteoric dissolution generates significant secondary porosity under a more humid climate. Slower subsidence enhances diagenetic overprinting during repeated exposure events. Single-phase streamline simulations show how early diagenesis develops more permeable fairways within the finer-grained condensed units that can act as thief zones for flow from the grainier but less diagenetically altered cyclic units.
The East African Plateau provides a spectacular example of geodynamic plateau uplift, active continental rifting, and associated climatic forcing. It is an integral part of the East African Rift System and has an average elevation of approximately 1,000 m. Its location coincides with a negative Bouguer gravity anomaly with a semi-circular shape, closely related to a mantle plume, which influences the Cenozoic crustal development since its impingement in Eocene-Oligocene time. The uplift of the East African Plateau, preceding volcanism, and rifting formed an important orographic barrier and tectonically controlled environment, which is profoundly influenced by climate driven processes. Its location within the equatorial realm supports recently proposed hypotheses, that topographic changes in this region must be considered as the dominant forcing factor influencing atmospheric circulation patterns and rainfall distribution. The uplift of this region has therefore often been associated with fundamental climatic and environmental changes in East Africa and adjacent regions. While the far-reaching influence of the plateau uplift is widely accepted, the timing and the magnitude of the uplift are ambiguous and are still subject to ongoing discussion. This dilemma stems from the lack of datable, geomorphically meaningful reference horizons that could record surface uplift. In order to quantify the amount of plateau uplift and to find evidence for the existence of significant relief along the East African Plateau prior to rifting, I analyzed and modeled one of the longest terrestrial lava flows; the 300-km-long Yatta phonolite flow in Kenya. This lava flow is 13.5 Ma old and originated in the region that now corresponds to the eastern rift shoulders. The phonolitic flow utilized an old riverbed that once drained the eastern flank of the plateau. Due to differential erosion this lava flow now forms a positive relief above the parallel-flowing Athi River, which is mimicking the course of the paleo-river. My approach is a lava-flow modeling, based on an improved composition and temperature dependent method to parameterize the flow of an arbitrary lava in a rectangular-shaped channel. The essential growth pattern is described by a one-dimensional model, in which Newtonian rheological flow advance is governed by the development of viscosity and/or velocity in the internal parts of the lava-flow front. Comparing assessments of different magma compositions reveal that length-dominated, channelized lava flows are characterized by high effusion rates, rapid emplacement under approximately isothermal conditions, and laminar flow. By integrating the Yatta lava flow dimensions and the covered paleo-topography (slope angle) into the model, I was able to determine the pre-rift topography of the East African Plateau. The modeling results yield a pre-rift slope of at least 0.2°, suggesting that the lava flow must have originated at a minimum elevation of 1,400 m. Hence, high topography in the region of the present-day Kenya Rift must have existed by at least 13.5 Ma. This inferred mid-Miocene uplift coincides with the two-step expansion of grasslands, as well as important radiation and speciation events in tropical Africa. Accordingly, the combination of my results regarding the Yatta lava flow emplacement history, its location, and its morphologic character, validates it as a suitable “paleo-tiltmeter” and has thus to be considered as an important topographic and volcanic feature for the topographic evolution in East Africa.
Cenozoic uplift of the East African Plateau has been associated with fundamental climatic and environmental changes in East Africa and adjacent regions. While this influence is widely accepted, the timing and the magnitude of plateau uplift have remained unclear. This uncertainty stems from the lack of datable, geomorphically meaningful reference horizons that could record surface uplift. Here, we document the existence of significant relief along the East African Plateau prior to rifting, as inferred from modeling the emplacement history of one of the longest terrestrial lava flows, the similar to 300-km-long Yatta phonolite flow in Kenya. This 13.5 Ma lava flow originated on the present-day eastern Kenya Rift flank, and utilized a riverbed that once routed runoff from the eastern rim of the plateau. Combining an empirical viscosity model with subsequent cooling and using the Yatta lava flow geometry and underlying paleotopography (slope angle), we found that the prerift slope was at least 0.2 degrees, suggesting that the lava flow originated at a minimum elevation of 1400 m. Hence, high paleotopography in the Kenya Rift region must have existed by at least 13.5 Ma. We infer from this that middle Miocene uplift occurred, which coincides with the two-step expansion of grasslands, as well as important radiation and speciation events in tropical Africa.
Timing and magnitude of surface uplift are key to understanding the impact of crustal deformation and topographic growth on atmospheric circulation, environmental conditions, and surface processes. Uplift of the East African Plateau is linked to mantle processes, but paleoaltimetry data are too scarce to constrain plateau evolution and subsequent vertical motions associated with rifting. Here, we assess the paleotopographic implications of a beaked whale fossil (Ziphiidae) from the Turkana region of Kenya found 740 km inland from the present-day coastline of the Indian Ocean at an elevation of 620 m. The specimen is similar to 17 My old and represents the oldest derived beaked whale known, consistent with molecular estimates of the emergence of modern straptoothed whales (Mesoplodon). The whale traveled from the Indian Ocean inland along an eastward-directed drainage system controlled by the Cretaceous Anza Graben and was stranded slightly above sea level. Surface uplift from near sea level coincides with paleoclimatic change from a humid environment to highly variable and much drier conditions, which altered biotic communities and drove evolution in east Africa, including that of primates.
Isostasy is one of the oldest and most widely applied concepts in the geosciences, but the geoscientific community lacks a coherent, easy-to-use tool to simulate flexure of a realistic (i.e., laterally heterogeneous) lithosphere under an arbitrary set of surface loads. Such a model is needed for studies of mountain building, sedimentary basin formation, glaciation, sea-level change, and other tectonic, geodynamic, and surface processes. Here I present gFlex (for GNU flexure), an open-source model that can produce analytical and finite difference solutions for lithospheric flexure in one (profile) and two (map view) dimensions. To simulate the flexural isostatic response to an imposed load, it can be used by itself or within GRASS GIS for better integration with field data. gFlex is also a component with the Community Surface Dynamics Modeling System (CSDMS) and Landlab modeling frameworks for coupling with a wide range of Earth-surface-related models, and can be coupled to additional models within Python scripts. As an example of this in-script coupling, I simulate the effects of spatially variable lithospheric thickness on a modeled Iceland ice cap. Finite difference solutions in gFlex can use any of five types of boundary conditions: 0-displacement, 0-slope (i.e., clamped); 0-slope, 0-shear; 0-moment, 0-shear (i.e., broken plate); mirror symmetry; and periodic. Typical calculations with gFlex require << 1 s to similar to 1 min on a personal laptop computer. These characteristics - multiple ways to run the model, multiple solution methods, multiple boundary conditions, and short compute time - make gFlex an effective tool for flexural isostatic modeling across the geosciences.
Alluvial and transport-limited bedrock rivers constitute the majority of fluvial systems on Earth. Their long profiles hold clues to their present state and past evolution. We currently possess first-principles-based governing equations for flow, sediment transport, and channel morphodynamics in these systems, which we lack for detachment-limited bedrock rivers. Here we formally couple these equations for transport-limited gravel-bed river long-profile evolution. The result is a new predictive relationship whose functional form and parameters are grounded in theory and defined through experimental data. From this, we produce a power-law analytical solution and a finite-difference numerical solution to long-profile evolution. Steady-state channel concavity and steepness are diagnostic of external drivers: concavity decreases with increasing uplift rate, and steepness increases with an increasing sediment-to-water supply ratio. Constraining free parameters explains common observations of river form: to match observed channel concavities, gravel-sized sediments must weather and fine - typically rapidly - and valleys typically should widen gradually. To match the empirical square-root width-discharge scaling in equilibrium-width gravel-bed rivers, downstream fining must occur. The ability to assign a cause to such observations is the direct result of a deductive approach to developing equations for landscape evolution.
Relative pollen productivity (RPP) estimates are fractionate values, often in relation to Poaceae, that allow vegetation cover to be estimated from pollen counts with the help of models. RPP estimates are especially used in the scientific community in Europe and China, with a few studies in North America. Here we present a comprehensive compilation of available northern hemispheric RPP studies and their results arising from 51 publications with 60 sites and 131 taxa. This compilation allows scientists to identify data gaps in need of further RPP analyses but can also aid them in finding an RPP set for their study region. We also present a taxonomically harmonised, unified RPP dataset for the Northern Hemisphere and subsets for North America (including Greenland), Europe (including arctic Russia), and China, which we generated from the available studies. The unified dataset gives the mean RPP for 55 harmonised taxa as well as fall speeds, which are necessary to reconstruct vegetation cover from pollen counts and RPP values. Data are openly available at https://doi.org/10.1594/PANGAEA.922661 (Wieczorek and Herzschuh, 2020).
Tree stands in the boreal treeline ecotone are, in addition to climate change, impacted by disturbances such as fire, water-related disturbances and logging. We aim to understand how these disturbances affect growth, age structure, and spatial patterns of larch stands in the north-eastern Siberian treeline ecotone (lower Kolyma River region), an insufficiently researched region. Stand structure of Larix cajanderi Mayr was studied at seven sites impacted by disturbances. Maximum tree age ranged from 44 to 300 years. Young to medium-aged stands had, independent of disturbance type, the highest stand densities with over 4000 larch trees per ha. These sites also had the highest growth rates for tree height and stem diameter. Overall lowest stand densities were found in a polygonal field at the northern end of the study area, with larches growing in distinct " tree islands". At all sites, saplings are significantly clustered. Differences in fire severity led to contrasting stand structures with respect to tree, recruit, and overall stand densities. While a low severity fire resulted in low-density stands with high proportions of small and young larches, high severity fires resulted in high-density stands with high proportions of big trees. At waterdisturbed sites, stand structure varied between waterlogged and drained sites and latitude. These mixed effects of climate and disturbance make it difficult to predict future stand characteristics and the treeline position.
Arctic and alpine treelines worldwide differ in their reactions to climate change. A northward advance of or densification within the treeline ecotone will likely influence climate-vegetation feedback mechanisms. In our study, which was conducted in the Taimyr Depression in the North Siberian Lowlands, w present a combined field-and model-based approach helping us to better understand the population processes involved in the responses of the whole treeline ecotone, spanning from closed forest to single-tree tundra, to climate warming. Using information on stand structure, tree age, and seed quality and quantity from seven sites, we investigate effects of intra-specific competition and seed availability on the specific impact of recent climate warming on larch stands. Field data show that tree density is highest in the forest-tundra, and average tree size decreases from closed forest to single-tree tundra. Age-structure analyses indicate that the trees in the closed forest and forest-tundra have been present for at least similar to 240 yr. At all sites except the most southerly ones, past establishment is positively correlated with regional temperature increase. In the single-tree tundra, however, a change in growth form from krummholz to erect trees, beginning similar to 130 yr ago, rather than establishment date has been recorded. Seed mass decreases from south to north, while seed quantity increases. Simulations with LAVESI (Larix Vegetation Simulator) further suggest that relative density changes strongly in response to a warming signal in the forest-tundra while intra-specific competition limits densification in the closed forest and seed limitation hinders densification in the single-tree tundra. We find striking differences in strength and timing of responses to recent climate warming. While forest-tundra stands recently densified, recruitment is almost non-existent at the southern and northern end of the ecotone due to autecological processes. Palaeo-treelines may therefore be inappropriate to infer past temperature changes at a fine scale. Moreover, a lagged treeline response to past warming will, via feedback mechanisms, influence climate change in the future.
This study monitors regional changes in the crystallinity of carbonaceous matter (CM) by applying Micro-Raman spectroscopy to a total of 214 metasediment samples (largely so-called Bundnerschiefer) dominantly metamorphosed under blueschist- to amphibolite-facies conditions. They were collected within the northeastern margin of the Lepontine dome and easterly adjacent areas of the Swiss Central Alps. Three-dimensional mapping of isotemperature contours in map and profile views shows that the isotemperature contours associated with the Miocene Barrow-type Lepontine metamorphic event cut across refolded nappe contacts, both along and across strike within the northeastern margin of the Lepontine dome and adjacent areas. Further to the northeast, the isotemperature contours reflect temperatures reached during the Late Eocene subduction-related blueschist-facies event and/or during subsequent near-isothermal decompression; these contours appear folded by younger, large-scale post-nappe-stacking folds. A substantial jump in the recorded maximum temperatures across the tectonic contact between the frontal Adula nappe complex and surrounding metasediments indicates that this contact accommodated differential tectonic movement of the Adula nappe with respect to the enveloping Bundnerschiefer after maximum temperatures were reached within the northern Adula nappe, i.e. after Late Eocene time.
The investigated HP/LT metasedimentary units of the Valaisan and adjacent European domains occupy a key position in the Alpine belt for understanding the transition from early subduction-related HP/LT metamorphism to collision-related Barrovian overprint and the evolution of mountain belts in general. The timing of high-pressure metamorphism, subsequent retrogression and following Barrow-type overprint was studied by Ar-40/Ar-39 dating of biotite and several white mica generations that are well characterized in terms of mineral chemistry, texture and associated mineral assemblages. Four distinct age populations of white mica record peak pressure conditions (42-40 Ma) and several stages of subsequent retrograde metamorphic evolution (36-25 Ma). Biotite isotopic analyses yield consistent apparent ages that cluster around 18-16 Ma for the Barrow-type thermal overprint. The recorded isotopic data reveal a significant time gap in the order of some 20 Ma between subduction-related HP/LT metamorphism and collision-related Barrovian overprint, supporting the notion of a polymetamorphic evolution associated with a bimodal P-T path.
A common problem in ecology is identifying the relationship between relief and site properties obtainable only by point measurements. The method of Multi-Scale Landscape Analysis (MSLA) identifies such correlations. MSLA combines frequency filtering of the digital elevation model (DEM) with an estimation of the optimum filter coefficients using an optimization procedure. Tested using point data of soil decarbonation from a German young moraine landscape, MSLA provided significant results. Implemented within open source software SAMT. MSLA is comfortable and flexible to use, offering applications for numerous other spatial analysis problems.
Challenges in understanding the hydrologic controls on the mobility of slow-moving landslides
(2011)
Slow-moving landslides are a wide-spread type of active mass movement, can cause severe damages to infrastructure, and may be a precursor of sudden catastrophic slope failures. Pore-water pressure is commonly regarded as the most important among a number of possible factors controlling landslide velocity. We used high-resolution monitoring data to explore the relations of landslide mobility and hydrologic processes at the Heumoser landslide in Austria, which is characterized by continuous slow movement along a shear zone. Movement rates showed a seasonality that was associated with elevated pore-water pressures. Pore pressure monitoring revealed a system of confined and separated aquifers with differing dynamics. Analysis of a simple infinite slope mobility model showed that small variations in parameters, along with measured pore pressure dynamics, provided a perfect match to our observations. Modeling showed a stabilizing effect of snow cover due to the additional load. This finding was supported by a multiple regression model, which further suggested that effective pore pressures at the slip surface were partially differing from the borehole observations and were related to preferential infiltration and subsurface flow in adjacent areas. It appears that in a setting like the Heumoser landslide, hydrologic processes delicately influence slope mobility through their control on pore pressure dynamics and the weight of the landslide body, which challenges observation and modeling. Moreover, it appears that their simplicity, and especially their high sensitivity to parameter variations, limits the conclusions that can be drawn from infinite slope models.
Agricultural soils have a high potential for sequestration of atmospheric carbon due to their volume and several promising management options. However, there is a remarkable lack of information about the status quo of organic carbon in agricultural soils. In this study a comprehensive data set of 384 cropland soils and 333 grassland soils within the state of Bavaria in southeast Germany was analyzed in order to provide representative information on total amount, regional distribution and driving parameters of soil organic carbon (SOC) and nitrogen (N) in agricultural soils of central Europe. The results showed that grassland soils stored higher amounts of SOC (11.8 kg m(-2)) and N (0.92 kg m(-2)) than cropland soils (9.0 and 0.66 kg m(-2), respectively) due to moisture-induced accumulation of soil organic matter (SOM) in B horizons. Surprisingly, no distinct differences were found for the A horizons since tillage led to a relocation of SOM with depth in cropland soils. Statistical analyses of driving factors for SOM storage revealed soil moisture, represented by the topographic wetness index (TWI), as the most important parameter for both cropland and grassland soils. Climate effects (mean annual temperature and precipitation) were of minor importance in agricultural soils because management options counteracted them to a certain extent, particularly in cropland soils. The distribution of SOC and N stocks within Bavaria based on agricultural regions confirmed the importance of soil moisture since the highest cropland SOC and N stocks were found for tertiary hills and loess regions, which exhibited large areas with potentially high soil moisture content in extant floodplains. Grassland soils showed the highest accumulation of SOC and N in the Alps and Pre-Alps as a result of low temperatures, high amounts of precipitation and high soil moisture content in areas of glacial denudation. Soil class was identified as a further driving parameter for SOC and N storage in cropland soils. In total, cropland and grassland soils in Bavaria store 242 and 134 Mt SOC as well as 19 and 12 Mt N down to a soil depth of 1 m or the parent material, respectively.
Carbon storage capacity of semi-arid grassland soils and sequestration potentials in northern China
(2015)
Organic carbon (OC) sequestration in degraded semi-arid environments by improved soil management is assumed to contribute substantially to climate change mitigation. However, information about the soil organic carbon (SOC) sequestration potential in steppe soils and their current saturation status remains unknown. In this study, we estimated the OC storage capacity of semi-arid grassland soils on the basis of remote, natural steppe fragments in northern China. Based on the maximum OC saturation of silt and clay particles <20m, OC sequestration potentials of degraded steppe soils (grazing land, arable land, eroded areas) were estimated. The analysis of natural grassland soils revealed a strong linear regression between the proportion of the fine fraction and its OC content, confirming the importance of silt and clay particles for OC stabilization in steppe soils. This relationship was similar to derived regressions in temperate and tropical soils but on a lower level, probably due to a lower C input and different clay mineralogy. In relation to the estimated OC storage capacity, degraded steppe soils showed a high OC saturation of 78-85% despite massive SOC losses due to unsustainable land use. As a result, the potential of degraded grassland soils to sequester additional OC was generally low. This can be related to a relatively high contribution of labile SOC, which is preferentially lost in the course of soil degradation. Moreover, wind erosion leads to substantial loss of silt and clay particles and consequently results in a direct loss of the ability to stabilize additional OC. Our findings indicate that the SOC loss in semi-arid environments induced by intensive land use is largely irreversible. Observed SOC increases after improved land management mainly result in an accumulation of labile SOC prone to land use/climate changes and therefore cannot be regarded as contribution to long-term OC sequestration.
Temperate forest soils of central Europe are regarded as important pools for soil organic carbon (SOC) and thought to have a high potential for carbon (C) sequestration. However, comprehensive data on total SOC storage, particularly under different forest types, and its drivers is limited. In this study, we analyzed a forest data set of 596 completely sampled soil profiles down to the parent material or to a depth of 1 m within Bavaria in southeast Germany in order to determine representative SOC stocks under different forest types in central Europe and the impact of different environmental parameters. We calculated a total median SOC stock of 9.8 kg m(-2) which is considerably lower compared with many other inventories within central Europe that used modelled instead of measured soil properties. Statistical analyses revealed climate as controlling parameter for the storage of SOC with increasing stocks in cool, humid mountainous regions and a strong decrease in areas with higher temperatures. No significant differences of total SOC storage were found between broadleaf, coniferous and mixed forests. However, coniferous forests stored around 35% of total SOC in the labile organic layer that is prone to human disturbance, forest fires and rising temperatures. In contrast, mixed and broadleaf forests stored the major part of SOC in the mineral soil. Moreover, these two forest types showed unchanged or even slightly increased mineral SOC stocks with higher temperatures, whereas SOC stocks in mineral soils under coniferous forest were distinctly lower. We conclude that mixed and broadleaf forests are more advantageous for C sequestration than coniferous forests. An intensified incorporation of broadleaf species in extent coniferous forests of Bavaria would prevent substantial SOC losses as a result of rising temperatures in the course of climate change.
Different upper tail indicators exist to characterize heavy tail phenomena, but no comparative study has been carried out so far. We evaluate the shape parameter (GEV), obesity index, Gini index and upper tail ratio (UTR) against a novel benchmark of tail heaviness - the surprise factor. Sensitivity analyses to sample size and changes in scale-to-location ratio are carried out in bootstrap experiments. The UTR replicates the surprise factor best but is most uncertain and only comparable between records of similar length. For samples with symmetric Lorenz curves, shape parameter, obesity and Gini indices provide consistent indications. For asymmetric Lorenz curves, however, the first two tend to overestimate, whereas Gini index tends to underestimate tail heaviness. We suggest the use of a combination of shape parameter, obesity and Gini index to characterize tail heaviness. These indicators should be supported with calculation of the Lorenz asymmetry coefficients and interpreted with caution.
Land-use concepts provide decision support for the most efficient usage options according to sustainable development and multifunctionality requirements. However, developments in landscape-related, agricultural production schemes are primarily driven by economic benefits. Therefore, most agricultural land-use concepts tackle particular problems or interests and lack a systemic perspective. As a result, we discuss a conceptual model for future site-specific agricultural land-use with an inbuilt requirement for adequate experimental sites to enable monitoring systems for a new generation of ecosystem models and for new approaches to address science-stakeholder interactions.
The "Bohemians" of New Zealand - an ethnic group? In 1982 a small group of New Zealanders established contacts with the region of origin of their about 200 German-speaking ancestors who had emigrated from Bohemia for economic reasons in the 1860s and 1870s. They had all come from about twenty villages situated west of plzen and founded a rural settlement and participated in the foundation of a second one in New Zealand. Since World War I there had been no further contacts between the emigrants and their descendants on the one hand and their relatives in Bohemia on the other hand. For two reasons new contacts were established after such a long time: (1) the back-to the-roots-movement had spread to New Zealand from the USA, Canada and Australia, and (2) the status of cultural diversity keeps being enhanced in New Zealand since about 1970. These processes also influenced those people in New Zealand who call themselves "Bohemians" because of their ancestors' region of origin. Their total number is estimated at 10, 000 to 15, 000 at present. Up to now hardly any attention was attached to them in New Zealand by academic research and/or the general public. This paper discusses the history and today's situation of the former immigrants' community as well as the New Zealand "Bohemians " in general, raising the questions to what extent they can be defined as an ethnic group now and whether they will retain their status as a specific group in future
Understanding the key factors influencing the water quality of large river systems forms an important basis for the assessment and protection of cross-regional ecosystems and the implementation of adapted water management concepts. However, identifying these factors requires in-depth comprehension of the unique environmental systems, which can only be achieved by detailed water quality monitoring.
Within the scope of the joint science and sports event "Elbschwimmstaffel" (swimming relay on the river Elbe) in June/July 2017 organized by the German Ministry of Education and Research, water quality data were acquired along a 550 km long stretch of the Elbe River in Germany. During the survey, eight physiochemical water quality parameters were recorded in high spatial and temporal resolution with the BIOFISH multisensor system. Multivariate statistical methods were applied to identify and delineate processes influencing the water quality.
The BIOFISH dataset revealed that phytoplankton activity has a major impact on the water quality of the Elbe River in the summer months. The results suggest that phytoplankton biomass constitutes a substantial proportion of the suspended particles and that photosynthetic activity of phytoplankton is closely related to significant temporal changes in pH and oxygen saturation.
An evaluation of the BIOFISH data based on the combination of statistical analysis with weather and discharge data shows that the hydrological and meteorological history of the sampled water body was the main driver of phytoplankton dynamics. This study demonstrates the capacity of longitudinal river surveys with the BIOFISH or similar systems for water quality assessment, the identification of pollution sources and their utilization for online in situ monitoring of rivers.
In 1915, Alfred Wegener published his hypotheses of plate tectonics that revolutionised the world for geologists. Since then, many scientists have studied the evolution of continents and especially the geologic structure of orogens: the most visible consequence of tectonic processes. Although the morphology and landscape evolution of mountain belts can be observed due to surface processes, the driving force and dynamics at lithosphere scale are less well understood despite the fact that rocks from deeper levels of orogenic belts are in places exposed at the surface. In this thesis, such formerly deeply-buried (ultra-) high-pressure rocks, in particular eclogite facies series, have been studied in order to reveal details about the formation and exhumation conditions and rates and thus provide insights into the geodynamics of the most spectacular orogenic belt in the world: the Himalaya. The specific area investigated was the Kaghan Valley in Pakistan (NW Himalaya). Following closure of the Tethyan Ocean by ca. 55-50 Ma, the northward subduction of the leading edge of India beneath the Eurasian Plate and subsequent collision initiated a long-lived process of intracrustal thrusting that continues today. The continental crust of India – granitic basement, Paleozoic and Mesozoic cover series and Permo-Triassic dykes, sills and lavas – has been buried partly to mantle depths. Today, these rocks crop out as eclogites, amphibolites and gneisses within the Higher Himalayan Crystalline between low-grade metamorphosed rocks (600-640°C/ ca. 5 kbar) of the Lesser Himalaya and Tethyan sediments. Beside tectonically driven exhumation mechanisms the channel flow model, that describes a denudation focused ductile extrusion of low viscosity material developed in the middle to lower crust beneath the Tibetan Plateau, has been postulated. To get insights into the lithospheric and crustal processes that have initiated and driven the exhumation of this (ultra-) high-pressure rocks, mineralogical, petrological and isotope-geochemical investigations have been performed. They provide insights into 1) the depths and temperatures to which these rocks were buried, 2) the pressures and temperatures the rocks have experienced during their exhumation, 3) the timing of these processes 4) and the velocity with which these rocks have been brought back to the surface. In detail, through microscopical studies, the identification of key minerals, microprobe analyses, standard geothermobarometry and modelling using an effective bulk rock composition it has been shown that published exhumation paths are incomplete. In particular, the eclogites of the northern Kaghan Valley were buried to depths of 140-100 km (36-30 kbar) at 790-640°C. Subsequently, cooling during decompression (exhumation) towards 40-35 km (17-10 kbar) and 630-580°C has been superseded by a phase of reheating to about 720-650°C at roughly the same depth before final exhumation has taken place. In the southern-most part of the study area, amphibolite facies assemblages with formation conditions similar to the deduced reheating phase indicate a juxtaposition of both areas after the eclogite facies stage and thus a stacking of Indian Plate units. Radiometric dating of zircon, titanite and rutile by U-Pb and amphibole and micas by Ar-Ar reveal peak pressure conditions at 47-48 Ma. With a maximum exhumation rate of 14 cm/a these rocks reached the crust-mantle boundary at 40-35 km within 1 Ma. Subsequent exhumation (46-41 Ma, 40-35 km) decelerated to ca. 1 mm/a at the base of the continental crust but rose again to about 2 mm/a in the period of 41-31 Ma, equivalent to 35-20 km. Apatite fission track (AFT) and (U-Th)/He ages from eclogites, amphibolites, micaschists and gneisses yielded moderate Oligocene to Miocene cooling rates of about 10°C/Ma in the high altitude northern parts of the Kaghan Valley using the mineral-pair method. AFT ages are of 24.5±3.8 to 15.6±2.1 Ma whereas apatite (U-Th)/He analyses yielded ages between 21.0±0.6 and 5.3±0.2 Ma. The southern-most part of the Valley is dominated by younger late Miocene to Pliocene apatite fission track ages of 7.6±2.1 and 4.0±0.5 Ma that support earlier tectonically and petrologically findings of a juxtaposition and stack of Indian Plate units. As this nappe is tectonically lowermost, a later distinct exhumation and uplift driven by thrusting along the Main Boundary Thrust is inferred. A multi-stage exhumation path is evident from petrological, isotope-geochemical and low temperature thermochronology investigations. Buoyancy driven exhumation caused an initial rapid exhumation: exhumation as fast as recent normal plate movements (ca. 10 cm/a). As the exhuming units reached the crust-mantle boundary the process slowed down due to changes in buoyancy. Most likely, this exhumation pause has initiated the reheating event that is petrologically evident (e.g. glaucophane rimmed by hornblende, ilmenite overgrowth of rutile). Late stage processes involved widespread thrusting and folding with accompanied regional greenschist facies metamorphism, whereby contemporaneous thrusting on the Batal Thrust (seen by some authors equivalent to the MCT) and back sliding of the Kohistan Arc along the inverse reactivated Main Mantle Thrust caused final exposure of these rocks. Similar circumstances have been seen at Tso Morari, Ladakh, India, 200 km further east where comparable rock assemblages occur. In conclusion, as exhumation was already done well before the initiation of the monsoonal system, climate dependent effects (erosion) appear negligible in comparison to far-field tectonic effects.