Refine
Year of publication
Document Type
- Article (2644)
- Doctoral Thesis (409)
- Postprint (138)
- Other (75)
- Review (51)
- Preprint (17)
- Monograph/Edited Volume (13)
- Conference Proceeding (12)
- Habilitation Thesis (11)
- Master's Thesis (6)
Language
- English (3378) (remove)
Keywords
- climate change (50)
- Holocene (44)
- erosion (27)
- permafrost (27)
- Himalaya (26)
- remote sensing (25)
- Climate change (23)
- Tibetan Plateau (22)
- Earthquake source observations (21)
- Pollen (21)
Institute
- Institut für Geowissenschaften (3378) (remove)
Natural gas can be temporarily stored in a variety of underground facilities, such as depleted gas and oil fields, natural aquifers and caverns in salt rocks. Being extensively monitored during operations, these systems provide a favourable opportunity to investigate how pressure varies in time and space and possibly induces/triggers earthquakes on nearby faults. Elaborate and detailed numerical modelling techniques are often applied to study gas reservoirs. Here we show the possibilities and discuss the limitations of a flexible and easily formulated tool that can be straightforwardly applied to simulate temporal pore-pressure variations and study the relation with recorded microseismic events. We use the software POEL (POroELastic diffusion and deformation) which computes the poroelastic response to fluid injection/extraction in a horizontally layered poroelastic structure. We further develop its application to address the presence of vertical impermeable faults bounding the reservoir and of multiple injection/extraction sources. Exploiting available information on the reservoir geometry and physical parameters, and records of injection/extraction rates for a gas reservoir in southern Europe, we perform an extensive parametric study considering different model configurations. Comparing modelled spatiotemporal pore-pressure variations with in situ measurements, we show that the inclusion of vertical impermeable faults provides an improvement in reproducing the observations and results in pore-pressure accumulation near the faults and in a variation of the temporal pore-pressure diffusion pattern. To study the relation between gas storage activity and recorded local microseismicity, we applied different seismicity models based on the estimated porepressure distribution. This analysis helps to understand the spatial distribution of seismicity and its temporal modulation. The results show that the observed microseismicity could be partly linked to the storage activity, but the contribution of tectonic background seismicity cannot be excluded.
Reservoir-triggered seismicity has been observed near dams during construction, impoundment, and cyclic filling in many parts of the earth. In Turkey, the number of dams has increased substantially over the last decade, with Ataturk Dam being the largest dam in Turkey with a total water capacity of 48.7 billion m(3). After the construction of the dam, the monitoring network has improved. Considering earthquakes above the long-term completeness magnitude of M-C = 3.5, the local seismicity rate has substantially increased after the filling of the reservoir. Recently, two damaging earthquakes of M-w 5.5 and M-w 5.1 occurred in the town of Samsat near the Ataturk Reservoir in 2017 and 2018, respectively. In this study, we analyze the spatio-temporal evolution of seismicity and its source properties in relation to the temporal water-level variations and the stresses resulting from surface loading and pore-pressure diffusion. We find that water-level and seismicity rate are anti-correlated, which is explained by the stabilization effect of the gravitational induced stress imposed by water loading on the local faults. On the other hand, we find that the overall effective stress in the seismogenic zone increased over decades due to pore-pressure diffusion, explaining the enhanced background seismicity during recent years. Additionally, we observe a progressive decrease of the Gutenberg-Richter b-value. Our results indicate that the stressing rate finally focused on the region where the two damaging earthquakes occurred in 2017 and 2018.
The accelerating climatic changes and new infrastructure development across the Arctic require more robust risk and environmental assessment, but thus far there is no consistent record of human impact. We provide a first panarctic satellite-based record of expanding infrastructure and anthropogenic impacts along all permafrost affected coasts (100 km buffer, approximate to 6.2 Mio km(2)), named the Sentinel-1/2 derived Arctic Coastal Human Impact (SACHI) dataset. The completeness and thematic content goes beyond traditional satellite based approaches as well as other publicly accessible data sources. Three classes are considered: linear transport infrastructure (roads and railways), buildings, and other impacted area. C-band synthetic aperture radar and multi-spectral information (2016-2020) is exploited within a machine learning framework (gradient boosting machines and deep learning) and combined for retrieval with 10 m nominal resolution. In total, an area of 1243 km(2) constitutes human-built infrastructure as of 2016-2020. Depending on region, SACHI contains 8%-48% more information (human presence) than in OpenStreetMap. 221 (78%) more settlements are identified than in a recently published dataset for this region. 47% is not covered in a global night-time light dataset from 2016. At least 15% (180 km(2)) correspond to new or increased detectable human impact since 2000 according to a Landsat-based normalized difference vegetation index trend comparison within the analysis extent. Most of the expanded presence occurred in Russia, but also some in Canada and US. 31% and 5% of impacted area associated predominantly with oil/gas and mining industry respectively has appeared after 2000. 55% of the identified human impacted area will be shifting to above 0 C-circle ground temperature at two meter depth by 2050 if current permafrost warming trends continue at the pace of the last two decades, highlighting the critical importance to better understand how much and where Arctic infrastructure may become threatened by permafrost thaw.
We revise the conceptual model of calcite varves and present, for the first time, a dual lake monitoring study in two alkaline lakes providing new insights into the seasonal sedimentation processes forming these varves. The study lakes, Tiefer See in NE Germany and Czechowskie in N Poland, have distinct morphology and bathymetry, and therefore, they are ideal to decipher local effects on seasonal deposition. The monitoring setup in both lakes is largely identical and includes instrumental observation of (i) meteorological parameters, (ii) chemical profiling of the lake water column including water sampling, and (iii) sediment trapping at both bi-weekly and monthly intervals. We then compare our monitoring data with varve micro-facies in the sediment record. One main finding is that calcite varves form complex laminae triplets rather than simple couplets as commonly thought. Sedimentation of varve sub-layers in both lakes is largely dependent on the lake mixing dynamics and results from the same seasonality, commencing with diatom blooms in spring turning into a pulse of calcite precipitation in summer and terminating with a re-suspension layer in autumn and winter, composed of calcite patches, plant fragments and benthic diatoms. Despite the common seasonal cycle, the share of each of these depositional phases in the total annual sediment yield is different between the lakes. In Lake Tiefer See calcite sedimentation has the highest yields, whereas in Lake Czechowskie, the so far underestimated re-suspension sub-layer dominates the sediment accumulation. Even in undisturbed varved sediments, re-suspended material becomes integrated in the sediment fabric and makes up an important share of calcite varves. Thus, while the biogeochemical lake cycle defines the varves' autochthonous components and micro-facies, the physical setting plays an important role in determining the varve sub-layers' proportion.
Large earthquakes can increase the amount of water feeding stream flows, raise groundwater levels, and thus grant plant roots more access to water in water-limited environments. We examine growth and photosynthetic responses of Pine plantations to the Maule M-w 8.8 earthquake in headwater catchments of Chile's Coastal Range. We combine high-resolution wood anatomic (lumen area) and biogeochemical (delta 13C of wood cellulose) proxies of daily to weekly tree growth sampled from trees on floodplains and close to ridge lines. We find that, immediately after the earthquake, at least two out of six tree trees on valley floors had increased lumen area and decreased delta 13C, while trees on hillslopes had a reverse trend. Our results indicate a control of soil water on this response, largely consistent with models that predict how enhanced postseismic vertical soil permeability causes groundwater levels to rise on valley floors, but fall along the ridges. Statistical analysis with boosted regression trees indicates that streamflow discharge gained predictive importance for photosynthetic activity on the ridges, but lost importance on the valley floor after the earthquake. We infer that earthquakes may stimulate ecohydrological conditions favoring tree growth over days to weeks by triggering stomatal opening. The weak and short-lived signals that we identified, however, show that such responses are only valid under water-limited, rather than energy-limited tree, growth. Hence, dendrochronological studies targeted at annual resolution may overlook some earthquake effects on tree vitality.
Seafloor spreading at slow rates can be accommodated on large-offset oceanic detachment faults (ODFs), that exhume lower crustal and mantle rocks in footwall domes termed oceanic core complexes (OCCs). Footwall rocks experience large rotation during exhumation, yet important aspects of the kinematics-particularly the relative roles of solid-block rotation and flexure-are not clearly understood. Using a high-resolution numerical model, we explore the exhumation kinematics in the footwall beneath an emergent ODF/OCC. A key feature of the models is that footwall motion is dominated by solid-block rotation, accommodated by the nonplanar, concave-down fault interface. A consequence is that curvature measured along the ODF is representative of a neutral stress configuration, rather than a "bent" one. Instead, it is in the subsequent process of "apparent unbending" that significant flexural stresses are developed in the model footwall. The brittle strain associated with apparent unbending is produced dominantly in extension, beneath the OCC, consistent with earthquake clustering observed in the Trans-Atlantic Geotraverse at the Mid-Atlantic Ridge.
The wetland cover is defined as the spatially homogenous region of a wetland attributed to the underlying biophysical conditions such as vegetation, turbidity, hydric soil, and the amount of water.
Here, we present a novel method to derive the wetland-cover types (WCTs) combining three commonly used multispectral indices, NDVI, MNDWI, and NDTI, in three large Ramsar wetlands located in different geomorphic and climatic settings across India. These wetlands include the Kaabar Tal, a floodplain wetland in east Ganga Plains, Chilika Lagoon, a coastal wetland in eastern India, and Nal Sarovar in semi-arid western India.
The novelty of our approach is that the derived WCTs are stable in space and time, and therefore, a given WCT across different wetlands or within different zones of a large wetland will imply similar underlying biophysical attributes.
The WCTs can therefore provide a novel tool for monitoring and change detection of wetland cover types.
We have automated the proposed WCT algorithm using the Google Earth Engine (GEE) environment and by developing ArcGIS tools. The method can be implemented on any wetland and using any multispectral imagery dataset with visible and NIR bands.
The proposed methodology is simple yet robust and easy to implement and, therefore, holds significant importance in wetland monitoring and management.
The Alpine mountains in central Europe are characterized by a heterogeneous crust accumulating different tectonic units and blocks in close proximity to sedimentary foreland basins. Centroid moment tensor inversion provides insight into the faulting mechanisms of earthquakes and related tectonic processes but is significantly aggravated in such an environment. Thanks to the dense AlpArray seismic network and our flexible bootstrap-based inversion tool Grond, we are able to test different setups with respect to the uncertainties of the obtained moment tensors and centroid locations. We evaluate the influence of frequency bands, azimuthal gaps, input data types, and distance ranges and study the occurrence and reliability of non-double-couple (DC) components. We infer that for most earthquakes (M-w >= 3.3) a combination of time domain full waveforms and frequency domain amplitude spectra in a frequency band of 0.02-0.07 Hz is suitable. Relying on the results of our methodological tests, we perform deviatoric moment tensor (MT) inversions for events with M-w > 3.0. Here, we present 75 solutions for earthquakes between January 2016 and December 2019 and analyze our results in the seismotectonic context of historical earthquakes, seismic activity of the last 3 decades, and GNSS deformation data. We study regions of comparably high seismic activity during the last decades, namely the Western Alps, the region around Lake Garda, and the eastern Southern Alps, as well as clusters further from the study region, i.e., in the northern Dinarides and the Apennines. Seismicity is particularly low in the Eastern Alps and in parts of the Central Alps. We apply a clustering algorithm to focal mechanisms, considering additional mechanisms from existing catalogs. Related to the N-S compressional regime, E-W-to-ENE-WSW-striking thrust faulting is mainly observed in the Friuli area in the eastern Southern Alps. Strike-slip faulting with a similarly oriented pressure axis is observed along the northern margin of the Central Alps and in the northern Dinarides. NW-SE-striking normal faulting is observed in the NW Alps, showing a similar strike direction to normal faulting earthquakes in the Apennines. Both our centroid depths and hypocentral depths in existing catalogs indicate that Alpine seismicity is predominantly very shallow; about 80% of the studied events have depths shallower than 10 km.
The magnitude of earthquakes on continental normal faults rarely exceeds 7.0 Mw. However, because of their vicinity to large population centers they can be highly destructive.
Long recurrence time, relatively small deformations, and limited observations hinder our understanding of the deformation patterns and mechanisms controlling the magnitude of events.
Here, this problem is addressed with 2D thermomechanical modeling of normal fault seismic cycles.
The 2020 Samos, Greece Mw7.0 earthquake is used as an example as it is one of the largest and most studied continental normal fault earthquakes. The modeling approach employs visco-elasto-plastic rheology, compressibility, free surface, and a rate-and-state friction law for the fault.
Modeling of the Samos earthquake suggests the pore fluid pressure ratio on the fault ranges from 0 to 0.7. The model demonstrates that most of the deformation during interseismic and coseismic periods, besides on the fault, occurs in the hanging wall and footwall below the seismogenic part of the fault. The largest vertical surface displacement during the earthquake is the subsidence of the hanging wall in the vicinity of the fault, while the uplift of the footwall and remote part of the hanging wall is significantly smaller.
Modeling of the seismic cycles on normal faults with different setups shows the dependency of the magnitude on the thermal profile and dipping angle of the fault; low heat flow and low dipping angle are favorable conditions for the largest events, while steep normal faults in the areas of high heat flow tend to have the smallest magnitudes.
In a warming Arctic, permafrost-related disturbances, such as retrogressive thaw slumps (RTS), are becoming more abundant and dynamic, with serious implications for permafrost stability and bio-geochemical cycles on local to regional scales. Despite recent advances in the field of earth observation, many of these have remained undetected as RTS are highly dynamic, small, and scattered across the remote permafrost region. Here, we assessed the potential strengths and limitations of using deep learning for the automatic segmentation of RTS using PlanetScope satellite imagery, ArcticDEM and auxiliary datasets. We analyzed the transferability and potential for pan-Arctic upscaling and regional cross-validation, with independent training and validation regions, in six different thaw slump-affected regions in Canada and Russia. We further tested state-of-the-art model architectures (UNet, UNet++, DeepLabv3) and encoder networks to find optimal model configurations for potential upscaling to continental scales. The best deep learning models achieved mixed results from good to very good agreement in four of the six regions (maxIoU: 0.39 to 0.58; Lena River, Horton Delta, Herschel Island, Kolguev Island), while they failed in two regions (Banks Island, Tuktoyaktuk). Of the tested architectures, UNet++ performed the best. The large variance in regional performance highlights the requirement for a sufficient quantity, quality and spatial variability in the training data used for segmenting RTS across diverse permafrost landscapes, in varying environmental conditions. With our highly automated and configurable workflow, we see great potential for the transfer to active RTS clusters (e.g., Peel Plateau) and upscaling to much larger regions.
The formation of the Central Andes dates back to similar to 50 Ma, but its most pronounced episode, including the growth of the Altiplano-Puna Plateau and pulsatile tectonic shortening phases, occurred within the last 25 Ma.
The reason for this evolution remains unexplained. Using geodynamic numerical modeling we infer that the primary cause of the pulses of tectonic shortening and growth of the Central Andes is the changing geometry of the subducted Nazca plate, and particularly the steepening of the mid-mantle slab segment which results in a slowing down of the trench retreat and subsequent increase in shortening of the advancing South America plate.
This steepening first happens after the end of the flat slab episode at similar to 25 Ma, and later during the buckling and stagnation of the slab in the mantle transition zone. Processes that mechanically weaken the lithosphere of the South America plate, as suggested in previous studies, enhance the intensity of the shortening events.
These processes include delamination of the mantle lithosphere and weakening of foreland sediments.
Our new modeling results are consistent with the timing and amplitude of the deformation from geological data in the Central Andes at the Altiplano latitude.
Plain Language Summary
The Central Andes is a subduction-type orogeny that formed as a result of the interaction between the Nazca oceanic plate and the South American continental plate over the last 50 million years. Growth of the Andes is primarily the result of crustal shortening. Nevertheless, "geological" data compiled from previous studies have shown that phases of drastic pulsatile shortening occur at 15 and 5 Ma.
In this study, we used high-resolution 2D numerical geodynamic simulations to investigate the link between oceanic and continental plate dynamics and their interaction. We find that when the oceanic plate steepens in the mantle transition zone, the trench retreat is hindered. Coupled with the weakening of the continental plate through the slab flattening and subsequent delamination of the lithospheric mantle, this leads to pulsatile shortening phases of a magnitude equivalent to that suggested by the data.
Magma-filled dikes may feed erupting fissures that lead to alignments of craters developing at the surface, yet the details of activity and migrating eruptions at the crater row are difficult to monitor and are hardly understood.
The 2021 Tajogaite eruption at the Cumbre Vieja, La Palma (Spain), lasted 85 days and developed a pronounced alignment of craters that may be related to changes within the volcano edifice.
Here, we use COSMO-SkyMed satellite radar data and ground-based time-lapse photographs, offering a high-resolution dataset to explore the locations and characteristics of evolving craters.
Our results show that the craters evolve both gradually and suddenly and can be divided into three main phases. Phase 1, lasting the first 6 weeks of the eruption, was characterized by a NW-SE linear evolution of up to seven craters emerging on the growing cone.
Following two partial collapses of the cone to the northwest and a seismicity increase at depth, Phase 2 started and caused a propagation of the main activity toward the southeastern side, together with the presence of up to 11 craters along this main NW-SE trend. Associated with strong deep and shallow earthquakes, Phase 3 was initiated and continued for the final 2 weeks of the eruption, expressed by the development of up to 18 craters, which became dominant and clustered in the southeastern sector in early December 2021. In Phase 3, a second and oblique alignment and surface fracture was identified.
Our findings that crater and eruption changes coincide together with an increase in seismic activity at depth point to a deep driver leading to crater and morphology changes at the surface.
These also suggest that crater distributions might allow for improved monitoring of changes occurring at depth, and vice versa, such that strong seismicity changes at depth may herald the migration and new formation of craters, which have major implications for the assessment of tephra and lava flow hazards on volcanoes.
The Sea of Marmara is a tectonically active basin that straddles the North Anatolian Fault Zone (NAFZ), a major strike-slip fault that separates the Eurasian and Anatolian tectonic plates. The Main Marmara Fault (MMF), which is part of the NAFZ, contains an approximately 150 km long seismotectonic segment that has not ruptured since 1766. A key question for seismic hazard and risk assessment is whether or not the next rupture along this segment is likely to produce one major earthquake or a series of smaller earthquakes. Geomechanical characteristics such as along-strike variations in rock strength may provide an important control on seismotectonic segmentation. We find that variations in lithospheric strength throughout the Marmara region control the mechanical segmentation of the MMF and help explain its long-term seismotectonic segmentation. In particular, a strong crust that is mechanically coupled to the upper mantle spatially correlates with aseismic patches, where the MMF bends and changes its strike in response to the presence of high-density lower crustal bodies. Between the bends, mechanically weaker crustal domains that are decoupled from the mantle indicate a predominance of creeping. These results are highly relevant for the ongoing debate regarding the characteristics of the Marmara seismic gap, especially in view of the seismic hazard (Mw > 7) in the densely populated Marmara region.
On 12 August 2021, a > 220 s lasting complex earthquake with Mw > 8.2 hit the South Sandwich Trench. Due to its remote location and short interevent times, reported earthquake parameters varied significantly between different international agencies. We studied the complex rupture by combining different seismic source characterization techniques sensitive to different frequency ranges based on teleseismic broadband recordings from 0.001 to 2 Hz, including point and finite fault inversions and the back-projection of high-frequency signals.
We also determined moment tensor solutions for 88 aftershocks. The rupture initiated simultaneously with a rupture equivalent to a M-w 7.6 thrust earthquake in the deep part of the seismogenic zone in the central subduction interface and a shallow megathrust rupture, which propagated unilaterally to the south with a very slow rupture velocity of 1.2 km/s and varying strike following the curvature of the trench.
The slow rupture covered nearly two-thirds of the entire subduction zone length, and with M-w 8.2 released the bulk of the total moment of the whole earthquake.
Tsunami modeling indicates the inferred shallow rupture can explain the tsunami records. The southern segment of the shallow rupture overlaps with another activation of the deeper part of the megathrust equivalent to M-w 7.6. The aftershock distribution confirms the extent and curvature of the rupture. Some mechanisms are consistent with the mainshocks, but many indicate also activation of secondary faults. Rupture velocities and radiated frequencies varied strongly between different stages of the rupture, which might explain the variability of published source parameters.
Plain Language Summary
The earthquake of 12 August 2021 along the deep-sea trench of the South Sandwich Islands in the South Atlantic reached a magnitude of 8.2 and triggered a tsunami. The automatic earthquake parameter determination of different agencies showed very different results shortly after the earthquake and partially underestimated the tsunami potential of the earthquake.
A possible reason was the complex rupture process and that the tsunami was generated by a long and shallow slow slip rupture sandwiched between more conventional fast slip subevents at its northern and southern ends. In addition, the fault surface, which extended over 450 km, was highly curved striking 150 degrees-220 degrees.
We investigated the different components of the seismic wavefields in different frequency ranges and with different methods.
The analysis shows how even complex earthquakes can be deciphered by combining analyzing methods. The comparison with aftershocks and the triggered tsunami waves confirms our model that explains the South Sandwich rupture by four subevents in the plate boundary along the curved deep-sea trench. Here, the depth, rupture velocities, and slip on each segment of the rupture vary considerably.
The method can also be applied to other megathrust earthquakes and help to further improve tsunami warnings in the future.
New middle Miocene to Pliocene (~14–3 Ma) apatite fission track (AFT) cooling ages combined with published K–Ar/Ar–Ar and zircon fission track (ZFT) ages from the Hazara and Swat regions of Pakistan are used to explain the Oligocene to Pliocene structural evolution in the Western Himalaya. The structural model explains the distribution of K–Ar/Ar–Ar ages in three distinct age groups (Proterozoic, Paleozoic-Mesozoic, and Eocene to Oligocene). The Proterozoic to Mesozoic sequence of northern Hazara and Swat experienced elevated temperature and pressure conditions, evident by reset Eocene to Oligocene K–Ar/Ar–Ar hornblende and Eocene to Miocene muscovite ages, caused by Kohistan overthrusting the Indian margin during and after the India–Asia collision. Samples from the Indus syntaxis with Paleo to Mesoproterozoic K–Ar/Ar–Ar hornblende ages and Eocene to Oligocene Ar–Ar muscovite ages show no signs of Cenozoic metamorphism; these samples were thermally imprinted up to the Ar–Ar muscovite closure temperature. Neoproterozoic to Lower Paleozoic rocks from the southern parts of Hazara and Swat show Mesozoic to Oligocene partially reset Ar–Ar muscovite ages and preservation of Ordovician metamorphism. The combined analysis of published K–Ar/Ar–Ar (muscovite), ZFT, and new AFT ages (~14–12 Ma) suggests that the Main Central thrust/Panjal thrust was active from Oligocene to early Miocene (~30–18 Ma), and the Nathia-Gali and Main Boundary thrusts were active from the middle to late Miocene (~14–9 Ma) in the Hazara area. New and published AFT ages (~6–3 Ma) from the Indus syntaxis suggest that early Pliocene tectonic thickening in the hinterland formed the N–S trending Indus anticline, creating an erosional half window in the Main Mantle thrust, forming the Indus syntaxis, and dividing the Main Central thrust sheet into the Hazara and Swat segments.
We propose lacunarity as a novel recurrence quantification measure and illustrate its efficacy to detect dynamical regime transitions which are exhibited by many complex real-world systems. We carry out a recurrence plot-based analysis for different paradigmatic systems and nonlinear empirical data in order to demonstrate the ability of our method to detect dynamical transitions ranging across different temporal scales. It succeeds to distinguish states of varying dynamical complexity in the presence of noise and non-stationarity, even when the time series is of short length. In contrast to traditional recurrence quantifiers, no specification of minimal line lengths is required and geometric features beyond linear structures in the recurrence plot can be accounted for. This makes lacunarity more broadly applicable as a recurrence quantification measure. Lacunarity is usually interpreted as a measure of heterogeneity or translational invariance of an arbitrary spatial pattern. In application to recurrence plots, it quantifies the degree of heterogeneity in the temporal recurrence patterns at all relevant time scales. We demonstrate the potential of the proposed method when applied to empirical data, namely time series of acoustic pressure fluctuations from a turbulent combustor. Recurrence lacunarity captures both the rich variability in dynamical complexity of acoustic pressure fluctuations and shifting time scales encoded in the recurrence plots. Furthermore, it contributes to a better distinction between stable operation and near blowout states of combustors.
The present-day structure of the Eastern Cordillera in NW Argentina is governed by structural and lithological heterogeneities inherited from preceding deformation phases, which influence the localization of newly-formed faults and the inversion of pre-existing structures.
The Salta Rift Basin formed during a Late Jurassic-Cretaceous extensional phase and created a dominant structural and stratigraphic imprint in NW Argentina that is partic-ularly evident within the Eastern Cordillera, where uplift and exhumation have exposed the Salta Group syn-rift succession.
Although in general, the Salta Group rests upon Paleozoic rocks, locally the Tacuru Group forms an intermediate succession, consisting of interfingering eolian sandstones and proximal fault-related conglomerates with a Jurassic maximum depositional age. This succession might be the key to unraveling the Mesozoic history of NW Argentina, prior to the deposition of the Salta Group.
The conglomerates represent the earliest deposits related to extension in the western Lomas de Olmedo sub-basin, which is also documented in predominantly Jurassic zircon (U-Th-Sm)/He cooling ages of the rift shoulders. The detrital zircon U-Pb age signature and sandstone provenance of the Tacuru Group conglomerates differs strongly from the Salta Group syn-rift strata, which show a more regional signal.
These variations and the angularity of the unconformity may be connected to a rotation of the extension direction in the western Lomas de Olmedo sub-basin.
In the Gasht-Masuleh area in the Alborz Mountains, gabbroic magma intruded Palaeozoic metasediments and Mesozoic sediments and crystallised as isotropic and cumulate gabbros. LREE enrichment points to relatively low degrees of mantle melting and depletion of Ti, Nb and Ta relative to primitive mantle points to an arc related component in the magma. Clinopyroxene compositions indicate MORB to arc signatures. U–Pb zircon crystallisation ages of 99.5 ± 0.6 Ma and 99.4 ± 0.6 Ma and phlogopite 40Ar/39Ar ages of 97.1 ± 0.4 Ma, 97.5 ± 0.4 Ma, 97.1 ± 0.1 Ma, within 2σ error, indicate that gabbro intrusion occurred in the (Albian-)Cenomanian (mid-Cretaceous). As active subduction did not take place in the Cretaceous in North Iran, the small volume mafic magmatism in the Gasht-Masuleh area must be due to local, extension-related mantle melting. Melting was most likely caused by far field effects triggered by roll-back of the Neo-Tethys subducting slab. As subduction took place at a distance of ~ 400 km (present distance) from the Alborz Mountains, the observed arc geochemical signatures must be inherited from a previous subduction event and concomitant mantle metasomatism, possibly in combination with contamination of the magma by crustal material.
In recurrence analysis, the tau-recurrence rate encodes the periods of the cycles of the underlying high-dimensional time series. It, thus, plays a similar role to the autocorrelation for scalar time-series in encoding temporal correlations.
However, its Fourier decomposition does not have a clean interpretation. Thus, there is no satisfactory analogue to the power spectrum in recurrence analysis.
We introduce a novel method to decompose the tau-recurrence rate using an over-complete basis of Dirac combs together with sparsity regularization.
We show that this decomposition, the inter-spike spectrum, naturally provides an analogue to the power spectrum for recurrence analysis in the sense that it reveals the dominant periodicities of the underlying time series.
We show that the inter-spike spectrum correctly identifies patterns and transitions in the underlying system in a wide variety of examples and is robust to measurement noise.
Late Miocene-Pliocene onset of fluvial incision of the Cauca River Canyon in the Northern Andes
(2022)
The incision of kilometer-scale canyons into high-standing topography is often used to constrain the surface uplift history of mountain ranges, controlled by tectonic and geodynamic processes.
However, changes in climate may also be responsible for canyon incision. This study deciphers the timing of incision of the similar to 2.5-km-deep Cauca River Canyon in the Central Cordillera of the Northern Andes using the cooling (exhumation) history of rocks from the canyon walls and a regional analysis of channel steepness in rivers.
Ten bedrock samples and one detrital sample were collected on the eastern border of the canyon between 300 m and 2300 m of elevation.
Bedrock and detrital AFT data yield ages from 50 to 38 Ma, while two bed-rock AHe ages from the valley bottom yield ages of 7-6 Ma.
The AHe ages and inverse thermal history models reveal a previously unidentified late Miocene (ca. 7-6 Ma) pulse of exhumation that we interpret as the age of a single incision event that formed the Cauca River Canyon.
We conclude that the Cauca River Canyon was carved as a response to rock uplift in the northern Central Cordillera and propagation of an erosion wave into the mountain range starting in the latest Miocene.
Strain NGK35T is a motile, Gram-stain-negative, rod-shaped (1.0-2.1 mu m long and 0.6-0.8 mu m wide), aerobic bacterium that was isolated from plastic-polluted landfill soil. The strain grew at temperatures between 6 and 37 degrees C (optimum, 28 degrees C), in 0-10 % NaCl (optimum, 1 %) and at pH 6.0-9.5 (optimum, pH 7.5-8.5).
It was positive for cytochrome c oxidase, catalase as well as H2S production, and hydrolysed casein and urea. It used a variety of different carbon sources including citrate, lactate and pyruvate.
The predominant membrane fatty acids were C-16:1 cis9 and C-16:0, followed by C-17:0 cyclo and C-18:1 cis11. The major polar lipids were phosphatidylglycerol and phosphatidylethanolamine, followed by diphosphatidyglycerol. The only quinone was ubiquinone Q-8. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain NGK35(T) belongs to the genus Paenalcaligenes (family Alcaligenaceae), appearing most closely related to Paenalcaligenes hominis CCUG 53761A(T) (96.90 %) and Paenalcaligenes suwonensis ABC02-12(T) (96.94 %).
The genomic DNA G+C content of strain NGK35(T) was 52.1 mol%. Genome-based calculations (genome-to-genome distance, average nucleotide identity and DNA G+C content) clearly indicated that the isolate represents a novel species within the genus Paenalcaligenes.
Based on phenotypic and molecular characterization, strain NGK35(T) can clearly be differentiated from its phylogenetic neighbours establishing a novel species, for which the name Paenalcaligenes niemegkensis sp. nov. is proposed.
The type strain is NGK35T (=DSM 113270(T)=NCCB 100854(T)).
A marine sediment record from the central Bering Sea, spanning the last 20 thousand years (ka), was studied to unravel the depositional history with regard to terrigenous sediment supply and biogenic sedimentation. Methodic approaches comprised the inference of accumulation rates of siliciclastic and biogenic components, grain-size analysis, and (clay) mineralogy, as well as paleoclimatic modelling. Changes in the depositional history provides insight into land-ocean linkages of paleoenvironmental changes. During the finale of the Last Glacial Maximum, the depositional environment was characterized by hemipelagic background sedimentation. A marked change in the terrigenous sediment provenance during the late Heinrich 1 Stadial (15.7-14.5 ka), indicated by increases in kaolinite and a high glaciofluvial influx of clay, gives evidence of the deglaciation of the Brooks Range in the hinterland of Alaska. This meltwater pulse also stimulated the postglacial onset of biological productivity. Glacial melt implies regional climate warming during a time of widespread cooling on the northern hemisphere. Our simulation experiment with a coupled climate model suggests atmospheric teleconnections to the North Atlantic, with impacts on the dynamics of the Aleutian Low system that gave rise to warmer winters and an early onset of spring during that time. The late deglacial period between 14.5 and 11.0 ka was characterized by enhanced fluvial runoff and biological productivity in the course of climate amelioration, sea-level rise, seasonal sea-ice retreat, and permafrost thaw in the hinterland. The latter processes temporarily stalled during the Younger Dryas stadial (12.9-11.7 ka) and commenced again during the Preboreal (earliest Holocene), after 11.7 ka. High river runoff might have fertilized the Bering Sea and contributed to enhanced upper ocean stratification. Since 11.0 ka, advanced transgression has shifted the coast line and fluvial influence of the Yukon River away from the study site. The opening of the Bering Strait strengthened contour currents along the continental slope, leaving behind winnowed sand-rich sediments through the early to mid-Holocene, with non-deposition occurring since about 6.0 ka.
Rapid humidity changes across the Northern South China Sea during the last similar to 40 kyrs
(2021)
A key aspect of East Asian climate is its summer monsoonal system which influences nearly one-third of the world's population. Recent results indicate that the primary response of the East Asian summer monsoon (EASM) to anthropogenic forced climate warming may be a shift in geographical range instead of an intensity change, which would lead to spatial coexistence of floods and droughts over southeastern Asia. The predicted EASM variability in the future has made it paramount to study its past changes and the associated tempo-spatial pattern of aridity and humidity in its purview. In order to decipher past changes in EASM, we applied a multi-proxy geochemical approach to the sediment core ORI-891-16-P1 located in the northern South China Sea. The position of this sediment core on top of a seamount makes it uniquely sensitive to changes in the terrigenous input into northern South China Sea unbiased by sea level-induced downslope transport processes. Utilizing the ln(Ti/Ca) ratio throughout the sediment sequence we trace terrigenous influx changes reflecting EASM prevalence during the last similar to 40 kyrs. Based on the comparison of our results to previous studies we infer that the Last Glacial Maximum (LGM; similar to 20 ka BP) was characterized by a steep N-S humidity gradient. This spatial pattern was in line with a southward shift or contraction of the summer monsoonal trough of 10-15 degrees from its current position toward the centre of the South China Sea. Superimposed on orbital time scale fluctuations we also find strong indication of millennial-scale variability related to Heinrich Stadials. The impact of Heinrich Stadials on the EASM seems amplified during insolation minima, while high summer insolation seems to buffer the monsoonal system to such perturbations. We infer that (i) the humidity-aridity distribution during the LGM mimics predictions of the proposed future EASM configuration, and (ii) that the sensitivity of the EASM to weakening in the Atlantic Meridional Overturning Circulation is the strongest since the last glacial.
With the present study, we introduce a fast and robust method to calculate the source displacement spectra of small earthquakes on a local to regional scale. The work is based on the publicly available Qopen method of full envelope inversion, which is further tuned for the given purpose. Important source parameters-seismic moment, moment magnitude, corner frequency, and high-frequency fall off-are determined from the source spectra by fitting a simple earthquake source model. The method is demonstrated by means of a data set comprising the 2018 West Bohemia earthquake swarm. We report moment magnitudes, corner frequencies, and centroid moment tensors inverted from short-period body waves with the Grond package for all earthquakes with a local magnitude larger than 1.8. Moment magnitudes calculated by envelope inversion show a very good agreement to moment magnitudes resulting from the probabilisitc moment tensor inversion. Furthermore, source displacement spectra from envelope inversion show a good agreement with spectra obtained by multiple taper analysis of the direct onsets of body waves but are not affected by the large scatter of the second. The seismic moments obtained with the envelope inversion scale with corner frequencies according to M-0 proportional to f(c)(-4.7). Earthquakes of the present data set result in a smaller stress drop for smaller magnitudes. Self-similarity of earthquake rupture is not observed. In addition, we report frequency-dependent site amplification at the used stations.
The Shanderman lamprophyre dykes crop out in the western part of the Alborz Mountains (Talesh).
These rocks are classified as camptonites, composed of primary olivine, Ti-rich diopside, kaersutite, biotite, plagioclase, K-feldspar, and minor Ti-rich spinels, magnetite, pentlandite-pyrrhotite/chalcopyrite, and powellite-scheelite. Secondary analcime-wairakite, serpentines, and prehnite are common minor minerals within the studied rocks.
Olivine, Ti-rich diopside, spinel, and amphibole show distinct chemical zoning. Spinels display a core-to-rim decrease in Cr2O3, MgO, and Al2O3 concentrations and an increase in TiO2 and FeOT (total Fe as FeO), reflecting the oxidation state increase due to hydrothermal fluid influx. Low SiO2 contents (< 42 wt%), high MgO (12.44 to 13.98 wt%), and Fe2O3T (12.76 to 13.43 wt%), Cr (318-537 mu g/g) and Ni (231-327 mu g/g) contents indicate the ultrabasic nature of the rocks.
The samples show potassic character (2.1-2.8 wt% K2O), along with elevated LREE and LILE, and also exhibit minor positive Eu anomalies (Eu/Eu* = 1.09 to 1.20).
Olivine-spinel geothermometry indicates a maximum crystallization temperature of 1227 degrees C (ave. 988 degrees C +/- 65 degrees C).
Exsolution of pentlandite-pyrrhotite/chalcopyrite solid solutions occurred during magma cooling and crystallization. At lower temperatures, analcime-wairakite and prehnite partially replaced plagioclases.
The geochemical modeling of the rocks indicates the Shanderman lamprophyre magmas were derived from low-grade melting (< 5%) of amphibole-bearing garnet lherzolite source without or with very few phlogopites.
The primary magma of Shanderman lamprophyres was derived from a depth of similar to 135 km by partial melting of a metasomatized mantle source in a post-collisional environment.
Ground subsidence caused by natural or anthropogenic processes affects major urban areas worldwide. Sinkhole formation and infrastructure fractures have intensified in the federal capital of Maceio (Alagoas, Brazil) since early 2018, forcing authorities to relocate affected residents and place buildings under demolition. In this study, we present a 16-year history (2004-2020) of surface displacement, which shows precursory deformations in 2004-2005, reaching a maximum cumulative subsidence of approximately 200 cm near the Mundau Lagoon coast in November 2020. By integrating the displacement observations with numerical source modelling, we suggest that extensive subsidence can be primarily associated with the removal of localized, deep-seated material at the location and depth where salt is mined. We discuss the accelerating subsidence rates, influence of severe precipitation events on the aforementioned geological instability, and related hazards. This study suggests that feedback destabilization mechanisms may arise in evaporite systems due to anthropogenic activities, fostering enhanced and complex superficial ground deformation.
Geochemical homogeneity in shale is often assumed when tracing subsurface fluids and characterizing sedimentary basins. This study presents measurements of the bulk gas composition, stable isotopes, and noble gas volume fraction and isotopes for shale gas samples collected from gas wells in the Wufeng-Longmaxi Shale, the southern Sichuan Basin, China. The dryness [C-1 /(C-2 + C-3)] ranging from 166.3 to 251.2, combined with delta C-13(1) and delta DC1 that vary from -28.8 to -27.3 parts per thousand and - 153 to -145 parts per thousand, respectively, point to a late mature thermogenic origin of hydrocarbon gas. He-3/He-4 ratios of gas samples are around 0.01 times the air value suggesting dominantly crust-derived He. Ne-21/Ne-22 and Ar-40/Ar-36 ratios of many gas samples are higher than the corresponding air values indicating the mixing of crustal and atmospheric noble gases. Multiple dichotomous patterns are observed in noble gas signatures of forelimb and backlimb samples, and depression and crest samples. Ne-20/Ne-22 ratios of some crest samples are higher than that of depression samples in the backlimb, pointing to the presence of diffusion-driven fractionation that is likely caused by the long-distance migration from depression to crest. Elemental ratios of air-derived noble gas isotopes - Ne-22/Ar-36, Kr-84/Ar-36, and Xe-132/Ar-36 are compared to the recharge water values, suggesting the interactions of oil, gas, and water phases in the shale over geologic time. Forelimb samples generally display older ages than backlimb samples, indicating a larger flux of external radiogenic He-4 due to the higher density of deep faults in the forelimb area caused by the basementinvolved deformation. The basement-involved deformation also causes pore collapse especially in the forelimb leading to a lower porosity that results in a more pristine noble gas signature in the forelimb due to the reduced impact of younger recharge water.
Equilibrium mass-dependent ("stable") isotopic fractionation of an element during magmatic processes is driven by a contrast in bonding environment between minerals and silicate melt, which is expressed as an isotopic fractionation factor.
A quantitative understanding of such isotopic fractionation factors is vital to interpret observed isotopic variations in magmatic rocks.
It is well known that the local environment and the bond strength of an element dictate the sign and magnitude of isotopic fractionation between minerals, but it is uncertain how the structure and chemical composition of a silicate melt can affect mineral-melt isotopic fractionation factors.
To explore this, we studied the coordination environment of nickel (Ni) in different silicate glasses using extended X-ray absorption fine structure (EXAFS) measurements at the German synchrotron X-ray source (DESY).
We determined -Ni-O bond lengths in a suite of synthetic but near-natural silicate glasses using EXAFS and found that the former vary systematically with melt alkalinity, which is best described by the parameter ln[1 + (Na + K)/Ca]. With increasing melt alkalinity, Ni occupies more IV-fold coordinated sites, which are associated with a shorter -Ni-O bond length. Next, we use the ionic model, which allows to predict isotopic fractionation factors based on the difference in bond length between two phases.
We find that more alkaline melts have a stronger preference for the heavier isotopes of Ni than less alkaline melts. This implies that the magnitude of mineral-melt Ni isotope fractionation factors, for instance between olivine and melt, will depend on the alkalinity of the melt.
At magmatic temperatures, however, the variation in fractionation factors caused by melt alkalinity will rarely exceed 0.05 parts per thousand and is thus mostly negligible, in particular in the realm of basaltic melt compositions. Nevertheless, the relationship between melt alkalinity and fractionation factor reported here can be used to extrapolate empirical data for mineral-melt Ni isotope fractionation factors, once such data become available, to the full range of magma compositions on Earth and other Solar System bodies.
Origin and migration of fluoride in the area of the Aluto Volcanic Complex (Main Ethiopian Rift)
(2022)
Fluoride-enriched ground and surface waters represent a major health risk for the local population in many areas along the East African Rift.
The present study investigates the origin of fluoride and the reason for its accumulation in the rift waters, following two hypotheses: (i) fluid-rock-interactions release fluoride from minerals into the water and (ii) magmatic-derived fluoride-containing liquids and gases migrate along permeable fault zones until they dissolve in ground-and surface water or be released to the atmosphere.
Rock-, gas, water-, soil-, and plant samples were collected from the area within and close by the Aluto Volcanic Complex, which is part of the Main Ethiopian Rift. Most analyzed waters showed fluoride concentrations above the drinking water limit (> 1.5 mg/L) with the highest values in hot springs (up to 70 mg/L) and the geothermal well (76 mg/L).
In the solid phase, a high fluoride content was found in acid volcanic rocks (ignimbrite: 4391 ppm; rhyolite: 3248 ppm) as well as in pumice (up to 1955 ppm). The fluoride content of soil samples collected within the volcanic complex varied between 82 and 1036 ppm, whereas former lake sediments from outside the Aluto Volcanic Complex showed higher fluoride contents ranging from 674 to 8171 ppm. Identified fluoride-rich minerals are various amphiboles (about 3 wt.-% F-), fluor-apatite (about 3 wt % F-), minerals of the fluorcaphite group (about 5 wt.-% F-), parisite (up to 9 wt.-% F-), and fluorite (CaF2).
Elevated fluoride concentrations were also measured in some gas samples from fumaroles (up to 50 ppm) and in plant samples collected next to the fumaroles (up to 65 ppm). Leaching experiments of solid samples with deionized water demonstrated that fluoride can easily be mobilized from pumice and sediments but hardly from consolidated rocks.
This fluoride release increased with temperature (up to 150 C) and correlated roughly with dis-solved silica indicating the binding of some fluoride to the amorphous or weakly crystalline silica fraction. Based on these results it was concluded that fluoride migrates via different pathways through the environment: At the depth of the magma chamber during magmatic differentiation fluoride enriches initially in the magmatic melt and accumulates in some late-crystallizing minerals of igneous rocks such as fluorite or mica. Upon volcanic eruption fluoride is predominantly incorporated in the glass (ignimbrite) and ash phase (pumice).
On the surface, these rocks are exposed to weathering and fluoride leaches partly out into the aqueous phase. Soft and porous rocks such as pumice release fluoride first whereas hard extrusive/igneous rocks are less prone to weathering and retain the fluoride. Pumice and (lake) sediments might act both, as source and as sink for fluoride in the area.
Although some fluoride might drain from the surface into the deeper subsoil, we conclude that magmatic fluids (liquid and gas) contribute more substantially to the overall fluoride ground water content because (i) of the much higher fluoride content in deep geothermal waters and hot springs as compared to surface-near waters; (ii) active geothermal surface manifestations located along fault zones indicate that fluids migrate from deep magmatic intrusions (as gas and liquid) towards the surface, where fluoride dissolves in groundwater; (iii) and the good correlation between bicarbonate (deriving from dissolution of magmatic CO2) and fluoride content in all analyzed water samples.
Forward stratigraphic modelling is a fast-developing modelling approach, used to test conceptual models, and predict stratigraphic architecture and depositional facies from basin to reservoir scales.
Published subsurface applications demonstrate its added value by integrating multidisciplinary data as well as geological concepts into its constraints.
When applied to carbonate depositional systems, composed of multiple sediment factories, the cooperating and interdependent production mechanisms remain poorly studied.
By applying the technique to a well-studied section of the Maldives carbonate platform, a specific model design-adapted to the geological age and setting, and constrained by available data-sheds light on the interaction of its carbonate producers.
The results yield a naturalistic depositional facies distribution and offer insight in the changing relationship between biotic communities during the platform evolution.
After calibration, the reference model unequivocally links the formerly proposed genetic model to the seismostratigraphic architecture.
Furthermore, the results show how environmental changes (seemingly of secondary impact compared to changes in physical accommodation in the stratigraphic record) can induce substantial fluctuations in carbonate production rates of biotic communities, affect the ecological accommodation, and thus impact the platform architecture. Therefore, it is crucial to treat carbonate production rates during periods of environmental change as variables with associated uncertainties in a forward stratigraphic model setup.
Floodplain wetlands are critical for sustaining various ecological and hydrological functions in a riverine environment. Severe anthropogenic alterations and human occupation of floodplains have threatened these wetlands in several parts of the world. A major handicap in designing sustainable restoration and monitoring strategies for these wetlands is the lack of scientific process-based understanding and information on the basin-scale controls of their degradation. Here, we offer a novel approach to integrate the connectivity of the wetlands with the surrounding landscape along with other attributes such as stream density, hydrometeorological parameters, and groundwater dynamics to explain their degradation and then to prioritise them for restoration and monitoring.
We hypothesise that the best possible connectivity scenario for the existence of a wetland would be if (a) the wetland has a high connectivity with its upslope area, and (b) the wetland has a low connectivity with its downslope region.
The first condition ensures the flow of water into the wetland and the second condition allows longer water residence time in the wetland. Accordingly, we define four connectivity-based wetland health scenarios-good, no impact, bad, and worst.
We have implemented the proposed method in 3226 wetlands in the Ramganga Basin in north India. Further, we have applied specific selection criteria, such as distance from the nearest stream and stream density, to prioritise the wetlands for restoration and monitoring.
We conclude that the connectivity analysis offers a quick process-based assessment of wetlands' health status and serves as an important criterion to prioritise the wetlands for developing appropriate management strategies.
This research demonstrated the application of hydrochemical data and stable water isotopes of delta O-18 and delta D (or delta 2H) for evaluating the relationship between surface water in Lake Bosumtwi and the underlying groundwater system.
It aimed at determining the presence or absence of a hydraulic relationship and for evaluating the possible direction of flow at the interface between the two reservoirs. The study also estimated evaporative losses of infiltrating rainwater as it transits the unsaturated zone and provided important information on the hydrological processes prevailing in the area.
The results of Q-Mode hierarchical cluster analysis (HCA) clearly differentiate the lake water from the groundwater based on their spatial relationship.
These results indicated that groundwater recharge occurs on the hilltops of the crater, where it is slightly acidic with low levels of dissolved ions, characterised by short residence time and rapid unrestricted vertical infiltration and recharge.
The groundwater becomes more mineralized with longer contact times and deeper circulation with the host rock, while it flows from the recharge areas towards the lake at lower elevations.
Analyses of delta O-18 and delta D showed a high evaporation rate on the lake surface (90%) with a significant evaporative enrichment, whereas groundwater showed no significant isotopic variations.
Thus suggesting that the aquifers have been recharged by recent meteoric water that has undergone some evaporative enrichment since the study established an evaporation rate of water infiltrating the unsaturated zone ranging from 54 to 60%. Both reservoirs do not appear to be hydraulically connected, and where such a connection exists, it is expected to favour the lake.
A circular, single-contig Methanobacterium sp. metagenome-assembled genome (MAG) was recovered from high-CO2 enrichments inoculated with drill core material from the tectonic Eger Rift terrestrial subsurface.
Annotation of the recovered MAG highlighted putative methanogenesis genes, providing valuable information on archaeal activity in the deep biosphere.
Refractive index provides fundamental insights into the electronic structure of materials. At high pressure, however, the determination of refractive index and its wavelength dispersion is challenging, which limits our understanding of how physical properties of even simple materials, such as MgO, evolve with pressure.
Here, we report on the measurement of room-temperature refractive index of MgO up to similar to 140 GPa.
The refractive index of MgO at 600 nm decreases by similar to 2.4% from similar to 1.737 at 1 atm to similar to 1.696 (+/- 0.017) at similar to 140 GPa.
Despite the index at 600 nm is essentially pressure independent, the absolute wavelength dispersion of the refractive index at 550-870 nm decreases by similar to 28% from similar to 0.015 at 1 atm to similar to 0.011 (+/- 8.04 x 10(-4)) at similar to 103 GPa.
Single-effective-oscillator analysis of our refractive index data suggests that the bandgap of MgO increases by similar to 1.1 eV from 7.4 eV at 1 atm to similar to 8.5 (+/- 0.6) eV at similar to 103 GPa.
The Gram-negative bacterium Paenalcaligenes niemegkensis NGK35(T) was isolated from plastic debris in an abandoned landfill.
It has the ability to grow on polyethylene and hexadecane as the sole carbon sources. Here, we report the corresponding draft genome, which contains 3.66 Mbp and is characterized by a G+C content of 52.1%.
International migration patterns, at the global level, can to a large extent be explained through economic factors in origin and destination countries.
On the other hand, it has been shown that global climate change is likely to affect economic development over the coming decades.
Here, we demonstrate how these future climate impacts on national income levels could alter the global migration landscape.
Using an empirically calibrated global migration model, we investigate two separate mechanisms. The first is through destination-country income, which has been shown consistently to have a positive effect on immigration.
As countries' income levels relative to each other are projected to change in the future both due to different rates of economic growth and due to different levels of climate change impacts, the relative distribution of immigration across destination countries also changes as a result, all else being equal.
Second, emigration rates have been found to have a complex, inverted U-shaped dependence on origin-country income.
Given the available migration flow data, it is unclear whether this dependence-found in spatio-temporal panel data-also pertains to changes in a given migration flow over time. If it does, then climate change will additionally affect migration patterns through origin countries' emigration rates, as the relative and absolute positions of countries on the migration "hump" change.
We illustrate these different possibilities, and the corresponding effects of 3 degrees C global warming (above pre-industrial) on global migration patterns, using climate model projections and two different methods for estimating climate change effects on macroeconomic development.
Seismicity models are probabilistic forecasts of earthquake rates to support seismic hazard assessment.
Physics-based models allow extrapolating previously unsampled parameter ranges and enable conclusions on underlying tectonic or human-induced processes.
The Coulomb Failure (CF) and the rate-and-state (RS) models are two widely used physics-based seismicity models both assuming pre-existing populations of faults responding to Coulomb stress changes.
The CF model depends on the absolute Coulomb stress and assumes instantaneous triggering if stress exceeds a threshold, while the RS model only depends on stress changes.
Both models can predict background earthquake rates and time-dependent stress effects, but the RS model with its three independent parameters can additionally explain delayed aftershock triggering.
This study introduces a modified CF model where the instantaneous triggering is replaced by a mean time-to-failure depending on the absolute stress value.
For the specific choice of an exponential dependence on stress and a stationary initial seismicity rate, we show that the model leads to identical results as the RS model and reproduces the Omori-Utsu relation for aftershock decays as well stress-shadowing effects.
Thus, both CF and RS models can be seen as special cases of the new model. However, the new stress response model can also account for subcritical initial stress conditions and alternative functions of the mean time-to-failure depending on the problem and fracture mode.
The most profound consequences of the presence of Ca-Mg carbonates (CaCO3-MgCO3) in the Earth's upper mantle may be to lower the melting temperatures of the mantle and control the melt composition.
Low-degree partial melting of a carbonate-bearing mantle produces CO2-rich, silica-poor melts compositionally imposed by the melting relations of carbonates.
Thus, understanding the melting relations in the CaCO3-MgCO3 system facilitates the interpretation of natural carbonate-bearing silicate systems.
We report the melting relations of the CaCO3-MgCO3 system and the partition coefficient of trace elements between carbonates and carbonate melt from experiments at high pressure (6 and 9 GPa) and temperature (1300-1800 degrees C) using a rocking multi-anvil press. In the absence of water, Ca-Mg carbonates are stable along geothermal gradients typical of subducting slabs.
Ca-Mg carbonates ( similar to Mg0.1-0.9Ca0.9-0.1CO3) partially melt beneath mid-ocean ridges and in plume settings. Ca-Mg carbonates melt incongruently, forming periclase crystals and carbonate melt between 4 and 9 GPa.
Furthermore, we show that the rare earth element (REE) signature of Group-I kimberlites, namely strong REE fractionation and depletion of heavy REE relative to the primitive mantle, is resembled by carbonate melt in equilibrium with Ca-bearing magnesite and periclase at 6 and 9 GPa.
This suggests that the dolomite-magnesite join of the CaCO3-MgCO3 system might be useful to approximate the REE signature of carbonate-rich melts parental to kimberlites.
A metagenome-assembled genome (MAG), named Methanosarcina sp. strain ERenArc_MAG2, was obtained from a 3-month-old H-2/CO2 atmosphere enrichment culture, originally inoculated with 60-m deep drill core sediment collected from the tectonic Eger Rift terrestrial subsurface.
Annotation of the recovered draft genome revealed putative archaeal methanogenesis genes in the deep biosphere.
A metagenome-assembled genome (MAG), named Methanosarcina sp. strain ERenArc_MAG2, was obtained from a 3-month-old H-2/CO2 atmosphere enrichment culture, originally inoculated with 60-m deep drill core sediment collected from the tectonic Eger Rift terrestrial subsurface.
Annotation of the recovered draft genome revealed putative archaeal methanogenesis genes in the deep biosphere.
Information regarding the spatial distribution of soil water content is key in many disciplines and applications including soil and atmospheric sciences, hydrology, and agricultural engineering.
Thus, within the past decades various experimental methods and strategies have been developed to map spatial variations in soil moisture distribution and to monitor temporal changes.
Our study examines the combination of electrical resistivity mapping and point observations of soil moisture to infer the spatial and the temporal variability of soil moisture.
Over a period of around two years, we performed field measurements on six days to collect repeated electrical resistivity mapping data for a nine-hectare test site south-east of Berlin, Germany.
Permanently installed TDR probes, temporary TDR measurements within permanently installed tubes, and gravimetric measurements using soil samples provided soil moisture data at various selected points.
In addition, soil analysis and classification results are available for 132 regularly distributed positions up to depths of 1.2 m.
We compare and link three-dimensional resistivity models obtained via data inversion to soil composition and soil moisture as provided by our point data.
Both the soil samples and the resistivity models indicate a two-layer medium characterized by a sandy top layer with varying thickness and a loamy bottom soil.
For all six field campaigns, we observe similar resistivity patterns reflecting the temporally stable influence of soil texture.
While the overall patterns are stable, the range of resistivity values changes with soil moisture. Finally, to estimate spatial models of soil moisture, we link our soil moisture and resistivity data using empirical petrophysical models relying on a second order polynomial function.
We observe a mean prediction error for soil moisture of +/-0.034 m3 & BULL; m? 3 using all observation points while we notice that point-specific models further reduce the error.
Thus, we conclude that our experimental and data analysis strategies represent a reliable approach to establish site-specific models and to estimate three-dimensional moisture distribution including temporal variations.
Plate tectonic boundaries constitute the suture zones between tectonic plates. They are shaped by a variety of distinct and interrelated processes and play a key role in geohazards and georesource formation. Many of these processes have been previously studied, while many others remain unaddressed or undiscovered. In this work, the geodynamic numerical modeling software ASPECT is applied to shed light on further process interactions at continental plate boundaries. In contrast to natural data, geodynamic modeling has the advantage that processes can be directly quantified and that all parameters can be analyzed over the entire evolution of a structure. Furthermore, processes and interactions can be singled out from complex settings because the modeler has full control over all of the parameters involved. To account for the simplifying character of models in general, I have chosen to study generic geological settings with a focus on the processes and interactions rather than precisely reconstructing a specific region of the Earth.
In Chapter 2, 2D models of continental rifts with different crustal thicknesses between 20 and 50 km and extension velocities in the range of 0.5-10 mm/yr are used to obtain a speed limit for the thermal steady-state assumption, commonly employed to address the temperature fields of continental rifts worldwide. Because the tectonic deformation from ongoing rifting outpaces heat conduction, the temperature field is not in equilibrium, but is characterized by a transient, tectonically-induced heat flow signal. As a result, I find that isotherm depths of the geodynamic evolution models are shallower than a temperature distribution in equilibrium would suggest. This is particularly important for deep isotherms and narrow rifts. In narrow rifts, the magnitude of the transient temperature signal limits a well-founded applicability of the thermal steady-state assumption to extension velocities of 0.5-2 mm/yr. Estimation of the crustal temperature field affects conclusions on all temperature-dependent processes ranging from mineral assemblages to the feasible exploitation of a geothermal reservoir.
In Chapter 3, I model the interactions of different rheologies with the kinematics of folding and faulting using the example of fault-propagation folds in the Andean fold-and-thrust belt. The evolution of the velocity fields from geodynamic models are compared with those from trishear models of the same structure. While the latter use only geometric and kinematic constraints of the main fault, the geodynamic models capture viscous, plastic, and elastic deformation in the entire model domain. I find that both models work equally well for early, and thus relatively simple stages of folding and faulting, while results differ for more complex situations where off-fault deformation and secondary faulting are present. As fault-propagation folds can play an important role in the formation of reservoirs, knowledge of fluid pathways, for example via fractures and faults, is crucial for their characterization.
Chapter 4 deals with a bending transform fault and the interconnections between tectonics and surface processes. In particular, the tectonic evolution of the Dead Sea Fault is addressed where a releasing bend forms the Dead Sea pull-apart basin, while a restraining bend further to the North resulted in the formation of the Lebanese mountains. I ran 3D coupled geodynamic and surface evolution models that included both types of bends in a single setup. I tested various randomized initial strain distributions, showing that basin asymmetry is a consequence of strain localization. Furthermore, by varying the surface process efficiency, I find that the deposition of sediment in the pull-apart basin not only controls basin depth, but also results in a crustal flow component that increases uplift at the restraining bend.
Finally, in Chapter 5, I present the computational basis for adding further complexity to plate boundary models in ASPECT with the implementation of earthquake-like behavior using the rate-and-state friction framework. Despite earthquakes happening on a relatively small time scale, there are many interactions between the seismic cycle and the long time spans of other geodynamic processes. Amongst others, the crustal state of stress as well as the presence of fluids or changes in temperature may alter the frictional behavior of a fault segment. My work provides the basis for a realistic setup of involved structures and processes, which is therefore important to obtain a meaningful estimate for earthquake hazards.
While these findings improve our understanding of continental plate boundaries, further development of geodynamic software may help to reveal even more processes and interactions in the future.
Tuning systems of traditional Georgian singing determined from a new corpus of field recordings
(2022)
In this study we examine the tonal organization of the 2016 GVM dataset, a newly-created corpus of high-quality multimedia field recordings of traditional Georgian singing with a focus on Svaneti. For this purpose, we developed a new processing pipeline for the computational analysis of non-western polyphonic music which was subsequently applied to the complete 2016 GVM dataset.
To evaluate under what conditions a single tuning system is representative of current Svan performance practice, we examined the stability of the obtained tuning systems from an ensemble-, a song-, and a corpus-related perspective.
Furthermore, we compared the resulting Svan tuning systems with the tuning systems obtained for the Erkomaishvili dataset (Rosenzweig et al., 2020) in the study by Scherbaum et al. (2020). In comparison to a 12-TET (12-tone-equal-temperament) system, the Erkomaishvili and the Svan tuning systems are surprisingly similar.
Both systems show a strong presence of pure fourths (500 cents) and fifths (700 cents), and 'neutral' thirds (peaking around 350 cents) as well as 'neutral' sixths.
In addition, the sizes of the melodic and the harmonic seconds in both tuning systems differ systematically from each other, with the size of the harmonic second being systematically larger than the melodic one.
We examine sub-caldera structure and seismicity in and around the Okataina Volcanic Centre (OVC), Taupo Volcanic Zone, New Zealand, using seismic data collected over 4+ years with a temporary array of broadband and short-period seismometers, supplementing data from the permanent New Zealand seismometer network.
We derive a new 3-D image of P-wave seismic velocity for the upper crust in the region, using double-difference seismic tomography and utilising waveform cross-correlations. We subsequently relocate 6989 earthquakes which occurred in the region over the 2010-2021 time period, using the 3D velocity model.
The seismicity distribution shows spatial clusters west of Lake Rotomahana, as well as beneath Haroharo and the Makatiti Dome, inside the Okataina caldera. Beneath Makatiti Dome 90% of the events are shallower than 7.7 +/- 0.5 km.
Outside of the Okataina caldera event relocations highlight short (similar to 3-4 km long) streaks of activity in the Ngakuru graben, part of the active Taupo Rift southwest of Okataina caldera. Inside the OVC the relocated seismicity beneath Makatiti appears closely associated with low (similar to 10%) P-wave velocity anomalies, which we resolve in the similar to 5-to-8-km depth range beneath the Okataina caldera, and which are likely related to partial melt and/or fluid-volatile pathways.
Moment tensor analyses for two larger-magnitude events (M(L)4.5 and M(L)4.9) near Haroharo indicate normal faulting, with NNE-SSW fault strike, but with positive CLVD and positive isotropic components when allowing for a full moment tensor, consistent with a magmatic environment with degassing and/or fluid migration.
An accurate estimation of crop yield under climate change scenarios is essential to quantify our ability to feed a growing population and develop agronomic adaptations to meet future food demand. A coordinated evaluation of yield simulations from process-based eco-physiological models for climate change impact assessment is still missing for soybean, the most widely grown grain legume and the main source of protein in our food chain. In this first soybean multi-model study, we used ten prominent models capable of simulating soybean yield under varying temperature and atmospheric CO2 concentration [CO2] to quantify the uncertainty in soybean yield simulations in response to these factors. Models were first parametrized with high quality measured data from five contrasting environments. We found considerable variability among models in simulated yield responses to increasing temperature and [CO2]. For example, under a + 3 degrees C temperature rise in our coolest location in Argentina, some models simulated that yield would reduce as much as 24%, while others simulated yield increases up to 29%. In our warmest location in Brazil, the models simulated a yield reduction ranging from a 38% decrease under + 3 degrees C temperature rise to no effect on yield. Similarly, when increasing [CO2] from 360 to 540 ppm, the models simulated a yield increase that ranged from 6% to 31%. Model calibration did not reduce variability across models but had an unexpected effect on modifying yield responses to temperature for some of the models. The high uncertainty in model responses indicates the limited applicability of individual models for climate change food projections. However, the ensemble mean of simulations across models was an effective tool to reduce the high uncertainty in soybean yield simulations associated with individual models and their parametrization. Ensemble mean yield responses to temperature and [CO2] were similar to those reported from the literature. Our study is the first demonstration of the benefits achieved from using an ensemble of grain legume models for climate change food projections, and highlights that further soybean model development with experiments under elevated [CO2] and temperature is needed to reduce the uncertainty from the individual models.
Geodetic studies of crustal deformation using Global Navigation Satellite System (GNSS, earlier commonly referred to as Global Positioning System, GPS) measurements at CSIR-NGRI started in 1995 with the installation of a permanent GNSS station at CSIR-NGRI Hyderabad which later became an International GNSS Service (IGS) site. The CSIR-NGRI started expanding its GNSS networks after 2003 with more focussed studies through installation in the NE India, Himalayan arc, Andaman subduction zone, stable and failed rift regions of India plate. In each instance, these measurements helped in unravelling the geodynamics of the region and seismic hazard assessment, e.g., the discovery of a plate boundary fault in the Indo-Burmese wedge, rate and mode of strain accumulation and its spatial variation in the Garhwal-Kumaun and Kashmir region of the Himalayan arc, the influence of non-tectonic deformation on tectonic deformation in the Himalayan arc, nature of crustal deformation through earthquake cycle in the Andaman Sumatra subduction zone, and localised deformation in the intraplate region and across the paleo rift regions. Besides these, GNSS measurements initiated in the Antarctica region have helped in understanding the plate motion and influence of seasonal variations on deformation. Another important by-product of the GNSS observations is the capabilities of these observations in understanding the ionospheric variations due to earthquake processes and also due to solar eclipse. We summarize these outcomes in this article.
Climatic and hydrological variability as a driver of the Lake Gościąż biota during the Younger Dryas
(2022)
The Younger Dryas (YD) is a roughly 1,100-year cold period marking the end of the last glaciation. Climate modelling for northern Europe indicates high summer temperatures and strong continentality. In eastern Europe, the scale of temperature variation and its influence on ecosystems is weakly recognised. Here, we present a multi-proxy reconstruction of YD conditions from Lake Gos ' ciaz (central Poland). The decadal-resolution analysis of its annually varved sediments indicates an initial decrease in Chironomidae-inferred mean July air temperature followed by steady warming. The pollen-inferred winter-to-summer temperature amplitude and annual precip-itation is highest at the Allerod/YD transition and the early YD (ca. 12.7-12.4 ky cal BP) and YD/Holocene (11.7-11.4 ka cal BP) transition. Temperature and precipitation were the main reasons for lake level fluctuations as reflected in the planktonic/littoral Cladocera ratio. The lake's diatom-inferred total phosphorus decreased with increasing summer temperature from about mid YD. Windy conditions in the early YD until ~12.3 ka cal BP caused water mixing and a short-lived/temporary increase in nutrient availability for phytoplankton. The Chironomidae-inferred summer temperature and pollen inferred summer temperature, winter temperature and annual precipitation herein are one of only a few in eastern Europe conducted with such high resolution.
Understanding eastern African paleoclimate is critical for contextualizing early human evolution, adaptation, and dispersal, yet Pleistocene climate of this region and its governing mechanisms remain poorly understood due to the lack of long, orbitally-resolved, terrestrial paleoclimate records. Here we present leaf wax hydrogen isotope records of rainfall from paleolake sediment cores from key time windows that resolve long-term trends, variations, and high-latitude effects on tropical African precipitation. Eastern African rainfall was dominantly controlled by variations in low-latitude summer insolation during most of the early and middle Pleistocene, with little evidence that glacial-interglacial cycles impacted rainfall until the late Pleistocene. We observe the influence of high-latitude-driven climate processes emerging from the last interglacial (Marine Isotope Stage 5) to the present, an interval when glacial-interglacial cycles were strong and insolation forcing was weak. Our results demonstrate a variable response of eastern African rainfall to low-latitude insolation forcing and high-latitude-driven climate change, likely related to the relative strengths of these forcings through time and a threshold in monsoon sensitivity. We observe little difference in mean rainfall between the early, middle, and late Pleistocene, which suggests that orbitally-driven climate variations likely played a more significant role than gradual change in the relationship between early humans and their environment.
The immense advances in computer power achieved in the last decades have had a significant impact in Earth science, providing valuable research outputs that allow the simulation of complex natural processes and systems, and generating improved forecasts. The development and implementation of innovative geoscientific software is currently evolving towards a sustainable and efficient development by integrating models of different aspects of the Earth system. This will set the foundation for a future digital twin of the Earth. The codification and update of this software require great effort from research groups and therefore, it needs to be preserved for its reuse by future generations of geoscientists. Here, we report on Geo-Soft-CoRe, a Geoscientific Software & Code Repository, hosted at the archive DIGITAL.CSIC. This is an open source, multidisciplinary and multiscale collection of software and code developed to analyze different aspects of the Earth system, encompassing tools to: 1) analyze climate variability; 2) assess hazards, and 3) characterize the structure and dynamics of the solid Earth. Due to the broad range of applications of these software packages, this collection is useful not only for basic research in Earth science, but also for applied research and educational purposes, reducing the gap between the geosciences and the society. By providing each software and code with a permanent identifier (DOI), we ensure its self-sustainability and accomplish the FAIR (Findable, Accessible, Interoperable and Reusable) principles. Therefore, we aim for a more transparent science, transferring knowledge in an easier way to the geoscience community, and encouraging an integrated use of computational infrastructure.
The martensitic transformation is a fundamental physical phenomenon at the origin of important industrial applications. However, the underlying microscopic mechanism, which is of critical importance to explain the outstanding mechanical properties of martensitic materials, is still not fully understood. This is because for most martensitic materials the transformation is a fast process that makes in situ studies extremely challenging. Noble solids krypton and xenon undergo a progressive pressure-induced face-centered cubic (fcc) to hexagonal close-packed (hcp) martensitic transition with a very wide coexistence domain. Here, we took advantage of this unique feature to study the detailed transformation progress at the atomic level by employing in situ x-ray diffraction and absorption spectroscopy. We evidenced a four-stage pathway and suggest that the lattice mismatch between the fcc and hcp forms plays a key role in the generation of strain. We also determined precisely the effect of the transformation on the compression behavior of these materials.