Refine
Year of publication
Document Type
- Article (2601) (remove)
Language
- English (2601) (remove)
Keywords
- Holocene (35)
- climate change (30)
- Earthquake source observations (21)
- Climate change (18)
- Pollen (18)
- permafrost (18)
- Himalaya (17)
- Tibetan Plateau (17)
- erosion (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2601) (remove)
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 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.
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 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.
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.
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.
Metabasites were sampled from rock series of the subducted margin of the Indian Plate, the so-called Higher Himalayan Crystalline, in the Upper Kaghan Valley, Pakistan. These vary from corona dolerites, cropping out around Saif- ul-Muluk in the south, to coesite-eclogite close to the suture zone against rocks of the Kohistan arc in the north. Bulk rock major- and trace-element chemistry reveals essentially a single protolith as the source for five different eclogite types, which differ in fabric, modal mineralogy as well as in mineral chemistry. The study of newly-collected samples reveals coesite (confirmed by in situ Raman spectroscopy) in both garnet and omphacite. All eclogites show growth of amphiboles during exhumation. Within some coesite-bearing eclogites the presence of glaucophane cores to barroisite is noted whereas in most samples porphyroblastic sodic-calcic amphiboles are rimmed by more aluminous calcic amphibole (pargasite, tschermakite, and edenite). Eclogite facies rutile is replaced by ilmenite which itself is commonly surrounded by titanite. In addition, some eclogite bodies show leucocratic segregations containing phengite, quartz, zoisite and/or kyanite. The important implication is that the complex exhumation path shows stages of initial cooling during decompression (formation of glaucophane) followed by reheating: a very similar situation to that reported for the coesite-bearing eclogite series of the Tso Morari massif, India, 450 km to the south-east.