Refine
Has Fulltext
- no (84)
Year of publication
- 2021 (84) (remove)
Document Type
- Article (84) (remove)
Is part of the Bibliography
- yes (84)
Keywords
- PHREEQC (5)
- reactive transport (4)
- Exhumation (2)
- heterogeneity (2)
- numerical modeling (2)
- permafrost (2)
- repository far-field (2)
- sorption (2)
- weathering (2)
- (compliant) cracks (1)
- 3D root (1)
- ALOS World 3D (1)
- ASTER GDEM (1)
- Afyon volcanics (1)
- Aluminium– silicates (1)
- Andean plateau (1)
- Anthropocene (1)
- Ar-39 (1)
- Ar-40 (1)
- Arctic Ocean (1)
- Balanced cross section (1)
- Baltic Sea coast (1)
- Big Naryn complex (1)
- Bushveld Complex (1)
- Central Asia (1)
- Contamination (1)
- Copernicus DEM (1)
- Corundum (1)
- Crystal nucleation (1)
- Cuesta de la chilca (1)
- Cyanobium (1)
- DEM noise (1)
- Dam engineering (1)
- Duplex (1)
- Dynamic bulk modulus (1)
- Décollement (1)
- Earthquake ground motions (1)
- Earthquake hazards (1)
- Earthquake source observations (1)
- Elastic response (1)
- Energy and society (1)
- Environmental sciences (1)
- Equant (stiff) pores (1)
- Europe (1)
- Experimental (1)
- Experimental petrology (1)
- Famatinian (1)
- Fold and thrust belt (1)
- Forschungsmethodik (1)
- Fourier analysis (1)
- Fracture sealing (1)
- Fracture wall offset (1)
- GNSS-integrated water vapour (1)
- Geochemical indices (1)
- Geochemical modelling (1)
- Geochemistry (1)
- Geology (1)
- Geophysics (1)
- Geothermal (1)
- Geothermobarometry (1)
- Germany (1)
- Glaciers (1)
- Global Navigation Satellite System (GNSS) (1)
- Granite (1)
- Grundwasser-Oberflächenwasser-Interaktion (1)
- HVSR (1)
- Heavy metals (1)
- Hierarchical clustering (1)
- High-pressure (1)
- High-temperature experiments (1)
- Hydraulic aperture (1)
- Hydro-meteorological indices (1)
- Inverse theory (1)
- Joint inversion (1)
- Kenya Rift (1)
- Kinematics of crustal and mantle deformation (1)
- Kuilyu complex (1)
- Kurpsai dam (1)
- Lake Issyk‐ Kul (1)
- Lake-surfaces variability (1)
- Landfill (1)
- Landslide inventory (1)
- Landslide susceptibility (1)
- Late Holocene (1)
- Leones Glacier (1)
- Logistic regression (1)
- Lower crust (1)
- MATLAB (1)
- Matrix diffusion (1)
- Mechanical aperture (1)
- Mechanical properties (1)
- Melainabacteria (1)
- Merensky Reef (1)
- Metastability (1)
- Methylheptadecanes (1)
- Micro-structure (1)
- Minor (1)
- Mont Terri (1)
- Multivariate statistical analysis (1)
- Nanolayers (1)
- Nanorocks (1)
- Narrow (1)
- Nearshore zone (1)
- Neogene (1)
- Normal (1)
- North (1)
- Northern Patagonian Icefield (1)
- Numerische Modellierung (1)
- Ocean color remote sensing (1)
- Ocean sciences (1)
- Operational and environmental effects (1)
- Orbital forcing (1)
- Ostseeküste (1)
- P-waves (1)
- PGE (1)
- PGM (1)
- PT path (1)
- Paleoclimate (1)
- Panasqueira (1)
- Peridotitic melts (1)
- Piedmont Zone (1)
- Pitzer (1)
- Plate tectonics (1)
- Porous sandstone (1)
- Poshtuk (1)
- Precambrain (1)
- Precipitation variability (1)
- Probability distributions (1)
- Quartz (1)
- Raoult's law (1)
- Reactive transport (1)
- Regional forcings (1)
- Remote sensing (1)
- Rheology: crust and (1)
- Risk assessment (1)
- S-waves (1)
- SRTM (1)
- Seismic cycle (1)
- Seismic noise (1)
- Self-propping fracture (1)
- Sericytochromatia (1)
- Shan (1)
- Sierras Pampeanas Orientales (1)
- Site effects (1)
- Slope (1)
- Snow-cover (1)
- Soil contamination (1)
- Solid Earth sciences (1)
- Source longevity (1)
- South American Monsoon (1)
- South Korea (1)
- Southern Kyrgyzstan (1)
- Soziale Medien (1)
- Space geodetic surveys (1)
- Spectral analysis (1)
- Static bulk modulus (1)
- Strain (1)
- Strong-motion (1)
- Structural health monitoring (1)
- Sturmhochwasser (1)
- Surface roughness (1)
- System (1)
- TRMM (1)
- TanDEM-X (1)
- Thermal treatment (1)
- Tian (1)
- Tianshan orogenic belt (1)
- Tonian–Cryogenian (1)
- Tungsten–tin deposits (1)
- Turbidity retrieval (1)
- Twitter (1)
- Unloading (1)
- Unterrichtsmethoden (1)
- Upper Indus Basin (1)
- Varves (1)
- Versalzung (1)
- Volatiles (1)
- Volcanology (1)
- Wave propagation (1)
- Waveform inversion (1)
- Western Alps (1)
- WorldDEM (1)
- Younger Dryas (1)
- ZHe tracer thermochronology (1)
- acoustic emissions (1)
- adaptation (1)
- algorithms (1)
- alteration geochemistry (1)
- amplicon sequencing (1)
- analog experiment (1)
- arctic (1)
- arctic water bodies (1)
- asymmetric rifting (1)
- attributes (1)
- automated monitoring (1)
- basin analysis (1)
- blueSeis-3A sensors (1)
- body waves (1)
- broadband seismometer (1)
- bubble trap (1)
- calc-alkaline (1)
- calderas (1)
- carbon stabilization (1)
- carnallite (1)
- cementation (1)
- central Poland (1)
- change detection (1)
- chloritoid stability (1)
- climate change (1)
- climate change adaptation (1)
- climate-change mitigation (1)
- climate-change policy (1)
- coexisting phases (1)
- coherency (1)
- coherent noise (1)
- components (1)
- continental rifts (1)
- convection (1)
- convective available potential energy (CAPE) (1)
- core-shell (1)
- correlation (1)
- damage (1)
- damage assessment (1)
- debris flows (1)
- density-driven (1)
- density-driven flow (1)
- diffusion (1)
- digital core reconstruction (1)
- digital grain sizing (1)
- dike (1)
- dike propagation (1)
- downstream fining (1)
- dryland (1)
- earthquake source (1)
- earthquakes (1)
- ecological succession (1)
- elastic rock properties (1)
- elastic waves (1)
- emergency (1)
- emergency preparedness (1)
- energy transfer (1)
- environmental seismology (1)
- equations (1)
- erosion (1)
- eruptions (1)
- eruptive cycle (1)
- extracellular DNA (1)
- extreme rainfall (1)
- facies (1)
- fatigue (1)
- flood events (1)
- flood risk management (1)
- fluvial geomorphology (1)
- formation damage (1)
- fracture growth (1)
- freeze-thaw cycling (1)
- frost‐ cracking (1)
- full (1)
- fundamental frequency (1)
- gas supply conditions (1)
- geodynamics (1)
- geologic hazards (1)
- geothermal energy (1)
- geyser (1)
- glacial erosion (1)
- governance (1)
- grain (1)
- grain-size distribution (1)
- ground-penetrating radar (1)
- groundwater-surface water (1)
- heuristic analysis (1)
- high pressure (1)
- high temperature (1)
- huddle test (1)
- hydrate formation process (1)
- hydraulic fracturing (1)
- hydro-geological system (1)
- hydrothermal systems (1)
- in (1)
- inheritance (1)
- injections (1)
- injectivity (1)
- innovation (1)
- insurance (1)
- interaction (1)
- interpretation (1)
- intracellular DNA (1)
- lake monitoring (1)
- lake sediments (1)
- lake stratification (1)
- laminography (1)
- landsat (1)
- lanthanides (1)
- lanthanoid migration (1)
- lithobiont (1)
- lithosphere (1)
- location (1)
- machine learning (1)
- magma ocean (1)
- magmatic arc (1)
- magmatic-hydrothermal systems (1)
- magnetostratigraphy (1)
- mass (1)
- meta-ophiolites (1)
- meta-pelites (1)
- micro-CT scan (1)
- mixed gas hydrates (1)
- modeling (1)
- moment tensor (1)
- mountain geomorphology (1)
- multi-spectral analysis (1)
- multi‐ temporal analyses (1)
- natural hazard (1)
- natural hazards (1)
- nature-based solutions (NBS) (1)
- neutron imaging (1)
- noise cross (1)
- numerical (1)
- numerical model (1)
- pebblecounts (1)
- permafrost lowlands (1)
- permafrost thaw (1)
- permeability-porosity relationship (1)
- phreeqc (1)
- picocyanobacteria diversity (1)
- pore-scale (1)
- porous media (1)
- preparedness (1)
- pseudosection (1)
- psychrotolerant (1)
- radial flow (1)
- rainfall events (1)
- recovery (1)
- redox processes (1)
- relief (1)
- response (1)
- retrogressive thaw slumps (1)
- rhizosphere biogeochemistry (1)
- rift variability (1)
- rifting (1)
- rock (1)
- rockfalls (1)
- root activity (1)
- root-soil interaction (1)
- rotational seismology (1)
- salinization (1)
- salt dissolution (1)
- sediment production (1)
- sedimentation pattern (1)
- sedimentology (1)
- seismic monitoring (1)
- seismic waves (1)
- seismology (1)
- sentinel-2 (1)
- shoreline movement direction (1)
- shoreline movement rate (1)
- signal-to-noise ratio (1)
- simulation (1)
- situ Raman spectroscopy (1)
- size fractions (1)
- slope failure (1)
- social (1)
- socioeconomic scenarios (1)
- solar radiation (1)
- source (1)
- source direction (1)
- storm flood (1)
- stress (1)
- stress field (1)
- structural (1)
- subarctic (1)
- sulfur solubility (1)
- supergene ores (1)
- sustainability (1)
- swarms (1)
- system architecture (1)
- temperature extremes (1)
- thaw (1)
- thermal (1)
- thermo-erosion (1)
- tilt (1)
- time series (1)
- topographic asymmetry (1)
- topography (1)
- tourmaline (1)
- trace elements (1)
- trajectories (1)
- transferability (1)
- tremor (1)
- upconversion nanoparticles (1)
- uranium speciation (1)
- urban infill development (1)
- vegetation cover (1)
- very high resolution imagery (1)
- volcanic earthquakes (1)
- wasting (1)
- water rock interaction (1)
- waveform (1)
- white mica (1)
Institute
- Institut für Geowissenschaften (84) (remove)
Ground-penetrating radar (GPR) is a standard geophysical technique used to image near-surface structures in sedimentary environments. In such environments, GPR data acquisition and processing are increasingly following 3D strategies. However, the processed GPR data volumes are typically still interpreted using selected 2D slices and manual concepts such as GPR facies analyses. In seismic volume interpretation, the application of (semi-)automated and reproducible approaches such as 3D attribute analyses as well as the production of attribute-based facies models are common practices today. In contrast, the field of 3D GPR attribute analyses and corresponding facies models is largely untapped. We have developed and applied a workflow to produce 3D attribute-based GPR facies models comprising the dominant sedimentary reflection patterns in a GPR volume, which images complex sandy structures on the dune island of Spiekeroog (Northern Germany). After presenting our field site and details regarding our data acquisition and processing, we calculate and filter 3D texture attributes to generate a database comprising the dominant texture features of our GPR data. Then, we perform a dimensionality reduction of this database to obtain meta texture attributes, which we analyze and integrate using composite imaging and (also considering additional geometric information) fuzzy c-means cluster analysis resulting in a classified GPR facies model. Considering our facies model and a corresponding GPR facies chart, we interpret our GPR data set in terms of near-surface sedimentary units, the corresponding depositional environments, and the recent formation history at our field site. Thus, we demonstrate the potential of our workflow, which represents a novel and clear strategy to perform a more objective and consistent interpretation of 3D GPR data collected across different sedimentary environments.
We present a new autoclave that enables in situ characterization of hydrothermal fluids at high pressures and high temperatures at synchrotron x-ray radiation sources. The autoclave has been specifically designed to enable x-ray absorption spectroscopy in fluids with applications to mineral solubility and element speciation analysis in hydrothermal fluids in complex compositions. However, other applications, such as Raman spectroscopy, in high-pressure fluids are also possible with the autoclave. First experiments were run at pressures between 100 and 600 bars and at temperatures between 25 degrees C and 550 degrees C, and preliminary results on scheelite dissolution in fluids of different compositions show that the autoclave is well suited to study the behavior of ore-forming metals at P-T conditions relevant to the Earth's crust.
Atmospheric water vapour content is a key variable that controls the development of deep convective storms and rainfall extremes over the central Andes. Direct measurements of water vapour are challenging; however, recent developments in microwave processing allow the use of phase delays from L-band radar to measure the water vapour content throughout the atmosphere: Global Navigation Satellite System (GNSS)-based integrated water vapour (IWV) monitoring shows promising results to measure vertically integrated water vapour at high temporal resolutions. Previous works also identified convective available potential energy (CAPE) as a key climatic variable for the formation of deep convective storms and rainfall in the central Andes. Our analysis relies on GNSS data from the Argentine Continuous Satellite Monitoring Network, Red Argentina de Monitoreo Satelital Continuo (RAMSAC) network from 1999 to 2013. CAPE is derived from version 2.0 of the ECMWF’s (European Centre for Medium-Range Weather Forecasts) Re-Analysis (ERA-interim) and rainfall from the TRMM (Tropical Rainfall Measuring Mission) product. In this study, we first analyse the rainfall characteristics of two GNSS-IWV stations by comparing their complementary cumulative distribution function (CCDF). Second, we separately derive the relation between rainfall vs. CAPE and GNSS-IWV. Based on our distribution fitting analysis, we observe an exponential relation of rainfall to GNSS-IWV. In contrast, we report a power-law relationship between the daily mean value of rainfall and CAPE at the GNSS-IWV station locations in the eastern central Andes that is close to the theoretical relationship based on parcel theory. Third, we generate a joint regression model through a multivariable regression analysis using CAPE and GNSS-IWV to explain the contribution of both variables in the presence of each other to extreme rainfall during the austral summer season. We found that rainfall can be characterised with a higher statistical significance for higher rainfall quantiles, e.g., the 0.9 quantile based on goodness-of-fit criterion for quantile regression. We observed different contributions of CAPE and GNSS-IWV to rainfall for each station for the 0.9 quantile. Fourth, we identify the temporal relation between extreme rainfall (the 90th, 95th, and 99th percentiles) and both GNSS-IWV and CAPE at 6 h time steps. We observed an increase before the rainfall event and at the time of peak rainfall—both for GNSS-integrated water vapour and CAPE. We show higher values of CAPE and GNSS-IWV for higher rainfall percentiles (99th and 95th percentiles) compared to the 90th percentile at a 6-h temporal scale. Based on our correlation analyses and the dynamics of the time series, we show that both GNSS-IWV and CAPE had comparable magnitudes, and we argue to consider both climatic variables when investigating their effect on rainfall extremes.
Mineral resource exploration and mining is an essential part of today's high-tech industry. Elements such as rare-earth elements (REEs) and copper are, therefore, in high demand. Modern exploration techniques from multiple platforms (e.g., spaceborne and airborne), to detect and map the spectral characteristics of the materials of interest, require spectral libraries as an essential reference. They include field and laboratory spectral information in combination with geochemical analyses for validation. Here, we present a collection of REE- and copper-related hyperspectral spectra with associated geochemical information. The libraries contain reflectance spectra from rare-earth element oxides, REE-bearing minerals, copper-bearing minerals and mine surface samples from the Apliki copper-gold-pyrite mine in the Republic of Cyprus. The samples were measured with the HySpex imaging spectrometers in the visible and near infrared (VNIR) and shortwave infrared (SWIR) range (400-2500 nm). The geochemical validation of each sample is provided with the reflectance spectra. The spectral libraries are openly available to assist future mineral mapping campaigns and laboratory spectroscopic analyses. The spectral libraries and corresponding geochemistry are published via GFZ Data Services with the following DOIs: https://doi.org/10.5880/GFZ.1.4.2019.004 (13 REE-bearing minerals and 16 oxide powders, Koerting et al., 2019a), https://doi.org/10.5880/GFZ.1.4.2019.003 (20 copper-bearing minerals, Koellner et al., 2019), and https://doi.org/10.5880/GFZ.1.4.2019.005 (37 copper-bearing surface material samples from the Apliki coppergold-pyrite mine in Cyprus, Koerting et al., 2019b). All spectral libraries are united and comparable by the internally consistent method of hyperspectral data acquisition in the laboratory.
Ambitious climate policies, as well as economic development, education, technological progress and less resource-intensive lifestyles, are crucial elements for progress towards the UN Sustainable Development Goals (SDGs). However, using an integrated modelling framework covering 56 indicators or proxies across all 17 SDGs, we show that they are insufficient to reach the targets. An additional sustainable development package, including international climate finance, progressive redistribution of carbon pricing revenues, sufficient and healthy nutrition and improved access to modern energy, enables a more comprehensive sustainable development pathway. We quantify climate and SDG outcomes, showing that these interventions substantially boost progress towards many aspects of the UN Agenda 2030 and simultaneously facilitate reaching ambitious climate targets. Nonetheless, several important gaps remain; for example, with respect to the eradication of extreme poverty (180 million people remaining in 2030). These gaps can be closed by 2050 for many SDGs while also respecting the 1.5 °C target and several other planetary boundaries.
A tale of shifting relations
(2021)
Understanding the dynamics between the East Asian summer (EASM) and winter monsoon (EAWM) is needed to predict their variability under future global warming scenarios. Here, we investigate the relationship between EASM and EAWM as well as the mechanisms driving their variability during the last 10,000 years by stacking marine and terrestrial (non-speleothem) proxy records from the East Asian realm. This provides a regional and proxy independent signal for both monsoonal systems. The respective signal was subsequently analysed using a linear regression model. We find that the phase relationship between EASM and EAWM is not time-constant and significantly depends on orbital configuration changes. In addition, changes in the Atlantic Meridional Overturning circulation, Arctic sea-ice coverage, El Niño-Southern Oscillation and Sun Spot numbers contributed to millennial scale changes in the EASM and EAWM during the Holocene. We also argue that the bulk signal of monsoonal activity captured by the stacked non-speleothem proxy records supports the previously argued bias of speleothem climatic archives to moisture source changes and/or seasonality.
Centroid moment tensor (CMT) parameters can be estimated from seismic waveforms. Since these data indirectly observe the deformation process, CMTs are inferred as solutions to inverse problems which are generally underdetermined and require significant assumptions, including assumptions about data noise. Broadly speaking, we consider noise to include both theory and measurement errors, where theory errors are due to assumptions in the inverse problem and measurement errors are caused by the measurement process. While data errors are routinely included in parameter estimation for full CMTs, less attention has been paid to theory errors related to velocity-model uncertainties and how these affect the resulting moment-tensor (MT) uncertainties. Therefore, rigorous uncertainty quantification for CMTs may require theory-error estimation which becomes a problem of specifying noise models. Various noise models have been proposed, and these rely on several assumptions. All approaches quantify theory errors by estimating the covariance matrix of data residuals. However, this estimation can be based on explicit modelling, empirical estimation and/or ignore or include covariances. We quantitatively compare several approaches by presenting parameter and uncertainty estimates in nonlinear full CMT estimation for several simulated data sets and regional field data of the M-1 4.4, 2015 June 13 Fox Creek, Canada, event. While our main focus is at regional distances, the tested approaches are general and implemented for arbitrary source model choice. These include known or unknown centroid locations, full MTs, deviatoric MTs and double-couple MTs. We demonstrate that velocity-model uncertainties can profoundly affect parameter estimation and that their inclusion leads to more realistic parameter uncertainty quantification. However, not all approaches perform equally well. Including theory errors by estimating non-stationary (non-Toeplitz) error covariance matrices via iterative schemes during Monte Carlo sampling performs best and is computationally most efficient. In general, including velocity-model uncertainties is most important in cases where velocity structure is poorly known.
Concentrations of lead (Pb), zinc (Zn), copper (Cu), mercury (Hg), and arsenic (As) in soils at the Kpone landfill site (Ghana) were determined using Atomic Absorption Spectrophotometry (AAS). Further analyses allowed establishing the degree of heavy metals (HMs) pollution, suitability of the soils for agriculture, sources of the HMs and their ecological and health risks. The site was divided into five zones, A, B, C, D, and E, and in all, seventeen (17) soil samples were collected. Average concentrations of Cu fell within the allowable range for agricultural soils in all the zones while average concentrations of Pb, Zn, Hg, and As exceeded the range in some or all the zones. Concentrations of the HMs generally exceeded their respective background value, with all zones showing very high degree of HMs contamination. The pollution load index (PLI) was 16.48, signifying extreme HMs pollution of the entire site. Multivariate statistical analyses revealed that Cu, Zn, and Pb in the soils originated from the deposited waste materials as well as traffic-related activities (e.g. wear and tear of tyres, brakes, and engines) at the site. Hg also originated from the deposited waste materials as well as cement production and oil and coal combustion activities in the study area, while As derived from industrial discharges and metal smelting activities. All the zones exhibited very high ecological risk. The carcinogenic and non-carcinogenic health risks posed by the HMs were also above acceptable levels, with children being more vulnerable than adults to these health risks.
Barite scalings are a common cause of permanent formation damage to deep geothermal reservoirs. Well injectivity can be impaired because the ooling of saline fluids reduces the solubility of barite, and the continuous re-injection of supersaturated fluids forces barite to precipitate in the host rock. Stimulated reservoirs in the Upper Rhine Graben often have multiple relevant flow paths in the porous matrix and fracture zones, sometimes spanning multiple stratigraphical units to achieve the economically necessary injectivity. While the influence of barite scaling on injectivity has been investigated for purely porous media, the role of fractures within reservoirs consisting of both fractured and porous sections is still not well understood. Here, we present hydro-chemical simulations of a dual-layer geothermal reservoir to study the long-term impact of barite scale formation on well injectivity. Our results show that, compared to purely porous reservoirs, fractured porous reservoirs have a significantly reduced scaling risk by up to 50%, depending on the flow rate ratio of fractures. Injectivity loss is doubled, however, if the amount of active fractures is increased by one order of magnitude, while the mean fracture aperture is decreased, provided the fractured aquifer dictates the injection rate. We conclude that fractured, and especially hydraulically stimulated, reservoirs are generally less affected by barite scaling and that large, but few, fractures are favourable. We present a scaling score for fractured-porous reservoirs, which is composed of easily derivable quantities such as the radial equilibrium length and precipitation potential. This score is suggested for use approximating the scaling potential and its impact on injectivity of a fractured-porous reservoir for geothermal exploitation.
Nature-based solutions (NBS) are seen as a promising adaptation measure that sustainably deals with diverse societal challenges, while simultaneously delivering multiple benefits. Nature-based solutions have been highlighted as a resilient and sustainable means of mitigating floods and other hazards globally. This study examined diverging conceptualizations of NBS, as well as the attitudinal (for example, emotions and beliefs) and contextual (for example, legal and political aspects) barriers and drivers of NBS for flood risks in South Korea. Semistructured interviews were conducted with 11 experts and focused on the topic of flood risk measures and NBS case studies. The analysis found 11 barriers and five drivers in the attitudinal domain, and 13 barriers and two drivers in the contextual domain. Most experts see direct monetary benefits as an important attitudinal factor for the public. Meanwhile, the cost-effectiveness of NBS and their capacity to cope with flood risks were deemed influential factors that could lead decision makers to opt for NBS. Among the contextual factors, insufficient systems to integrate NBS in practice and the ideologicalization of NBS policy were found to be peculiar barriers, which hinder consistent realization of initiatives and a long-term national plan for NBS. Understanding the barriers and drivers related to the mainstreaming of NBS is critical if we are to make the most of such solutions for society and nature. It is also essential that we have a shared definition, expectation, and vision of NBS.
Quantitative geomorphic research depends on accurate topographic data often collected via remote sensing. Lidar, and photogrammetric methods like structure-from-motion, provide the highest quality data for generating digital elevation models (DEMs). Unfortunately, these data are restricted to relatively small areas, and may be expensive or time-consuming to collect. Global and near-global DEMs with 1 arcsec (∼30 m) ground sampling from spaceborne radar and optical sensors offer an alternative gridded, continuous surface at the cost of resolution and accuracy. Accuracy is typically defined with respect to external datasets, often, but not always, in the form of point or profile measurements from sources like differential Global Navigation Satellite System (GNSS), spaceborne lidar (e.g., ICESat), and other geodetic measurements. Vertical point or profile accuracy metrics can miss the pixel-to-pixel variability (sometimes called DEM noise) that is unrelated to true topographic signal, but rather sensor-, orbital-, and/or processing-related artifacts. This is most concerning in selecting a DEM for geomorphic analysis, as this variability can affect derivatives of elevation (e.g., slope and curvature) and impact flow routing. We use (near) global DEMs at 1 arcsec resolution (SRTM, ASTER, ALOS, TanDEM-X, and the recently released Copernicus) and develop new internal accuracy metrics to assess inter-pixel variability without reference data. Our study area is in the arid, steep Central Andes, and is nearly vegetation-free, creating ideal conditions for remote sensing of the bare-earth surface. We use a novel hillshade-filtering approach to detrend long-wavelength topographic signals and accentuate short-wavelength variability. Fourier transformations of the spatial signal to the frequency domain allows us to quantify: 1) artifacts in the un-projected 1 arcsec DEMs at wavelengths greater than the Nyquist (twice the nominal resolution, so > 2 arcsec); and 2) the relative variance of adjacent pixels in DEMs resampled to 30-m resolution (UTM projected). We translate results into their impact on hillslope and channel slope calculations, and we highlight the quality of the five DEMs. We find that the Copernicus DEM, which is based on a carefully edited commercial version of the TanDEM-X, provides the highest quality landscape representation, and should become the preferred DEM for topographic analysis in areas without sufficient coverage of higher-quality local DEMs.
The Big Naryn Complex (BNC) in the East Djetim-Too Range of the Kyrgyz Middle Tianshan block is a tectonized, at least 2 km thick sequence of predominantly felsic to intermediate volcanic rocks intruded by porphyric rhyolite sills. It overlies a basement of metamorphic rocks and is overlain by late Neoproterozoic Djetim-Too Formation sediments; these also occur as tectonic intercalations in the BNC. The up to ca. 1100 m thick Lower Member is composed of predominantly rhyolites-to-dacites and minor basalts, while the at least 900 m thick pyroclastic Upper Member is dominated by rhyolitic-to-dacitic ignimbrites. Porphyric rhyolite sills are concentrated at the top of the Lower Member. A Lower Member rhyolite and a sill sample have LA-ICP-MS U-Pb zircon crystallization ages of 726.1 +/- 2.2 Ma and 720.3 +/- 6.5 Ma, respectively, showing that most of the magmatism occurred within a short time span in the late Tonian-early Cryogenian. Inherited zircons in the sill sample have Neoarchean (2.63, 2.64 Ga), Paleo- (2.33-1.81 Ga), Meso- (1.55 Ga), and Neoproterozoic (ca. 815 Ma) ages, and were derived from a heterogeneous Kuilyu Complex basement. A 1751 +/- 7 Ma Ar-40/Ar-39 age for amphibole from metagabbro is the age of cooling subsequent to Paleoproterozoic metamorphism of the Kuilyu Complex. The large amount of pyroclastic rocks, and their major and trace element compositions, the presence of Neoarchean to Neoproterozoic inherited zircons and a depositional basement of metamorphic rocks point to formation of the BNC in a continental magmatic arc setting.
Insolation differences play a primary role in controlling microclimate and vegetation cover, which together influence the development of topography. Topographic asymmetry (TA), or slope differences between terrain aspects, has been well documented in small-scale, field-based, and modeling studies. Here we combine a suite of environmental (e.g., vegetation, temperature, solar insolation) and topographic (e.g., elevation, drainage network) data to explore the driving mechanisms and markers of TA on a global scale. Using a novel empirical TA analysis method, we find that (1) steeper terrain has higher TA magnitudes, (2) globally, pole-facing terrain is on average steeper than equator-facing terrain, especially in mid-latitude, tectonically quiescent, and vegetated landscapes, and (3) high-elevation and low-temperature regions tend to have terrain steepened toward the equator. We further show that there are distinct differences in climate and vegetation cover across terrain aspects, and that TA is reflected in the size and form of fluvial drainage networks. Our work supports the argument that insolation asymmetries engender differences in local microclimates and vegetation on opposing terrain aspects, which broadly encourage the development of asymmetric topography across a range of lithologic, tectonic, geomorphic, and climatic settings.
The Eastern Mediterranean is the most seismically active region in Europe due to the complex interactions of the Arabian, African, and Eurasian tectonic plates. Deformation is achieved by faulting in the brittle crust, distributed flow in the viscoelastic lower-crust and mantle, and Hellenic subduction, but the long-term partitioning of these mechanisms is still unknown. We exploit an extensive suite of geodetic observations to build a kinematic model connecting strike-slip deformation, extension, subduction, and shear localization across Anatolia and the Aegean Sea by mapping the distribution of slip and strain accumulation on major active geological structures. We find that tectonic escape is facilitated by a plate-boundary-like, translithospheric shear zone extending from the Gulf of Evia to the Turkish-Iranian Plateau that underlies the surface trace of the North Anatolian Fault. Additional deformation in Anatolia is taken up by a series of smaller-scale conjugate shear zones that reach the upper mantle, the largest of which is located beneath the East Anatolian Fault. Rapid north-south extension in the western part of the system, driven primarily by Hellenic Trench retreat, is accommodated by rotation and broadening of the North Anatolian mantle shear zone from the Sea of Marmara across the north Aegean Sea, and by a system of distributed transform faults and rifts including the rapidly extending Gulf of Corinth in central Greece and the active grabens of western Turkey. Africa-Eurasia convergence along the Hellenic Arc occurs at a median rate of 49.8mm yr(-1) in a largely trench-normal direction except near eastern Crete where variably oriented slip on the megathrust coincides with mixed-mode and strike-slip deformation in the overlying accretionary wedge near the Ptolemy-Pliny-Strabo trenches. Our kinematic model illustrates the competing roles the North Anatolian mantle shear zone, Hellenic Trench, overlying mantle wedge, and active crustal faults play in accommodating tectonic indentation, slab rollback and associated Aegean extension. Viscoelastic flow in the lower crust and upper mantle dominate the surface velocity field across much of Anatolia and a clear transition to megathrust-related slab pull occurs in western Turkey, the Aegean Sea and Greece. Crustal scale faults and the Hellenic wedge contribute only a minor amount to the large-scale, regional pattern of Eastern Mediterranean interseismic surface deformation.
Complex, time-dependent, and asymmetric rift geometries are observed throughout the East African Rift System (EARS) and are well documented, for instance, in the Kenya Rift. To unravel asymmetric rifting processes in this region, we conduct 2D geodynamic models. We use the finite element software ASPECT employing visco-plastic rheologies, mesh-refinement, distributed random noise seeding, and a free surface. In contrast to many previous numerical modeling studies that aimed at understanding final rifted margin symmetry, we explicitly focus on initial rifting stages to assess geodynamic controls on strain localization and fault evolution. We thereby link to geological and geophysical observations from the Southern and Central Kenya Rift. Our models suggest a three-stage early rift evolution that dynamically bridges previously inferred fault-configuration phases of the eastern EARS branch: (1) accommodation of initial strain localization by a single border fault and flexure of the hanging-wall crust, (2) faulting in the hanging-wall and increasing upper-crustal faulting in the rift-basin center, and (3) loss of pronounced early stage asymmetry prior to basinward localization of deformation. This evolution may provide a template for understanding early extensional faulting in other branches of the East African Rift and in asymmetric rifts worldwide. By modifying the initial random noise distribution that approximates small-scale tectonic inheritance, we show that a spectrum of first-order fault configurations with variable symmetry can be produced in models with an otherwise identical setup. This approach sheds new light on along-strike rift variability controls in active asymmetric rifts and proximal rifted margins.
The metastable paragenesis of corundum and quartz is rare in nature but common in laboratory experiments where according to thermodynamic predictions aluminum-silicate polymorphs should form. We demonstrate here that the existence of a hydrous, silicon-bearing, nanometer-thick layer (called "HSNL") on the corundum surface can explain this metastability in experimental studies without invoking unspecific kinetic inhibition. We investigated experimentally formed corundum reaction products synthesized during hydrothermal and piston-cylinder experiments at 500-800 degrees C and 0.25-1.8 GPa and found that this HSNL formed inside and on the corundum crystals, thereby controlling the growth behavior of its host. The HSNL represents a substitution of Al with Si and H along the basal plane of corundum. Along the interface of corundum and quartz, the HSNL effectively isolates the bulk phases corundum and quartz from each other, thus apparently preventing their reaction to the stable aluminum silicate. High temperatures and prolonged experimental duration lead to recrystallization of corundum including the HSNL and to the formation of quartz + fluid inclusions inside the host crystal. This process reduces the phase boundary area between the bulk phases, thereby providing further opportunity to expand their coexistence. In addition to its small size, its transient nature makes it difficult to detect the HSNL in experiments and even more so in natural samples. Our findings emphasize the potential impact of nanometer-sized phases on geochemical reaction pathways and kinetics under metamorphic conditions in one of the most important chemical systems of the Earth's crust.
The computational costs associated with coupled reactive transport simulations are mostly due to the chemical subsystem: replacing it with a pre-trained statistical surrogate is a promising strategy to achieve decisive speedups at the price of small accuracy losses and thus to extend the scale of problems which can be handled. We introduce a hierarchical coupling scheme in which "full-physics" equation-based geochemical simulations are partially replaced by surrogates. Errors in mass balance resulting from multivariate surrogate predictions effectively assess the accuracy of multivariate regressions at runtime: inaccurate surrogate predictions are rejected and the more expensive equation-based simulations are run instead. Gradient boosting regressors such as XGBoost, not requiring data standardization and being able to handle Tweedie distributions, proved to be a suitable emulator. Finally, we devise a surrogate approach based on geochemical knowledge, which overcomes the issue of robustness when encountering previously unseen data and which can serve as a basis for further development of hybrid physics-AI modelling.
Diagenetic trends of synthetic reservoir sandstone properties assessed by digital rock physics
(2021)
Quantifying interactions and dependencies among geometric, hydraulic and mechanical properties of reservoir sandstones is of particular importance for the exploration and utilisation of the geological subsurface and can be assessed by synthetic sandstones comprising the microstructural complexity of natural rocks. In the present study, three highly resolved samples of the Fontainebleau, Berea and Bentheim sandstones are generated by means of a process-based approach, which combines the gravity-driven deposition of irregularly shaped grains and their diagenetic cementation by three different schemes. The resulting evolution in porosity, permeability and rock stiffness is examined and compared to the respective micro-computer tomographic (micro-CT) scans. The grain contact-preferential scheme implies a progressive clogging of small throats and consequently produces considerably less connected and stiffer samples than the two other schemes. By contrast, uniform quartz overgrowth continuously alters the pore space and leads to the lowest elastic properties. The proposed stress-dependent cementation scheme combines both approaches of contact-cement and quartz overgrowth, resulting in granulometric, hydraulic and elastic properties equivalent to those of the respective micro-CT scans, where bulk moduli slightly deviate by 0.8%, 4.9% and 2.5% for the Fontainebleau, Berea and Bentheim sandstone, respectively. The synthetic samples can be further altered to examine the impact of mineral dissolution or precipitation as well as fracturing on various petrophysical correlations, which is of particular relevance for numerous aspects of a sustainable subsurface utilisation.
The eruption frequency of geysers can be studied easily on the surface. However, details of the internal structure including possible water and gas filled chambers feeding eruptions and the driving mechanisms often remain elusive. We used a multidisciplinary network of seismometers, video cameras, water pressure sensors and one tiltmeter to study the eruptive cycle, internal structure, and mechanisms driving the eruptive cycle of Strokkur geyser in June 2018. An eruptive cycle at Strokkur always consists of four phases: (1) Eruption, (2) post-eruptive conduit refilling, (3) gas filling of the bubble trap, and (4) regular bubble collapse at shallow depth in the conduit. For a typical single eruption 19 +/- 4 bubble collapses occur in Phase 3 and 8 +/- 2 collapses in Phase 4 at a mean spacing of 1.52 +/- 0.29 and 24.5 +/- 5.9 s, respectively. These collapses release latent heat to the fluid in the bubble trap (Phase 3) and later to the fluid in the conduit (Phase 4). The latter eventually reaches thermodynamic conditions for an eruption. Single to sextuple eruptions have similar spacings between bubble collapses and are likely fed from the same bubble trap at 23.7 +/- 4.4 m depth, 13-23 m west of the conduit. However, the duration of the eruption and recharging phase linearly increases likely due to a larger water, gas and heat loss from the system. Our tremor data provides documented evidence for a bubble trap beneath a pool geyser.