Refine
Year of publication
Document Type
- Article (2809) (remove)
Keywords
- Holocene (35)
- climate change (31)
- Earthquake source observations (21)
- Climate change (18)
- erosion (18)
- permafrost (18)
- Himalaya (17)
- Pollen (17)
- Tibetan Plateau (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2809) (remove)
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.
Raman spectroscopic quantification of tetrahedral boron in synthetic aluminum-rich tourmaline
(2021)
The Raman spectra of five B-[4]-bearing tourmalines of different composition synthesized at 700 degrees C/4.0 GPa (including first-time synthesis of Na-Li-B-[4]-tourmaline, Ca-Li-B-[4]-tourmaline, and Ca-bearing square-B-[4]-tourmaline) reveal a strong correlation between the tetrahedral boron content and the summed relative intensity of all OH-stretching bands between 3300-3430 cm(-1). The band shift to low wavenumbers is explained by strong O3-H center dot center dot center dot O5 hydrogen bridge bonding. Applying the regression equation to natural B-[4]-bearing tourmaline from the Koralpe (Austria) reproduces the EMPA-derived value perfectly [EMPA: 0.67(12) B-[4] pfu vs. Raman: 0.66(13) B-[4] pfu]. This demonstrates that Raman spectroscopy provides a fast and easy-to-use tool for the quantification of tetrahedral boron in tourmaline. The knowledge of the amount of tetrahedral boron in tourmaline has important implications for the better understanding and modeling of B-isotope fractionation between tourmaline and fluid/melt, widely used as a tracer of mass transfer processes.
In this study, we investigate numerically the hydro-mechanical behavior of fractured crystalline rock due to one of the five hydraulic stimulations at the Pohang Enhanced Geothermal site in South Korea. We use the commercial code FracMan (Golder Associates) that enables studying hydro-mechanical coupled processes in fractured media in three dimensions combining the finite element method with a discrete fracture network. The software is used to simulate fluid pressure perturbation at fractures during hydraulic stimulation. Our numerical simulation shows that pressure history matching can be obtained by partitioning the treatment into separate phases. This results in adjusted stress-aperture relationships. The evolution of aperture adjustment implies that the stimulation mechanism could be a combination of hydraulic fracturing and shearing. The simulated extent of the 0.01 MPa overpressure contour at the end of the treatment equals to similar to 180 m around the injection point.
Pedogenic carbonate is widespread at mid latitudes where warm and dry conditions favor soil carbonate growth from spring to fall. The mechanisms and timing of pedogenic carbonate formation are more ambiguous in the tropical domain, where long periods of soil water saturation and high soil respiration enhance calcite dissolution. This paper provides stable carbon, oxygen and clumped isotope values from Quaternary and Miocene pedogenic carbonates in the tropical domain of Myanmar, in areas characterized by warm (>18°C) winters and annual rainfall up to 1,700 mm. We show that carbonate growth in Myanmar is delayed to the driest and coldest months of the year by sustained monsoonal rainfall from mid spring to late fall. The range of isotopic variability in Quaternary pedogenic carbonates can be solely explained by temporal changes of carbonate growth within the dry season, from winter to early spring. We propose that high soil moisture year-round in the tropical domain narrows carbonate growth to the driest months and makes it particularly sensitive to the seasonal distribution of rainfall. This sensitivity is also enabled by high winter temperatures, allowing carbonate growth to occur outside the warmest months of the year. This high sensitivity is expected to be more prominent in the geological record during times with higher temperatures and greater expansion of the tropical realm. Clumped isotope temperatures, δ13C and δ18O values of tropical pedogenic carbonates are impacted by changes of both rainfall seasonality and surface temperatures; this sensitivity can potentially be used to track past tropical rainfall distribution.
The Gutenberg-Richter (GR) and the Omori-Utsu (OU) law describe the earthquakes' energy release and temporal clustering and are thus of great importance for seismic hazard assessment. Motivated by experimental results, which indicate stress-dependent parameters, we consider a combined global data set of 127 main shock-aftershock sequences and perform a systematic study of the relationship between main shock-induced stress changes and associated seismicity patterns. For this purpose, we calculate space-dependent Coulomb Stress (& UDelta;CFS) and alternative receiver-independent stress metrics in the surrounding of the main shocks. Our results indicate a clear positive correlation between the GR b-value and the induced stress, contrasting expectations from laboratory experiments and suggesting a crucial role of structural heterogeneity and strength variations. Furthermore, we demonstrate that the aftershock productivity increases nonlinearly with stress, while the OU parameters c and p systematically decrease for increasing stress changes. Our partly unexpected findings can have an important impact on future estimations of the aftershock hazard.
Large earthquakes are usually modeled with simple planar fault surfaces or a combination of several planar fault segments. However, in general, earthquakes occur on faults that are non-planar and exhibit significant geometrical variations in both the along-strike and down-dip directions at all spatial scales. Mapping of surface fault ruptures and high-resolution geodetic observations are increasingly revealing complex fault geometries near the surface and accurate locations of aftershocks often indicate geometrical complexities at depth. With better geodetic data and observations of fault ruptures, more details of complex fault geometries can be estimated resulting in more realistic fault models of large earthquakes. To address this topic, we here parametrize non-planar fault geometries with a set of polynomial parameters that allow for both along-strike and down-dip variations in the fault geometry. Our methodology uses Bayesian inference to estimate the non-planar fault parameters from geodetic data, yielding an ensemble of plausible models that characterize the uncertainties of the non-planar fault geometry and the fault slip. The method is demonstrated using synthetic tests considering slip spatially distributed on a single continuous finite non-planar fault surface with varying dip and strike angles both in the down-dip and along-strike directions. The results show that fault-slip estimations can be biased when a simple planar fault geometry is assumed in presence of significant non-planar geometrical variations. Our method can help to model earthquake fault sources in a more realistic way and may be extended to include multiple non-planar fault segments or other geometrical fault complexities.
Earthquake source parameters such as seismic stress drop and corner frequency are observed to vary widely, leading to persistent discussion on potential scaling of stress drop and event size. Physical mechanisms that govern stress drop variations arc difficult to evaluate in nature and are more readily studied in controlled laboratory experiments. We perform two stick-slip experiments on fractured (rough) and cut (smooth) Westerly granite samples to explore fault roughness effects on acoustic emission (AE) source parameters. We separate large stick-slip events that generally saturate the seismic recording system from populations of smaller AE events which are sensitive to fault stresses prior to slip. AE event populations show many similarities to natural seismicity and may be interpreted as laboratory equivalent of natural microseismic events. We then compare the temporal evolution of mechanical data such as measured stress release during slip to temporal changes in stress drops derived from Alis using the spectral ratio technique. We report on two primary observations: (1) In contrast to most case studies for natural earthquakes, we observe a strong increase in seismic stress drop with AE size. (2) The scaling of stress drop with magnitude is governed by fault roughness, whereby the rough fault shows a more rapid increase of the stress drop magnitude relation with progressing large stick-slip events than the smooth fault. The overall range of AE sizes on the rough surface is influenced by both the average grain size and the width of the fault core. The magnitudes of the smallest AE events on smooth faults may also be governed by grain size. However, AEs significantly grow beyond peak roughness and the width of the fault core. Our laboratory tests highlight that source parameters vary substantially in the presence of fault zone heterogeneity (i.e. roughness and narrow grain size distribution), which may affect seismic energy partitioning and static stress drops of small and large AE events.
The West Burma Terrane (WBT) is a small terrane bounded to the east by the Asian Sibumasu Block and to the west by the Indo-Burman Ranges (IBR), the latter being an exhumed accretionary prism that formed during subduction of Indian oceanic lithosphere beneath Asia. Understanding the geological history of the WBT is important for reconstruction of the closure history of the Tethys Ocean and India-Asia collision. Currently there are major discrepancies in the proposed timings of collision between the WBT with both India and Asia; whether the WBT collided with India or Asia first is debated, and proposed timings of collisions stretch from the Mesozoic to the Cenozoic. We undertook a multi-technique provenance study involving petrography, detrital zircon U-Pb and Hf analyses, rutile U-Pb analyses and Sr-Nd bulk rock analyses on sediments of the Central Myanmar Basins of the WBT. We determined that the first arrival of Asian material into the basin occurred after the earliest late Eocene and by the early Oligocene, thus placing a minimum constraint on the timing of WBT-Asia collision. Our low temperature thermochronological study of the IBR records two periods of exhumation, in the early-middle Eocene, and at the Oligo-Miocene boundary. The Eocene event may be associated with the collision of the WBT with India. The later event at the Oligo-Miocene boundary may be associated with changes in wedge dynamics resulting from increased sediment supply to the system; however a number of other possible causes provide equally plausible explanations for both events.
At the junction of greenhouse and icehouse climate states, the Eocene-Oligocene Transition (EOT) is a key moment in Cenozoic climate history. While it is associated with severe extinctions and biodiversity turnovers on land, the role of terrestrial climate evolution remains poorly resolved, especially the associated changes in seasonality. Some paleobotanical and geochemical continental records in parts of the Northern Hemisphere suggest the EOT is associated with a marked cooling in winter, leading to the development of more pronounced seasons (i.e., an increase in the mean annual range of temperature, MATR). However, the MATR increase has been barely studied by climate models and large uncertainties remain on its origin, geographical extent and impact. In order to better understand and describe temperature seasonality changes between the middle Eocene and the early Oligocene, we use the Earth system model IPSL-CM5A2 and a set of simulations reconstructing the EOT through three major climate forcings: pCO(2) decrease (1120, 840 and 560 ppm), the Antarctic ice-sheet (AIS) formation and the associated sea-level decrease. Our simulations suggest that pCO(2) lowering alone is not sufficient to explain the seasonality evolution described by the data through the EOT but rather that the combined effects of pCO(2) , AIS formation and increased continentality provide the best data-model agreement.pCO(2) decrease induces a zonal pattern with alternating increasing and decreasing seasonality bands particularly strong in the northern high latitudes (up to 8 degrees C MATR increase) due to sea-ice and surface albedo feedback. Conversely, the onset of the AIS is responsible for a more constant surface albedo yearly, which leads to a strong decrease in seasonality in the southern midlatitudes to high latitudes (> 40 degrees S). Finally, continental areas that emerged due to the sea-level lowering cause the largest increase in seasonality and explain most of the global heterogeneity in MATR changes (1MATR) patterns. The Delta MATR patterns we reconstruct are generally consistent with the variability of the EOT biotic crisis intensity across the Northern Hemisphere and provide insights on their underlying mechanisms.