Refine
Year of publication
- 2024 (5)
- 2023 (34)
- 2022 (114)
- 2021 (126)
- 2020 (205)
- 2019 (276)
- 2018 (314)
- 2017 (298)
- 2016 (292)
- 2015 (256)
- 2014 (228)
- 2013 (187)
- 2012 (171)
- 2011 (141)
- 2010 (90)
- 2009 (120)
- 2008 (22)
- 2007 (15)
- 2006 (78)
- 2005 (73)
- 2004 (66)
- 2003 (24)
- 2002 (19)
- 2001 (30)
- 2000 (46)
- 1999 (64)
- 1998 (75)
- 1997 (74)
- 1996 (49)
- 1995 (29)
- 1994 (33)
- 1993 (3)
- 1992 (4)
- 1991 (1)
Document Type
- Article (2712)
- Doctoral Thesis (504)
- Postprint (139)
- Other (73)
- Review (51)
- Monograph/Edited Volume (34)
- Preprint (17)
- Conference Proceeding (12)
- Habilitation Thesis (12)
- Master's Thesis (6)
Keywords
- climate change (49)
- Holocene (44)
- erosion (28)
- Himalaya (26)
- remote sensing (25)
- permafrost (23)
- Climate change (21)
- Tibetan Plateau (21)
- Pollen (20)
- climate (20)
Institute
- Institut für Geowissenschaften (3567) (remove)
The new in situ geodynamic laboratory established in the framework of the ICDP Eger project aims to develop the most modern, comprehensive, multiparameter laboratory at depth for studying earthquake swarms, crustal fluid flow, mantle-derived CO2 and helium degassing, and processes of the deep biosphere. In order to reach a new level of high-frequency, near-source and multiparameter observation of earthquake swarms and related phenomena, such a laboratory comprises a set of shallow boreholes with high-frequency 3-D seismic arrays as well as modern continuous real-time fluid monitoring at depth and the study of the deep biosphere.
This laboratory is located in the western part of the Eger Rift at the border of the Czech Republic and Germany (in the West Bohemia–Vogtland geodynamic region) and comprises a set of five boreholes around the seismoactive zone. To date, all monitoring boreholes have been drilled. This includes the seismic monitoring boreholes S1, S2 and S3 in the crystalline units north and east of the major Nový Kostel seismogenic zone, borehole F3 in the Hartoušov mofette field and borehole S4 in the newly discovered Bažina maar near Libá. Supplementary borehole P1 is being prepared in the Neualbenreuth maar for paleoclimate and biological research. At each of these sites, a borehole broadband seismometer will be installed, and sites S1, S2 and S3 will also host a 3-D seismic array composed of a vertical geophone chain and surface seismic array. Seismic instrumenting has been completed in the S1 borehole and is in preparation in the remaining four monitoring boreholes. The continuous fluid monitoring site of Hartoušov includes three boreholes, F1, F2 and F3, and a pilot monitoring phase is underway. The laboratory also enables one to analyze microbial activity at CO2 mofettes and maar structures in the context of changes in habitats. The drillings into the maar volcanoes contribute to a better understanding of the Quaternary paleoclimate and volcanic activity.
The spatial pattern of extreme precipitation from 40 years of gauge data in the central Himalaya
(2022)
The topography of the Himalaya exerts a substantial control on the spatial distribution of monsoonal rainfall, which is a vital water source for the regional economy and population. But the occurrence of short-lived and high-intensity precipitation results in socio-economic losses. This study relies on 40 years of daily data from 204 ground stations in Nepal to derive extreme precipitation thresholds, amounts, and days at the 95th percentile. We additionally determine the precipitation magnitude-frequency relation. We observe that extreme precipitation amounts follow an almost uniform band parallel to topographic contour lines in the southern Himalaya mountains in central and eastern Nepal but not in western Nepal. The relationship of extreme precipitation indices with topographic relief shows that extreme precipitation thresholds decrease with increasing elevation, but extreme precipitation days increase in higher elevation areas. Furthermore, stations above 1 km elevation exhibit a power-law relation in the rainfall magnitude-frequency framework. Stations at higher elevations generally have lower values of power-law exponents than low elevation areas. This suggests a fundamentally different behaviour of the rainfall distribution and an increased occurrence of extreme rainfall storms in the high elevation areas of Nepal.
A mid-aperture broad-band test array (OBS array DOCTAR) was deployed from June 2011 to April 2012 about 100 km north of the Gloria fault in the Eastern North Atlantic in about 5000 m water depth. In addition arrays were installed on Madeira Island and in western Portugal mainland. For the first time in the Eastern North Atlantic, we recorded a large number of high frequency Po and So waves from local and regional small and moderate earthquakes (M-L < 4). An incoherent beamforming method was adapted to scan continuous data for such Po and So arrivals applying a sliding window waveform migration and frequency-wavenumber technique. We identify about 320 Po and 1550 So arrivals and compare the phase onsets with the ISC catalogue (ISC 2015) for the same time span. Up to a distance of 6 degrees to the DOCTAR stations all events listed in the ISC catalogue could be associated to Po and So phases. Arrivals from events in more than 10 degrees distance could be identified only in some cases. Only few Po and/or So arrivals were detected for earthquakes from the European and African continental area, the continental shelf regions and for earthquakes within or northwest of the Azores plateau. Unexpectedly, earthquake clusters are detected within the oceanic plates north and south of the Gloria fault and far from plate boundaries, indicating active intraplate structures. We also observe and locate numerous small magnitude earthquakes on the segment of the Gloria fault directly south of DOCTAR, which likely coincides with the rupture of the 25 November 1941 event. Local small magnitude earthquakes located beneath DOCTAR show hypocentres up to 30 km depth and strike-slip focal mechanisms. A comparison with detections at temporary mid-aperture arrays on Madeira and in western Portugal shows that the deep ocean array performs much better than the island and the continental array regarding the detection threshold for events in the oceanic plates.
We conclude that sparsely distributed mid-aperture seismic arrays in the deep ocean could decrease the detection and location threshold for seismicity with M-L < 4 in the oceanic plate and might constitute a valuable tool to monitor oceanic plate seismicity.
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.
The Gofa Province and the Chew Bahir Basin of southern Ethiopia constitute tectonically active regions, where the Southern Main Ethiopian Rift converges with the Northern Kenya Rift through a wide zone of extensional deformation with several north to northeast-trending, left-stepping en-e & PRIME;chelon basins. This sector of the Southern Main Ethiopian Rift is characterized by a semi-arid climate and a largely uniform lithology, and thus provides ideal conditions for studying the different parameters that define the tectonic and geomorphic features of this complex kinematic transfer zone. In this study, the degree of tectonic activity, spatiotemporal variations in extension, and the nature of kinematic linkage between different fault systems of the transfer zone are constrained by detailed quantitative geomorphic analysis of river catchments and focused field work. We analyzed fluvial and landscape morphometric characteristics in combination with structural, seismicity, and climatic data to better evaluate the tectono-geomorphic history of this transfer zone. Our data reveal significant north-south variations in the degree of extension from the Sawula Basin in the north (mature) to the Chew Bahir Basin in the south (juvenile). First, normalized channel-steepness indices and the spatial arrangement of knickpoints in footwall-draining streams suggest a gradual, southward shift in extensional deformation and recent tectonic activity. Second, based on 1-k(m) radius local relief and mean-hillslope maximum values that are consistent with ksn anomalies, we confirm strain localization within zones of fault interaction. Third, morphometric indices such as hypsometry, basin asymmetry factor, and valley floor width to valley height ratio also indicate a north to south gradient in tectonic activity, highlighting the importance of such a wide transfer zone with diffuse extension linking different rift segments during the break-up of continental crust.
On 7 January 2020, an M-w 6.4 earthquake occurred in the northeastern Caribbean, a few kilometers offshore of the island of Puerto Rico. It was the mainshock of a complex seismic sequence, characterized by a large number of energetic earthquakes illuminating an east-west elongated area along the southwestern coast of Puerto Rico. Deformation fields constrained by Interferometric Synthetic Aperture Radar and Global Navigation Satellite System data indicate that the coseismic movements affected only the western part of the island. To assess the mainshock's source fault parameters, we combined the geodetically derived coseismic deformation with teleseismic waveforms using Bayesian inference. The results indicate a roughly east-west oriented fault, dipping northward and accommodating similar to 1.4 m of transtensional motion. Besides, the determined location and orientation parameters suggest an offshore continuation of the recently mapped North Boqueron Bay-Punta Montalva fault in southwest Puerto Rico. This highlights the existence of unmapped faults with moderate-to-large earthquake potential within the Puerto Rico region.
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.
The new in situ geodynamic laboratory established in the framework of the ICDP Eger project aims to develop the most modern, comprehensive, multiparameter laboratory at depth for studying earthquake swarms, crustal fluid flow, mantle-derived CO2 and helium degassing, and processes of the deep biosphere. In order to reach a new level of high-frequency, near-source and multiparameter observation of earthquake swarms and related phenomena, such a laboratory comprises a set of shallow boreholes with high-frequency 3-D seismic arrays as well as modern continuous real-time fluid monitoring at depth and the study of the deep biosphere.
This laboratory is located in the western part of the Eger Rift at the border of the Czech Republic and Germany (in the West Bohemia-Vogtland geodynamic region) and comprises a set of five boreholes around the seismoactive zone. To date, all monitoring boreholes have been drilled. This includes the seismic monitoring boreholes S1, S2 and S3 in the crystalline units north and east of the major Novy Kostel seismogenic zone, borehole F3 in the Hartousov mofette field and borehole S4 in the newly discovered Bazina maar near Liba. Supplementary borehole P1 is being prepared in the Neualbenreuth maar for paleoclimate and biological research. At each of these sites, a borehole broadband seismometer will be installed, and sites S1, S2 and S3 will also host a 3-D seismic array composed of a vertical geophone chain and surface seismic array. Seismic instrumenting has been completed in the S1 borehole and is in preparation in the remaining four monitoring boreholes. The continuous fluid monitoring site of Hartousov includes three boreholes, F1, F2 and F3, and a pilot monitoring phase is underway. The laboratory also enables one to analyze microbial activity at CO2 mofettes and maar structures in the context of changes in habitats. The drillings into the maar volcanoes contribute to a better understanding of the Quaternary paleoclimate and volcanic activity.
In this article, we address the question of how observed ground-motion data can most effectively be modeled for engineering seismological purposes. Toward this goal, we use a data-driven method, based on a deep-learning autoencoder with a variable number of nodes in the bottleneck layer, to determine how many parameters are needed to reconstruct synthetic and observed ground-motion data in terms of their median values and scatter. The reconstruction error as a function of the number of nodes in the bottleneck is used as an indicator of the underlying dimensionality of ground-motion data, that is, the minimum number of predictor variables needed in a ground-motion model. Two synthetic and one observed datasets are studied to prove the performance of the proposed method. We find that mapping ground-motion data to a 2D manifold primarily captures magnitude and distance information and is suited for an approximate data reconstruction. The data reconstruction improves with an increasing number of bottleneck nodes of up to three and four, but it saturates if more nodes are added to the bottleneck.
Detailed organic geochemical and carbon isotopic (delta C-13 and Delta C-14) analyses are performed on permafrost deposits affected by coastal erosion (Herschel Island, Canadian Beaufort Sea) and adjacent marine sediments (Herschel Basin) to understand the fate of organic carbon in Arctic nearshore environments. We use an end-member model based on the carbon isotopic composition of bulk organic matter to identify sources of organic carbon. Monte Carlo simulations are applied to quantify the contribution of coastal permafrost erosion to the sedimentary carbon budget. The models suggest that similar to 40% of all carbon released by local coastal permafrost erosion is efficiently trapped and sequestered in the nearshore zone. This highlights the importance of sedimentary traps in environments such as basins, lagoons, troughs, and canyons for the carbon sequestration in previously poorly investigated, nearshore areas.
Plain Language Summary Increasing air and sea surface temperatures at high latitudes leads to accelerated thaw, destabilization, and erosion of perennially frozen soils (i.e., permafrost), which are often rich in organic carbon. Coastal erosion leads to an increased mobilization of organic carbon into the Arctic Ocean, which there can be converted into greenhouse gases and may therefore contribute to further warming. Carbon decomposition can be limited if organic matter is efficiently deposited on the seafloor, buried in marine sediments, and thus removed from the short-term carbon cycle. Basins, canyons, and troughs near the coastline can serve as sediment traps and potentially accommodate large quantities of organic carbon along the Arctic coast. Here we use biomarkers (source-specific molecules), stable carbon isotopes, and radiocarbon to identify the sources of organic carbon in the nearshore zone of the southern Canadian Beaufort Sea near Herschel Island. We quantify the contribution of coastal permafrost erosion to the sedimentary carbon budget of the area and estimate that more than a third of all carbon released by local permafrost erosion is efficiently trapped in marine sediments. This highlights the importance of regional sediment traps for carbon sequestration.
The current awareness of the high importance of urban green leads to a stronger need for tools to comprehensively represent urban green and its benefits. A common scientific approach is the development of urban ecosystem services (UES) based on remote sensing methods at the city or district level. Urban planning, however, requires fine-grained data that match local management practices. Hence, this study linked local biotope and tree mapping methods to the concept of ecosystem services. The methodology was tested in an inner-city district in SW Germany, comparing publicly accessible areas and non-accessible courtyards. The results provide area-specific [m(2)] information on the green inventory at the microscale, whereas derived stock and UES indicators form the basis for comparative analyses regarding climate adaptation and biodiversity. In the case study, there are ten times more micro-scale green spaces in private courtyards than in the public space, as well as twice as many trees. The approach transfers a scientific concept into municipal planning practice, enables the quantitative assessment of urban green at the microscale and illustrates the importance for green stock data in private areas to enhance decision support in urban development. Different aspects concerning data collection and data availability are critically discussed.
Sulfate reduction is the quantitatively most important process to degrade organic matter in anoxic marine sediment and has been studied intensively in a variety of settings. Guaymas Basin, a young marginal ocean basin, offers the unique opportunity to study sulfate reduction in an environment characterized by organic-rich sediment, high sedimentation rates, and high geothermal gradients (100-958 degrees C km(-1)). We measured sulfate reduction rates (SRR) in samples taken during the International Ocean Discovery Program (IODP) Expedition 385 using incubation experiments with radiolabeled (SO42-)-S-35 carried out at in situ pressure and temperature. The highest SRR (387 nmol cm(-3) d(-1)) was recorded in near-surface sediments from Site U1548C, which had the steepest geothermal gradient (958 degrees C km(-1)). At this site, SRR were generally over an order of magnitude higher than at similar depths at other sites (e.g., 387-157 nmol cm(-3) d(-1) at 1.9 mbsf from Site U1548C vs. 46-1.0 nmol cm(-3) d(-1) at 2.1 mbsf from Site U1552B). Site U1546D is characterized by a sill intrusion, but it had already reached thermal equilibrium and SRR were in the same range as nearby Site U1545C, which is minimally affected by sills. The wide temperature range observed at each drill site suggests major shifts in microbial community composition with very different temperature optima but awaits confirmation by molecular biological analyses. At the transition between the mesophilic and thermophilic range around 40 degrees C-60 degrees C, sulfate-reducing activity appears to be decreased, particularly in more oligotrophic settings, but shows a slight recovery at higher temperatures.
We study the source properties of the 2005 Kashmir earthquake and its aftershocks to unravel the seismotectonics of the NW Himalayan syntaxis. The mainshock and larger aftershocks have been simultaneously relocated using phase data. We use back-projection of high-frequency energy from multiple teleseismic arrays to model the spatio-temporal evolution of the mainshock rupture. Our analysis reveal a bilateral rupture, which initially propagated SE and then NW of the epicenter, with an average rupture velocity of similar to 2 km s(-1). The area of maximum energy release is parallel to and bound by the surface rupture. Incorporating rupture propagation and velocity, we model the mainshock as a line source using P- and SH-waveform inversion. Our result confirms that the mainshock occurred on a NE dipping (similar to 35 degrees) fault plane, with centroid depth of similar to 10 km. Integrated source time function show that majority of the energy was released in the first similar to 20 s, and was confined above the hypocenter. From waveform inverted fault dimension and seismic moment, we argue that the mainshock had an additional similar to 25 km blind rupture beyond the NW Himalayan syntaxis. Combining this with findings from previous studies, we conjecture that the blind rupture propagated NW of the syntaxis underneath a weak detachment overlain by infra-Cambrian salt layer, and terminated in a wedge thrust. All moderate-to-large aftershocks, NW of the mainshock rupture, are concentrated at the edge of the blind rupture termination. Source modeling of these aftershocks reveal thrust mechanism with centroid depths of 2-10 km, and fault planes oriented subparallel to the mainshock rupture. To study the influence of mainshock rupture on aftershock occurrence, we compute Coulomb failure stress on aftershock faults. All these aftershocks lie in the positive Coulomb stress change region. This suggest that the aftershocks have been triggered by either co-seismic or post-seismic slip on the mainshock fault.
Geysers are hot springs whose frequency of water eruptions remain poorly understood. We set up a local broadband seismic network for 1 year at Strokkur geyser, Iceland, and developed an unprecedented catalog of 73,466 eruptions. We detected 50,135 single eruptions but find that the geyser is also characterized by sets of up to six eruptions in quick succession. The number of single to sextuple eruptions exponentially decreased, while the mean waiting time after an eruption linearly increased (3.7 to 16.4 min). While secondary eruptions within double to sextuple eruptions have a smaller mean seismic amplitude, the amplitude of the first eruption is comparable for all eruption types. We statistically model the eruption frequency assuming discharges proportional to the eruption multiplicity and a constant probability for subsequent events within a multituple eruption. The waiting time after an eruption is predictable but not the type or amplitude of the next one. <br /> Plain Language Summary Geysers are springs that often erupt in hot water fountains. They erupt more often than volcanoes but are quite similar. Nevertheless, it is poorly understood how often volcanoes and also geysers erupt. We created a list of 73,466 eruption times of Strokkur geyser, Iceland, from 1 year of seismic data. The geyser erupted one to six times in quick succession. We found 50,135 single eruptions but only 1 sextuple eruption, while the mean waiting time increased from 3.7 min after single eruptions to 16.4 min after sextuple eruptions. Mean amplitudes of each eruption type were higher for single eruptions, but all first eruptions in a succession were similar in height. Assuming a constant heat inflow at depth, we can predict the waiting time after an eruption but not the type or amplitude of the next one.
Multiple P-T-d-t paths reveal the evolution of the final Nuna assembly in northeast Australia
(2020)
The final assembly of the Mesoproterozoic supercontinent Nuna was marked by the collision of Laurentia and Australia at 1.60 Ga, which is recorded in the Georgetown Inlier of NE Australia. Here, we decipher the metamorphic evolution of this final Nuna collisional event using petrostructural analysis, major and trace element compositions of key minerals, thermodynamic modelling, and multi-method geochronology. The Georgetown Inlier is characterised by deformed and metamorphosed 1.70-1.62 Ga sedimentary and mafic rocks, which were intruded byc. 1.56 Ga old S-type granites. Garnet Lu-Hf and monazite U-Pb isotopic analyses distinguish two major metamorphic events (M1 atc. 1.60 Ga and M2 atc. 1.55 Ga), which allows at least two composite fabrics to be identified at the regional scale-c. 1.60 Ga S1 (consisting in fabrics S1a and S1b) andc. 1.55 Ga S2 (including fabrics S2a and S2b). Also, three tectono-metamorphic domains are distinguished: (a) the western domain, with S1 defined by low-P(LP) greenschist facies assemblages; (b) the central domain, where S1 fabric is preserved as medium-P(MP) amphibolite facies relicts, and locally as inclusion trails in garnet wrapped by the regionally dominant low-Pamphibolite facies S2 fabric; and (c) the eastern domain dominated by upper amphibolite to granulite facies S2 foliation. In the central domain, 1.60 GaMP-medium-T(MT) metamorphism (M1) developed within the staurolite-garnet stability field, with conditions ranging from 530-550 degrees C at 6-7 kbar (garnet cores) to 620-650 degrees C at 8-9 kbar (garnet rims), and it is associated with S1 fabric. The onset of 1.55 GaLP-high-T(HT) metamorphism (M2) is marked by replacement of staurolite by andalusite (M2a/D2a), which was subsequently pseudomorphed by sillimanite (M2b/D2b) where granite and migmatite are abundant.P-Tconditions ranged from 600 to 680 degrees C and 4-6 kbar for the M2b sillimanite stage. 1.60 Ga garnet relicts within the S2 foliation highlight the progressive obliteration of the S1 fabric by regional S2 in the central zone during peak M2 metamorphism. In the eastern migmatitic complex, partial melting of paragneiss and amphibolite occurred syn- to post-S2, at 730-770 degrees C and 6-8 kbar, and at 750-790 degrees C and 6 kbar, respectively. The pressure-temperature-deformation-time paths reconstructed for the Georgetown Inlier suggest ac. 1.60 Ga M1/D1 event recorded under greenschist facies conditions in the western domain and under medium-Pand medium-Tconditions in the central domain. This event was followed by the regional 1.56-1.54 Ga low-Pand high-Tphase (M2/D2), extensively recorded in the central and eastern domains. Decompression between these two metamorphic events is ascribed to an episode of exhumation. The two-stage evolution supports the previous hypothesis that the Georgetown Inlier preserves continental collisional and subsequent thermal perturbation associated with granite emplacement.
The main Marmara fault (MMF) extends for 150 km through the Sea of Marmara and forms the only portion of the North Anatolian fault zone that has not ruptured in a large event (Mw >7) for the last 250 yr. Accordingly, this portion is potentially a major source contributing to the seismic hazard of the Istanbul region. On 26 September 2019, a sequence of moderate-sized events started along the MMF only 20 km south of Istanbul and were widely felt by the population. The largest three events, 26 September Mw 5.8 (10:59 UTC), 26 September 2019 Mw 4.1 (11:26 UTC), and 20 January 2020 Mw 4.7 were recorded by numerous strong-motion seismic stations and the resulting ground motions were compared to the predicted means resulting from a set of the most recent ground-motion prediction equations (GMPEs). The estimated residuals were used to investigate the spatial variation of ground motion across the Marmara region. Our results show a strong azimuthal trend in ground-motion residuals, which might indicate systematically repeating directivity effects toward the eastern Marmara region.
Diffusive transport and sorption processes of uranium in the Swiss Opalinus Clay were investigated as a function of partial pressure of carbon dioxide pCO(2), varying mineralogy in the facies and associated changes in porewater composition. Simulations were conducted in one-dimensional diffusion models on the 100 m-scale for a time of one million years using a bottom-up approach based on mechanistic surface complexation models as well as cation exchange to quantify sorption. Speciation calculations have shown, uranium is mainly present as U(VI) and must therefore be considered as mobile for in-situ conditions. Uranium migrated up to 26 m in both, the sandy and the carbonate-rich facies, whereas in the shaly facies 16 m was the maximum. The main species was the anionic complex CaUO2(CO3)(3)(2-) . Hence, anion exclusion was taken into account and further reduced the migration distances by 30 %. The concentrations of calcium and carbonates reflected by the set pCO(2) determine speciation and activity of uranium and consequently the sorption behaviour. Our simulation results allow for the first time to prioritize on the far-field scale the governing parameters for diffusion and sorption of uranium and hence outline the sensitivity of the system. Sorption processes are controlled in descending priority by the carbonate and calcium concentrations, pH, pe and the clay mineral content. Therefore, the variation in porewater composition resulting from the heterogeneity of the facies in the Opalinus Clay formation needs to be considered in the assessment of uranium migration in the far field of a potential repository.
Barite scales in geothermal installations are a highly unwanted effect of circulating deep saline fluids. They build up in the reservoir if supersaturated fluids are re-injected, leading to irreversible loss of injectivity. A model is presented for calculating the total expected barite precipitation. To determine the related injectivity decline over time, the spatial precipitation distribution in the subsurface near the injection well is assessed by modelling barite growth kinetics in a radially diverging Darcy flow domain. Flow and reservoir properties as well as fluid chemistry are chosen to represent reservoirs subject to geothermal exploration located in the North German Basin (NGB) and the Upper Rhine Graben (URG) in Germany. Fluids encountered at similar depths are hotter in the URG, while they are more saline in the NGB. The associated scaling amount normalised to flow rate is similar for both regions. The predicted injectivity decline after 10 years, on the other hand, is far greater for the NGB (64%) compared to the URG (24%), due to the temperature- and salinity-dependent precipitation rate. The systems in the NGB are at higher risk. Finally, a lightweight score is developed for approximating the injectivity loss using the Damkohler number, flow rate and total barite scaling potential. This formula can be easily applied to geothermal installations without running complex reactive transport simulations.
Despite more than half a century of hominin fossil discoveries in eastern Africa, the regional environmental context of hominin evolution and dispersal is not well established due to the lack of continuous palaeoenvironmental records from one of the proven habitats of early human populations, particularly for the Pleistocene epoch. Here we present a 620,000-year environmental record from Chew Bahir, southern Ethiopia, which is proximal to key fossil sites. Our record documents the potential influence of different episodes of climatic variability on hominin biological and cultural transformation. The appearance of high anatomical diversity in hominin groups coincides with long-lasting and relatively stable humid conditions from similar to 620,000 to 275,000 years bp (episodes 1-6), interrupted by several abrupt and extreme hydroclimate perturbations. A pattern of pronounced climatic cyclicity transformed habitats during episodes 7-9 (similar to 275,000-60,000 years bp), a crucial phase encompassing the gradual transition from Acheulean to Middle Stone Age technologies, the emergence of Homo sapiens in eastern Africa and key human social and cultural innovations. Those accumulative innovations plus the alignment of humid pulses between northeastern Africa and the eastern Mediterranean during high-frequency climate oscillations of episodes 10-12 (similar to 60,000-10,000 years bp) could have facilitated the global dispersal of H. sapiens.
Atmospheric dynamics of extreme discharge events from 1979 to 2016 in the southern Central Andes
(2020)
During the South-American Monsoon season, deep convective systems occur at the eastern flank of the Central Andes leading to heavy rainfall and flooding. We investigate the large- and meso-scale atmospheric dynamics associated with extreme discharge events (> 99.9th percentile) observed in two major river catchments meridionally stretching from humid to semi-arid conditions in the southern Central Andes. Based on daily gauge time series and ERA-Interim reanalysis, we made the following three key observations: (1) for the period 1940-2016 daily discharge exhibits more pronounced variability in the southern, semi-arid than in the northern, humid catchments. This is due to a smaller ratio of discharge magnitudes between intermediate (0.2 year return period) and rare events (20 year return period) in the semi-arid compared to the humid areas; (2) The climatological composites of the 40 largest discharge events showed characteristic atmospheric features of cold surges based on 5-day time-lagged sequences of geopotential height at different levels in the troposphere; (3) A subjective classification revealed that 80% of the 40 largest discharge events are mainly associated with the north-northeastward migration of frontal systems and 2/3 of these are cold fronts, i.e. cold surges. This work highlights the importance of cold surges and their related atmospheric processes for the generation of heavy rainfall events and floods in the southern Central Andes.
This study, based on new and high quality in situ observations, quantifies for the first time, the individual contributions of light-absorbing aerosols (black carbon (BC), brown carbon (BrC) and dust) to aerosol absorption over the Indo-Gangetic Plain (IGP) and the Himalayan foothill region, a relatively poorly studied region with several sensitive ecosystems of global importance, as well as highly vulnerable populations. The annual and seasonal average single scattering albedo (SSA) over Kathmandu is the lowest of all the locations. The SSA over Kathmandu is < 0.89 during all seasons, which confirms the dominance of light-absorbing carbonaceous aerosols from local and regional sources over Kathmandu. It is observed here that the SSA decreases with increasing elevation, confirming the dominance of light absorbing carbonaceous aerosols at higher elevations. In contrast, the SSA over the IGP does not exhibit a pronounced spatial variation. BC dominates (>= 75%) the aerosol absorption over the IGP and the Himalayan foothills throughout the year. Higher BC concentration at elevated locations in the Himalayas leads to lower SSA at elevated locations in the Himalayas. The contribution of dust to aerosol absorption is higher throughout the year over the IGP than over the Himalayan foothills. The aerosol absorption over South Asia is very high, exceeding available observations over East Asia, and also exceeds previous model estimates. This quantification will be valuable as observational constraints to help improve regional simulations of climate change, impacts on the glaciers and the hydrological cycle, and will help to direct the focus towards BC as the main contributor to aerosol-induced warming in the region.
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.
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.
Mountain ranges can fundamentally influence the physical and and chemical processes that shape Earths’ surface. With elevations of up to several kilometers they create climatic enclaves by interacting with atmospheric circulation and hydrologic systems, thus leading to a specific distribution of flora and fauna. As a result, the interiors of many Cenozoic mountain ranges are characterized by an arid climate, internally drained and sediment-filled basins, as well as unique ecosystems that are isolated from the adjacent humid, low-elevation regions along their flanks and forelands. These high-altitude interiors of orogens are often characterized by low relief and coalesced sedimentary basins, commonly referred to as plateaus, tectono-geomorphic entities that result from the complex interactions between mantle-driven geological and tectonic conditions and superposed atmospheric and hydrological processes. The efficiency of these processes and the fate of orogenic plateaus is therefore closely tied to the balance of constructive and destructive processes – tectonic uplift and erosion, respectively. In numerous geological studies it has been shown that mountain ranges are delicate systems that can be obliterated by an imbalance of these underlying forces. As such, Cenozoic mountain ranges might not persist on long geological timescales and will be destroyed by erosion or tectonic collapse. Advancing headward erosion of river systems that drain the flanks of the orogen may ultimately sever the internal drainage conditions and the maintenance of storage of sediments within the plateau, leading to destruction of plateau morphology and connectivity with the foreland. Orogenic collapse may be associated with the changeover from a compressional stress field with regional shortening and topographic growth, to a tensional stress field with regional extensional deformation and ensuing incision of the plateau. While the latter case is well-expressed by active extensional faults in the interior parts of the Tibetan Plateau and the Himalaya, for example, the former has been attributed to have breached the internally drained areas of the high-elevation sectors of the Iranian Plateau.
In the case of the Andes of South America and their internally drained Altiplano-Puna Plateau, signs of both processes have been previously described. However, in the orogenic collapse scenario the nature of the extensional structures had been primarily investigated in the northern and southern terminations of the plateau; in some cases, the extensional faults were even regarded to be inactive. After a shallow earthquake in 2020 within the Eastern Cordillera of Argentina that was associated with extensional deformation, the state of active deformation and the character of the stress field in the central parts of the plateau received renewed interest to explain a series of extensional structures in the northernmost sectors of the plateau in north-western Argentina. This study addresses (1) the issue of tectonic orogenic collapse of the Andes and the destruction of plateau morphology by studying the fill and erosion history of the central eastern Andean Plateau using sedimentological and geochronological data and (2) the kinematics, timing and magnitude of extensional structures that form well-expressed fault scarps in sediments of the regional San Juan del Oro surface, which is an integral part of the Andean Plateau and adjacent morphotectonic provinces to the east.
Importantly, sediment properties and depositional ages document that the San Juan del Oro Surface was not part of the internally-drained Andean Plateau, but rather associated with a foreland-directed drainage system, which was modified by the Andean orogeny and that became successively incorporated into the orogen by the eastward-migration of the Andean deformation front during late Miocene – Pliocene time. Structural and geomorphic observations within the plateau indicate that extensional processes must have been repeatedly active between the late Miocene and Holocene supporting the notion of plateau-wide extensional processes, potentially associated with Mw ~ 7 earthquakes. The close relationship between extensional joints and fault orientations underscores that 3 was oriented horizontally in NW-SE direction and 1 was vertical. This unambiguously documents that the observed deformation is related to gravitational forces that drive the orogenic collapse of the plateau. Applied geochronological analyses suggest that normal faulting in the northern Puna was active at about 3 Ma, based on paired cosmogenic nuclide dating of sediment fill units. Possibly due to regional normal faulting the drainage system within the plateau was modified, promoting fluvial incision.
The SiDroForest (Siberian drone-mapped forest inventory) data collection is an attempt to remedy the scarcity of forest structure data in the circumboreal region by providing adjusted and labeled tree-level and vegetation plot-level data for machine learning and upscaling purposes. We present datasets of vegetation composition and tree and plot level forest structure for two important vegetation transition zones in Siberia, Russia; the summergreen-evergreen transition zone in Central Yakutia and the tundra-taiga transition zone in Chukotka (NE Siberia). The SiDroForest data collection consists of four datasets that contain different complementary data types that together support in-depth analyses from different perspectives of Siberian Forest plot data for multi-purpose applications. i. Dataset 1 provides unmanned aerial vehicle (UAV)-borne data products covering the vegetation plots surveyed during fieldwork (Kruse et al., 2021, ). The dataset includes structure-from-motion (SfM) point clouds and red-green-blue (RGB) and red-green-near-infrared (RGN) orthomosaics. From the orthomosaics, point-cloud products were created such as the digital elevation model (DEM), canopy height model (CHM), digital surface model (DSM) and the digital terrain model (DTM). The point-cloud products provide information on the three-dimensional (3D) structure of the forest at each plot. Dataset 2 contains spatial data in the form of point and polygon shapefiles of 872 individually labeled trees and shrubs that were recorded during fieldwork at the same vegetation plots (van Geffen et al., 2021c, ). The dataset contains information on tree height, crown diameter, and species type. These tree and shrub individually labeled point and polygon shapefiles were generated on top of the RGB UVA orthoimages. The individual tree information collected during the expedition such as tree height, crown diameter, and vitality are provided in table format. This dataset can be used to link individual information on trees to the location of the specific tree in the SfM point clouds, providing for example, opportunity to validate the extracted tree height from the first dataset. The dataset provides unique insights into the current state of individual trees and shrubs and allows for monitoring the effects of climate change on these individuals in the future. Dataset 3 contains a synthesis of 10 000 generated images and masks that have the tree crowns of two species of larch ( and ) automatically extracted from the RGB UAV images in the common objects in context (COCO) format (van Geffen et al., 2021a, ). As machine-learning algorithms need a large dataset to train on, the synthetic dataset was specifically created to be used for machine-learning algorithms to detect Siberian larch species. Larix gmeliniiLarix cajanderiDataset 4 contains Sentinel-2 (S-2) Level-2 bottom-of-atmosphere processed labeled image patches with seasonal information and annotated vegetation categories covering the vegetation plots (van Geffen et al., 2021b, ). The dataset is created with the aim of providing a small ready-to-use validation and training dataset to be used in various vegetation-related machine-learning tasks. It enhances the data collection as it allows classification of a larger area with the provided vegetation classes. The SiDroForest data collection serves a variety of user communities. <br /> The detailed vegetation cover and structure information in the first two datasets are of use for ecological applications, on one hand for summergreen and evergreen needle-leaf forests and also for tundra-taiga ecotones. Datasets 1 and 2 further support the generation and validation of land cover remote-sensing products in radar and optical remote sensing. In addition to providing information on forest structure and vegetation composition of the vegetation plots, the third and fourth datasets are prepared as training and validation data for machine-learning purposes. For example, the synthetic tree-crown dataset is generated from the raw UAV images and optimized to be used in neural networks. Furthermore, the fourth SiDroForest dataset contains S-2 labeled image patches processed to a high standard that provide training data on vegetation class categories for machine-learning classification with JavaScript Object Notation (JSON) labels provided. The SiDroForest data collection adds unique insights into remote hard-to-reach circumboreal forest regions.
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.
Tectonically active coasts are dynamic environments characterized by the presence of multiple marine terraces formed by the combined effects of wave erosion, tectonic uplift, and sea-level oscillations at glacialcycle timescales. Well-preserved erosional terraces from the last interglacial sea-level highstand are ideal marker horizons for reconstructing past sea-level positions and calculating vertical displacement rates. We carried out an almost continuous mapping of the last interglacial marine terrace along similar to 5000 km of the western coast of South America between 1 degrees N and 40 degrees S. We used quantitatively replicable approaches constrained by published terrace-age estimates to ultimately compare elevations and patterns of uplifted terraces with tectonic and climatic parameters in order to evaluate the controlling mechanisms for the formation and preservation of marine terraces and crustal deformation. Uncertainties were estimated on the basis of measurement errors and the distance from referencing points. Overall, our results indicate a median elevation of 30.1 m, which would imply a median uplift rate of 0.22 m kyr(-1) averaged over the past similar to 125 kyr. The patterns of terrace elevation and uplift rate display high-amplitude (similar to 100-200 m) and long-wavelength (similar to 10(2) km) structures at the Manta Peninsula (Ecuador), the San Juan de Marcona area (central Peru), and the Arauco Peninsula (south-central Chile). Medium-wavelength structures occur at the Mejillones Peninsula and Topocalma in Chile, while short-wavelength (< 10 km) features are for instance located near Los Vilos, Valparaiso, and Carranza, Chile. We interpret the long-wavelength deformation to be controlled by deep-seated processes at the plate interface such as the subduction of major bathymetric anomalies like the Nazca and Carnegie ridges. In contrast, short-wavelength deformation may be primarily controlled by sources in the upper plate such as crustal faulting, which, however, may also be associated with the subduction of topographically less pronounced bathymetric anomalies. Latitudinal differences in climate additionally control the formation and preservation of marine terraces. Based on our synopsis we propose that increasing wave height and tidal range result in enhanced erosion and morphologically well-defined marine terraces in south-central Chile. Our study emphasizes the importance of using systematic measurements and uniform, quantitative methodologies to characterize and correctly interpret marine terraces at regional scales, especially if they are used to unravel the tectonic and climatic forcing mechanisms of their formation. This database is an integral part of the World Atlas of Last Interglacial Shorelines (WALIS), published online at https://doi.org/10.5281/zenodo.4309748 (Freisleben et al., 2020).
We analyze the spatiotemporal evolution of seismicity during a sequence of moderate (an M-w 4.7 foreshock and M-w 5.8 mainshock) earthquakes occurring in September 2019 at the transition between a creeping and a locked segment of the North Anatolian fault in the central Sea of Marmara, northwest Turkey. To investigate in detail the seismicity evolution, we apply a matched-filter technique to continuous waveforms, thus reducing the magnitude threshold for detection. Sequences of foreshocks preceding the two largest events are clearly seen, exhibiting two different behaviors: a long-term activation of the seismicity along the entire fault segment and a short-term concentration around the epicenters of the large events. We suggest a two-scale preparation phase, with aseismic slip preparing the mainshock final rupture a few days before, and a cascade mechanism leading to the nucleation of the mainshock. Thus, our study shows a combination of seismic and aseismic slip during the foreshock sequence changing the strength of the fault, bringing it closer to failure.
The geometry of carbonate platforms reflects the interaction of several factors. However, the impact of carbonate-producing organisms has been poorly investigated so far. This study applies stratigraphic forward modelling (SFM) and sensitivity analysis to examine, referenced to the Miocene Llucmajor Platform, the effect of changes of dominant biotic production in the oligophotic and euphotic zones on platform geometry. Our results show that the complex interplay of carbonate production rates, bathymetry and variations in accommodation space control the platform geometry. The main driver of progradation is the oligophotic production of rhodalgal sediments during the lowstands. This study demonstrates that platform geometry and internal architecture varies significantly according to the interaction of the predominant carbonate-producing biotas. The input parameters for this study are based on well-understood Miocene carbonate biotas with characteristic euphotic, oligophotic and photo-independent carbonate production in which it is crucial that each carbonate-producing class is modelled explicitly within the simulation run and not averaged with a single carbonate production-depth profile. This is important in subsurface exploration studies based on stratigraphic forward models where the overall platform geometry may be approximated through calibration runs, and constrained by seismic surveys and wellbores. However, the internal architecture is likely to be oversimplified without an in-depth understanding of the target carbonate system and a transfer to forward modelling parameters.
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.
The LArge-scale Reservoir Simulator (LARS) has been previously developed to study hydrate dissociation in hydrate-bearing systems under in-situ conditions. In the present study, a numerical framework of equations of state describing hydrate formation at equilibrium conditions has been elaborated and integrated with a numerical flow and transport simulator to investigate a multi-stage hydrate formation experiment undertaken in LARS. A verification of the implemented modeling framework has been carried out by benchmarking it against another established numerical code. Three-dimensional (3D) model calibration has been performed based on laboratory data available from temperature sensors, fluid sampling, and electrical resistivity tomography. The simulation results demonstrate that temperature profiles, spatial hydrate distribution, and bulk hydrate saturation are consistent with the observations. Furthermore, our numerical framework can be applied to calibrate geophysical measurements, optimize post-processing workflows for monitoring data, improve the design of hydrate formation experiments, and investigate the temporal evolution of sub-permafrost methane hydrate reservoirs.
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>
El plateau Andino es el segundo plateau orogénico más grande del mundo y se ubica en los Andes Centrales, desarrollado en un sistema orogénico no colisional. Se extiende desde el sur del Perú (15°S), hasta el norte de Argentina y Chile (27°30´S). A partir de los 24°S y prologándose hacia el sur, el plateau Andino se denomina Puna y está caracterizado por un sistema de cuencas endorreicas y salares delimitados por cordones montañosos. Entre los 26° y 27°30´S, la Puna encuentra su límite austral en una zona de transición entre una zona de subducción normal y una zona de subducción plana o “flat slab” que se prolonga hasta los 33°S. Diversos estudios documentan la ocurrencia de un aumento del espesor cortical, y levantamiento episódico y diacrónico del relieve, alcanzando su configuración actual durante el Mioceno tardío. Posteriormente, el plateau habría experimentado un cambio en el estilo de deformación dominado por procesos extensionales evidenciado por fallas y terremotos de cinemática normal. Sin embargo, en el borde sur del plateau de la Puna y en las áreas delimitadas con el resto del orógeno, la variación del campo de esfuerzo no está del todo comprendida, reflejando una excelente oportunidad para evaluar cómo el campo de esfuerzo puede evolucionar durante el desarrollo del orógeno y cómo puede verse afectado por la presencia/ausencia de un plateau orogénico, así como también por la existencia de anisotropías estructurales propias de cada unidad morfotectónica.
Esta Tesis investiga la relación entre la deformación cortical somera y la evolución en tiempo y espacio del campo de esfuerzos en el sector sur del plateau Andino, durante el cenozoico tardío. Para realizar esta investigación, se utilizaron técnicas de obtención de edades radiométricas con el método Uranio-Plomo (U-Pb), análisis de fallas mesoscópicas para la obtención de tensores de esfuerzos y delimitación de la orientación de los ejes principales de esfuerzos, análisis de anisotropía de susceptibilidad magnética en rocas sedimentarias y volcanoclásticas para estimar direcciones de acortamiento o direcciones de transporte sedimentario, técnicas de modelado cinemático para llegar a una aproximación de las estructuras corticales profundas asociadas a la deformación allí registrada, y un análisis morfométrico para la identificación de indicadores geomorfológicos asociados a deformación producto de la actividad tectónica cuaternaria.
Combinando estos resultados con los antecedentes previamente documentados, el estudio revela una compleja variación del campo de esfuerzo caracterizado por cambios en la orientación y permutaciones verticales de los ejes principales de esfuerzos, durante cada régimen de deformación, durante los últimos ~24 Ma. La evolución del campo de esfuerzos puede ser asociada temporalmente a tres fases orogénicas involucradas con la evolución de los Andes Centrales en esta latitud: (1) una primera fase con un régimen de esfuerzos compresivos de acortamiento E-O documentado desde el Eoceno, Oligoceno tardío hasta el Mioceno medio en el área, coincide con la fase de construcción andina, engrosamiento y crecimiento de la corteza y levantamiento topográfico; (2) una segunda fase caracterizada por un régimen de esfuerzos de transcurrencia, a partir de los ~11 Ma en el borde occidental y compresión y transcurrencia a los~5 Ma en el borde oriental del plateau de la Puna, y un régimen de esfuerzo compresivos en Famatina y las Sierras Pampeanas interpretado como una transición entre la construcción orogénica del Neógeno y la máxima acumulación de deformación y el alzamiento topográfico del plateau de la Puna, y (3) una tercera fase donde el régimen se caracteriza por la transcurrencia en la Puna y en su borde occidental y en su borde oriental con las Sierras Pampeanas, después de ~5-4 Ma, interpretado como un régimen de esfuerzos controlados por el engrosamiento cortical desarrollado a lo largo del borde sur del plateau Altiplano/Puna, previo a un colapso orogénico. Los resultados dejan en evidencia que el borde del plateau experimentó el paso desde un régimen compresivo hacia uno transcurrente, que se diferencia de la extensión documentada hacia el norte en el plateau Andino para el mismo período. Cambios en los esfuerzos similares han sido documentado durante la construcción del plateau Tibetano, en donde un régimen de esfuerzo predominantemente compresivo cambió a un régimen de transcurrente cuando el plateau habría alcanzado la mitad de su elevación actual, y que posteriormente derivó en un régimen extensional, entre 14 y 4 Ma, cuando la altitud del plateau fue superior al 80% respecto a su actitud actual, lo que podría estar indicando que los regímenes transcurrentes representan etapas transicionales entre las zonas externas del plateau bajo compresión y las zonas internas, en las que los regímenes extensionales son más viables de ocurrir.
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 creation of building exposure models for seismic risk assessment is frequently challenging due to the lack of availability of detailed information on building structures. Different strategies have been developed in recent years to overcome this, including the use of census data, remote sensing imagery and volunteered graphic information (VGI). This paper presents the development of a building-by-building exposure model based exclusively on openly available datasets, including both VGI and census statistics, which are defined at different levels of spatial resolution and for different moments in time. The initial model stemming purely from building-level data is enriched with statistics aggregated at the neighbourhood and city level by means of a Monte Carlo simulation that enables the generation of full realisations of damage estimates when using the exposure model in the context of an earthquake scenario calculation. Though applicable to any other region of interest where analogous datasets are available, the workflow and approach followed are explained by focusing on the case of the German city of Cologne, for which a scenario earthquake is defined and the potential damage is calculated. The resulting exposure model and damage estimates are presented, and it is shown that the latter are broadly consistent with damage data from the 1978 Albstadt earthquake, notwithstanding the differences in the scenario. Through this real-world application we demonstrate the potential of VGI and open data to be used for exposure modelling for natural risk assessment, when combined with suitable knowledge on building fragility and accounting for the inherent uncertainties.
Neoarchean (similar to 2.73-2.70 Ga) accretionary history of the eastern Dharwar Craton, India
(2022)
Cratonic mid-crustal plutons may contain supracrustal enclaves that preserve evidence of an earlier growth history. The Eastern Dharwar craton records Neoarchean two-stage accretionary sequential growth (2.70 and 2.55 Ga) and a chronology of their enclaves could refine orogenic models. To test whether the metamorphic history of their enclaves was related to any of these stages, phase equilibria modelling and combined Lu-Hf and Sm-Nd geochronology on garnet were conducted on metapsammite, now preserved as garnet-orthopyroxene-cordierite gneiss. Phase equilibria modelling indicates peak metamorphic conditions, similar to 850 degrees C and similar to 8.5 kbar (M1a), were followed by near isothermal decompression to 5-6 kbar (M1b) and isobaric cooling to similar to 800 degrees C (M1c). The thermobaric gradient related to peak metamorphic conditions, similar to 30 degrees C kbar(-1), is typical of collisional orogens. Regression of the whole-rock and garnet, for sample S17b, yield Lu-Hf isochron ages of 2733 +/- 29 Ma, and for sample S18, 2724 +/- 13 Ma. A Lu-Hf weighted mean age for the porphyroblastic garnet suggests growth at 2725.5 +/- 11.9 Ma during the M1a-M1b stages. In contrast, the whole-rock sample S17b and the garnet fractions yield a Sm-Nd isochron age of 2696 +/- 10 Ma. From sample S18 the whole rock, garnet fractions, and orthopyroxene yield an isochron age of 2683 +/- 15 Ma. The garnet Sm-Nd weighted mean age at 2692.0 +/- 8.3 Ma constrains the M1b-M1c stages. We suggest that the protoliths to these supracrustal enclaves were deposited in an arc tectonic setting and underwent thickening followed by heating during peeled-back lithospheric convergence. Therefore, the earliest of the craton-forming accretionary stages is preserved as the similar to 2.73 Ga granulite-facies enclaves, marginally older than the 2.70-2.65 Ga cratonic greenstone volcanism. Tectonic exhumation of these mid-crustal granulite enclaves was in response to the late-Proterozoic (similar to 1.7 Ga) Bhopalpatnam orogeny.
Magmatic continental rifts often constitute nascent plate boundaries, yet long-term extension rates and transient rate changes associated with these early stages of continental breakup remain difficult to determine. Here, we derive a time-averaged minimum extension rate for the inner graben of the Northern Kenya Rift (NKR) of the East African Rift System for the last 0.5 m.y. We use the TanDEM-X science digital elevation model to evaluate fault-scarp geometries and determine fault throws across the volcano-tectonic axis of the inner graben of the NKR. Along rift-perpendicular profiles, amounts of cumulative extension are determined, and by integrating four new Ar-40/Ar-39 radiometric dates for the Silali volcano into the existing geochronology of the faulted volcanic units, time-averaged extension rates are calculated. This study reveals that in the inner graben of the NKR, the long-term extension rate based on mid-Pleistocene to recent brittle deformation has minimum values of 1.0-1.6 mm yr(-1), locally with values up to 2.0 mm yr(-1). A comparison with the decadal, geodetically determined extension rate reveals that at least 65% of the extension must be accommodated within a narrow, 20-km-wide zone of the inner rift. In light of virtually inactive border faults of the NKR, we show that extension is focused in the region of the active volcano-tectonic axis in the inner graben, thus highlighting the maturing of continental rifting in the NKR.
Here I present a comparison between two of the most widely used reduced-complexity models for the representation of sediment transport and deposition processes, namely the transport-limited (or TL) model and the under-capacity (or xi-q) model more recently developed by Davy and Lague (2009). Using both models, I investigate the behavior of a sedimentary continental system of length L fed by a fixed sedimentary flux from a catchment of size A(0) in a nearby active orogen through which sediments transit to a fixed base level representing a large river, a lake or an ocean. This comparison shows that the two models share the same steady-state solution, for which I derive a simple 1D analytical expression that reproduces the major features of such sedimentary systems: a steep fan that connects to a shallower alluvial plain. The resulting fan geometry obeys basic observational constraints on fan size and slope with respect to the upstream drainage area, A(0). The solution is strongly dependent on the size of the system, L, in comparison to a distance L-0, which is determined by the size of A(0), and gives rise to two fundamentally different types of sedimentary systems: a constrained system where L < L-0 and open systems where L > L-0. I derive simple expressions that show the dependence of the system response time on the system characteristics, such as its length, the size of the upstream catchment area, the amplitude of the incoming sedimentary flux and the respective rate parameters (diffusivity or erodibility) for each of the two models. I show that the xi-q model predicts longer response times. I demonstrate that although the manner in which signals propagates through the sedimentary system differs greatly between the two models, they both predict that perturbations that last longer than the response time of the system can be recorded in the stratigraphy of the sedimentary system and in particular of the fan. Interestingly, the xi-q model predicts that all perturbations in the incoming sedimentary flux will be transmitted through the system, whereas the TL model predicts that rapid perturbations cannot. I finally discuss why and under which conditions these differences are important and propose observational ways to determine which of the two models is most appropriate to represent natural systems.
Earthquake modeling is the key to a profound understanding of a rupture. Its kinematics or dynamics are derived from advanced rupture models that allow, for example, to reconstruct the direction and velocity of the rupture front or the evolving slip distribution behind the rupture front. Such models are often parameterized by a lattice of interacting sub-faults with many degrees of freedom, where, for example, the time history of the slip and rake on each sub-fault are inverted. To avoid overfitting or other numerical instabilities during a finite-fault estimation, most models are stabilized by geometric rather than physical constraints such as smoothing.
As a basis for the inversion approach of this study, we build on a new pseudo-dynamic rupture model (PDR) with only a few free parameters and a simple geometry as a physics-based solution of an earthquake rupture. The PDR derives the instantaneous slip from a given stress drop on the fault plane, with boundary conditions on the developing crack surface guaranteed at all times via a boundary element approach. As a side product, the source time function on each point on the rupture plane is not constraint and develops by itself without additional parametrization. The code was made publicly available as part of the Pyrocko and Grond Python packages. The approach was compared with conventional modeling for different earthquakes. For example, for the Mw 7.1 2016 Kumamoto, Japan, earthquake, the effects of geometric changes in the rupture surface on the slip and slip rate distributions could be reproduced by simply projecting stress vectors. For the Mw 7.5 2018 Palu, Indonesia, strike-slip earthquake, we also modelled rupture propagation using the 2D Eikonal equation and assuming a linear relationship between rupture and shear wave velocity. This allowed us to give a deeper and faster propagating rupture front and the resulting upward refraction as a new possible explanation for the apparent supershear observed at the Earth's surface.
The thesis investigates three aspects of earthquake inversion using PDR: (1) to test whether implementing a simplified rupture model with few parameters into a probabilistic Bayesian scheme without constraining geometric parameters is feasible, and whether this leads to fast and robust results that can be used for subsequent fast information systems (e.g., ground motion predictions). (2) To investigate whether combining broadband and strong-motion seismic records together with near-field ground deformation data improves the reliability of estimated rupture models in a Bayesian inversion. (3) To investigate whether a complex rupture can be represented by the inversion of multiple PDR sources and for what type of earthquakes this is recommended.
I developed the PDR inversion approach and applied the joint data inversions to two seismic sequences in different tectonic settings. Using multiple frequency bands and a multiple source inversion approach, I captured the multi-modal behaviour of the Mw 8.2 2021 South Sandwich subduction earthquake with a large, curved and slow rupturing shallow earthquake bounded by two faster and deeper smaller events. I could cross-validate the results with other methods, i.e., P-wave energy back-projection, a clustering analysis of aftershocks and a simple tsunami forward model.
The joint analysis of ground deformation and seismic data within a multiple source inversion also shed light on an earthquake triplet, which occurred in July 2022 in SE Iran. From the inversion and aftershock relocalization, I found indications for a vertical separation between the shallower mainshocks within the sedimentary cover and deeper aftershocks at the sediment-basement interface. The vertical offset could be caused by the ductile response of the evident salt layer to stress perturbations from the mainshocks.
The applications highlight the versatility of the simple PDR in probabilistic seismic source inversion capturing features of rather different, complex earthquakes. Limitations, as the evident focus on the major slip patches of the rupture are discussed as well as differences to other finite fault modeling methods.
We examined the relationship between the mechanical strength of the lithosphere and the distribution of seismicity within the overriding continental plate of the southern Central Andes (SCA, 29 degrees-39 degrees S), where the oceanic Nazca Plate changes its subduction angle between 33 degrees S and 35 degrees S, from subhorizontal in the north (<5 degrees) to steep in the south (similar to 30 degrees). We computed the long-term lithospheric strength based on an existing 3D model describing variations in thickness, density, and temperature of the main geological units forming the lithosphere of the SCA and adjacent forearc and foreland regions. The comparison between our results and seismicity within the overriding plate (upper-plate seismicity) shows that most of the events occur within the modeled brittle domain of the lithosphere. The depth where the deformation mode switches from brittle frictional to thermally activated ductile creep provides a conservative lower bound to the seismogenic zone in the overriding plate of the study area. We also found that the majority of upper-plate earthquakes occurs within the realm of first-order contrasts in integrated strength (12.7-13.3 log Pam in the Andean orogen vs. 13.5-13.9 log Pam in the forearc and the foreland). Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the effects of slab steepening.
Draft Genome Sequence of Nocardioides alcanivorans NGK65(T), a Hexadecane-Degrading Bacterium
(2022)
The Gram-positive bacterium Nocardioides alcanivorans NGK65(T) was isolated from plastic-polluted soil and cultivated on medium with polyethylene as the single carbon source. Nanopore sequencing revealed the presence of candidate enzymes for the biodegradation of polyethylene. Here, we report the draft genome of this newly described member of the terrestrial plastisphere.
Data recorded by distributed acoustic sensing (DAS) along an optical fibre sample the spatial and temporal properties of seismic wavefields at high spatial density. Often leading to massive amount of data when collected for seismic monitoring along many kilometre long cables. The spatially coherent signals from weak seismic arrivals within the data are often obscured by incoherent noise. We present a flexible and computationally efficient filtering technique, which makes use of the dense spatial and temporal sampling of the data and that can handle the large amount of data. The presented adaptive frequency-wavenumber filter suppresses the incoherent seismic noise while amplifying the coherent wavefield. We analyse the response of the filter in time and spectral domain, and we demonstrate its performance on a noisy data set that was recorded in a vertical borehole observatory showing active and passive seismic phase arrivals. Lastly, we present a performant open-source software implementation enabling real-time filtering of large DAS data sets.
The deformation style of mountain belts is greatly influenced by the upper plate architecture created during preceding deformation phases. The Mesozoic Salta Rift extensional phase has created a dominant structural and lithological framework that controls Cenozoic deformation and exhumation patterns in the Central Andes. Studying the nature of these pre-existing anisotropies is a key to understanding the spatiotemporal distribution of exhumation and its controlling factors. The Eastern Cordillera in particular, has a structural grain that is in part controlled by Salta Rift structures and their orientation relative to Andean shortening. As a result, there are areas in which Andean deformation prevails and areas where the influence of the Salta Rift is the main control on deformation patterns.
Between 23 and 24°S, lithological and structural heterogeneities imposed by the Lomas de Olmedo sub-basin (Salta Rift basin) affect the development of the Eastern Cordillera fold-and-thrust belt. The inverted northern margin of the sub-basin now forms the southern boundary of the intermontane Cianzo basin. The former western margin of the sub-basin is located at the confluence of the Subandean Zone, the Santa Barbara System and the Eastern Cordillera. Here, the Salta Rift basin architecture is responsible for the distribution of these morphotectonic provinces. In this study we use a multi-method approach consisting of low-temperature (U-Th-Sm)/He and apatite fission track thermochronology, detrital geochronology, structural and sedimentological analyses to investigate the Mesozoic structural inheritance of the Lomas de Olmedo sub-basin and Cenozoic exhumation patterns.
Characterization of the extension-related Tacurú Group as an intermediate succession between Paleozoic basement and the syn-rift infill of the Lomas de Olmedo sub-basin reveals a Jurassic maximum depositional age. Zircon (U-Th-Sm)/He cooling ages record a pre-Cretaceous onset of exhumation for the rift shoulders in the northern part of the sub-basin, whereas the western shoulder shows a more recent onset (140–115 Ma). Variations in the sedimentary thickness of syn- and post-rift strata document the evolution of accommodation space in the sub-basin. While the thickness of syn-rift strata increases rapidly toward the northern basin margin, the post-rift strata thickness decreases toward the margin and forms a condensed section on the rift shoulder.
Inversion of Salta Rift structures commenced between the late Oligocene and Miocene (24–15 Ma) in the ranges surrounding the Cianzo basin. The eastern and western limbs of the Cianzo syncline, located in the hanging wall of the basin-bounding Hornocal fault, show diachronous exhumation. At the same time, western fault blocks of Tilcara Range, south of the Cianzo basin, began exhuming in the late Oligocene to early Miocene (26–16 Ma). Eastward propagation to the frontal thrust and to the Paleozoic strata east of the Tilcara Range occurred in the middle Miocene (22–10 Ma) and the late Miocene–early Pliocene (10–4 Ma), respectively.
During the period 750-600 Ma ago, prior to the final break-up of the supercontinent Rodinia, the crust of both the North American Craton and Baltica was intruded by significant amounts of rift-related magmas originating from the mantle. In the Proterozoic crust of Southern Norway, the 580 Ma old Fen carbonatite-ultramafic complex is a representative of this type of rocks. In this paper, we report the occurrence of an ultramafic lamprophyre dyke which possibly is linked to the Fen complex, although Ar-40/Ar-39 data from phenocrystic phlogopite from the dyke gave an age of 686 +/- 9 Ma. The lamprophyre dyke was recently discovered in one of the Kongsberg silver mines at Vinoren, Norway. Whole rock geochemistry, geochronological and mineralogical data from the ultramafic lamprophyre dyke are presented aiming to elucidate its origin and possible geodynamic setting. From the whole-rock composition of the Vinoren dyke, the rock could be recognized as transitional between carbonatite and kimberlite-II (orangeite). From its diagnostic mineralogy, the rock is classified as aillikite. The compositions and xenocrystic nature of several of the major and accessory minerals from the Vinoren aillikite are characteristic for diamondiferous rocks (kimberlites/lamproites/UML): Phlogopite with kinoshitalite-rich rims, chromite-spinel-ulvospinel series, Mg- and Mn-rich ilmenites, rutile and lucasite-(Ce). We suggest that the aillikite melt formed during partial melting of a MARID (mica-amphibole-rutile-ilmenite-diopside)-like source under CO2 fluxing. The pre-rifting geodynamic setting of the Vinoren aillikite before the Rodinia supercontinent breakup suggests a relatively thick SCLM (Subcontinental Lithospheric Mantle) during this stage and might indicate a diamond-bearing source for the parental melt. This is in contrast to the about 100 Ma younger Fen complex, which were derived from a thin SCLM.
Fifteen N-butylpyridinium salts - five monometallic [C4Py](2)[MBr4] and ten bimetallic [C4Py](2)[(M0.5M0.5Br4)-M-a-Br-b] (M=Co, Cu, Mn, Ni, Zn) - were synthesized, and their structures and thermal and electrochemical properties were studied. All the compounds are ionic liquids (ILs) with melting points between 64 and 101 degrees C. Powder and single-crystal X-ray diffraction show that all ILs are isostructural. The electrochemical stability windows of the ILs are between 2 and 3 V. The conductivities at room temperature are between 10(-5) and 10(-6) S cm(-1). At elevated temperatures, the conductivities reach up to 10(-4) S cm(-1) at 70 degrees C. The structures and properties of the current bromide-based ILs were also compared with those of previous examples using chloride ligands, which illustrated differences and similarities between the two groups of ILs.
The steady increase of ground-motion data not only allows new possibilities but also comes with new challenges in the development of ground-motion models (GMMs). Data classification techniques (e.g., cluster analysis) do not only produce deterministic classifications but also probabilistic classifications (e.g., probabilities for each datum to belong to a given class or cluster). One challenge is the integration of such continuous classification in regressions for GMM development such as the widely used mixed-effects model. We address this issue by introducing an extension of the mixed-effects model to incorporate data weighting. The parameter estimation of the mixed-effects model, that is, fixed-effects coefficients of the GMMs and the random-effects variances, are based on the weighted likelihood function, which also provides analytic uncertainty estimates. The data weighting permits for earthquake classification beyond the classical, expert-driven, binary classification based, for example, on event depth, distance to trench, style of faulting, and fault dip angle. We apply Angular Classification with Expectation-maximization, an algorithm to identify clusters of nodal planes from focal mechanisms to differentiate between, for example, interface- and intraslab-type events. Classification is continuous, that is, no event belongs completely to one class, which is taken into account in the ground-motion modeling. The theoretical framework described in this article allows for a fully automatic calibration of ground-motion models using large databases with automated classification and processing of earthquake and ground-motion data. As an example, we developed a GMM on the basis of the GMM by Montalva et al. (2017) with data from the strong-motion flat file of Bastias and Montalva (2016) with similar to 2400 records from 319 events in the Chilean subduction zone. Our GMM with the data-driven classification is comparable to the expert-classification-based model. Furthermore, the model shows temporal variations of the between-event residuals before and after large earthquakes in the region.
High-pressure experiments were performed to investigate the effectiveness, rate and mechanism of carbonation of serpentinites by a carbon-saturated COH fluid at 1.5-2.5 GPa and 375-700 degrees C. This allows a better understanding of the fate and redistribution of slab-derived carbonic fluids when they react with the partially hydrated mantle within and above the subducting slab under pressure and temperature conditions corresponding to the forearc mantle. Interactions between carbon-saturated CO2-H2O-CH4 fluids and serpentinite were investigated using natural serpentinite cylinders with natural grain sizes and shapes in piston-cylinder experiments. The volatile composition of post-run fluids was quantified by gas chromatography. Solid phases were examined by Raman spectroscopy, electron microscopy and laser ablation inductively coupled plasma mass spectrometry. Textures, porosity and phase abundances of recovered rock cores were visualized and quantified by three-dimensional, high-resolution computed tomography. We find that carbonation of serpentinites is efficient at sequestering CO2 from the interacting fluid into newly formed magnesite. Time-series experiments demonstrate that carbonation is completed within similar to 96 h at 2 GPa and 600 degrees C. With decreasing CO2, aq antigorite is replaced first by magnesite + quartz followed by magnesite + talc + chlorite in distinct, metasomatic fronts. Above antigorite stability magnesite + enstatite + talc + chlorite occur additionally. The formation of fluid-permeable reaction zones enhances the reaction rate and efficiency of carbonation. Carbonation probably occurs via an interface-coupled replacement process, whereby interconnected porosity is present within reaction zones after the experiment. Consequently, carbonation of serpentinites is self-promoting and efficient even if fluid flow is channelized into veins. We conclude that significant amounts of carbonates may accumulate, over time, in the hydrated forearc mantle.
The Mackenzie Delta (MD) is a permafrost-bearing region along the coasts of the Canadian Arctic which exhibits high sub-permafrost gas hydrate (GH) reserves. The GH occurring at the Mallik site in the MD is dominated by thermogenic methane (CH4), which migrated from deep conventional hydrocarbon reservoirs, very likely through the present fault systems. Therefore, it is assumed that fluid flow transports dissolved CH4 upward and out of the deeper overpressurized reservoirs via the existing polygonal fault system and then forms the GH accumulations in the Kugmallit-Mackenzie Bay Sequences. We investigate the feasibility of this mechanism with a thermo-hydraulic-chemical numerical model, representing a cross section of the Mallik site. We present the first simulations that consider permafrost formation and thawing, as well as the formation of GH accumulations sourced from the upward migrating CH4-rich formation fluid. The simulation results show that temperature distribution, as well as the thickness and base of the ice-bearing permafrost are consistent with corresponding field observations. The primary driver for the spatial GH distribution is the permeability of the host sediments. Thus, the hypothesis on GH formation by dissolved CH4 originating from deeper geological reservoirs is successfully validated. Furthermore, our results demonstrate that the permafrost has been substantially heated to 0.8-1.3 degrees C, triggered by the global temperature increase of about 0.44 degrees C and further enhanced by the Arctic Amplification effect at the Mallik site from the early 1970s to the mid-2000s.
Evolution of chemical bonding and spin-pairing energy in ferropericlase across Its spin transition
(2022)
The evolution of chemical bonding in ferropericlase, (Mg,Fe)O, with pressure may affect the physical and chemical properties of the Earth's lower mantle. Here, we report high-pressure optical absorption spectra of single-crystalline ferropericlase ((Mg0.87Fe0.13)O) up to 135 GPa. Combined with a re-evaluation of published partial fluorescence yield X-ray absorption spectroscopy data, we show that the covalency of the Fe-O bond increases with pressure, but the iron spin transition at 57-76.5 GPa reverses this trend. The qualitative crossover in chemical bonding suggests that the spin-pairing transition weakens the Fe-O bond in ferropericlase. We find, that the spin transition in ferropericlase is caused by both the increase of the ligand field-splitting energy and the decrease in the spin-pairing energy of high-spin Fe2+.
Forested areas are assumed not to be influenced by erosion processes. However, forest soils of Northern Germany in a hummocky ground moraine landscape can sometimes exhibit a very shallow thickness on crest positions and buried soils on slope positions. The question consequently is: Are these on-going or ancient erosional and depositional processes? Plutonium isotopes act as soil erosion/deposition tracers for recent (last few decades) processes. Here, we quantified the 239+240PU inventories in a small, forested catchment (ancient forest "Melzower Forst", deciduous trees), which is characterised by a hummocky terrain including a kettle hole. Soil development depths (depth to C horizon) and 239+240PU inventories along a catena of sixteen different profiles were determined and correlated to relief parameters. Moreover, we compared different modelling approaches to derive erosion rates from Pu data. <br /> We find a strong relationship between soil development depths, distance-to-sink and topography along the catena. Fully developed Retisols (thicknesses > 1 m) in the colluvium overlay old land surfaces as documented by fossil Ah horizons. However, we found no relationship of Pu-based erosion rates to any relief parameter. Instead, 239+240PU inventories showed a very high local, spatial variability (36-70 Bq m(-2)). Low annual rainfall, spatially distributed interception and stem flow might explain the high variability of the 239+240PU inventories, giving rise to a patchy input pattern. Different models resulted in quite similar erosion and deposition rates (max: -5 t ha(-1) yr(-1) to +7.3 t ha(-1) yr(-1)). Although some rates are rather high, the magnitude of soil erosion and deposition - in terms of soil thickness change - is negligible during the last 55 years. The partially high values are an effect of the patchy Pu deposition on the forest floor. This forest has been protected for at least 240 years. Therefore rather natural events and anthropogenic activities during medieval times or even earlier must have caused the observed soil pattern, which documents strong erosion and deposition processes.