Filtern
Erscheinungsjahr
- 2021 (130) (entfernen)
Dokumenttyp
- Wissenschaftlicher Artikel (85)
- Dissertation (26)
- Postprint (16)
- Monographie/Sammelband (1)
- Habilitation (1)
- Masterarbeit (1)
Gehört zur Bibliographie
- ja (130)
Schlagworte
- PHREEQC (5)
- reactive transport (4)
- ALOS World 3D (2)
- ASTER GDEM (2)
- Alpen (2)
- Alps (2)
- Arctic Ocean (2)
- Argentina (2)
- Argentinien (2)
- Big Naryn complex (2)
- Copernicus DEM (2)
- DEM noise (2)
- Environmental sciences (2)
- Exhumation (2)
- Forschungsmethodik (2)
- Fourier analysis (2)
- GNSS-integrated water vapour (2)
- Germany (2)
- Global Navigation Satellite System (GNSS) (2)
- HVSR (2)
- Kuilyu complex (2)
- Nearshore zone (2)
- Ocean color remote sensing (2)
- Ocean sciences (2)
- Paläoklimadynamik (2)
- Rheologie (2)
- SRTM (2)
- Solid Earth sciences (2)
- Soziale Medien (2)
- TRMM (2)
- TanDEM-X (2)
- Tianshan orogenic belt (2)
- Tonian–Cryogenian (2)
- Turbidity retrieval (2)
- Twitter (2)
- Unterrichtsmethoden (2)
- WorldDEM (2)
- Younger Dryas (2)
- blueSeis-3A sensors (2)
- calc-alkaline (2)
- central Poland (2)
- climate change (2)
- climate change adaptation (2)
- coherency (2)
- coherent noise (2)
- convective available potential energy (CAPE) (2)
- correlation (2)
- damage (2)
- elastic rock properties (2)
- environmental seismology (2)
- extreme rainfall (2)
- fatigue (2)
- fundamental frequency (2)
- governance (2)
- heterogeneity (2)
- heuristic analysis (2)
- huddle test (2)
- innovation (2)
- insurance (2)
- lake monitoring (2)
- lake sediments (2)
- magmatic arc (2)
- mass (2)
- mountain geomorphology (2)
- natural hazard (2)
- noise cross (2)
- numerical modeling (2)
- permafrost (2)
- preparedness (2)
- recovery (2)
- repository far-field (2)
- rotational seismology (2)
- sedimentation pattern (2)
- seismic monitoring (2)
- slope failure (2)
- social (2)
- sorption (2)
- source direction (2)
- time series (2)
- transferability (2)
- urban infill development (2)
- wasting (2)
- weathering (2)
- (compliant) cracks (1)
- 3D root (1)
- Abbaufrontkartierung (1)
- Abbildende Spektroskopie (1)
- African climate (1)
- Afrikanisches Klima (1)
- Afyon volcanics (1)
- AlpArray (1)
- Altersmodelierung (1)
- Aluminium– silicates (1)
- Andean plateau (1)
- Anden (1)
- Andes (1)
- Anthropocene (1)
- Antibiotikaresistenz (1)
- Aphanizomenon (1)
- Ar-39 (1)
- Ar-40 (1)
- Arctic (1)
- Arctic climate change (1)
- Arktis (1)
- Bakterien (1)
- Balanced cross section (1)
- Baltic Sea coast (1)
- Biodiversität (1)
- Biogeowissenschaften (1)
- Biomarker (1)
- Boundary element method (1)
- Bruchmechanik (1)
- Bushveld Complex (1)
- Caimancito oil field (1)
- Caimancito-Ölfeld (1)
- Campo petrolero Caimancito (1)
- Central Asia (1)
- Chirete (1)
- Cinética de fases (1)
- Cinética del querógeno (1)
- Contamination (1)
- Corundum (1)
- Cretaceous basin (1)
- Crystal nucleation (1)
- Cuenca Cretácica (1)
- Cuenca intramontana (1)
- Cuesta de la chilca (1)
- Cyanobium (1)
- Cyclostratigraphy (1)
- Dam engineering (1)
- Deformación cuaternaria (1)
- Deformation (1)
- Deuterium Exzesses (1)
- Diaguita (1)
- Dichtemodellierung (1)
- Digitale Gesteinsphysik (1)
- Drucklösungsprozesse (1)
- Duplex (1)
- Dynamic bulk modulus (1)
- Décollement (1)
- EXAFS (1)
- Earthquake forecasting (1)
- Earthquake ground motions (1)
- Earthquake hazards (1)
- Earthquake source observations (1)
- Earthquakes (1)
- East African Rift (1)
- Einengungsraten (1)
- Elastic response (1)
- Elastische Gesteinseigenschaften (1)
- Energy and society (1)
- Equant (stiff) pores (1)
- Erdbeben (1)
- Erdbebenvorhersage (1)
- Erdbeobachtung (1)
- Europe (1)
- Experimental (1)
- Experimental petrology (1)
- Famatinian (1)
- Fernerkundung (1)
- Fluid (1)
- Fluid-Gesteinswechselwirkungen (1)
- Fluide (1)
- Fluidströmung (1)
- Flüsse (1)
- Fold and thrust belt (1)
- Formación Yacoraite (1)
- Fracture mechanics (1)
- Fracture sealing (1)
- Fracture wall offset (1)
- GDGT (1)
- Gasgeochemie (1)
- Geochemical indices (1)
- Geochemical modelling (1)
- Geochemistry (1)
- Geology (1)
- Geomorphologie (1)
- Geophysics (1)
- Geothermal (1)
- Geothermal monitoring (1)
- Geothermisches Monitoring (1)
- Geothermobarometry (1)
- Glaciers (1)
- Global earthquake data (1)
- Globale Erdbebenkatalogdaten (1)
- Gläser (1)
- Granite (1)
- Gravity (1)
- Grundwasser-Oberflächenwasser-Interaktion (1)
- Heavy metals (1)
- Hierarchical clustering (1)
- High-pressure (1)
- High-temperature experiments (1)
- Himalaya (1)
- Hochdruck (1)
- Holocene (1)
- Hydraulic aperture (1)
- Hydro-meteorological indices (1)
- Hämolyse (1)
- IHPV (1)
- Interseismic strain rates (1)
- Interseismiche Dehnungsraten (1)
- Intramontane basin (1)
- Intramontanes Becken (1)
- Inverse theory (1)
- Joint inversion (1)
- Karbonatschmelze (1)
- Kenia (1)
- Kenya (1)
- Kenya Rift (1)
- Kerogenkinetik (1)
- Kinematics of crustal and mantle deformation (1)
- Klimawandel (1)
- Krankheitserreger (1)
- Kreidebecken (1)
- Kurpsai dam (1)
- LGM (1)
- Lake Issyk‐ Kul (1)
- Lake-surfaces variability (1)
- Landfill (1)
- Landslide inventory (1)
- Landslide susceptibility (1)
- Late Holocene (1)
- Leones Glacier (1)
- Lithospheric strength (1)
- Lithosphärenfestigkeit (1)
- Logistic regression (1)
- Lower crust (1)
- MATLAB (1)
- MC-ICP-MS (1)
- Marmarameer (1)
- Martinez del Tineo (1)
- Matrix diffusion (1)
- Mechanical aperture (1)
- Mechanical properties (1)
- Melainabacteria (1)
- Melilitit (1)
- Merensky Reef (1)
- Metastability (1)
- Methylheptadecanes (1)
- Micro-structure (1)
- Mikro-CT (1)
- Minor (1)
- Model (1)
- Modell (1)
- Moment tensor inversion (1)
- Momententensorinversion (1)
- Monsoonal variability (1)
- Mont Terri (1)
- Multivariate statistical analysis (1)
- Nanolayers (1)
- Nanorocks (1)
- Narrow (1)
- Neogene (1)
- Neotektonik (1)
- Nordanatolische Störungszone (1)
- Normal (1)
- North (1)
- North Anatolian Fault Zone (1)
- North Atlantic climate change (1)
- Northern Patagonian Icefield (1)
- Numerische Modellierung (1)
- Operational and environmental effects (1)
- Orbital forcing (1)
- Organofazies (1)
- Ostafrikanisches Rift (1)
- Ostseeküste (1)
- P-waves (1)
- PGE (1)
- PGM (1)
- PT path (1)
- Paleoclimate (1)
- Paleoclimate dynamics (1)
- Paleoclimatic dynamics (1)
- Paläohydrologie (1)
- Paläoökologie (1)
- Panasqueira (1)
- Passive Kontinentalränder (1)
- Passive margins (1)
- Peridotitic melts (1)
- Permafrost-Taustörungen (1)
- Permafrostlandschaften (1)
- Permeabilität (1)
- Permeabilitätsentwicklung (1)
- PhaseKinetics (1)
- Phasenkinetik (1)
- Piedmont Zone (1)
- Pilze (1)
- Pitzer (1)
- Planktothrix (1)
- Plate tectonics (1)
- Plio-Pleistocene timeframe (1)
- Plio-Pleistozän Zeitfenster (1)
- Pollen (1)
- Pollenproduktivitätsschätzungen (1)
- Porous sandstone (1)
- Porphyrische Kupferlagerstätte (1)
- Poshtuk (1)
- Precambrain (1)
- Precipitation variability (1)
- Probability distributions (1)
- Quartz (1)
- Quartäre Deformation (1)
- Quaternary deformation (1)
- REVEALS (1)
- Ramos X-11 (1)
- Randelementmethode (1)
- Raoult's law (1)
- Reactive transport (1)
- Regional forcings (1)
- Reifegradbezogene Biomarker (1)
- Remote Sensing (1)
- Remote sensing (1)
- Rheological modeling (1)
- Rheologische Modellierung (1)
- Rheology (1)
- Rheology: crust and (1)
- Risk assessment (1)
- Riss (1)
- Río Pescado (1)
- S-waves (1)
- Salt Range (1)
- Salta basin (1)
- Salta-Becken (1)
- Sandstein (1)
- Savannen-Ökologie (1)
- Schwerkraft (1)
- Sea of Marmara (1)
- Seismic cycle (1)
- Seismic noise (1)
- Seismotektonik (1)
- Self-propping fracture (1)
- Sericytochromatia (1)
- Shan (1)
- Shortening rates (1)
- Sierras Pampeanas Orientales (1)
- Silikatschmelze (1)
- Site effects (1)
- Slope (1)
- Snow-cover (1)
- Snowella (1)
- Soil contamination (1)
- Source longevity (1)
- South American Monsoon (1)
- South Korea (1)
- Southern Kyrgyzstan (1)
- Space geodetic surveys (1)
- Spectral analysis (1)
- Spektroskopie (1)
- Static bulk modulus (1)
- Statistical seismology (1)
- Statistische Seismologie (1)
- Strahlung und chemische Eigenschaften (1)
- Strain (1)
- Strong-motion (1)
- Structural health monitoring (1)
- Strukturgeologie (1)
- Strukturmodellierung (1)
- Sturmhochwasser (1)
- Superconducting gravimetry (1)
- Supraleit-Gravimetrie (1)
- Surface roughness (1)
- Synthetische Sandsteine (1)
- System (1)
- Tasas de acortamiento (1)
- Tektonik (1)
- Thermal (1)
- Thermal modeling (1)
- Thermal treatment (1)
- Thermisch (1)
- Thermische Modellierung (1)
- Thermotoleranz (1)
- Tian (1)
- Time series analysis (1)
- Tungsten–tin deposits (1)
- UV fs laser ablation (1)
- Uk’37 (1)
- Unloading (1)
- Upper Indus Basin (1)
- Varves (1)
- Verdichtung (1)
- Versalzung (1)
- Virulenz (1)
- Volatiles (1)
- Volcanology (1)
- Wave propagation (1)
- Waveform inversion (1)
- Western Alps (1)
- XANES (1)
- XAS (1)
- Yacoraite Formation (1)
- ZHe tracer thermochronology (1)
- Zeitreihenanalyse (1)
- Zeitserie (1)
- Zustands-Übergangs-Modelle (1)
- Zyklostratigraphie (1)
- acoustic emissions (1)
- adaptation (1)
- age modeling (1)
- age-related biomarkers (1)
- algorithms (1)
- alteration geochemistry (1)
- amplicon sequencing (1)
- analog experiment (1)
- antibiotic resistance (1)
- antimicrobial resistance (1)
- arctic (1)
- arctic water bodies (1)
- arktische Klimaänderung (1)
- asymmetric rifting (1)
- attributes (1)
- automated monitoring (1)
- bacteria (1)
- basin analysis (1)
- biodegradación (1)
- biodegradation (1)
- biodiversity (1)
- biogeoscience (1)
- biologischer Abbau (1)
- biomarcadores diagnósticos de edad (1)
- biomarcadores diagnósticos de madurez (1)
- biomarker (1)
- body waves (1)
- broadband seismometer (1)
- bubble trap (1)
- calderas (1)
- carbon stabilization (1)
- carbonate melt (1)
- carnallite (1)
- cementation (1)
- change detection (1)
- chloritoid stability (1)
- climate-change mitigation (1)
- climate-change policy (1)
- coexisting phases (1)
- components (1)
- compression (1)
- continental rifts (1)
- convection (1)
- convolutional neural network (1)
- core-shell (1)
- cyanobacteria sedimentation (1)
- d-excess (1)
- damage assessment (1)
- debris flows (1)
- deep learning (1)
- deformation (1)
- density modeling (1)
- density-driven (1)
- density-driven flow (1)
- diffusion (1)
- digital core reconstruction (1)
- digital grain sizing (1)
- digital rock physics (1)
- dike (1)
- dike propagation (1)
- downstream fining (1)
- dryland (1)
- earthquake source (1)
- earthquakes (1)
- ecological succession (1)
- elastic waves (1)
- emergency (1)
- emergency preparedness (1)
- energy transfer (1)
- environmental reconstruction (1)
- equations (1)
- erosion (1)
- eruptions (1)
- eruptive cycle (1)
- exhumation (1)
- extracellular DNA (1)
- facies (1)
- faltendes neuronales Netzwerk (1)
- fault bend fold (1)
- flood events (1)
- flood risk management (1)
- fluid flow (1)
- fluid-rock interactions (1)
- fluvial geomorphology (1)
- formation damage (1)
- fracture growth (1)
- freeze-thaw cycling (1)
- frost‐ cracking (1)
- full (1)
- fungi (1)
- gas geochemistry (1)
- gas supply conditions (1)
- geodynamics (1)
- geologic hazards (1)
- geological hyperspectral image classification (1)
- geologische hyperspektrale Bildklassifikation (1)
- geomorphology (1)
- geothermal energy (1)
- geothermal exploration (1)
- geothermal monitoring (1)
- geothermische Exploration (1)
- geothermische Überwachung (1)
- geyser (1)
- glacial erosion (1)
- glasses (1)
- grain (1)
- grain-size distribution (1)
- ground-penetrating radar (1)
- groundwater-surface water (1)
- hemolysis (1)
- high pressure (1)
- high temperature (1)
- high-pressure (1)
- hydrate formation process (1)
- hydraulic fracturing (1)
- hydro-geological system (1)
- hydrothermal systems (1)
- hyperspectral (1)
- hyperspektral (1)
- imaging spectroscopy (1)
- in (1)
- inheritance (1)
- injections (1)
- injectivity (1)
- interaction (1)
- interpretation (1)
- intracellular DNA (1)
- kerogen kinetics (1)
- lake stratification (1)
- laminography (1)
- landsat (1)
- lanthanides (1)
- lanthanoid migration (1)
- lithobiont (1)
- lithosphere (1)
- local structure (1)
- location (1)
- lokale Struktur (1)
- machine learning (1)
- magma ocean (1)
- magmatic-hydrothermal systems (1)
- magnetostratigraphy (1)
- manaderos de petróleo (1)
- matríz mineral (1)
- maturity-related biomarkers (1)
- melilitite (1)
- meta-ophiolites (1)
- meta-pelites (1)
- micro-CT (1)
- micro-CT scan (1)
- mine face mapping (1)
- mineral matrix (1)
- mineralische Matrix (1)
- mixed gas hydrates (1)
- modeling (1)
- moment tensor (1)
- monsunale Variabilität (1)
- multi-spectral analysis (1)
- multivariate Statistik (1)
- multivariate statistics (1)
- multi‐ temporal analyses (1)
- n-alkane (1)
- n-alkanes (1)
- natural hazards (1)
- nature-based solutions (NBS) (1)
- natürlichen Ölaustritten (1)
- neotectonics (1)
- neutron imaging (1)
- nordatlantische Klimaänderung (1)
- northern high latitudes (1)
- numerical (1)
- numerical model (1)
- nördliche hohe Breitengrade (1)
- oberflächennahe Geophysík (1)
- offener Tagebau (1)
- oil seeps (1)
- open pit mining (1)
- organofacies (1)
- paleoecology (1)
- paleohydrology (1)
- pathogens (1)
- pebblecounts (1)
- permafrost landscapes (1)
- permafrost lowlands (1)
- permafrost thaw (1)
- permafrost thaw disturbances (1)
- permeability (1)
- permeability evolution (1)
- permeability-porosity relationship (1)
- phreeqc (1)
- picocyanobacteria diversity (1)
- pollen (1)
- pore-scale (1)
- porous media (1)
- porphyry copper deposit (1)
- pressure solution (1)
- pseudosection (1)
- psychrotolerant (1)
- radiación y propiedades químicas (1)
- radial flow (1)
- radiation and chemical properties (1)
- radiolyse (1)
- radiolysis (1)
- radiólisis (1)
- rainfall events (1)
- ramp (1)
- redox processes (1)
- relief (1)
- response (1)
- retrogressive thaw slumps (1)
- rheology (1)
- rhizosphere biogeochemistry (1)
- rift variability (1)
- rifting (1)
- rivers (1)
- rock (1)
- rock fracture (1)
- rockfalls (1)
- root activity (1)
- root-soil interaction (1)
- salinization (1)
- salt dissolution (1)
- sandstones (1)
- savanna ecology (1)
- sedaDNA (1)
- sediment production (1)
- sediment traps (1)
- sedimentary ancient DNA (1)
- sedimentology (1)
- sedimentäre alte DNA (1)
- seismic interpretation (1)
- seismic waves (1)
- seismische Interpretation (1)
- seismology (1)
- seismotectonic (1)
- sentinel-2 (1)
- shallow geophysics (1)
- shoreline movement direction (1)
- shoreline movement rate (1)
- signal-to-noise ratio (1)
- silicate melt (1)
- silicate melts (1)
- silikatische Schmelzen (1)
- simulation (1)
- situ Raman spectroscopy (1)
- size fractions (1)
- socioeconomic scenarios (1)
- solar radiation (1)
- source (1)
- spatial correlation (1)
- spectroscopy (1)
- stabilen Isotopen (1)
- stable isotope (1)
- state-transition models (1)
- storm flood (1)
- stress (1)
- stress field (1)
- structural (1)
- structural geology (1)
- structural modelling (1)
- subarctic (1)
- sulfur solubility (1)
- supergene ores (1)
- sustainability (1)
- swarms (1)
- synthetic sandstone (1)
- system architecture (1)
- tectonics (1)
- temperature (1)
- temperature extremes (1)
- thaw (1)
- thermal (1)
- thermal field (1)
- thermisches Feld (1)
- thermo-erosion (1)
- thermotolerance (1)
- tilt (1)
- topographic asymmetry (1)
- topography (1)
- tourmaline (1)
- trace elements (1)
- trajectories (1)
- tremor (1)
- upconversion nanoparticles (1)
- uranium speciation (1)
- vegetation cover (1)
- very high resolution imagery (1)
- virulence (1)
- volcanic earthquakes (1)
- water rock interaction (1)
- waveform (1)
- white mica (1)
- Þeistareykir Iceland (1)
- Þeistareykir Island (1)
Institut
- Institut für Geowissenschaften (130) (entfernen)
Volcanic tremor extraction and earthquake detection using music information retrieval algorithms
(2021)
Volcanic tremor signals are usually observed before or during volcanic eruptions and must be monitored to evaluate the volcanic activity. A challenge in studying seismic signals of volcanic origin is the coexistence of transient signal swarms and long-lasting volcanic tremor signals. Separating transient events from volcanic tremors can, therefore, contrib-ute to improving upon our understanding of the underlying physical processes. Exploiting the idea of harmonic-percussive separation in musical signal processing, we develop a method to extract the harmonic volcanic tremor signals and to detect tran-sient events from seismic recordings. Based on the similarity properties of spectrogram frames in the time-frequency domain, we decompose the signal into two separate spec-trograms representing repeating (harmonic) and nonrepeating (transient) patterns, which correspond to volcanic tremor signals and earthquake signals, respectively. We reconstruct the harmonic tremor signal in the time domain from the complex spectrogram of the repeating pattern by only considering the phase components for the frequency range in which the tremor amplitude spectrum is significantly contribut-ing to the energy of the signal. The reconstructed signal is, therefore, clean tremor signal without transient events. Furthermore, we derive a characteristic function suitable for the detection of tran-sient events (e.g., earthquakes) by integrating amplitudes of the nonrepeating spectro-gram over frequency at each time frame. Considering transient events like earthquakes, 78% of the events are detected for signal-to-noise ratio = 0.1 in our semisynthetic tests. In addition, we compared the number of detected earthquakes using our method for one month of continuous data recorded during the Holuhraun 2014-2015 eruption in Iceland with the bulletin presented in Agustsdottir et al. (2019). Our single station event detection algorithm identified 84% of the bulletin events. Moreover, we detected a total of 12,619 events, which is more than twice the number of the bulletin events.
Rheology describes the flow of matter under the influence of stress, and - related to solids- it investigates how solids subjected to stresses deform. As the deformation of the Earth’s outer layers, the lithosphere and the crust, is a major focus of rheological studies, rheology in the geosciences describes how strain evolves in rocks of variable composition and temperature under tectonic stresses. It is here where deformation processes shape the form of ocean basins and mountain belts that ultimately result from the complex interplay between lithospheric plate motion and the susceptibility of rocks to the influence of plate-tectonic forces. A rigorous study of the strength of the lithosphere and deformation phenomena thus requires in-depth studies of the rheological characteristics of the involved materials and the temporal framework of deformation processes.
This dissertation aims at analyzing the influence of the physical configuration of the lithosphere on the present-day thermal field and the overall rheological characteristics of the lithosphere to better understand variable expressions in the formation of passive continental margins and the behavior of strike-slip fault zones. The main methodological approach chosen is to estimate the present-day thermal field and the strength of the lithosphere by 3-D numerical modeling. The distribution of rock properties is provided by 3-D structural models, which are used as the basis for the thermal and rheological modeling. The structural models are based on geophysical and geological data integration, additionally constrained by 3-D density modeling. More specifically, to decipher the thermal and rheological characteristics of the lithosphere in both oceanic and continental domains, sedimentary basins in the Sea of Marmara (continental transform setting), the SW African passive margin (old oceanic crust), and the Norwegian passive margin (young oceanic crust) were selected for this study.
The Sea of Marmara, in northwestern Turkey, is located where the dextral North Anatolian Fault zone (NAFZ) accommodates the westward escape of the Anatolian Plate toward the Aegean. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the lateral crustal heterogeneities is presented for the first time in this study. Here, I use different gravity datasets and the general non-uniqueness in potential field modeling, to propose three possible end-member scenarios of crustal configuration. The models suggest that pronounced gravitational anomalies in the basin originate from significant density heterogeneities within the crust. The rheological modeling reveals that associated variations in lithospheric strength control the mechanical segmentation of the NAFZ. Importantly, a strong crust that is mechanically coupled to the upper mantle spatially correlates with aseismic patches where the fault bends and changes its strike in response to the presence of high-density lower crustal bodies. Between the bends, mechanically weaker crustal domains that are decoupled from the mantle are characterized by creep.
For the passive margins of SW Africa and Norway, two previously published 3-D conductive and lithospheric-scale thermal models were analyzed. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data, such as seismic observations and the gravity field. Here, the rheological modeling suggests that the present-day lithospheric strength across the oceanic domain is ultimately affected by the age and past thermal and tectonic processes as well as the depth of the thermal lithosphere-asthenosphere boundary, while the configuration of the crystalline crust dominantly controls the rheological behavior of the lithosphere beneath the continental domains of both passive margins.
The thermal and rheological models show that the variations of lithospheric strength are fundamentally influenced by the temperature distribution within the lithosphere. Moreover, as the composition of the lithosphere significantly influences the present-day thermal field, it therefore also affects the rheological characteristics of the lithosphere. Overall my studies add to our understanding of regional tectonic deformation processes and the long-term behavior of sedimentary basins; they confirm other analyses that have pointed out that crustal heterogeneities in the continents result in diverse lithospheric thermal characteristics, which in turn results in higher complexity and variations of rheological behavior compared to oceanic domains with a thinner, more homogeneous crust.
Transient permeability in porous and fractured sandstones mediated by fluid-rock interactions
(2021)
Understanding the fluid transport properties of subsurface rocks is essential for a large number of geotechnical applications, such as hydrocarbon (oil/gas) exploitation, geological storage (CO2/fluids), and geothermal reservoir utilization. To date, the hydromechanically-dependent fluid flow patterns in porous media and single macroscopic rock fractures have received numerous investigations and are relatively well understood. In contrast, fluid-rock interactions, which may permanently affect rock permeability by reshaping the structure and changing connectivity of pore throats or fracture apertures, need to be further elaborated. This is of significant importance for improving the knowledge of the long-term evolution of rock transport properties and evaluating a reservoir’ sustainability. The thesis focuses on geothermal energy utilization, e.g., seasonal heat storage in aquifers and enhanced geothermal systems, where single fluid flow in porous rocks and rock fracture networks under various pressure and temperature conditions dominates.
In this experimental study, outcrop samples (i.e., Flechtinger sandstone, an illite-bearing Lower Permian rock, and Fontainebleau sandstone, consisting of pure quartz) were used for flow-through experiments under simulated hydrothermal conditions. The themes of the thesis are (1) the investigation of clay particle migration in intact Flechtinger sandstone and the coincident permeability damage upon cyclic temperature and fluid salinity variations; (2) the determination of hydro-mechanical properties of self-propping fractures in Flechtinger and Fontainebleau sandstones with different fracture features and contrasting mechanical properties; and (3) the investigation of the time-dependent fracture aperture evolution of Fontainebleau sandstone induced by fluid-rock interactions (i.e., predominantly pressure solution). Overall, the thesis aims to unravel the mechanisms of the instantaneous reduction (i.e., direct responses to thermo-hydro-mechanical-chemical (THMC) conditions) and progressively-cumulative changes (i.e., time-dependence) of rock transport properties.
Permeability of intact Flechtinger sandstone samples was measured under each constant condition, where temperature (room temperature up to 145 °C) and fluid salinity (NaCl: 0 ~ 2 mol/l) were stepwise changed. Mercury intrusion porosimetry (MIP), electron microprobe analysis (EMPA), and scanning electron microscopy (SEM) were performed to investigate the changes of local porosity, microstructures, and clay element contents before and after the experiments. The results indicate that the permeability of illite-bearing Flechtinger sandstones will be impaired by heating and exposure to low salinity pore fluids. The chemically induced permeability variations prove to be path-dependent concerning the applied succession of fluid salinity changes. The permeability decay induced by a temperature increase and a fluid salinity reduction operates by relatively independent mechanisms, i.e., thermo-mechanical and thermo-chemical effects.
Further, the hydro-mechanical investigations of single macroscopic fractures (aligned, mismatched tensile fractures, and smooth saw-cut fractures) illustrate that a relative fracture wall offset could significantly increase fracture aperture and permeability, but the degree of increase depends on fracture surface roughness. X-ray computed tomography (CT) demonstrates that the contact area ratio after the pressure cycles is inversely correlated to the fracture offset. Moreover, rock mechanical properties, determining the strength of contact asperities, are crucial so that relatively harder rock (i.e., Fontainebleau sandstone) would have a higher self-propping potential for sustainable permeability during pressurization. This implies that self-propping rough fractures with a sufficient displacement are efficient pathways for fluid flow if the rock matrix is mechanically strong.
Finally, two long-term flow-through experiments with Fontainebleau sandstone samples containing single fractures were conducted with an intermittent flow (~140 days) and continuous flow (~120 days), respectively. Permeability and fluid element concentrations were measured throughout the experiments. Permeability reduction occurred at the beginning stage when the stress was applied, while it converged at later stages, even under stressed conditions. Fluid chemistry and microstructure observations demonstrate that pressure solution governs the long-term fracture aperture deformation, with remarkable effects of the pore fluid (Si) concentration and the structure of contact grain boundaries. The retardation and the cessation of rock fracture deformation are mainly induced by the contact stress decrease due to contact area enlargement and a dissolved mass accumulation within the contact boundaries. This work implies that fracture closure under constant (pressure/stress and temperature) conditions is likely a spontaneous process, especially at the beginning stage after pressurization when the contact area is relatively small. In contrast, a contact area growth yields changes of fracture closure behavior due to the evolution of contact boundaries and concurrent changes in their diffusive properties. Fracture aperture and thus permeability will likely be sustainable in the long term if no other processes (e.g., mineral precipitations in the open void space) occur.
Mixed sand- and gravel-bed rivers record erosion, transport, and fining signals in their bedload size distributions. Thus, grain-size data are imperative for studying these processes. However, collecting hundreds to thousands of pebble measurements in steep and dynamic high-mountain river settings remains challenging. Using the recently published digital grain-sizing algorithm PebbleCounts, we were able to survey seven large (>= 1,000 m2) channel cross-sections and measure thousands to tens-of-thousands of grains per survey along a 100-km stretch of the trunk stream of the Toro Basin in Northwest Argentina. The study region traverses a steep topographic and environmental gradient on the eastern margin of the Central Andean Plateau. Careful counting and validation allows us to identify measurement errors and constrain percentile uncertainties using large sample sizes. In the coarse >= 2.5 cm fraction of bedload, only the uppermost size percentiles (>= 95th) vary significantly downstream, whereas the 50th and 84th percentiles show less variability. We note a relation between increases in these upper percentiles and along-channel junctions with large, steep tributaries. This signal is strongly influenced by lithology and geologic structures, and mixed with local hillslope input. In steep catchments like the Toro Basin, we suggest nonlinear relationships between geomorphic metrics and grain size, whereby the steepest parts of the landscape exert primary control on the upper grain-size percentiles. Thus, average or median metrics that do not apply weights or thresholds to steeper topography may be less predictive of grain-size distributions in such settings.
The Upper Indus Basin (UIB), which covers a wide range of climatic and topographic settings, provides an ideal venue to explore the relationship between climate and topography. While the distribution of snow and glaciers is spatially and temporally heterogeneous, there exist regions with similar elevation-snow relationships. In this work, we construct elevation-binned snow-cover statistics to analyze 3415 watersheds and 7357 glaciers in the UIB region. We group both glaciers and watersheds using a hierarchical clustering approach and find that (1) watershed clusters mirror large-scale moisture transport patterns and (2) are highly dependent on median watershed elevation. (3) Glacier clusters are spatially heterogeneous and are less strongly controlled by elevation, but rather by local topographic parameters that modify solar insolation. Our clustering approach allows us to clearly define self-similar snow-topographic regions. Eastern watersheds in the UIB show a steep snow cover-elevation relationship whereas watersheds in the central and western UIB have moderately sloped relationships, but cluster in distinct groups. We highlight this snow-cover-topographic transition zone and argue that these watersheds have different hydrologic responses than other regions. Our hierarchical clustering approach provides a potential new framework to use in defining climatic zones in the cyrosphere based on empirical data.
A multidisciplinary approach to the study of collisional orogenic belts can improve our knowledge of their geodynamic evolution and may suggest new tectonic models, especially for (U)HP rocks inside the accretionary wedge. In the Western Alps, wherein nappes of different origin are stacked, having recorded different metamorphic peaks at different stages of the orogenic evolution. This study focuses on the External (EPZ) and Internal (IPZ) ophiolitic units of the Piedmont Zone (Susa Valley, Western Alps), which were deformed throughout four tectonometamorphic phases (D1 to D4), developing different foliations and cleavages (S1 to S4) at different metamorphic conditions. The IPZ and EPZ are separated by a shear zone (i.e. the Susa Shear Zone (SSZ)) during which a related mylonitic foliation (SM) developed. S1 developed at high pressure conditions (Epidote-eclogite vs. Lawsonite-blueschist facies conditions for IPZ and EPZ, respectively), as suggested by the composition of white mica (i.e. phengite), whereas S2 developed at low pressure conditions (Epidote-greenschist facies conditions in both IPZ and EPZ) and is defined by muscovite. White mica defining the SM mylonitic foliation (T1) is mostly defined by phengite, while the T2-related disjunctive cleavage is defined by fine-grained muscovite. The relative chronology inferred from meso-and micro-structural observations suggests that T1 was near-coeval with respect to the D2, while T2 developed during D4. A new set of radiometric ages of the main metamorphic foliations were obtained by in situ Ar/Ar dating on white mica. Different generations of white mica defining S1 and S2 foliations in both the IPZ and EPZ and SM in the SSZ, were dated and two main groups of ages were obtained. In both IPZ and EPZ, S1 foliation developed at-46-41 Ma, while S2 foliation developed at-40-36 Ma and was nearly coeval with the SM mylonitic foliation (-39-36 Ma). Comparison between structural, petrological and geochronological data allows to define time of coupling of the different units and consequently to infer new tectonic implications for the exhumation of meta-ophiolites of the Piedmont Zone within axial sector of the Western Alps.
The Central Andes region in South America is characterized by a complex and heterogeneous deformation system. Recorded seismic activity and mapped neotectonic structures indicate that most of the intraplate deformation is located along the margins of the orogen, in the transitions to the foreland and the forearc. Furthermore, the actively deforming provinces of the foreland exhibit distinct deformation styles that vary along strike, as well as characteristic distributions of seismicity with depth. The style of deformation transitions from thin-skinned in the north to thick-skinned in the south, and the thickness of the seismogenic layer increases to the south. Based on geological/geophysical observations and numerical modelling, the most commonly invoked causes for the observed heterogeneity are the variations in sediment thickness and composition, the presence of inherited structures, and changes in the dip of the subducting Nazca plate. However, there are still no comprehensive investigations on the relationship between the lithospheric composition of the Central Andes, its rheological state and the observed deformation processes. The central aim of this dissertation is therefore to explore the link between the nature of the lithosphere in the region and the location of active deformation. The study of the lithospheric composition by means of independent-data integration establishes a strong base to assess the thermal and rheological state of the Central Andes and its adjacent lowlands, which alternatively provide new foundations to understand the complex deformation of the region. In this line, the general workflow of the dissertation consists in the construction of a 3D data-derived and gravity-constrained density model of the Central Andean lithosphere, followed by the simulation of the steady-state conductive thermal field and the calculation of strength distribution. Additionally, the dynamic response of the orogen-foreland system to intraplate compression is evaluated by means of 3D geodynamic modelling.
The results of the modelling approach suggest that the inherited heterogeneous composition of the lithosphere controls the present-day thermal and rheological state of the Central Andes, which in turn influence the location and depth of active deformation processes. Most of the seismic activity and neo--tectonic structures are spatially correlated to regions of modelled high strength gradients, in the transition from the felsic, hot and weak orogenic lithosphere to the more mafic, cooler and stronger lithosphere beneath the forearc and the foreland. Moreover, the results of the dynamic simulation show a strong localization of deviatoric strain rate second invariants in the same region suggesting that shortening is accommodated at the transition zones between weak and strong domains. The vertical distribution of seismic activity appears to be influenced by the rheological state of the lithosphere as well. The depth at which the frequency distribution of hypocenters starts to decrease in the different morphotectonic units correlates with the position of the modelled brittle-ductile transitions; accordingly, a fraction of the seismic activity is located within the ductile part of the crust. An exhaustive analysis shows that practically all the seismicity in the region is restricted above the 600°C isotherm, in coincidence with the upper temperature limit for brittle behavior of olivine. Therefore, the occurrence of earthquakes below the modelled brittle-ductile could be explained by the presence of strong residual mafic rocks from past tectonic events. Another potential cause of deep earthquakes is the existence of inherited shear zones in which brittle behavior is favored through a decrease in the friction coefficient. This hypothesis is particularly suitable for the broken foreland provinces of the Santa Barbara System and the Pampean Ranges, where geological studies indicate successive reactivation of structures through time. Particularly in the Santa Barbara System, the results indicate that both mafic rocks and a reduction in friction are required to account for the observed deep seismic events.
White mica and tourmaline are the dominant hydrothermal alteration minerals at the world-class Panasqueira W-Sn-Cu deposit in Portugal. Thus, understanding the controls on their chemical composition helps to constrain ore formation processes at this deposit and determine their usefulness as pathfinder minerals for mineralization in general. We combine whole-rock geochemistry of altered and unaltered metasedimentary host rocks with in situ LA-ICP-MS measurements of tourmaline and white mica from the alteration halo. Principal component analysis (PCA) is used to better identify geochemical patterns and trends of hydrothermal alteration in the datasets. The hydrothermally altered metasediments are enriched in As, Sn, Cs, Li, W, F, Cu, Rb, Zn, Tl, and Pb relative to unaltered samples. In situ mineral analyses show that most of these elements preferentially partition into white mica over tourmaline (Li, Rb, Cs, Tl, W, and Sn), whereas Zn is enriched in tourmaline. White mica has distinct compositions in different settings within the deposit (greisen, vein selvages, wall rock alteration zone, late fault zone), indicating a compositional evolution with time. In contrast, tourmaline from different settings overlaps in composition, which is ascribed to a stronger dependence on host rock composition and also to the effects of chemical zoning and microinclusions affecting the LA-ICP-MS analyses. Hence, in this deposit, white mica is the better recorder of the fluid composition. The calculated trace-element contents of the Panasqueira mineralizing fluid based on the mica data and estimates of mica-fluid partition coefficients are in good agreement with previous fluid-inclusion analyses. A compilation of mica and tourmaline trace-element compositions from Panasqueira and other W-Sn deposits shows that white mica has good potential as a pathfinder mineral, with characteristically high Li, Cs, Rb, Sn, and W contents. The trace-element contents of hydrothermal tourmaline are more variable. Nevertheless, the compiled data suggest that high Sn and Li contents are distinctive for tourmaline from W-Sn deposits.
The Arctic is greatly impacted by climate change. The increase in air temperature drives the thawing of permafrost and an increase in coastal erosion and river discharge. This leads to a greater input of sediment and organic matter into coastal waters, which substantially impacts the ecosystems by reducing light transmission through the water column and altering the biogeochemistry, but also the subsistence economy of local people, and changes in climate because of the transformation of organic matter into greenhouse gases. Yet, the quantification of suspended sediment in Arctic coastal and nearshore waters remains unsatisfactory due to the absence of dedicated algorithms to resolve the high loads occurring in the close vicinity of the shoreline. In this study we present the Arctic Nearshore Turbidity Algorithm (ANTA), the first reflectance-turbidity relationship specifically targeted towards Arctic nearshore waters that is tuned with in-situ measurements from the nearshore waters of Herschel Island Qikiqtaruk in the western Canadian Arctic. A semi-empirical model was calibrated for several relevant sensors in ocean color remote sensing, including MODIS, Sentinel 3 (OLCI), Landsat 8 (OLI), and Sentinel 2 (MSI), as well as the older Landsat sensors TM and ETM+. The ANTA performed better with Landsat 8 than with Sentinel 2 and Sentinel 3. The application of the ANTA to Sentinel 2 imagery that matches in-situ turbidity samples taken in Adventfjorden, Svalbard, shows transferability to nearshore areas beyond Herschel Island Qikiqtaruk.
The Arctic is greatly impacted by climate change. The increase in air temperature drives the thawing of permafrost and an increase in coastal erosion and river discharge. This leads to a greater input of sediment and organic matter into coastal waters, which substantially impacts the ecosystems by reducing light transmission through the water column and altering the biogeochemistry, but also the subsistence economy of local people, and changes in climate because of the transformation of organic matter into greenhouse gases. Yet, the quantification of suspended sediment in Arctic coastal and nearshore waters remains unsatisfactory due to the absence of dedicated algorithms to resolve the high loads occurring in the close vicinity of the shoreline. In this study we present the Arctic Nearshore Turbidity Algorithm (ANTA), the first reflectance-turbidity relationship specifically targeted towards Arctic nearshore waters that is tuned with in-situ measurements from the nearshore waters of Herschel Island Qikiqtaruk in the western Canadian Arctic. A semi-empirical model was calibrated for several relevant sensors in ocean color remote sensing, including MODIS, Sentinel 3 (OLCI), Landsat 8 (OLI), and Sentinel 2 (MSI), as well as the older Landsat sensors TM and ETM+. The ANTA performed better with Landsat 8 than with Sentinel 2 and Sentinel 3. The application of the ANTA to Sentinel 2 imagery that matches in-situ turbidity samples taken in Adventfjorden, Svalbard, shows transferability to nearshore areas beyond Herschel Island Qikiqtaruk.
Holocene temperature proxy records are commonly used in quantitative synthesis and model-data comparisons. However, comparing correlations between time series from records collected in proximity to one another with the expected correlations based on climate model simulations indicates either regional or noisy climate signals in Holocene temperature proxy records. In this study, we evaluate the consistency of spatial correlations present in Holocene proxy records with those found in data from the Last Glacial Maximum (LGM). Specifically, we predict correlations expected in LGM proxy records if the only difference to Holocene correlations would be due to more time uncertainty and more climate variability in the LGM. We compare this simple prediction to the actual correlation structure in the LGM proxy records. We found that time series data of ice-core stable isotope records and planktonic foraminifera Mg/Ca ratios were consistent between the Holocene and LGM periods, while time series of Uk'37 proxy records were not as we found no correlation between nearby LGM records. Our results support the finding of highly regional or noisy marine proxy records in the compilation analysed here and suggest the need for further studies on the role of climate proxies and the processes of climate signal recording and preservation.
The present work gives a detailed analysis of the metamorphic and structural evolution of the back-arc portion of the Famatinian Orogen exposed in the southern Sierra de Aconquija (Cuesta de La Chilca segment) in the Sierras Pampeanas Orientales (Eastern Pampean Sierras). The Pampeanas Orientales include from north to south the Aconquija, Ambato and Ancasti mountains. They are mainly composed of middle to high grade metasedimentary units and magmatic rocks.
At the south end of the Sierra de Aconquija, along an east to west segment extending over nearly 10 km (Cuesta de La Chilca), large volumes of metasedimentary rocks crop out. The eastern metasediments were defined as members of the El Portezuelo Metamorphic-Igneous Complex (EPMIC) or Eastern block and the western ones relate to the Quebrada del Molle Metamorphic Complex (QMMC) or Western block. The two blocks are divided by the La Chilca Shear Zone, which is reactivated as the Rio Chanarito fault.
The EPMIC, forming the hanging wall, is composed of schists, gneisses and rare amphibolites, calc- silicate schists, marbles and migmatites. The rocks underwent multiple episodes of deformation and a late high strain-rate episode with gradually increasing mylonitization to the west. Metamorphism progrades from a M-1 phase to the peak M-3, characterized by the reactions: Qtz + Pl + Bt +/- Ms -> Grt + Bt(2) + Pl(2) +/- Sil +/- Kfs, Qtz + Bt + Sil -> Crd + Kfs and Qtz + Grt + Sil -> Crd. The M-3 assemblage is coeval with the dominant foliation related to a third deformational phase (D-3).
The QMMC, forming the foot wall, is made up of fine-grained banded quartz - biotite schists with quartz veins and quartz-feldspar-rich pegmatites. To the east, schists are also overprinted by mylonitization. The M-3 peak assemblage is quartz + biotite + plagioclase +/- garnet +/- sillimanite +/- muscovite +/- ilmenite +/- magnetite +/- apatite.
The studied segment suffered multiphase deformation and metamorphism. Some of these phases can be correlated between both blocks. D-1 is locally preserved in scarce outcrops in the EPMIC but is the dominant in the QMMC, where S-1 is nearly parallel to S-0. In the EPMIC, D-2 is represented by the S-2 foliation, related to the F-2 folding that overprints S-1, with dominant strike NNW - SSE and high angles dip to the E. D-3 in the EPMIC have F-3 folds with axis oblique to S-2; the S-3 foliation has striking NW - SE dipping steeply to the E or W and develops interference patterns. In the QMMC, S-2 (D-2) is a discontinuous cleavage oblique to S-1 and transposed by S-3 (D-3), subparallel to S-1. Such structures in the QMMC developed at subsolidus conditions and could be correlated to those of the EPMIC, which formed under higher P-T conditions. The penetrative deformation D-2 in the EPMIC occurred during a prograde path with syntectonic growth of garnet reaching P-T conditions of 640 degrees C and 0.54 GPa in the EPMIC. This stage was followed by a penetrative deformation D-3 with syn-kinematic growth of garnet, cordierite and plagioclase. Peak P-T conditions calculated for M-3 are 710 degrees C and 0.60 GPa, preserved in the western part of the EPMIC, west of the unnamed fault.
The schists from the QMMC suffered the early low grade M-1 metamorphism with minimum PT conditions of ca 400 degrees C and 0.35 GPa, comparable to the fine schists (M-1) outcropping to the east. The D-2 deformation is associated with the prograde M-2 metamorphism. The penetrative D-3 stage is related to a medium grade metamorphism M-3, with peak conditions at ca 590 degrees C and 0.55 GPa.
The superimposed stages of deformation and metamorphism reaching high P-T conditions followed by isothermal decompression, defining a clockwise orogenic P-T path. During the Lower Paleozoic, folds were superimposed and recrystallization as well as partial melting at peak conditions occurred. Similar characteristics were described from the basement from other Famatinian-dominated locations of the Sierra de Aconquija and other ranges of the Sierras Pampeanas Orientales.
Multi-component (MC) diffusion simulations enable a process based and more precise approach to calculate transport and sorption compared to the commonly used single-component (SC) models following Fick's law. The MC approach takes into account the interaction of chemical species in the porewater with the diffuse double layer (DDL) adhering clay mineral surfaces. We studied the shaly, sandy and carbonate-rich facies of the Opalinus Clay. High clay contents dominate diffusion and sorption of uranium. The MC simulations show shorter diffusion lengths than the SC models due to anion exclusion from the DDL. This hampers diffusion of the predominant species CaUO2(CO3)32-. On the one side, species concentrations and ionic strengths of the porewater and on the other side surface charge of the clay minerals control the composition and behaviour of the DDL. For some instances, it amplifies the diffusion of uranium. We developed a workflow to transfer computationally intensive MC simulations to SC models via calibrated effective diffusion and distribution coefficients. Simulations for one million years depict maximum uranium diffusion lengths between 10 m and 35 m. With respect to the minimum requirement of a thickness of 100 m, the Opalinus Clay seems to be a suitable host rock for nuclear waste repositories.
Near-surface supergene ores of the Merensky Reef in the Bushveld Complex, South Africa, contain economic grades of platinum-group elements, however, these are currently uneconomic due to low recovery rates. This is the first study that investigates the variation in platinum-group elements in pristine and supergene samples of the Merensky Reef from five drill cores from the eastern Bushveld. The samples from the Richmond and Twickenham farms show different degrees of weathering. The whole-rock platinum-group element distribution was studied by inductively coupled plasma-mass spectrometry and the platinum-group minerals were investigated by reflected-light microscopy, scanning electron microscopy, and electron microprobe analysis. <br /> In pristine ("fresh") Merensky Reef samples, platinum-group elements occur mainly as discrete platinum-group minerals, such as platinum-group element-sulfides (cooperite-braggite) and laurite as well as subordinate platinum-group elementbismuthotellurides and platinum-group element-arsenides, and also in solid solution in sulfides (especially Pd in pentlandite). During weathering, Pd and S were removed, resulting in a platinum-group mineral mineralogy in the supergene Merensky Reef that mainly consists of relict platinum-group minerals, Pt-Fe alloys, and Pt-oxides/hydroxides. Additional proportions of platinum-group elements are hosted by Fe-hydroxides and secondary hydrosilicates (e.g., serpentine group minerals and chlorite). <br /> In supergene ores, only low recovery rates (ca. 40%) are achieved due to the polymodal and complex platinum-group element distribution. To achieve higher recovery rates for the platinum-group elements, hydrometallurgical or pyrometallurgical processing of the bulk ore would be required, which is not economically viable with existing technology.
The concentration of sulfur that can be dissolved in a silicate liquid is of fundamental importance because it is closely associated with several major Earth-related processes. Considerable effort has been made to understand the interplay between the effects of silicate melt composition and its capac-ity to retain sulfur, but the dependence on pressure and temperature is mostly based on experiments performed at pressures and temperatures below 6 GPa and 2073 K. Here we present a study of the effects of pressure and temperature on sulfur content at sulfide saturation of a peridotitic liquid. We performed 14 multi-anvil experiments using a peridotitic starting composition, and we produced 25 new measurements at conditions ranging from 7 to 23 GPa and 2173 to 2623 K. We analyzed the recovered samples using both electron microprobe and laser ablation ICP-MS. We compiled our data together with previously published data that were obtained at lower P-T conditions and with various silicate melt compositions. We present a new model based on this combined data set that encompasses the entire range of upper mantle pressure-temperature conditions, along with the effect of a wide range of silicate melt compositions. Our findings are consistent with earlier work based on extrapolation from lower-pressure and lower-temperature experiments and show a decrease of sulfur content at sulfide saturation (SCSS) with increasing pressure and an increase of SCSS with increasing temperature. We have extrapolated our results to pressure-temperature conditions of the Earth's primitive magma ocean, and show that FeS will exsolve from the molten silicate and can effectively be extracted to the core by a process that has been termed the "Hadean Matte." We also discuss briefly the implications of our results for the lunar magma ocean.
Assessing volcanic hazard in regions of distributed volcanism is challenging because of the uncertain location of future vents. A statistical-mechanical strategy to forecast such locations was recently proposed: here, we further develop and test it with analog models. We stress a gelatin block laterally and with surface excavations, and observe air-filled crack trajectories. We use the observed surface arrivals to sample the distributions of parameters describing the stress state of the gelatin block, combining deterministic crack trajectory simulations with a Monte Carlo approach. While the individual stress parameters remain unconstrained, we effectively retrieve their ratio and successfully forecast the arrival points of subsequent cracks.
Spectral analysis is a technique of time-series analysis that decomposes signals into linear combinations of harmonic components. Rooted in the 19th century, spectral analysis gained popularity in palaeoclimatology since the early 1980s. This was partly due to the availability of long time series of past climates, but also the development of new, partly adapted methods and the increasing spread of affordable personal computers. This paper reviews the most important methods of spectral analysis for palaeoclimate time series and discusses the prerequisites for their application as well as advantages and disadvantages. The paper also offers an overview of suitable software, as well as computer code for using the methods on synthetic examples.
Cyanobacteria are important primary producers in temperate freshwater ecosystems. However, studies on the seasonal and spatial distribution of cyanobacteria in deep lakes based on high-throughput DNA sequencing are still rare. In this study, we combined monthly water sampling and monitoring in 2019, amplicon sequence variants analysis (ASVs; a proxy for different species) and quantitative PCR targeting overall cyanobacteria abundance to describe the seasonal and spatial dynamics of cyanobacteria in the deep hard-water oligo-mesotrophic Lake Tiefer See, NE Germany. We observed significant seasonal variation in the cyanobacterial community composition (p < 0.05) in the epi- and metalimnion layers, but not in the hypolimnion. In winter-when the water column is mixed-picocyanobacteria (Synechococcus and Cyanobium) were dominant. With the onset of stratification in late spring, we observed potential niche specialization and coexistence among the cyanobacteria taxa driven mainly by light and nutrient dynamics. Specifically, ASVs assigned to picocyanobacteria and the genus Planktothrix were the main contributors to the formation of deep chlorophyll maxima along a light gradient. While Synechococcus and different Cyanobium ASVs were abundant in the epilimnion up to the base of the euphotic zone from spring to fall, Planktothrix mainly occurred in the metalimnetic layer below the euphotic zone where also overall cyanobacteria abundance was highest in summer. Our data revealed two potentially psychrotolerant (cold-adapted) Cyanobium species that appear to cope well under conditions of lower hypolimnetic water temperature and light as well as increasing sediment-released phosphate in the deeper waters in summer. The potential cold-adapted Cyanobium species were also dominant throughout the water column in fall and winter. Furthermore, Snowella and Microcystis-related ASVs were abundant in the water column during the onset of fall turnover. Altogether, these findings suggest previously unascertained and considerable spatiotemporal changes in the community of cyanobacteria on the species level especially within the genus Cyanobium in deep hard-water temperate lakes.
The Kohat fold and thrust belt in Pakistan shows a significantly different structural style due to the structural evolution on the double décollement compared to the rest of the Subhimalaya. In order to better understand the spatio-temporal structural evolution of the Kohat fold and thrust belt, we combine balanced cross sections with apatite (U?Th-Sm)/He (AHe) and apatite fission track (AFT) dating. The AHe and AFT ages appear to be totally reset, allowing us to date exhumation above structural ramps. The results suggest that deformation began on the frontal Surghar thrust at-15 Ma, predating or coeval with the development of the Main Boundary thrust at-12 Ma. Deformation propagated southward from the Main Boundary thrust on double de?collements between 10 Ma and 2 Ma, resulting in a disharmonic structural style inside the Kohat fold and thrust belt. Thermal modeling of the thermochronologic data suggest that samples inside Kohat fold and thrust belt experienced cooling due to formation of the duplexes; this deformation facilitated tectonic thickening of the wedge and erosion of the Miocene to Pliocene foreland strata. The spatial distribution of AHe and AFT ages in combination with the structural forward model suggest that, in the Kohat fold and thrust belt, the wedge deformed in-sequence as a supercritical wedge (-15-12 Ma), then readjusted by out-sequence deformation (-12-0 Ma) within the Kohat fold and thrust belt into a sub-critical wedge.
Leaching zones within potash seams generally represent a significant risk to subsurface mining operations and the construction of technical caverns in salt rocks, but their temporal and spatial formation has been investigated only rudimentarily to date. To the knowledge of the authors, current reactive transport simulation implementations are not capable to address hydraulic-chemical interactions within potash salt. For this reason, a reactive transport model has been developed and complemented by an innovative approach to calculate the interchange of minerals and solution at the water-rock interface. Using this model, a scenario analysis was carried out based on a carnallite-bearing potash seam. The results show that the evolution of leaching zones depends on the mineral composition and dissolution rate of the original salt rock, and that the formation can be classified by the dimensionless parameters of Peclet (Pe) and Damkohler (Da). For Pe > 2 and Da > 1, a funnel-shaped leaching zone is formed, otherwise the dissolution front is planar. Additionally, Da > 1 results in the formation of a sylvinitic zone and a flow barrier. Most scenarios represent hybrid forms of these cases. The simulated shapes and mineralogies are confirmed by literature data and can be used to assess the hazard potential.
Centroid moment tensor inversion can provide insight into ongoing tectonic processes and active faults. In the Alpine mountains (central Europe), challenges result from low signal-to-noise ratios of earthquakes with small to moderate magnitudes and complex wave propagation effects through the heterogeneous crustal structure of the mountain belt. In this thesis, I make use of the temporary installation of the dense AlpArray seismic network (AASN) to establish a work flow to study seismic source processes and enhance the knowledge of the Alpine seismicity. The cumulative thesis comprises four publications on the topics of large seismic networks, seismic source processes in the Alps, their link to tectonics and stress field, and the inclusion of small magnitude earthquakes into studies of active faults.
Dealing with hundreds of stations of the dense AASN requires the automated assessment of data and metadata quality. I developed the open source toolbox AutoStatsQ to perform an automated data quality control. Its first application to the AlpArray seismic network has revealed significant errors of amplitude gains and sensor orientations. A second application of the orientation test to the Turkish KOERI network, based on Rayleigh wave polarization, further illustrated the potential in comparison to a P wave polarization method. Taking advantage of the gain and orientation results of the AASN, I tested different inversion settings and input data types to approach the specific challenges of centroid moment tensor (CMT) inversions in the Alps. A comparative study was carried out to define the best fitting procedures.
The application to 4 years of seismicity in the Alps (2016-2019) substantially enhanced the amount of moment tensor solutions in the region. We provide a list of moment tensors solutions down to magnitude Mw 3.1. Spatial patterns of typical focal mechanisms were analyzed in the seismotectonic context, by comparing them to long-term seismicity, historical earthquakes and observations of strain rates. Additionally, we use our MT solutions to investigate stress regimes and orientations along the Alpine chain. Finally, I addressed the challenge of including smaller magnitude events into the study of active faults and source processes. The open-source toolbox Clusty was developed for the clustering of earthquakes based on waveforms recorded across a network of seismic stations. The similarity of waveforms reflects both, the location and the similarity of source mechanisms. Therefore the clustering bears the opportunity to identify earthquakes of similar faulting styles, even when centroid moment tensor inversion is not possible due to low signal-to-noise ratios of surface waves or oversimplified velocity models. The toolbox is described through an application to the Zakynthos 2018 aftershock sequence and I subsequently discuss its potential application to weak earthquakes (Mw<3.1) in the Alps.
Major challenges during geothermal exploration and exploitation include the structural-geological characterization of the geothermal system and the application of sustainable monitoring concepts to explain changes in a geothermal reservoir during production and/or reinjection of fluids. In the absence of sufficiently permeable reservoir rocks, faults and fracture networks are preferred drilling targets because they can facilitate the migration of hot and/or cold fluids. In volcanic-geothermal systems considerable amounts of gas emissions can be released at the earth surface, often related to these fluid-releasing structures.
In this thesis, I developed and evaluated different methodological approaches and measurement concepts to determine the spatial and temporal variation of several soil gas parameters to understand the structural control on fluid flow. In order to validate their potential as innovative geothermal exploration and monitoring tools, these methodological approaches were applied to three different volcanic-geothermal systems. At each site an individual survey design was developed regarding the site-specific questions.
The first study presents results of the combined measurement of CO2 flux, ground temperatures, and the analysis of isotope ratios (δ13CCO2, 3He/4He) across the main production area of the Los Humeros geothermal field, to identify locations with a connection to its supercritical (T > 374◦C and P > 221 bar) geothermal reservoir. The results of the systematic and large-scale (25 x 200 m) CO2 flux scouting survey proved to be a fast and flexible way to identify areas of anomalous degassing. Subsequent sampling with high resolution surveys revealed the actual extent and heterogenous pattern of anomalous degassing areas. They have been related to the internal fault hydraulic architecture and allowed to assess favourable structural settings for fluid flow such as fault intersections. Finally, areas of unknown structurally controlled permeability with a connection to the superhot geothermal reservoir have been determined, which represent promising targets for future geothermal exploration and development.
In the second study, I introduce a novel monitoring approach by examining the variation of CO2 flux to monitor changes in the reservoir induced by fluid reinjection. For that reason, an automated, multi-chamber CO2 flux system was deployed across the damage zone of a major normal fault crossing the Los Humeros geothermal field. Based on the results of the CO2 flux scouting survey, a suitable site was selected that had a connection to the geothermal reservoir, as identified by hydrothermal CO2 degassing and hot ground temperatures (> 50 °C). The results revealed a response of gas emissions to changes in reinjection rates within 24 h, proving an active hydraulic communication between the geothermal reservoir and the earth surface. This is a promising monitoring strategy that provides nearly real-time and in-situ data about changes in the reservoir and allows to timely react to unwanted changes (e.g., pressure decline, seismicity).
The third study presents results from the Aluto geothermal field in Ethiopia where an area-wide and multi-parameter analysis, consisting of measurements of CO2 flux, 222Rn, and 220Rn activity concentrations and ground temperatures was conducted to detect hidden permeable structures. 222Rn and 220Rn activity concentrations are evaluated as a complementary soil gas parameter to CO2 flux, to investigate their potential to understand tectono-volcanic degassing. The combined measurement of all parameters enabled to develop soil gas fingerprints, a novel visualization approach. Depending on the magnitude of gas emissions and their migration velocities the study area was divided in volcanic (heat), tectonic (structures), and volcano-tectonic dominated areas. Based on these concepts, volcano-tectonic dominated areas, where hot hydrothermal fluids migrate along permeable faults, present the most promising targets for future geothermal exploration and development in this geothermal field. Two of these areas have been identified in the south and south-east which have not yet been targeted for geothermal exploitation. Furthermore, two unknown areas of structural related permeability could be identified by 222Rn and 220Rn activity concentrations.
Eventually, the fourth study presents a novel measurement approach to detect structural controlled CO2 degassing, in Ngapouri geothermal area, New Zealand. For the first time, the tunable diode laser (TDL) method was applied in a low-degassing geothermal area, to evaluate its potential as a geothermal exploration method. Although the sampling approach is based on profile measurements, which leads to low spatial resolution, the results showed a link between known/inferred faults and increased CO2 concentrations. Thus, the TDL method proved to be a successful in the determination of structural related permeability, also in areas where no obvious geothermal activity is present. Once an area of anomalous CO2 concentrations has been identified, it can be easily complemented by CO2 flux grid measurements to determine the extent and orientation of the degassing segment.
With the results of this work, I was able to demonstrate the applicability of systematic and area-wide soil gas measurements for geothermal exploration and monitoring purposes. In particular, the combination of different soil gases using different measurement networks enables the identification and characterization of fluid-bearing structures and has not yet been used and/or tested as standard practice. The different studies present efficient and cost-effective workflows and demonstrate a hands-on approach to a successful and sustainable exploration and monitoring of geothermal resources. This minimizes the resource risk during geothermal project development. Finally, to advance the understanding of the complex structure and dynamics of geothermal systems, a combination of comprehensive and cutting-edge geological, geochemical, and geophysical exploration methods is essential.
Large rock slope failures play a pivotal role in long-term landscape evolution and are a major concern in land use planning and hazard aspects. While the failure phase and the time immediately prior to failure are increasingly well studied, the nature of the preparation phase remains enigmatic. This knowledge gap is due, to a large degree, to difficulties associated with instrumenting high mountain terrain and the local nature of classic monitoring methods, which does not allow integral observation of large rock volumes. Here, we analyse data from a small network of up to seven seismic sensors installed during July-October 2018 (with 43 days of data loss) at the summit of the Hochvogel, a 2592 m high Alpine peak. We develop proxy time series indicative of cyclic and progressive changes of the summit. Modal analysis, horizontal-to-vertical spectral ratio data and end-member modelling analysis reveal diurnal cycles of increasing and decreasing coupling stiffness of a 260,000 m(3) large, instable rock volume, due to thermal forcing. Relative seismic wave velocity changes also indicate diurnal accumulation and release of stress within the rock mass. At longer time scales, there is a systematic superimposed pattern of stress increased over multiple days and episodic stress release within a few days, expressed in an increased emission of short seismic pulses indicative of rock cracking. Our data provide essential first order information on the development of large-scale slope instabilities towards catastrophic failure. (c) 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd
Large rock slope failures play a pivotal role in long-term landscape evolution and are a major concern in land use planning and hazard aspects. While the failure phase and the time immediately prior to failure are increasingly well studied, the nature of the preparation phase remains enigmatic. This knowledge gap is due, to a large degree, to difficulties associated with instrumenting high mountain terrain and the local nature of classic monitoring methods, which does not allow integral observation of large rock volumes. Here, we analyse data from a small network of up to seven seismic sensors installed during July-October 2018 (with 43 days of data loss) at the summit of the Hochvogel, a 2592 m high Alpine peak. We develop proxy time series indicative of cyclic and progressive changes of the summit. Modal analysis, horizontal-to-vertical spectral ratio data and end-member modelling analysis reveal diurnal cycles of increasing and decreasing coupling stiffness of a 260,000 m(3) large, instable rock volume, due to thermal forcing. Relative seismic wave velocity changes also indicate diurnal accumulation and release of stress within the rock mass. At longer time scales, there is a systematic superimposed pattern of stress increased over multiple days and episodic stress release within a few days, expressed in an increased emission of short seismic pulses indicative of rock cracking. Our data provide essential first order information on the development of large-scale slope instabilities towards catastrophic failure. (c) 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd
The imagination of clearly separated core-shell structures is already outdated by the fact, that the nanoparticle core-shell structures remain in terms of efficiency behind their respective bulk material due to intermixing between core and shell dopant ions. In order to optimize the photoluminescence of core-shell UCNP the intermixing should be as small as possible and therefore, key parameters of this process need to be identified. In the present work the Ln(III) ion migration in the host lattices NaYF4 and NaGdF4 was monitored. These investigations have been performed by laser spectroscopy with help of lanthanide resonance energy transfer (LRET) between Eu(III) as donor and Pr(III) or Nd(III) as acceptor. The LRET is evaluated based on the Forster theory. The findings corroborate the literature and point out the migration of ions in the host lattices. Based on the introduced LRET model, the acceptor concentration in the surrounding of one donor depends clearly on the design of the applied core-shell-shell nanoparticles. In general, thinner intermediate insulating shells lead to higher acceptor concentration, stronger quenching of the Eu(III) donor and subsequently stronger sensitization of the Pr(III) or the Nd(III) acceptors. The choice of the host lattice as well as of the synthesis temperature are parameters to be considered for the intermixing process.
Models for the predictions of monetary losses from floods mainly blend data deemed to represent a single flood type and region. Moreover, these approaches largely ignore indicators of preparedness and how predictors may vary between regions and events, challenging the transferability of flood loss models. We use a flood loss database of 1812 German flood-affected households to explore how Bayesian multilevel models can estimate normalised flood damage stratified by event, region, or flood process type. Multilevel models acknowledge natural groups in the data and allow each group to learn from others. We obtain posterior estimates that differ between flood types, with credibly varying influences of water depth, contamination, duration, implementation of property-level precautionary measures, insurance, and previous flood experience; these influences overlap across most events or regions, however. We infer that the underlying damaging processes of distinct flood types deserve further attention. Each reported flood loss and affected region involved mixed flood types, likely explaining the uncertainty in the coefficients. Our results emphasise the need to consider flood types as an important step towards applying flood loss models elsewhere. We argue that failing to do so may unduly generalise the model and systematically bias loss estimations from empirical data.
Models for the predictions of monetary losses from floods mainly blend data deemed to represent a single flood type and region. Moreover, these approaches largely ignore indicators of preparedness and how predictors may vary between regions and events, challenging the transferability of flood loss models. We use a flood loss database of 1812 German flood-affected households to explore how Bayesian multilevel models can estimate normalised flood damage stratified by event, region, or flood process type. Multilevel models acknowledge natural groups in the data and allow each group to learn from others. We obtain posterior estimates that differ between flood types, with credibly varying influences of water depth, contamination, duration, implementation of property-level precautionary measures, insurance, and previous flood experience; these influences overlap across most events or regions, however. We infer that the underlying damaging processes of distinct flood types deserve further attention. Each reported flood loss and affected region involved mixed flood types, likely explaining the uncertainty in the coefficients. Our results emphasise the need to consider flood types as an important step towards applying flood loss models elsewhere. We argue that failing to do so may unduly generalise the model and systematically bias loss estimations from empirical data.
Permafrost is warming globally which leads to widespread permafrost thaw. Particularly ice-rich permafrost is vulnerable to rapid thaw and erosion, impacting whole landscapes and ecosystems. Retrogressive thaw slumps (RTS) are abrupt permafrost disturbances that expand by several meters each year and lead to an increased soil organic carbon release. Local Remote Sensing studies identified increasing RTS activity in the last two decades by increasing number of RTS or heightened RTS growth rates. However, a large-scale assessment across diverse permafrost regions and at high temporal resolution allowing to further determine RTS thaw dynamics and its main drivers is still lacking. In this study we apply the disturbance detection algorithm LandTrendr for automated large-scale RTS mapping and high temporal thaw dynamic assessment to North Siberia (8.1 x 106km2). We adapted and parametrised the temporal segmentation algorithm for abrupt disturbance detection to incorporate Landsat+Sentinel-2 mosaics, conducted spectral filtering, spatial masking and filtering, and a binary machine-learning object classification of the disturbance output to separate between RTS and false positives (F1 score: 0.609). Ground truth data for calibration and validation of the workflow was collected from 9 known RTS cluster sites using very highresolution RapidEye and PlanetScope imagery. Our study presents the first automated detection and assessment of RTS and their temporal dynamics at largescale for 2001-2019. We identified 50,895 RTS and a steady increase in RTS-affected area from 2001 to 2019 across North Siberia, with a more abrupt increase from 2016 onward. Overall the RTS-affected area increased by 331% compared to 2000 (2000: 20,158 ha, 2001-2019: 66,699 ha). Contrary to this, 5 focus sites show spatiotemporal variability in their annual RTS dynamics, with alternating periods of increased and decreased RTS development, indicating a close relationship to thaw drivers. The majority of identified RTS was active from 2000 onward and only a small proportion initiated during the assessment period, indicating that the increase in RTS-affected area was mainly caused by enlarging existing RTS and not by new RTS. The detected increase in RTS dynamics suggests advancing permafrost thaw and underlines the importance of assessing abrupt permafrost disturbances with high spatial and temporal resolution at large-scales. Obtaining such consistent disturbance products will help to parametrise regional and global climate change models.
In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D V-P and V-P/V-S models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1-4.2 M-L were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500m in their epicenter and similar to 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0-20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie-Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).
Precambrian meta-pelites of the Poshtuk area in northwest Iran contain the prograde mineral assemblage staurolite-garnet-chloritoid-muscovite-biotite that was replaced by the assemblage garnet-staurolite-chlorite-muscovite-biotite at peak metamorphic condition.
Whole-rock compositions reveal that high Fe, Al and Mn contents of their protolith rendered them prone to form these assemblages. Pseudosections calculated in KFMASH, MnKFMASH, and MnNCKFMASHO systems were used to investigate the P-T evolution of the samples. They clearly show the significant effect of MnO on the stability of the chloritoid-bearing assemblages and the formation of garnet through consumption of chlorite and chloritoid. The pseudosection in a T- aH(2)O diagram shows that the studied assemblage could be stable only at a(H2O)>0.8. X-Mg isopleths for garnet and biotite point to peak P-T conditions of about 3.75 kbar and 575 degrees C. Chloritoid stability is overstepped with such conditions. This can be attributed to thermal perturbation due to plutonism. It is concluded, metamorphism was primarily controlled by advective heat from magmatic intrusions in the Poshtuk area. The Precambrian basement complexes were extensively overprinted by the Pan-African Orogeny as well as younger magmatic and metamorphic activities associated to Alpine Orogeny during convergence of Arabian and Eurasian plate.
In recent decades, slope instability in high-mountain regions has often been linked to increase in temperature and the associated permafrost degradation and/or the increase in frequency/intensity of rainstorm events. In this context we analyzed the spatiotemporal evolution and potential controlling mechanisms of small- to medium-sized mass movements in a high-elevation catchment of the Italian Alps (Sulden/Solda basin). We found that slope-failure events (mostly in the form of rockfalls) have increased since the 2000s, whereas the occurrence of debris flows has increased only since 2010. The current climate-warming trend registered in the study area apparently increases the elevation of rockfall-detachment areas by approximately 300 m, mostly controlled by the combined effects of frost-cracking and permafrost thawing. In contrast, the occurrence of debris flows does not exhibit such an altitudinal shift, as it is primarily driven by extreme precipitation events exceeding the 75th percentile of the intensity-duration rainfall distribution. Potential debris-flow events in this environment may additionally be influenced by the accumulation of unconsolidated debris over time, which is then released during extreme rainfall events. Overall, there is evidence that the upper Sulden/Solda basin (above ca. 2500 m above sea level [a.s.l.]), and especially the areas in the proximity of glaciers, have experienced a significant decrease in slope stability since the 2000s, and that an increase in rockfalls and debris flows during spring and summer can be inferred. Our study thus confirms that "forward-looking" hazard mapping should be undertaken in these increasingly frequented, high-elevation areas of the Alps, as environmental change has elevated the overall hazard level in these regions.
Climatic change alters the frequency and intensity of natural hazards. In order to assess potential future changes in flood seasonality in the Rhine River Basin, we analyse changes in streamflow, snowmelt, precipitation, and evapotranspiration at 1.5, 2.0 and 3.0 ◦C global warming levels. The mesoscale Hydrological Model (mHM) forced with an ensemble of climate projection scenarios (five general circulation models under three representative concentration pathways) is used to simulate the present and future climate conditions of both, pluvial and nival hydrological regimes. Our results indicate that the interplay between changes in snowmelt- and rainfall-driven runoff is crucial to understand changes in streamflow maxima in the Rhine River. Climate projections suggest that future changes in flood characteristics in the entire Rhine River are controlled by both, more intense precipitation events and diminishing snow packs. The nature of this interplay defines the type of change in runoff peaks. On the sub-basin level (the Moselle River), more intense rainfall during winter is mostly counterbalanced by reduced snowmelt contribution to the streamflow. In the High Rhine (gauge at Basel), the strongest increases in streamflow maxima show up during winter, when strong increases in liquid precipitation intensity encounter almost unchanged snowmelt-driven runoff. The analysis of snowmelt events suggests that at no point in time during the snowmelt season, a warming climate results in an increase in the risk of snowmelt-driven flooding. We do not find indications of a transient merging of pluvial and nival floods due to climate warming.
Mehrphasige DNAPL-Pools zählen zu den häufigsten Ursachen für Grundwasserkontaminationen und sind bekannt für ihre Langlebigkeit. Obwohl Untersuchungen bereits gezeigt haben, dass die Phasen sich in ihrer Wasserlöslichkeit gegenseitig beeinflussen, werden diese Interaktionen von bisherigen Modellen vernachlässigt. Aus diesem Grund wurde ein semi-analytisches Berechnungsmodell entwickelt, welches die Poolzusammensetzung als zeitlich variabel behandelt. Basierend auf dem Raoult’schen Gesetz werden für jede Komponente Molanteil, effektive Wasserlöslichkeit und schließlich der Schadstoffaustrag infolge Advektion, Dispersion und Diffusion bestimmt. Die Ergebnisse unterscheiden sich deutlich von Studien an einphasigen Pools. So wird gezeigt, dass Schadstofffrachten über die Zeit sowohl zu- als auch abnehmen können und dass ohne Berücksichtigung des Raoult’schen Gesetzes sowohl die Langlebigkeit von DNAPL-Pools als auch die Dauer bis zur Unterschreitung von Grenzwerten teils deutlich unterschätzt wird. Eine Sensitivitätsanalyse zeigt zudem, dass schwer lösliche Nebenbestandteile nicht vernachlässigt werden dürfen, leicht lösliche hingegen schon.
This study is trying to understand the pre-eruptive magma storage and crystallization conditions of the Middle Miocene aged, silica-saturated trachytic rocks of the Afyon Volcanic Complex (AVC) in Western Anatolia, Turkey. Those rocks can be divided by their high K2O, K2O/Na2O ratio and Mg# into two groups, namely the intermediate-potassic (IPG) and the ultrapotassic (UPG). Here we are comparing calculated pressure (P) - temperature (T) conditions derived from geothermobarometric calculations of natural samples with results of high-pressure, high-temperature phase equilibria experiments. IPG samples are richer in silica (57-64 wt% SiO2), whereas UPG samples show intermediate SiO2 contents of 56-58 wt%. UPG are having high K2O contents ((>)9 wt %), K2O/Na2O ratios ((>)10 wt%) and Mg# values (75-77). IPG phenocrysts comprise plagioclase + biotite + amphibole + clinopyroxene +/- orthopyroxene +/- sanidine +/- phlogopite and oxides, while UPG mineralogical assemblage consists of amphibole + phlogopite + clinopyroxene + olivine + sanidine and oxides. IPG and UPG are enriched in Large-Ion Lithophile Elements (LILE), and both have negative anomalies in Nb, Sr, Zr and Ti elements. Additionally, IPG shows positive anomalies in Pb. Both IPG and UPG display enrichment in Light Rare Earth Elements (LREE), while IPG shows a more significant negative anomaly in Eu when compared to UPG. Plagioclase fractionation may play a role in magma generation. In IPG samples Ni and Cr values range between (3.3-18.8 ppm) and (2.6-27.8 ppm), respectively; whereas UPG samples have (119.1-120.7 ppm) Ni and (212.1-219.9 ppm) Cr. Dy/Yb ratios of IPG and UPG are higher than 2 and may indicate that garnet was present in the source. Geothermobarometric calculations for natural IPG clinopyroxene-melt pairs imply higher PT conditions (Dogan-Kulahci et al., 2015), while in this study high-pressure/high-temperature (HP/HT) phase equilibria experiments recreated the natural mineral assemblage at 2-4 kbar, 6-9 km and c. 900 degrees C. New plagioclase-melt calculations have confirmed lower mean magma storage temperatures, which are closer to the experimental results but still slightly elevated. Thus, trace element results of the natural rocks and experimental data may imply that a deep garnet-bearing magma source mixed with shallower magmas (IPG) was feeding the volcanic eruption.
Transport properties of potential host rocks for nuclear waste disposal are typically determined in laboratory or in-situ experiments under geochemically controlled and constant conditions. Such a homogeneous assumption is no longer applicable on the host rock scale as can be seen from the pore water profiles of the potential host rock Opalinus Clay at Mont Terri (Switzerland). The embedding aquifers are the hydro-geological boundaries, that established gradients in the 210 m thick low permeable section through diffusive exchange over millions of years. Present-day pore water profiles were confirmed by a data-driven as well as by a conceptual scenario. Based on the modelled profiles, the influence of the geochemical gradient on uranium migration was quantified by comparing the distances after one million years with results of common homogeneous models. Considering the heterogeneous system, uranium migrated up to 24 m farther through the formation depending on the source term position within the gradient and on the partial pressure of carbon dioxide pCO2 of the system. Migration lengths were almost equal for single- and multicomponent diffusion. Differences can predominantly be attributed to changes in the sorption capacity, whereby pCO2 governs how strong uranium migration is affected by the geochemical gradient. Thus, the governing parameters for uranium migration in the Opalinus Clay can be ordered in descending priority: pCO2, geochemical gradients, mineralogical heterogeneity.</p>
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.
POET (v0.1): speedup of many-core parallel reactive transport simulations with fast DHT lookups
(2021)
Coupled reactive transport simulations are extremely demanding in terms of required computational power, which hampers their application and leads to coarsened and oversimplified domains. The chemical sub-process represents the major bottleneck: its acceleration is an urgent challenge which gathers increasing interdisciplinary interest along with pressing requirements for subsurface utilization such as spent nuclear fuel storage, geothermal energy and CO2 storage. In this context we developed POET (POtsdam rEactive Transport), a research parallel reactive transport simulator integrating algorithmic improvements which decisively speed up coupled simulations. In particular, POET is designed with a master/worker architecture, which ensures computational efficiency in both multicore and cluster compute environments. POET does not rely on contiguous grid partitions for the parallelization of chemistry but forms work packages composed of grid cells distant from each other. Such scattering prevents particularly expensive geochemical simulations, usually concentrated in the vicinity of a reactive front, from generating load imbalance between the available CPUs (central processing units), as is often the case with classical partitions. Furthermore, POET leverages an original implementation of the distributed hash table (DHT) mechanism to cache the results of geochemical simulations for further reuse in subsequent time steps during the coupled simulation. The caching is hence particularly advantageous for initially chemically homogeneous simulations and for smooth reaction fronts. We tune the rounding employed in the DHT on a 2D benchmark to validate the caching approach, and we evaluate the performance gain of POET's master/worker architecture and the DHT speedup on a 3D benchmark comprising around 650 000 grid elements. The runtime for 200 coupling iterations, corresponding to 960 simulation days, reduced from about 24 h on 11 workers to 29 min on 719 workers. Activating the DHT reduces the runtime further to 2 h and 8 min respectively. Only with these kinds of reduced hardware requirements and computational costs is it possible to realistically perform the longterm complex reactive transport simulations, as well as perform the uncertainty analyses required by pressing societal challenges connected with subsurface utilization.
Rotational motions play a key role in measuring seismic wavefield properties. Using newly developed portable rotational instruments, it is now possible to directly measure rotational motions in a broad frequency range. Here, we investigated the instrumental self-noise and data quality in a huddle test in Fürstenfeldbruck, Germany, in August 2019. We compare the data from six rotational and three translational sensors. We studied the recorded signals using correlation, coherence analysis, and probabilistic power spectral densities. We sorted the coherent noise into five groups with respect to the similarities in frequency content and shape of the signals. These coherent noises were most likely caused by electrical devices, the dehumidifier system in the building, humans, and natural sources such as wind. We calculated self-noise levels through probabilistic power spectral densities and by applying the Sleeman method, a three-sensor method. Our results from both methods indicate that self-noise levels are stable between 0.5 and 40 Hz. Furthermore, we recorded the 29 August 2019 ML 3.4 Dettingen earthquake. The calculated source directions are found to be realistic for all sensors in comparison to the real back azimuth. We conclude that the five tested blueSeis-3A rotational sensors, when compared with respect to coherent noise, self-noise, and source direction, provide reliable and consistent results. Hence, field experiments with single rotational sensors can be undertaken.
Rotational motions play a key role in measuring seismic wavefield properties. Using newly developed portable rotational instruments, it is now possible to directly measure rotational motions in a broad frequency range. Here, we investigated the instrumental self-noise and data quality in a huddle test in Fürstenfeldbruck, Germany, in August 2019. We compare the data from six rotational and three translational sensors. We studied the recorded signals using correlation, coherence analysis, and probabilistic power spectral densities. We sorted the coherent noise into five groups with respect to the similarities in frequency content and shape of the signals. These coherent noises were most likely caused by electrical devices, the dehumidifier system in the building, humans, and natural sources such as wind. We calculated self-noise levels through probabilistic power spectral densities and by applying the Sleeman method, a three-sensor method. Our results from both methods indicate that self-noise levels are stable between 0.5 and 40 Hz. Furthermore, we recorded the 29 August 2019 ML 3.4 Dettingen earthquake. The calculated source directions are found to be realistic for all sensors in comparison to the real back azimuth. We conclude that the five tested blueSeis-3A rotational sensors, when compared with respect to coherent noise, self-noise, and source direction, provide reliable and consistent results. Hence, field experiments with single rotational sensors can be undertaken.
En el presente trabajo se realizó una investigación multidisciplinaria combinando métodos de geomorfología tectónica con estudios geofisicos y estructurales, focalizados principalmente en la caracterización neotectónica de ambos faldeos de la sierra de La Candelaria y del extremo sur de la cuenca de Metán. La zona de estudio se encuentra ubicada en la región limítrofe entre las provincias de Salta y Tucumán y pertenece a la provincia geológica del Sistema Santa Bárbara.
El principal objetivo consistió en contextualizar las evidencias de actividad tectónica cuaternaria de la región mediante la propuesta de un modelo estructural novedoso, con el propósito de incrementar la información disponible sobre estructuras neotectónicas y su potencial sismogénico. Con este fin, se aplicaron e integraron diversas técnicas tales como la interpretación de líneas sísmicas de reflexión, construcción de secciones estructurales balanceadas, y métodos geofísicos someros, para constatar el comportamiento en profundidad tanto de las estructuras geológicas identificadas en superficie como de las posibles fallas ciegas corticales involucradas.
En primer lugar, se realizó un relevamiento regional del área de estudio empleando imágenes satelitales multiespectrales LANDSAT y SENTINEL 2, que permitieron reconocer diferentes niveles de abanicos aluviales y terrazas fluviales cuaternarios. Mediante la determinación de diferentes indicadores morfométricos en modelos de elevación digital (MED), junto con observaciones de campo, fue posible identificar evidencias de deformación sobre dichos niveles cuaternarios que han sido relacionadas genéticamente con cuatro fallas neotectónicas. Tres de ellas (fallas Arias, El Quemado y Copo Quile) fueron seleccionadas para efectuar estudios de mayor detalle por medio de la aplicación de métodos de geofísica somera (tomografía eléctrica resistiva (ERT) y tomografía sísmica de refracción Sísmica (SRT)), que permitieron corroborar su existencia en profundidad, realizar inferencias geométricas y cinemáticas, y estimar la magnitud de la deformación reciente. Las fallas Arias y El Quemado fueron interpretadas como fallas inversas relacionadas con deslizamiento flexural interstratal, mientras que la falla Copo Quile se interpretó como una falla inversa ciega de bajo ángulo.También se realizó una interpretación conjunta de líneas sísmicas de reflexión y pozos exploratorios pertenecientes a áreas hidrocarburíferas de las cuencas de Choromoro y Metán con el fin de contextualizar las principales estructuras reconocidas en el marco estratigráfico y tectónico regional. Toda la información fue integrada en una sección estructural balanceada mediante técnicas de modelado cinemático. Dicho modelo permite inferir que la deformación cuaternaria reconocida está relacionada al desplazamiento del basamento a lo largo de un corrimiento ciego, responsable del levantamiento de la sierra de La Candelaria y el cerr Cantero. Asimismo, el modelo cinemático permite interpretar la ubicación aproximada de los principales niveles de despegue que controlan el estilo de deformación. El nivel de despegue más somero, que controla la deformación de la cobertura sedimentaria se encuentra a 4 km de profundidad, a 21 km se estima la presencia de otra zona de cizalla subhorizontal dentro del basamento.
Finalmente, a partir de la integración de todos los resultados obtenidos, se evaluó el potencial sismogénico de las fallas en la zona de estudio. Las fallas de primer orden que controlan la deformación en la zona son las responsables de los grandes terremotos. Mientras, las fallas Cuaternarias flexodeslizantes e inversas afectan solamente a la cobertura sedimentaria y serían estructuras de segundo orden que acomodan la deformación y fueron activadas durante el cuaternario con movimientos asísmicos y/o sísmicos de muy baja magnitud.
Estos resultados permiten inferir que el corrimiento La Candelaria constituye una fuente sismogénica potencial de importancia para la región, donde se ubican numerosas poblaciones y obras civiles de envergadura. Por otra parte, la sección estructural balanceada implica la presencia de otras fallas ciegas de distinto orden de magnitud que podrían ser posibles fuentes sismogénicas profundas adicionales, marcando la necesidad de continuar con el desarrollo de este tipo de estudios en esta región tectónicamente activa.
El flanco oriental de los Andes Centrales en el noroeste argentino es una zona caracterizada por serranías limitadas por fallas inversas que conforman un orógeno de piel gruesa activo con un patrón espacio-temporal no sistemático de deformación contraccional. Este patrón queda representado tanto por la dispersión de la actividad sísmica cortical como de la localización de las estructuras cuaternarias a través de la Cordillera Oriental y el Sistema de Santa Bárbara, configurando un frente orogénico difuso de más de 200 km de extensión. El estudio de la actividad neotectónica en esta región ha tomado más relevancia en los últimos años, mediante la aplicación de herramientas variadas, incluyendo técnicas de geomorfología tectónica, herramientas de teledetección, geodesia y estudios de campo convencionales. Los depósitos lacustres han demostrado ser, en numerosos ejemplos, excelentes marcadores de la actividad tectónica, dadas la horizontalidad original de sus capas y la susceptibilidad a los cambios del entorno. Es por ello que en este trabajo se analizaron los depósitos lacustres que afloran en el sector central de los valles Calchaquíes (región de Cafayate), para comprender cómo se acomoda la deformación cuaternaria en una de las cuencas intermontanas de la cuña orogénica activa.
El rumbo de las estructuras cuaternarias en el área de estudio es subparalelo al de las fallas que exhuman los cordones serranos circundantes. A partir del estudio estratigráfico, morfotectónico y estructural de los depósitos lacustres, se identificó un mínimo de cinco episodios de deformación afectando a la columna estratigráfica cuaternaria. Integrando perfiles estructurales balanceados con edades obtenidas en este trabajo y recopiladas de la bibliografía, se calcularon para el Pleistoceno mediotardío, tasas mínimas y máximas de acortamiento que varían entre 0,19-2,80 y 0,21-4,47 mm/a, respectivamente. Para comparar estos resultados con mediciones de la tectónica activa a escala regional se recopilaron datos de estaciones geodésicas del noroeste argentino, con los cuales se elaboró un perfil de velocidades horizontales. El perfil obtenido muestra un decrecimiento gradual de los vectores hacia el este, indicando actividad interna del orógeno en congruencia con los registros de actividad sísmica y compilación regional de las estructuras cuaternarias.
Además de la caracterización neotectónica de este sector de la Cordillera Oriental, el análisis estratigráfico de los depósitos lacustres ha permitido refinar la evolución geológica del sector central de los valles Calchaquíes durante el Cuaternario. De esta manera se han identificado al menos siete episodios de inundación lacustre relacionados con la desconexión del sistema fluvial con su nivel de base, dando lugar a sucesivos eventos de agradación y erosión. Las cotas máximas alcanzadas por los paleolagos, en conjunto con un modelo hidrológico previamente publicado para esta región, permitieron asimismo efectuar una comparación con el registro paleoclimático regional.
Los resultados de esta tesis representan un aporte significativo al conocimiento de la evolución tectónica y estratigráfica del sector central de los valles Calchaquíes durante el Cuaternario. Por otra parte, su integración a escala regional contribuye a comprender mejor la dinámica de la deformación en la cuña orogénica de piel gruesa del noroeste argentino.
The Salt Range in Pakistan exposes Precambrian to Pleistocene strata outcropping along the Salt Range Thrust (SRT). To better understand the in-situ Cambrian and Pliocene tectonic evolution of the Pakistan Subhimalaya, we have conducted low-temperature thermochronological analysis using apatite (U-Th-Sm)/He and fission track dating. We combine cooling ages from different samples located along the thrust front of the SRT into a thermal model that shows two major cooling events associated with rifting and regional erosion in the Late Palaeozoic and SRT activity since the Pliocene. Our results suggest that the SRT maintained a long-term average shortening rate of similar to 5-6 mm/yr and a high exhumation rate above the SRT ramp since similar to 4 Ma.
The correct orientation of seismic sensors is critical for studies such as full moment tensor inversion, receiver function analysis, and shear-wave splitting. Therefore, the orientation of horizontal components needs to be checked and verified systematically. This study relies on two different waveform-based approaches, to assess the sensor orientations of the broadband network of the Kandilli Observatory and Earthquake Research Institute (KOERI). The network is an important backbone for seismological research in the Eastern Mediterranean Region and provides a comprehensive seismic data set for the North Anatolian fault. In recent years, this region became a worldwide field laboratory for continental transform faults. A systematic survey of the sensor orientations of the entire network, as presented here, facilitates related seismic studies. We apply two independent orientation tests, based on the polarization of P waves and Rayleigh waves to 123 broadband seismic stations, covering a period of 15 yr (2004-2018). For 114 stations, we obtain stable results with both methods. Approximately, 80% of the results agree with each other within 10 degrees. Both methods indicate that about 40% of the stations are misoriented by more than 10 degrees. Among these, 20 stations are misoriented by more than 20 degrees. We observe temporal changes of sensor orientation that coincide with maintenance work or instrument replacement. We provide time-dependent sensor misorientation correction values for the KOERI network in the supplemental material.
This dissertation was carried out as part of the international and interdisciplinary graduate school StRATEGy. This group has set itself the goal of investigating geological processes that take place on different temporal and spatial scales and have shaped the southern central Andes. This study focuses on claystones and carbonates of the Yacoraite Fm. that were deposited between Maastricht and Dan in the Cretaceous Salta Rift Basin. The former rift basin is located in northwest Argentina and is divided into the sub-basins Tres Cruces, Metán-Alemanía and Lomas de Olmedo. The overall motivation for this study was to gain new knowledge about the evolution of marine and lacustrine conditions during the Yacoraite Fm. Deposit in the Tres Cruces and Metán-Alemanía sub-basins. Other important aspects that were examined within the scope of this dissertation are the conversion of organic matter from Yacoraite Fm. into oil and its genetic relationship to selected oils produced and natural oil spills. The results of my study show that the Yacoraite Fm. began to be deposited under marine conditions and that a lacustrine environment developed by the end of the deposition in the Tres Cruces and Metán-Alemanía Basins. In general, the kerogen of Yacoraite Fm. consists mainly of the kerogen types II, III and II / III mixtures. Kerogen type III is mainly found in samples from the Yacoraite Fm., whose TOC values are low. Due to the adsorption of hydrocarbons on the mineral surfaces (mineral matrix effect), the content of type III kerogen with Rock-Eval pyrolysis in these samples could be overestimated. Investigations using organic petrography show that the organic particles of Yacoraite Fm. mainly consist of alginites and some vitrinite-like particles. The pyrolysis GC of the rock samples showed that the Yacoraite Fm. generates low-sulfur oils with a predominantly low-wax, paraffinic-naphthenic-aromatic composition and paraffinic wax-rich oils. Small proportions of paraffinic, low-wax oils and a gas condensate-generating facies are also predicted. Here, too, mineral matrix effects were taken into account, which can lead to a quantitative overestimation of the gas-forming character.
The results of an additional 1D tank modeling carried out show that the beginning (10% TR) of the oil genesis took place between ≈10 Ma and ≈4 Ma. Most of the oil (from ≈50% to 65%) was generated prior to the development of structural traps formed during the Plio-Pleistocene Diaguita deformation phase. Only ≈10% of the total oil generated was formed and potentially trapped after the formation of structural traps. Important factors in the risk assessment of this petroleum system, which can determine the small amounts of generated and migrated oil, are the generally low TOC contents and the variable thickness of the Yacoraite Fm. Additional risks are associated with a low density of information about potentially existing reservoir structures and the quality of the overburden.
Barite stands out as one of the most ubiquitous scaling agents in deep geothermal systems, responsible for irreversible efficiency loss. Due to complex parameter interplay, it is imperative to utilise numerical simulations to investigate temporal and spatial precipitation effects. A one-dimensional reactive transport model is set up with heterogeneous nucleation and crystal growth kinetics. In line with geothermal systems in the North German Basin, the following parameters are considered in a sensitivity analysis: temperature (25 to 150 degrees C), pore pressure (10 to 50 MPa), fracture aperture (10(-4) to 10(-2) m), flow velocity (10(-3) to 10(0) m s(-1)), molar volume (50.3 to 55.6 cm(3) mol(-1)), contact angle for heterogeneous nucleation (0 degrees to 180 degrees), interfacial tension (0.07 to 0.134 J m(-2)), salinity (0.1 to 1.5 mol kgw(-1) NaCl), pH (5 to 7), and supersaturation ratio (1 to 30). Nucleation and consequently crystal growth can only begin if the threshold supersaturation is exceeded, therefore contact angle and interfacial tension are the most sensitive in terms of precipitation kinetics. If nucleation has occurred, crystal growth becomes the dominant process, which is mainly controlled by fracture aperture. Results show that fracture sealing takes place within months (median 33 days) and the affected range can be on the order of tens of metres (median 10 m). The presented models suggest that barite scaling must be recognised as a serious threat if the supersaturation threshold is exceeded, in which case, large fracture apertures could help to minimise kinetic rates. The models further are of use for adjusting the fluid injection temperature.
The sediment profile from Lake Goscia(z) over dot in central Poland comprises a continuous, seasonally resolved and exceptionally well-preserved archive of the Younger Dryas (YD) climate variation. This provides a unique opportunity for detailed investigation of lake system responses during periods of rapid climate cooling (YD onset) and warming (YD termination). The new varve record of Lake Goscia(z) over dot presented here spans 1662 years from the late Allerod (AL) to the early Preboreal (PB). Microscopic varve counting provides an independent chronology with a YD duration of 1149+14/-22 years, which confirms previous results of 1140 +/- 40 years. We link stable oxygen isotopes and chironomid-based air temperature reconstructions with the response of various geochemical and varve microfacies proxies especially focusing on the onset and termination of the YD. Cooling at the YD onset lasted similar to 180 years, which is about a century longer than the terminal warming that was completed in similar to 70 years. During the AL/YD transition, environmental proxy data lagged the onset of cooling by similar to 90 years and revealed an increase of lake productivity and internal lake re-suspension as well as slightly higher detrital sediment input. In contrast, rapid warming and environmental changes during the YD/PB transition occurred simultaneously. However, initial changes such as declining diatom deposition and detrital input occurred already a few centuries before the rapid warming at the YD/PB transition. These environmental changes likely reflect a gradual increase in summer air temperatures already during the YD. Our data indicate complex and differing environmental responses to the major climate changes related to the YD, which involve different proxy sensitivities and threshold processes.