ALBERT

Beschreibung: Carbonatites host some unique ore deposits, especially REE, and fractional crystallization might be a potentially powerful mechanism for control enrichment of carbonatitic magmas by these metals to economically significant levels. At present, data on distribution coefficients of REE during fractional crystallization of carbonatitic melts at volcanic conditions are extremely scarce. Here we present an experimental study of REE partitioning between carbonatitic melts and calcite in the system CaCO3-Na2CO3 with varying amounts of P2O5, F, Cl, SiO2, SO3 at 650–900 °C and 100 MPa using cold-seal pressure vessels and LA-ICP-MS. The presence of phosphorus in the system generally increases the distribution coefficients but its effect decreases with increasing concentration. The temperature factor is high: at 770–900 °C DREE ≥ 1, while at lower temperatures DREE become below unity. Silicon also promotes the fractionation of REE into calcite, while sulfur contributes to retention of REE in the melt. Our results imply that calcite may impose significant control upon REE fractionation at the early stages of crystallization of carbonatitic magmas and might be a closest proxy for monitoring the REE content in initial melt.

Beschreibung: Deep hydrothermal Mo, W and base metal mineralization found in the Detroit City portal of the Sweet Home mine in the Alma district (Colorado Mineral Belt) was deposited in response to magmatic activity and the formation of Climax-type Mo deposits during the Oligocene. This study presents extensive geochemical analyses of fluid inclusions inminerals from early greisen-like vein mineralization to better understand the fluid system responsible for ore formation. Quartz and fluorite, which are associated with molybdenite, huebnerite and/or pyrite mineralization, precipitated from low- to medium-salinity (1.5-11.5 wt.% equiv. NaCl), CO2-bearing fluids at temperatures between 360 and 415°C and probably under a fluctuating pressure regime at depths of at least 3.5 km. The formation of greisen-like and base metal mineralization at the Detroit City portal of the Sweet Home mine is related to fluids of different origin. Early magmatic fluids were the principal source for mantle-derived volatiles (CO2, H2S/SO2, noble gases) and mixed with significant amounts of heated meteoric water. Fluid mixing of magmatic fluids with meteoric water is constrained by δ 2H-δ 18O relationships of fluid inclusion water in different minerals (Fig. 1). Whether molybdenum was derived from magmatic fluids remains unclear. Fluid inclusions in huebnerite suggest that W originated from source rocks that are enriched in organic matter rather than from magmatic fluids. The deep hydrothermal mineralization at the Detroit City portal of the Sweet Home mine shows features similar to deep hydrothermal vein mineralization found in Climax-type Mo deposits and their periphery, suggesting that fluid migration and the deposition of ore and gangue minerals in the Sweet Home mine was triggered by a deep-seated magmatic intrusion. The findings of this study are in good agreement with the results of previous fluid inclusion studies of the mineralization of the Sweet Home mine (Lüders et al., 2009) and from Climaxtype Mo porphyry deposits in the Colorado Mineral Belt (Hall et al., 1974).

Beschreibung: The birth and expansion of continental plateaus exert a strong control on our planet's climate and the distribution and evolution of its biodiversity. It has been proposed that the Tibetan Plateau has been steadily growing by southward expansion. Here we demonstrate that the shape of the southeastern margin of the plateau has remained unchanged for the last 10 Myr despite vast amounts of exhumation. Our finding is based on a new, high-resolution thermochronological dataset from the deep gorges of the Salween and Mekong rivers, which we interpret using a physics-based model combined with an optimization method. We show that our scenario also agrees with a wide range of other, independent geological and geophysical data. This finding demonstrates that plateau margins can reach large-scale topographic steady-state between outward growth and surface erosion, which has important implications for our understanding of the evolution of Earth's climate and biodiversity in the recent geological past.

Beschreibung: By means of 4 years of Challenging Minisatellite Payload (CHAMP) zonal wind observations and a Thermosphere-Ionosphere Electrodynamics General Circulation Model simulation, longitudinal and seasonal variations of thermospheric superrotation at magnetic equator are investigated first. The superrotation shows four- and three-peaked longitudinal patterns in March and September equinoxes, and a two-peaked variation during solstices. The superrotation is stronger in December than in June solstice, and stronger in March than in September equinox. The mean correlation between the zonal variation of superrotation and nighttime eastward wind is about 0.8, while it is 0.6 with daytime westward zonal wind. The interaction between the neutral wind and geomagnetic field plays a more important role in the superrotation, rather than the F-region electron density. The lower atmospheric tides tend to suppress the superrotation, but contribute to the longitudinal patterns of the superrotation. The viscous force is also important for the longitudinal modulation of the superrotation.

Beschreibung: We report the concurrent observations of F-region plasma changes and field-aligned currents (FACs) above isolated proton auroras (IPAs) associated with electromagnetic ion cyclotron Pc1 waves. Key events on March 19, 2020 and September 12, 2018 show that ground magnetometers and all-sky imagers detected concurrent Pc1 wave and IPA, during which NOAA POES observed precipitating energetic protons. In the ionospheric F-layer above the IPA zone, the Swarm satellites observed transverse Pc1 waves, which span wider latitudes than IPA. Around IPA, Swarm also detected the bipolar FAC and localized plasma density enhancement, which is occasionally surrounded by wide/shallow depletion. This indicates that wave-induced proton precipitation contributes to the energy transfer from the magnetosphere to the ionosphere.

Beschreibung: The LArge-scale Reservoir Simulator (LARS) has been previously developed to study hydrate dissociation in hydrate-bearing systems under in-situ conditions. In the present study, a numerical framework of equations of state describing hydrate formation at equilibrium conditions has been elaborated and integrated with a numerical flow and transport simulator to investigate a multi-stage hydrate formation experiment undertaken in LARS. A verification of the implemented modeling framework has been carried out by benchmarking it against another established numerical code. Three-dimensional (3D) model calibration has been performed based on laboratory data available from temperature sensors, fluid sampling, and electrical resistivity tomography. The simulation results demonstrate that temperature profiles, spatial hydrate distribution, and bulk hydrate saturation are consistent with the observations. Furthermore, our numerical framework can be applied to calibrate geophysical measurements, optimize post-processing workflows for monitoring data, improve the design of hydrate formation experiments, and investigate the temporal evolution of sub-permafrost methane hydrate reservoirs.

Beschreibung: Tropical forests store 40–50 per cent of terrestrial vegetation carbon1. However, spatial variations in aboveground live tree biomass carbon (AGC) stocks remain poorly understood, in particular in tropical montane forests2. Owing to climatic and soil changes with increasing elevation3, AGC stocks are lower in tropical montane forests compared with lowland forests2. Here we assemble and analyse a dataset of structurally intact old-growth forests (AfriMont) spanning 44 montane sites in 12 African countries. We find that montane sites in the AfriMont plot network have a mean AGC stock of 149.4 megagrams of carbon per hectare (95% confidence interval 137.1–164.2), which is comparable to lowland forests in the African Tropical Rainforest Observation Network4 and about 70 per cent and 32 per cent higher than averages from plot networks in montane2,5,6 and lowland7 forests in the Neotropics, respectively. Notably, our results are two-thirds higher than the Intergovernmental Panel on Climate Change default values for these forests in Africa8. We find that the low stem density and high abundance of large trees of African lowland forests4 is mirrored in the montane forests sampled. This carbon store is endangered: we estimate that 0.8 million hectares of old-growth African montane forest have been lost since 2000. We provide country-specific montane forest AGC stock estimates modelled from our plot network to help to guide forest conservation and reforestation interventions. Our findings highlight the need for conserving these biodiverse9,10 and carbon-rich ecosystems.

Beschreibung: In this study, the effect of disturbance dynamo electric field (DDEF) induced by subauroral polarization streams (SAPS) on the variations of the equatorial electrojet (EEJ) and its counter electrojet (CEJ) during the geomagnetic storm on June 1, 2013 is analyzed in detail for the first time. Observations from ground-based magnetometers showed that the SAPS-induced EEJ flows westward and eastward in the daytime and dawn/dusk sectors, respectively. The effects of SAPS on EEJ are mainly associated with the changes of zonal ionospheric electric field, while the changes in the ionospheric conductivity play a secondary role. By using Thermosphere Ionosphere Electrodynamic General Circulation Model simulations, the zonal electric field induced by SAPS associated with the DDEF is examined. The results of the simulations show that the DDEF has a significant impact on the EEJ variability. The daytime westward EEJ at the dip equator is mainly driven by disturbance zonal wind, with secondary contributions from disturbance meridional wind. A similar mechanism can be observed in the dawn/dusk sector when the eastward EEJ is produced; however, it has a much weaker intensity than that during the daytime.

Beschreibung: CO2 release from particulate organic carbon (POC) oxidation during fluvial transit can influence climate over a range of timescales. Identifying the mechanistic controls on such carbon fluxes requires determining where POC oxidation occurs in river systems. While field data show POC oxidation and replacement moving downstream in lowland rivers, flume studies show that oxidation during active fluvial transport is limited. This suggests that most fluvial POC oxidation occurs during transient floodplain storage, but this idea has yet to be tested. Here, we isolate the influence of floodplain storage time on POC oxidation by exploiting a chronosequence of floodplain deposits above the modern groundwater table in the Rio Bermejo, Argentina. Measurements from 15 floodplain cores with depositional ages from 1 y to 20 ky show a progressive POC concentration decrease and 13C-enrichment with increasing time spent in floodplain storage. These results from the Rio Bermejo indicate that over 80% of fluvially-deposited POC can be oxidized over millennial timescales in aerated floodplains. Furthermore, POC in the oldest floodplain cores is more 14C-enriched than expected based on the independently-dated floodplain ages, indicating that a portion of this oxidized POC is replaced by autochthonous POC produced primarily by floodplain vegetation. We suggest floodplain storage timescales control the extent of oxidation of fluvially-deposited POC, and may play a prominent role in determining if rivers are significant atmospheric CO2 sources.

Beschreibung: Wuhan (China) is facing severe consolidation subsidence of soft soil and karst collapse hazards. To quantitatively explore the extent and causes of land subsidence in Wuhan, we performed multitemporal interferometry (MTI) analysis using synthetic aperture radar (SAR) data from the TerraSAR-X satellite from 2013 to 2017 and the Sentinel-1A satellite from 2015 to 2017. MTI results reveal four major subsidence zones in Wuhan, namely, Hankou (exceeding −6 cm/yr), Xudong-Qingshan (−3 cm/yr), Baishazhou-Jiangdi (−3 cm/yr), and Jianshe-Yangluo (−2 cm/yr). Accuracy assessment using 106 levelling benchmarks and cross-validation between the two InSAR-based results indicate an overall root-mean-square error (RMSE) of 2.5 and 3.1 mm/yr, respectively. Geophysical and geological analyses suggest that among the four major subsiding zones, Hankou, Xudong-Qingshan, and Jianshe-Yangluo are located in non-karstic soft soil areas, where shallow groundwater (〈 30 m) declines driven by engineering dewatering and industrial water depletion contribute directly to soft soil compaction. Subsidence in the Baishazhou-Jiangdi zone develops in the karst terrain with abundant underground caves and fissures, which are major natural factors for gradual subsidence and karst collapse. Spatial variation analysis of the geological conditions indicates that the stage of karst development plays the most important role in influencing kart subsidence, followed by municipal construction, proximity to major rivers, and overlying soil structure. Moreover, land subsidence in this zone is affected more via coupling effects from multiple factors. Risk zoning analysis integrating subsidence horizontal gradient, InSAR deformation rates, and municipal construction density show that the high-risk areas in Wuhan are mainly distributed in the Tianxingzhou and Baishazhou-Jiangdi zone, and generally spread along the metro lines.

Beschreibung: Geochemical processes such as mineral dissolution and precipitation alter the microstructure of rocks, and thereby affect their hydraulic and mechanical behaviour. Quantifying these property changes and considering them in reservoir simulations is essential for a sustainable utilisation of the geological subsurface. Due to the lack of alternatives, analytical methods and empirical relations are currently applied to estimate evolving hydraulic and mechanical rock properties associated with chemical reactions. However, the predictive capabilities of analytical approaches remain limited, since they assume idealised microstructures, and thus are not able to reflect property evolution for dynamic processes. Hence, aim of the present thesis is to improve the prediction of permeability and stiffness changes resulting from pore space alterations of reservoir sandstones. A detailed representation of rock microstructure, including the morphology and connectivity of pores, is essential to accurately determine physical rock properties. For that purpose, three-dimensional pore-scale models of typical reservoir sandstones, obtained from highly resolved micro-computed tomography (micro-CT), are used to numerically calculate permeability and stiffness. In order to adequately depict characteristic distributions of secondary minerals, the virtual samples are systematically altered and resulting trends among the geometric, hydraulic, and mechanical rock properties are quantified. It is demonstrated that the geochemical reaction regime controls the location of mineral precipitation within the pore space, and thereby crucially affects the permeability evolution. This emphasises the requirement of determining distinctive porosity-permeability relationships by means of digital pore-scale models. By contrast, a substantial impact of spatial alterations patterns on the stiffness evolution of reservoir sandstones are only observed in case of certain microstructures, such as highly porous granular rocks or sandstones comprising framework-supporting cementations. In order to construct synthetic granular samples a process-based approach is proposed including grain deposition and diagenetic cementation. It is demonstrated that the generated samples reliably represent the microstructural complexity of natural sandstones. Thereby, general limitations of imaging techniques can be overcome and various realisations of granular rocks can be flexibly produced. These can be further altered by virtual experiments, offering a fast and cost-effective way to examine the impact of precipitation, dissolution or fracturing on various petrophysical correlations. The presented research work provides methodological principles to quantify trends in permeability and stiffness resulting from geochemical processes. The calculated physical property relations are directly linked to pore-scale alterations, and thus have a higher accuracy than commonly applied analytical approaches. This will considerably improve the predictive capabilities of reservoir models, and is further relevant to assess and reduce potential risks, such as productivity or injectivity losses as well as reservoir compaction or fault reactivation. Hence, the proposed method is of paramount importance for a wide range of natural and engineered subsurface applications, including geothermal energy systems, hydrocarbon reservoirs, CO2 and energy storage as well as hydrothermal deposit exploration

Beschreibung: Geochemische Lösungs- und Fällungsprozesse verändern die Struktur des Porenraums und können dadurch die hydraulischen und mechanischen Gesteinseigenschaften erheblich beeinflussen. Die Quantifizierung dieser Parameteränderung und ihre Berücksichtigung in Reservoirmodellen ist entscheidend für eine nachhaltige Nutzung des geologischen Untergrunds. Aufgrund fehlender Alternativen werden dafür bisher analytische Methoden genutzt. Da diese Ansätze eine idealisierte Mikrostruktur annehmen, können insbesondere Änderungen der Gesteinseigenschaften infolge von dynamischen Prozessen nicht zuverlässig abgebildet werden. Ziel der vorliegenden Doktorarbeit ist es deshalb, die Entwicklung von Gesteinspermeabilitäten und -steifigkeiten aufgrund von Porenraumveränderungen genauer vorherzusagen. Für die möglichst exakte Bestimmung physikalischer Gesteinsparameter ist eine detaillierte Darstellung der Mikrostruktur notwendig. Basierend auf mikro-computertomographischen Scans werden daher hochaufgelöste, dreidimensionale Modelle typischer Reservoirsandsteine erstellt und Gesteinspermeabilität und -steifigkeit numerisch berechnet. Um charakteristische Verteilungen von Sekundärmineralen abzubilden, wird der Porenraum dieser virtuellen Sandsteinproben systematisch verändert und die resultierenden Auswirkungen auf die granulometrischen, hydraulischen und elastischen Gesteinseigenschaften bestimmt. Die Ergebnisse zeigen deutlich, dass charakteristische Fällungsmuster unterschiedlicher geochemischer Reaktionsregime die Permeabilität erheblich beeinflussen. Folglich ist die Nutzung von porenskaligen Modellen zur Bestimmung der Porosität-Permeabilitätsbeziehungen unbedingt notwendig. Im Gegensatz dazu ist die Verteilung von Sekundärmineralen für die Gesteinssteifigkeit nur bei bestimmten Mikrostrukturen von Bedeutung, hierzu zählen hochporöse Sandsteine oder solche mit Korngerüst-stützenden Zementierungen. In der Arbeit wird außerdem ein Ansatz zur Konstruktion granularer Gesteine vorgestellt, welcher sowohl die Kornsedimentation als auch die diagenetische Verfestigung umfasst. Es wird gezeigt, dass die synthetischen Proben die mikrostrukturelle Komplexität natürlicher Reservoirsandsteine gut abbilden. Dadurch können generelle Limitationen von bildgebenden Verfahren überwunden und unterschiedlichste virtuelle Repräsentationen von granularen Gesteinen generiert werden. Die synthetischen Proben können zukünftig in virtuellen Experimenten verwendet werden, um die Auswirkungen von Lösungs- und Fällungsreaktionen auf verschiedene petrophysikalische Korrelationen zu untersuchen. Die vorgestellte Arbeit liefert methodische Grundlagen zur Quantifizierung von Permeabilitäts- und Steifigkeitsänderungen infolge geochemischer Prozesse. Die berechneten petrophysikalischen Beziehungen basieren direkt auf mikrostrukturellen Veränderungen des Porenraums. Daher bieten sie eine genauere Vorhersage der Gesteinseigenschaften als herkömmliche analytische Methoden, wodurch sich die Aussagekraft von Reservoirmodellen erheblich verbessert. Somit können Risiken, wie Produktivitäts- oder Injektivitätsverluste sowie Reservoirkompaktion oder Störungsreaktivierung, verringert werden. Die präsentierten Ergebnisse sind daher relevant für verschiedenste Bereiche der geologischen Untergrundnutzung wie CO2- oder Energiespeicherung, Geothermie, Kohlenwasserstoffgewinnung sowie die Erkundung hydrothermaler Lagerstätten.

Beschreibung: Monitoring and controlling of mineral precipitation is a major challenge in several industry segments such as oil and gas or geothermal energy. Especially at elevated pressures, temperatures, and salinities, these reactions are poorly understood and difficult to predict. Real-time monitoring under such conditions, however, is crucial for reduction or prevention of scaling as well as for process modelling. In this study, a novel fibre-optical method for monitoring the barite precipitation process at high pressures and high temperatures was developed and tested. The applicability of a fibre-optical sensor was demonstrated and kinetic data derived from experiments using synthetic fluids that were oversaturated with barite. The ionic strength was varied between 0 and 5 mol/l NaCl to simulate geothermal brines typically found in deep sedimentary reservoirs. Experiments were performed at 25 and 150 °C. For the studied variations it was found that the ionic strength had a strong impact on the barite precipitation rate whereas temperature had only a minor influence. The kinetic rate constants were decreasing with increasing salinity ranging from 9 to 0.4 l‧mol−1‧s−1 at 25 °C and from 8.6 to 0.66 l‧mol−1‧s−1 at 150 °C.

Beschreibung: Earth’s lowermost mantle displays complex geological phenomena that likely result from its heterogeneous physical interaction with the core. Geophysical models of core-mantle interaction rely on the thermal and electrical conductivities of appropriate geomaterials which, however, have never been probed at representative pressure and temperature (P-T) conditions. Here we report on the opacity of single crystalline bridgmanite and ferropericlase and link it to their radiative and electrical conductivities. Our results show that light absorption in the visible spectral range is enhanced upon heating in both minerals but the rate of change in opacity with temperature is a factor of six higher in ferropericlase. As a result, bridgmanite in the lowermost mantle is moderately transparent while ferropericlase is highly opaque. Our measurements support previous indirect estimates of low (〈 1 W/m/K) and largely temperature-independent radiative conductivity in the lowermost mantle. This implies that the radiative mechanism has not contributed significantly to cooling the Earth’s core throughout the geologic time. Opaque ferropericlase is electrically conducting and mediates strong core-mantle electromagnetic coupling, explaining the intradecadal oscillations in the length of day, low secular geomagnetic variations in Central Pacific, and the preferred paths of geomagnetic pole reversals.

Beschreibung: The SWATH‐D experiment involved the deployment of a dense temporary broadband seismic network in the Eastern Alps. Its primary purpose was enhanced seismic imaging of the crust and crust–mantle transition, as well as improved constraints on local event locations and focal mechanisms in a complex part of the Alpine orogen. The study region is a key area of the Alps, where European crust in the north is juxtaposed and partially interwoven with Adriatic crust in the south, and a significant jump in the Moho depth was observed by the 2002 TRANSALP north–south profile. Here, a flip in subduction polarity has been suggested to occur. This dense network encompasses 163 stations and complements the larger‐scale sparser AlpArray seismic network. The nominal station spacing in SWATH‐D is 15 km in a high alpine, yet densely populated and industrialized region. We present here the challenges resulting from operating a large broadband network under these conditions and summarize how we addressed them, including the way we planned, deployed, maintained, and operated the stations in the field. Finally, we present some recommendations based on our experiences.

Beschreibung: The Pan-Eurasian Experiment (PEEX) Science Plan, released in 2015, addressed a need for a holistic system understanding and outlined the most urgent research needs for the rapidly changing Arctic-boreal region. Air quality in China, together with the long-range transport of atmospheric pollutants, was also indicated as one of the most crucial topics of the research agenda. These two geographical regions, the northern Eurasian Arctic-boreal region and China, especially the megacities in China, were identified as a “PEEX region”. It is also important to recognize that the PEEX geographical region is an area where science-based policy actions would have significant impacts on the global climate. This paper summarizes results obtained during the last 5 years in the northern Eurasian region, together with recent observations of the air quality in the urban environments in China, in the context of the PEEX programme. The main regions of interest are the Russian Arctic, northern Eurasian boreal forests (Siberia) and peatlands, and the megacities in China. We frame our analysis against research themes introduced in the PEEX Science Plan in 2015. We summarize recent progress towards an enhanced holistic understanding of the land–atmosphere–ocean systems feedbacks. We conclude that although the scientific knowledge in these regions has increased, the new results are in many cases insufficient, and there are still gaps in our understanding of large-scale climate–Earth surface interactions and feedbacks. This arises from limitations in research infrastructures, especially the lack of coordinated, continuous and comprehensive in situ observations of the study region as well as integrative data analyses, hindering a comprehensive system analysis. The fast-changing environment and ecosystem changes driven by climate change, socio-economic activities like the China Silk Road Initiative, and the global trends like urbanization further complicate such analyses. We recognize new topics with an increasing importance in the near future, especially “the enhancing biological sequestration capacity of greenhouse gases into forests and soils to mitigate climate change” and the “socio-economic development to tackle air quality issues”.

Beschreibung: Mantle metasomatism is an important process in subduction zones in which fluids from the dehydrating oceanic slab interact with the overlying upper mantle resulting in a chemical alteration of the mantle. Consequently, this fluid-rock interaction may influence the mantle rock's physical properties such as the deformation behavior. In order to study element redistribution during mantle metasomatism in the laboratory, we used the simplified model reaction olivine + quartz = orthopyroxene, where olivine acts as representative for the upper mantle and quartz as proxy for the metasomatizing agent. We conducted piston-cylinder experiments at 1.5 GPa and 950 to 1400 °C, lasting between 48 and 288 h, on samples containing a mixture of quartz and one set of synthesized forsterite samples doped with either Co, Ni, Mn, or Zn. Additionally, we tested the influence of either nominally anhydrous or hydrous experimental conditions on the chemical distribution of the respective dopant element by using either crushable alumina or natural CaF2 as pressure medium. Results of the chemical analyses of the recovered samples show dopant specific partitioning between doped forsterite and orthopyroxene independent of the confining pressure medium; except for the runs in which Ni-doped forsterite samples were used. The observed Ni- and Co-enrichment in forsterite samples may be used to identify mantle rocks that underwent mantle metasomatism in nature.

Beschreibung: The project DARE proposes an integrated study of seismic site effects on the deep and elongated Messinian Rhône Canyon (French Rhône Valley). Lithological information from boreholes reaching the bedrock and preliminary geophysical campaigns indicate that the canyon can reach locally 〉500 m and may be deeply incised. The strong material contrast between the sedimentary filling and the substratum, as well as its expected confined geometry make this canyon a good candidate for generating various kinds of multi-dimensional site effects. Waveform data are available from the GEOFON data centre, under network code Y7.

Beschreibung: The temporary seismic array of MySCOLAR in northern Myanmar consists of 30 broadband stations. The overall scientific goals are to understand the transition from continental collision to oceanic subduction, to quantify the partitioning of deformation in the accretionary prism, in the Burma Plate and along the strike-slip Sagaing fault system and to image the subducting Indian Plate beneath Myanmar and southwest China. The seismological analysis methods applied to this dataset will include location of local earthquakes and determining their focal mechanisms, surface wave tomography from ambient noise and earthquake data, body wave tomography from local and teleseismic earthquakes, full waveform inversion for Earth structure, receiver functions, and shear wave splitting. A subset of the stations was transmitting data in real time, and these stations contributed to real-time earthquake analysis by the Department of Meteorology and Hydrology (DMH) in Myanmar and the GEOFON earthquake monitoring service. Waveform data are available from the GEOFON data centre, under network code 6C.

Beschreibung: Near Reykjavik/Iceland, a "soft stimulation” geothermal experiment was performed in the frame of the DESTRESS project in 2019. The installed seismic stations consist of short period, and borehole stations in and around Geldinganes, NE of Reykjavik. The task of this network is the monitoring of the seismic events in the area around the stimulation site. The installation started in late 2018 with 6 short period stations (Reykjavik Energy). Since July 2019 additional seismic stations were integrated as a small scale array on the island Geldinganes and additional short period stations. A borehole geophone chain was installed with 17 short period 3-component geophones with a vertical spacing of 10 meter in the depth interval 1040m to 1200 m. Waveform data are available from the GEOFON data centre, under network code YG, and are embargoed until January 2025.

Beschreibung: The Villarrica Volcano is one of the most active volcanoes in South America and is located in a major tourism region. A dense seismological network is used to investigate the seismic characteristics of the volcano and its seismic structure tomographically with high spatial resolution. The network was in operation for 2 week from 01.03.2012 to 14.03.2012. It consisted of 30 3-component and 45 1-component short period seismographs covering an area of 2000 km*2. The covered area has a diameter of 45 km and includes the volcanic building.

Beschreibung: The Bransfield Strait is a seismically active extensional rift located between the Antarctic Peninsula and the South Shetland Islands. The Strait is partly located on continental crust including areas within the transition to seafloor spreading. The amphibious seismic network BRAVOSEIS is an international effort focused on the seismological research of submarine volcanoes and rift dynamics in the Bransfield Strait. This network is the onshore component of the entire network consisting of 15 broadband land stations deployed in the South Shetland Islands and Antarctic Peninsula between January 2018 and February 2020. The offshore components (network code ZX) include 9 broadband ocean bottom seismometers (OBS) across the Central Bransfield Basin and a group of 6 hydrophone moorings spanning the rift area of 200 x 100 km2, with inter-station distance of ~30 km. Additionally, a smaller offshore array consisting of 15 short-period OBSs with an aperture of 20 km and a narrow inter-station distance of ~4 km was deployed around the Orca submarine volcanic edifice south of King George Island. The data will be used to study the geodynamics of the Bransfield Strait and the evolution of the incipient rifting zone in the domain where extension has been suggested. Seismological methods will include earthquake location, source mechanism, surface wave analysis with ambient noise and earthquake data, receiver function and shear wave splitting. The results may shed light on the crustal structure and tectonic regime in the region and image the location and extent of magma accumulations related to submarine volcanic structures. Finally, the results should provide clues to assess the internal processes that occur in the submarine volcanoes of the area undergoing rifting. Waveform data are available from the GEOFON data centre, under network code 5M and are embargoed until Mar 2024. Acknowledgments: We thank all participants in the BRAVOSEIS 2018, 2019, and 2020 cruises, with a special acknowledgement to Capt. Jose Emilio Regodon and his crew at R/V Hesperides; Capt. Juan Carlos Hernandez and his crew at Sarmiento de Gamboa; Miki Ojeda, Ezequiel Gonzalez, and all the UTM staff involved in the planification and realization of the surveys. We also thank the Spanish Polar Committee and institutions involved in the management of the Spanish Antarctic campaigns and the development of the Spanish Polar Program. We are grateful for the help and support that we always find in the personnel of the Antarctic Bases, especially the Spanish Bases Juan Carlos I and Gabriel de Castilla.

Beschreibung: The main aim of this project is to investigate the crustal and mantle structure beneath the Longmenshan fault zone in China, based on a very dense passive seismology profile. The Longmenshan fault zone hosted the Wenchuan earthquake of May 2008 with a magnitude (Mw) of 7.9 and the Lushan earthquake of June 2013 with a magnitude (Mw) of 6.6. It is planned to mainly use the receiver-function method, to investigate the crustal and mantle structure beneath the Longmenshan fault zone. Waveform data are available from the GEOFON data center, under network code 4O under license CC BY-NC-ND 4.0, and are embargoed until February 2024.

Beschreibung: Oldoinyo Lengai in the North Tanzanian Divergence is the only active natrocarbonatite volcano world-wide and presents an important endmember magmatic system in a young rift segment of the East African Rift System. This volcano typically experiences long-duration episodes of natrocarbonatitic effusions with intermittent short-duration explosive eruptions. To better understand the role of the stress interactions and magma plumbing on the eruptions dynamics, this study aims to constrain the subsurface magmatic architecture of Oldoinyo Lengai, and a zone of ongoing seismicity and intrusive activity below the extinct 1 Ma-old Gelai shield volcano and active Naibor Soito monogenetic cone field. Waveform data are available from the GEOFON data centre, under network code 9J, and are embargoed until November 2025.

Beschreibung: On 26th of November 2019 an Mw 6.4 earthquake ruptured near the port town of Durres, only 25 km from Tirana, the capital of Albania. The earthquake caused major damage and killed 51 people, making it the deadliest earthquake in 2019 worldwide. The mainshock was relatively deep (~25 km) and of thrust type. In December 2019, a Hazard and Risk team (HART) from German Center for Geosciences (GFZ), Karslruhe Institute of Technology (KIT) in cooperation with the Institute of Geosciences, Energy, Water and Environment (IGEWE) of the Polytechnic University of Tirana, Albania installed a 30-station seismic network in the epicentral region to record aftershocks. Stations were equipped with short-period (1 Hz or 4.5 Hz) 3-component seismometers and CUBE data loggers recording continuously 100 sps. Waveform data are available from the GEOFON data centre, under network code 9K under CC-BY 4.0 license and are embargoed until January 2024.

Beschreibung: Irpinia seismic Array is part of the DEnse mulTi-paramEtriC observations and 4D high resoluTion imaging (DETECT) project focused on the acquisition of a unique multiparametric dataset and fosters collaboration among various institutions. The University of Naples Federico II (UniNa) and the German Research Centre for Geosciences (GFZ) are leading this effort carried out in collaboration with various local institutions and supported by the local municipalities. The DETECT project aims at exploiting very dense seismic networks deployed across a segmented fault system (Irpinia and Pergola-Melandro) to foster the development of scientific integrated methodologies for monitoring and imaging the fault behavior during the inter-seismic phase. The Irpinia seismic Array consists of a dense constellation of seismic antennas using more than 200 seismic stations deployed for one year. Each seismic antenna, with maximum aperture of ~2 km, uses one broad-band sensor, one short period sensor with 1 Hz and 8 with 4.5 Hz natural frequency. The antennas are deployed above and near the fault segments that generated during the last centuries many strong earthquakes in the southern Apennines. Waveform data are available from the GEOFON data centre, under network code ZK.

Beschreibung: Raw, SEGY and other supplementary data are presented from the seismic refraction / wide-angle reflection profile, TTZ-South, in Poland and Ukraine. The purpose of this 550 km long seismic profile was to reveal the lithospheric structure along the Teisseyre-Tornquist Zone (TTZ), a major geophysical boundary in Europe.

Beschreibung: The Villarrica Volcano is one of the most active volcanoes in South America and is located in a major tourism region. A dense temporal seismological network was installed to investigate the volcanic seismicity and the seismic structure of the edifice with seismic traveltime tomography at high spatial resolution. The network was in operation for 2 weeks from 01.03.2012 to 14.03.2012. It consisted of 30 three-component and 45 one-component short period seismographs covering an area of about 2000 km2. The covered area has a diameter of 45 km and includes the volcanic building.

Beschreibung: Study Region Vietnamese Mekong Delta. Study focus This study investigates the trends of groundwater levels (GWLs), the memory effect of alluvial aquifers, and the response times between surface water and groundwater across the Vietnamese Mekong Delta (VMD). Trend analysis, auto- and cross-correlation, and time-series decomposition were applied within a moving window approach to examine non-stationary behavior. New hydrological insights Our study revealed an effective connection between the shallowest aquifer unit (Holocene) and surface water, and a high potential for shallow groundwater recharge. However, low-permeable aquicludes separating the aquifers behave as low-pass filters that reduce the high‐frequency signals in the GWL variations, and limit the recharge to the deep groundwater. Declining GWLs (0.01−0.55 m/year) were detected for all aquifers throughout the 22 years of observation, indicating that the groundwater abstraction exceeds groundwater recharge. Stronger declining trends were detected for deeper groundwater. The dynamic trend analysis indicates that the decrease of GWLs accelerated continuously. The groundwater memory effect varied according to the geographical location, being shorter in shallow aquifers and flood-prone areas and longer in deep aquifers and coastal areas. Variation of the response time between the river and alluvial aquifers was controlled by groundwater depth and season. The response time was shorter during the flood season, indicating that the bulk of groundwater recharge occurred in the late flood season, particularly in the deep aquifers.

Beschreibung: Important aspects of geothermal exploration and exploitation are the assessment and mitigation of natural and/or induced seismicity, the imaging and resource assessment of a geothermal reservoir and the monitoring of the effects of the exploitation activities. With this thesis the analysis and application of various seismological tools for the planning, exploration, and monitoring of geothermal fields are given. The methods explored and further developed include survey design for microseismic network construction and assessment, local earthquake tomography, ambient noise tomography, and coda wave interferometry. These techniques were applied to study Los Humeros (Mexico), Theistareykir (Iceland) and Reykjanes (Iceland) geothermal fields. The geometry of seismological arrays is essential for high quality seismic event retrieval and minimal location errors. A sequential survey design algorithm that uses a quality measure based on the D-criterion was applied to extend the seismic network at Theistareykir and to qualify the geometry of the Reykjanes network. Assuming mean picking errors of tp =0.2 s and ts = 0.4 s, the extended Theistareykir network presented an improvement of ∼0.2 km for the computed hypocentral components of seismic events located within the new network. Conversely, we estimated that the Reykjanes network could spare up to 18 of its station locations and obtain comparable location errors nonetheless. This study showed thei mportance of prior survey design experiments to optimize the expenses for a geothermal project (required number of sensors) while obtaining good location estimates of expected seismic events (benefit/cost relations). To characterize the seismic structures at Los Humeros and Theistareykir geothermal fields, a local earthquake tomography and an ambient noise tomography were computed at both locations, respectively. A local earthquake tomography is feasible in areas with high seismicity and good ray coverage (earthquake/station geometries). On the other hand, an ambient noise tomography depends on a good and sufficiently dense station distribution. With the results of these studies, the seismic structures and the dynamics of these two producing fields were obtained for the first time. The seismicity distribution at Los Humeros was used to characterize structures and potential permeability enhancements in some of the existing faults. The retrieved Vp model was combined with available well log data and ultrasonic pulse measurements of collected rock samples to estimate the boundaries of different geologic units. The Vp/Vs model was then used in combination with resistivity data and surface CO2 measurements to deduce the geometry of the conductive clay cap (Vp/Vs ≤ 1.65 and resistivities ≤ 10 Ωm), to identify fluid (Vp reduction, Vp/Vs ≥ 1.71, and resistivities between ∼ 10-60 Ωm), and to locate gas bearing regions (Vp/Vs ≤ 1.55 and high surface CO2 concentrations). A similar study was carried out at Theistareykir, where the Vs model was combined with resistivity data to identify magmatic and/or hydrothermal bodies (Vs ≤ -7 %, resistivities ≤ 30 Ωm). An important conclusion from these studies is that the combination of seismic properties with additional geological and/or geophysical data avoids ambiguities and provides robust interpretations of the dynamics and structure of a geothermal reservoir. Finally, a coda wave interferometry technique (stretching method) was applied to two years of ambient noise records at the Theistareykir geothermal field with the aim to monitor possible velocity changes due to the exploitation activities. Here, the effects of the injection and production changes were very small on the computed Δv/v ratio and only a small long-term velocity reduction (possibly due to production) was detected (-0.05 %/year at the producing field compared to a regional -0.04 %/year). Such observations are also very relevant for the safe long-term continuation of exploitation activities. Although not yet a standard practice, the computation of these changes is very useful to control aseismic processes prior to potentially triggered/induced large seismic events and is complementary to microseismic monitoring. With these results, this thesis contributes to the efforts of the International Energy Agency to develop and increase the use of geothermal energy.

Beschreibung: Wichtige Aspekte der geothermischen Exploration und Nutzung sind die Bewertung und Reduzierung der natürlichen und/oder induzierten Seismizität, die Abbildung und Ressourcenbewertung eines geothermischen Reservoirs, sowie die Überwachung der Auswirkungen der Explorationsaktivitäten. In dieser Arbeit werden die Analysen und Anwendungen verschiedener seismologischer Methoden zur Planung, Erkundung und Überwachung geothermischer Felder vorgestellt. Die untersuchten und hier weiterentwickelten Methoden umfassen das Design für den Aufbau und die Bewertung mikroseismischer Netzwerke, die lokale Erdbebentomographie, die Ambient-Noise-Tomographie und die Coda-Wellen-Interferometrie. Diese Techniken werden angewendet, um die geothermischen Felder Los Humeros (Mexiko), Theistareykir (Island) und Reykjanes (Island) zu untersuchen. Die Geometrie seismologischer Arrays ist essentiell für eine gute Bestimmung seismischer Ereignisse mit kleinen Lokationsungenauigkeiten. Ein sequenzieller Algorithmus zum Design des Arrays, der eine Qualitätskennzahl auf der Grundlage des D-Kriteriums verwendet, wurde benutzt, um das seismische Netzwerk in Theistareykir zu erweitern und die Geometrie des Reykjanes-Netzwerks zu testen. Unter der Annahme von mittleren Ablesefehlern von tp = 0.2 s und ts = 0.4 s für P- und S- Wellen verbessert das erweiterte Theistareykir-Netzwerk die berechneten Hypo-Zentren um 0.2 km innerhalb des neuen Netzes. Das Reykjanes-Netz könnte andererseits um bis zu 18 Stationsstandorte reduziert werden und dennoch vergleichbare Lokationsgenauigkeiten erzielen. Diese Studie zeigte die Wichtigkeit vor den eigentlichen Feld-Experimenten Tests möglicher Array-Designs durchzuführen um die Kosten für ein geothermisches Projekt (erforderliche Anzahl von Sensoren) zu optimieren und gleichzeitig gute Lokationen für erwartete seismische Ereignisse zu erhalten (Nutzen/Kostenverhältniss). Um die seismischen Strukturen der geothermischen Felder Los Humeros und Theistareykir zu charakterisieren, wurden an beiden Standorten eine lokale Erdbebentomographie und eine Ambient-Noise-Tomographie berechnet. Eine lokale Erdbebentomographie ist in Gebieten mit hoher Seismizität und guter Strahlenabdeckung (Erdbeben/Stationsgeometrie) möglich. Eine Ambient-Noise-Tomographie hängt nur von einer guten und ausreichend dichten Stationsverteilung ab. Mit den Ergebnissen dieser Studien wurden dann erstmals die seismischen Strukturen und die Dynamik dieser beiden produzierenden Felder ermittelt. Die Seismizitätsverteilung in Los Humeros wurde verwendet, um Strukturen und potenzielle Verbesserungen in der Durchlässigkeit einiger Störungszonen zu charakterisieren. Das abgeleitete Vp-Modell wurde mit Bohrloch-Daten und Ultraschall-Messungen an Gesteinsproben kombiniert, um die Grenzen verschiedener geologischer Einheiten abzuschätzen. Das Vp/Vs-Modell wurde dann in Kombination mit Widerstandsdaten und Oberflächen-CO2-Messungen verwendet, um die Geometrie der leitfähigen Tonkappe abzuleiten (Vp/Vs ≤ 1.65 und Widerstand ≤ 10 Ωm), um Fluide zu identifizieren (reduzierte Vp Werte, Vp/Vs ≥ 1.71 und Widerstände zwischen 10-60 Ωm) und um gasführende Bereiche zu lokalisieren (Vp/Vs ≤ 1.55 und hohe CO2-Konzentrationen). Eine ähnliche Studie wurde in Theistareykir durchgeführt, wo das Vs-Modell mit Widerstandsdaten kombiniert wurde, um magmatische und/oder hydrothermale Körper zu identifizieren (Vs ≤ -7 %, Widerstände ≤ 30 Ωm). Eine wichtige Schlussfolgerung aus diesen Studien ist, dass die Kombination von seismischen Eigenschaften mit zusätzlichen geologischen und/oder geophysikalischen Daten Mehrdeutigkeiten vermeidet und robuste Interpretationen der Dynamik und Struktur eines geothermischen Reservoirs liefert. Weiterhin wurde eine Coda-Wellen-Interferometrie-Technik (Dehnungsmethode) auf zwei Jahre Ambient-Noise-Daten im Geothermiefeld Theistareykir angewendet, um mögliche Geschwindigkeitsänderungen aufgrund der Nutzung des Feldes zu überwachen. Hier waren die Auswirkungen der Injektions- und Produktionsveränderungen auf das Δv/v-Verhältnis sehr gering und nur eine kleine, möglicherweise produktionsbedingte, langfristige Geschwindigkeitsreduktion wurde festgestellt (-0.05 %/Jahr innerhalb des Produktionsbereichs im Vergleich zum regionalen Wert von -0.04 %/Jahr). Solche Beobachtungen sind für die sichere langfristige Ausbeutung von geothermischen Feldern von großer Bedeutung. Obwohl es noch keine Standardpraxis ist, ist die Berechnung dieser Änderungen weiterhin sehr nützlich, um aseismische Prozesse vor potenziell ausgelösten/induzierten großen seismischen Ereignissen zu steuern, sie ergänzen weiterhin die mikroseismische Überwachung. Mit diesen Ergebnissen trägt die vorgelegte Arbeit zu den Bemühungen der Internationalen Energieagentur bei, die Nutzung von Geothermie zu entwickeln und zu erhöhen.

Beschreibung: The Caribbean and South American tectonic plates bound the north-eastwards expulsion of the North Andean Block in western Venezuela. This complicated geodynamic setting resulted in the formation of major strike-slip fault systems and sizeable mountain chains. The 100 km wide Mérida Andes extend from the Colombian/Venezuelan border to the Coastal Cordillera. To the north and south, the Mérida Andes are bound by hydrocarbon-rich sedimentary basins. Knowledge of lithospheric structures, related to the formation of the Mérida Andes, is limited though, due to a lack of deep geophysical data. This thesis presents the results of the first broadband magnetotelluric profile crossing the Mérida Andes and the Maracaibo and Barinas-Apure foreland basins spanning a distance of 240 km. The MT dataset consists of 72 stations installed during March and April 2015 with a minimum recording period of 3 days per station. Geoelectrical strike and dimensionality analyses are consistent with one- or two- dimensional subsurface structures for the sedimentary basins yet also indicate a strong three- dimensional setting for the Mérida Andes. Even more significantly, these analyses showed the presence of off-profile features that influenced the data considerably, particularly at long periods. Therefore, a combination of 2D and 3D modelling was necessary for analysing the geoelectrical structures associated with this dataset. Off-profile structures can significantly affect the outcome of a 2D inversion. Thus, the systematic examination of the influence of 3D structures on 2D inversions was necessary to support the obtained result. Synthetic data sets derived from 3D modelling allowed identification and quantification of spurious off-profile features as well as smoothing artefacts due to limited areal station coverage of data collected along a profile. In general, structures in the 2D inversion are affected by the projection and rotation of the data resulting in sub-horizontal anomalies to reproduce the oblique extent of the fault systems and sedimentary basins. Moreover, a profile distributed dataset can limit the lateral resolution of a 3D inversion considerably. Hence, the effect of data distribution on a 3D inversion was carefully studied to determine the areas of the models that can be confidently explained by the data. To this end, several synthetic datasets were derived from 3D models with varying levels of complexity. The analysis of the synthetic datasets allows determination of the lateral resolution of the 3D models and identification of spurious shallow and deep features considered artefacts related to off-profile features. Furthermore, the inversion of synthetic models provided support to the geological interpretation of the recovered anomalies for the 2D and 3D modelling. The 2D and 3D inversion models were similar above the sedimentary basins and showed marked differences above the Mérida Andes, due to the 3D nature of this section. The inversion models show electrically conductive basins with depths of 2 to 5 km for the Barinas-Apure and 2 to 9 km for the Maracaibo basins. Many resistive bodies within the Maracaibo basin could be related to active deformation causing juxtaposition of older geological formations and younger basin sediments. A conductive zone under the Maracaibo Basin correlates spatially with the location of a Bouguer anomaly low and seem to describe the SE tilt of the Maracaibo Triangular Block under the Mérida Andes. This conductive zone is limited towards the mountain by the north-western thrust system, whose fault plane may function as a detachment surface reaching depth larger than 30 km in the 3D inversion models. The most prominent fault systems of the area, the Boconó and Valera Faults, cross-cut the Mérida Andes in NE-SW direction along its strike with a length 400 km and N-S direction at its centre with a length 60 km, respectively. Both faults are associated with sub-vertical zones of high electrical conductivity and sensitivity tests suggest that the Valera fault reach depths of up to 12 km. The Boconó fault can be considered a crustal structure with a depth up to 35 km. The observed anomalies seem to show a deep connection of the fault planes, possibly related to the formation of the fault systems in a transpressive regime. Conductive anomalies to the south of the Boconó Fault seem to represent a considerable back thrust structure well constrained between 3 and 10 km depth. The high conductivity of these structures is possibly related to weathering water from the surface and the accumulation of clay minerals in the fault gauges. However, fluids related to the flat and shallow subduction of the Caribbean Plate in north-western Venezuela could better explain the low resistivity of the deep structures (〉 15 km). A sizeable conductor at 50 km depth, which appears consistently in the 2D sections, could be identified as an inversion artefact caused by a conductor east of the profile. The 3D inversion places this structure 10 km to the east at 15 km deep. This model also shows depth connection (12 km depth) of the anomalies related to the Valera and Boconó faults with the off-profile conductor. The observed anomalies in the 2D and 3D inversion related to these conductors were tested and reproduced employing synthetic datasets, leading to the speculation that the high conductivity associated with the off-profile conductor may be related to the detachment of the Trujillo Block. The models obtained confirm the shape and distribution of the known geological structures related to the complicated geodynamic settings responsible for the formation of the Mérida Andes. These results partially support the "floating orogen hypothesis" developed to explain the geodynamic evolution of western Venezuela, and they highlight the relevance of the Trujillo Block in this process. However, they also show that features of known structures such as the Boconó fault system maximum depth, the back-thrusting in the Mérida Andes, and the relevance of the escape of the Trujillo Block in the tectonic processes need to be adjusted to the current knowledge.

Beschreibung: Die Karibische und Südamerikanische tektonischen Platten begrenzten die nordöstliche Abschiebung des Nord-Anden-Blocks im Westen Venezuelas. Diese komplizierte geodynamische Umgebung führte zur Bildung großer Blattverschiebungen-Verwerfungssysteme und beträchtlicher Gebirgsketten. Die 100 km breiten Mérida-Andes erstrecken sich von der Grenze zwischen Kolumbien und Venezuela bis zum Coastal Cordillera. Im Norden und Süden sind die Mérida-Anden von kohlenwasserstoffreichen Sedimentbecken umgeben. Das Wissen über lithosphärische Strukturen im Zusammenhang mit der Bildung der Mérida-Anden ist jedoch aufgrund des Mangels an tiefen geophysikalischen Daten begrenzt. Diese Arbeit präsentiert die Ergebnisse des ersten breitbandigen magnetotellurischen Profils, das die Mérida-Anden und die Vorlandbecken Maracaibo und Barinas-Apure über eine Entfernung von 240 km quert. Der MT-Datensatz besteht aus 72 Stationen, die im März und April 2015 mit einer Mindestaufzeichnungsdauer von 3 Tagen pro Station installiert wurden. Geoelektrische Streich- und Dimensionalitätsanalysen stimmen mit ein- oder zweidimensionalen Untergrundstrukturen für die Sedimentbecken überein, weisen jedoch auch auf eine starke dreidimensionale Strukturen in der Umgebung der Mérida-Andes hin.. Noch wichtiger ist, dass diese Analysen das Vorhandensein von Merkmalen außerhalb des Profils zeigten, die die Daten insbesondere für lange Perioden erheblich beeinflussen. Daher war eine Kombination aus 2D- und 3D-Modellierung erforderlich, um Störungssysteme und Sedimentbecken quer zum Profil zu reproduzieren. Off-Profile-Strukturen können das Ergebnis einer 2D-Inversion erheblich beeinflussen. Daher war die systematische Untersuchung des Einflusses von 3D-Strukturen auf 2D-Inversionen erforderlich, um das erhaltene Ergebnis zu verifizieren. Synthetische Datensätze, die aus der 3D-Modellierung abgeleitet wurden, ermöglichten die Identifizierung und Quantifizierung von störenden Strukturen außerhalb des Profils sowie die Glättung von Artefakten aufgrund der begrenzten Stationsüberdeckung der entlang eines Profils gesammelten Daten. Im Allgemeinen werden Strukturen in der 2D-Inversion durch die Projektion und Rotation der Daten beeinflusst, was zu flach stehenden Anomalien führt, um die schräge Ausdehnung der Verwerfungssysteme und Sedimentbecken zu reproduzieren. Darüber hinaus kann ein Datensatz entlang eines Profils die laterale Auflösung einer 3D-Inversion erheblich einschränken. Daher wurde die Auswirkung der Datenverteilung auf eine 3D-Inversion sorgfältig untersucht, um die Bereiche der Modelle zu bestimmen, die durch die Daten sicher erklärt werden können. Zu diesem Zweck wurden mehrere synthetische Datensätze aus 3D-Modellen mit unterschiedlicher Komplexität abgeleitet. Die Analyse der synthetischen Datensätze ermöglicht die Bestimmung der lateralen Auflösung der 3D-Modelle und die Identifizierung von störenden oberflächennahen und tiefen Merkmalen, die als Artefakte im Zusammenhang mit Strukturen außerhalb des Profils betrachtet werden. Darüber hinaus unterstützte die Inversion synthetischer Modelle die geologische Interpretation der reproduzierten Anomalien für die 2D- und 3D-Modellierung. Die 2D- und 3D-Inversionsmodelle stimmen über den Sedimentbecken überein. Aufgrund der 3D Strukturen über den MA ergaben sich jedoch deutliche unterschiede. Die Inversionsmodelle zeigen elektrisch leitende Becken mit Tiefen von 2 bis 5 km für das Barinas-Apure und 2 bis 9 km für das Maracaibo-Becken. Viele Gebiete höheren Widerstands im Maracaibo-Becken könnten mit einer aktiven Deformation zusammenhängen, die ein Nebeneinander älterer geologischer Formationen und jüngerer Beckensedimente verursacht. Eine besserleitende Zone unter dem Maracaibo-Becken korreliert räumlich mit der Lage einer Bouguer-Anomalie und scheint die SE-Neigung des Maracaibo-Dreiecksblocks unter den Mérida-Anden zu markieren. Diese leitende Zone ist in Richtung des Gebirges durch das nordwestliche Schubsystem begrenzt, dessen Störungsebene als Ablösefläche fungieren kann, die in den 3D-Inversionsmodellen eine Tiefe von mehr als 30 km erreicht. Die bekanntesten Störungssysteme des Gebiets, die Verwerfungen Boconó und Valera, kreuzen die Mérida-Anden in Nordost-Südwest-Richtung entlang ihres Streichens mit einer Länge von 400 km und die N-S-Richtung in ihrer Mitte mit einer Länge von 60 km. Beide Störungen sind durch steil stehende Zonen hoher elektrischer Leitfähigkeit verbunden. Sensitivitätsstudien legen nahe, dass die Valera-Störung Tiefen von bis zu 12 km erreicht. Die Boconó-Verwerfung kann als Krustenstruktur mit einer Tiefe von bis zu 35 km angesehen werden. Die beobachteten Anomalien scheinen eine tiefe Verbindung der Verwerfungsebenen zu zeigen, möglicherweise im Zusammenhang mit der Bildung der Verwerfungssysteme in einem transpressiven Regime. Leitfähige Anomalien südlich der Boconó-Verwerfung scheinen eine beträchtliche Rückschubstruktur darzustellen, die zwischen 3 und 10 km Tiefe gut lokalisiert ist. Die hohe Leitfähigkeit dieser Strukturen hängt möglicherweise mitWasser aus Verwitterungsprozessen nahe der Erdoberfläche und der Ansammlung von Tonmineralien in den Störungszonen zusammen. Alte Fluide im Zusammenhang mit flach stehenden und oberflächenahen Subduktion der Karibikplatte im Nordwesten Venezuelas könnten jedoch den geringen spezifischen Widerstand der tiefen Strukturen (〉 15 km) besser erklären. Ein beträchtlicher Leiter in einer Tiefe von 50 km, der in den 2D-Schnitten konsistent erscheint, konnte als Inversionsartefakt identifiziert werden, der durch einen Leiter östlich des Profils verursacht wird. Durch die 3D-Inversion wird diese Struktur 10 km östlich in 15 km Tiefe platziert. Dieses Modell zeigt auch die Tiefenverbindung (12 km Tiefe) der Anomalien im Zusammenhang mit den Störungen von Valera und Boconó mit dem Leiter außerhalb des Profils. Die beobachteten Anomalien in der 2D- und 3D-Inversion in Bezug auf diese Leiter wurden unter Verwendung synthetischer Datensätze getestet und reproduziert. Daher kann man annehmen das die mit dem Leiter abseits des Profils verbundene Leitfähigkeit mit der Ablösung des Trujilo Blocks zusammenhängt. Die erhaltenen Modelle bestätigen die Form und Verteilung der bekannten geologischen Strukturen im Zusammenhang mit dem komplizierten geodynamischen Millieu, welches für die Bildung der Mérida-Andes verantwortlich ist. Diese Ergebnisse stützen teilweise die "schwimmende Orogenhypothese", die entwickelt wurde, um die geodynamische Entwicklung West-Venezuelas zu erklären, und sie unterstreichen die Relevanz des Trujillo-Blocks in diesem Prozess. Sie zeigen jedoch auch, dass Merkmale bekannter Strukturen wie die maximale Tiefe des Boconó-Verwerfungssystems, das Zurückschieben in den Mérida-Anden und die Relevanz des Entweichens des Trujillo-Blocks in den tektonischen Prozessen an den aktuellen Kenntnisstand angepasst werden müssen.

Beschreibung: The Pacific Northwest (PNW) has substantial earthquake risk, both due to the offshore Cascadia megathrust fault but also other fault systems that produce earthquakes under the region's population centers. Forecasts of aftershocks following large earthquakes are thus highly desirable and require statistical models of a catalog of the PNW’s past earthquakes and aftershock sequences. This is complicated by the fact that the PNW contains multiple tectonic regimes hypothesized to have different aftershock dynamics as well as two types of earthquake clustering (aftershock sequences and swarms). The Epidemic-Type Aftershock Sequence (ETAS) model is a top-performing spatiotemporal point process model which describes the dynamics of earthquakes and aftershocks in a seismic region using a set of parameters. Typically, maximum likelihood estimation is used to fit ETAS to an earthquake catalog; however, the ETAS likelihood suffers from flatness near its optima, parameter correlation and numerical instability, making likelihood-based estimates less reliable. We present a Bayesian procedure for ETAS estimation, such that parameter estimates and uncertainty can be robustly quantified, even for small and complex catalogs like the PNW. The procedure is conditional on knowing which earthquakes triggered which aftershocks; this latent structure and the ETAS parameters are estimated iteratively. The procedure uses a Gibbs sampler to conditionally estimate the posterior distributions of each part of the model. We simulate several synthetic catalogs and test the modelling procedure, showing well-mixed posterior distributions centered on true parameter values. We also use the procedure to model the continental PNW, using a new catalog formed by algorthmically combining US and Canadian data sources and then, identifying and removing earthquake swarms. While MLEs are unstable and depend on both the optimization procedure and its initial values, Bayesian estimates are insensitive to these choices. Bayesian estimates also fit the catalog better than do MLEs. We use the Bayesian method to quantify the uncertainty in ETAS estimates when including swarms in the model or modelling across different tectonic regimes, as well as from catalog measurement error. Seismicity rate estimates and the earthquake forecasts they yield vary spatially and are usually represented as heat maps. While the visualization literature suggests that displaying forecast uncertainty improves understanding in users of forecast maps, research on uncertainty visualization (UV) is missing from earthquake science. In a pre-registered online experiment, we test the effectiveness of three UV techniques for displaying uncertainty in aftershock forecasts. Participants completed two map-reading tasks and a comparative judgment task, which demonstrated how successful a visualization was in reaching two key communication goals: indicating where many aftershocks and no aftershocks are likely (sure bets) and where the forecast is low but the uncertainty is high enough to imply potential risk (surprises). All visualizations performed equally well in the goal of communicating sure bet situations. But the visualization mapping the lower and upper bounds of an uncertainty interval was substantially better than the other map designs at communicating potential surprises. We discuss the implications of these experimental results for the communication of uncertainty in aftershock forecast maps.

Beschreibung: Carbonate has often been identified in aqueous carbonic inclusions in spodumene-bearing and other pegmatites, but its origin remains unclear. Here, the conditions at which carbonate and hydrogen carbonate can be generated from spodumene, CO2 and H2O, were studied using a hydrothermal diamond-anvil cell (HDAC) and Raman spectroscopy. In all experiments, spodumene persisted in aqueous carbonic solution up to the maximum temperature (600 to 800 °C). Heating of hydrogen carbonate/oxalate solutions produced CO2- and HCO−3-rich peralkaline fluids, which resulted in strong corrosion of spodumene (and polylithionite-trilithionite) and, in one run, formation of zabuyelite [Li2(CO3)] crystals at low temperatures. The experiments indicate that the reaction of spodumene with CO2 and H2O requires a peralkaline fluid to proceed rapidly. In addition, they show that spodumene crystallizes upon the heating of quartz, muscovite, and aqueous lithium carbonate solution. We conclude that if the aqueous fluid was rich in alkali hydrogen carbonate, zabuyelite in fluid inclusions in pegmatites can form both via a subsolidus reaction of CO2-bearing fluid inclusion with the spodumene host or by trapping a peralkaline fluid early in the evolution of simple or complex pegmatites. The results of our experimental study strengthen the conclusion that, although counterintuitive, hydrogen carbonate-rich peralkaline fluids may be involved in the evolution of peraluminous granitic pegmatites, in which peralkaline minerals are normally absent or very rare.

Beschreibung: Knowledge of groundwater flow is of high relevance for groundwater management and the planning of different subsurface utilizations such as deep geothermal facilities. While numerical models can help to understand the hydrodynamics of the targeted reservoir, their predictive capabilities are limited by the assumptions made in their set up. Among others, the choice of appropriate hydraulic boundary conditions, adopted to represent the regional to local flow dynamics in the simulation run, is of crucial importance for the final modelling result. In this publication we present the hydrogeological models to obtain results to quantify how and to which degree different upper hydraulic boundary conditions and vertical cross boundary fluid movement influence the calculated deep fluid conditions Therefore, we take the central Upper Rhine Graben area as a natural laboratory. The presented three models are set up with different sets of boundary conditions. The Reference Model uses the topography as upper hydraulic pressure surface of 0 kPa. In model S1, a subdued replica of the topography, which was built on the base of hydraulic head measurements is applied as the upper hydraulic boundary condition and in model S2 vertical cross boundary flow is implemented. Based on our results, we illustrate in the landing paper that for the Upper Rhine Graben specific characteristics of the flow field are robust and insensitive to the choice of imposed hydraulic boundary conditions, while specific local characteristics are more sensitive. Accordingly, these robust features characterizing the first order groundwater flow dynamics in the Upper Rhine Graben include: (i) a regional groundwater flow component descending from the graben shoulders to rise at its centre; (ii) infiltration of fluids across the graben shoulders, which locally rise along the main border faults; (iii) the presence of heterogeneous hydraulic potentials at the rift shoulders. The configuration of the adopted boundary conditions influence primarily calculated flow velocities and the absolute position of the upflow axis within the graben sediments. In addition, the choice of upper hydraulic boundary conditions exerts a direct control on the evolving local flow dynamics, with the degree of influence gradually decreasing with increasing depth. With respect to regional flow modelling of basin hosted, deep water resources, the main conclusions derived from this study are: (i) the often considered water table as an exact replica of the model topography (Reference Model) likely introduces a source of error in the simulations in regional hydraulic modelling approaches. Here, we show that these errors can be minimized by making use of a water table as upper boundary condition derived from available hydraulic head measurements (model S1). If the study area is part of a supra-regional flow system - like the central Upper Rhine Graben is part of the whole Upper Rhine Graben - the in- and outflow across vertical boundaries need to be considered (model S2).

Beschreibung: This dataset provides friction and elasticity data from ring shear and axial tests, respectively, on rock analogue materials used at the University Roma Tre (Rome, IT) in “Foamquake”, a novel seismotectonic analog model mimicking the megathrust seismic cycle (Mastella et al., under review). Two granular materials (quartz sand and Jasmine rice) have been characterized by means of internal friction coefficients µ and cohesions C. An elastic material (foam rubber) have been characterized by means of Young’s modulus E and Poisson’s ratio v. According to our analysis the granular materials show Mohr-Coulomb behaviour characterized by linear failure envelopes in the shear stress vs. normal load Mohr space. Peak, dynamic and reactivation friction coefficients of the quartz sand are µP = 0.69, µD = 0.56 and µR = 0.64, respectively. Cohesion ranges between 50 and 100 Pa. Rate-dependency of friction in quartz sand seems insignificant. Peak, dynamic and reactivation friction coefficients of the Jasmine rice are µP = 0.70, µD = 0.59 and µR = 0.61, respectively. Cohesion ranges between 30 and 50 Pa. Rate-weakening of Jasmine rice is c. 6% per tenfold change in shear velocity v. The Young’s modulus of the foam rubber has been constrained to 30 kPa, its Poisson’s ratio is v=0.1.

Beschreibung: This dataset includes particle image correlation data from 26 experiments performed with Foamquake, a novel analog seismotectonic model reproducing the megathrust seismic cycle. The seismotectonic model has been monitored by the means of a high-resolution top-view monitoring camera. The dataset presented here represents the particle image velocimetry surface velocity field extracted during the experimental model through the cross-correlation between consecutive images. This dataset is supplementary to Mastella et al. (2021) where detailed descriptions of models and experimental results can be found.

Beschreibung: Fluid inclusions in pegmatite minerals were studied using Raman spectroscopy to determine the carbon species. Carbon dioxide is very abundant in the aqueous liquid and vapor phases. Occasionally, CH4 was found in the vapor. In the aqueous liquid, HCO−3 was detected in fluid inclusions in tantalite-(Mn) from the Morrua Mine and in late-stage quartz from the Muiâne pegmatite and the Naipa Mine, all in the Alto Ligonha District, Mozambique. Moreover, we observed a carbonate (calcite group) in fluid inclusions in garnet from the Naipa Mine and in beryl from the Morrua Mine, both in the Alto Ligonha District, Mozambique, and a calcite-group carbonate and whewellite [CaC2O4∙H2O] in fluid inclusions in topaz from Khoroshiv, Ukraine. The occurrence of oxalate is interpreted to be due to a reaction of some form of carbon (possibly CO or bitumen) with a peralkaline fluid. Our results support the hypothesis that, although counterintuitive, hydrogen carbonate-rich peralkaline fluids may locally be involved in the evolution of peraluminous granitic pegmatites, in which peralkaline minerals are normally absent or very rare.

Beschreibung: One of the most fundamental questions in ecology is how many species inhabit the Earth. However, due to massive logistical and financial challenges and taxonomic difficulties connected to the species concept definition, the global numbers of species, including those of important and well-studied life forms such as trees, still remain largely unknown. Here, based on global ground-sourced data, we estimate the total tree species richness at global, continental, and biome levels. Our results indicate that there are ∼73,000 tree species globally, among which ∼9,000 tree species are yet to be discovered. Roughly 40% of undiscovered tree species are in South America. Moreover, almost one-third of all tree species to be discovered may be rare, with very low populations and limited spatial distribution (likely in remote tropical lowlands and mountains). These findings highlight the vulnerability of global forest biodiversity to anthropogenic changes in land use and climate, which disproportionately threaten rare species and thus, global tree richness.

Beschreibung: We report a new improvement in self-organized maps for geological interpretation of geophysical data. By using a multi-geophysical dataset recorded in the mining area of Thuringia, Germany, we show the results of replacing the typical feature analysis by a principal component analysis. By performing a transformation of the dataset according to a few of the principal components, we obtain a more detailed representation of the local geology than previous works. Results also show a significant improvement in processing time, while also minimizing influence of user´s interpretation.

Beschreibung: The Songliao Block is located in the eastern part of the Central Asian Orogenic Belt. Its creation and evolution are believed to be related to the closure of the Paleo-Asian and Mongol-Okhotsk oceans, and to the subduction of the Paleo-Pacific Ocean. The deep seismic reflection profiles showed that there are sloping mantle reflections below the Songliao Block, which are suspected to be the result of the convergence of three tectonic domains. However, it is still not clear the current structural form of the Songliao Block is caused by the direct action or not of the tectonic systems. This work used 138 broadband magnetotelluric stations to obtain a three-dimensional electrical structural model of the northern Songliao Block. The results showed there are orthogonal network fault systems, faulted basins, igneous rocks. And the Lindian fault depression is the center of the asthenospheric upwelling, the shallowest up to 45 km. Combined with evidence from seismic studies, we proposed that the superposition of tectonic systems may have produced weak tectonic zones. These zones provided channels for the later upward movement of fluids and melt, likely due to hydrous upwellings caused by the subduction of the Paleo-Pacific system.

Beschreibung: Understanding how resilient rangelands are to climatic disturbances such as drought is of major importance to land managers. The resilience of ecosystems can be reduced by livestock grazing and by environmental conditions. Most studies quantifying resilience are based on model simulations. However, natural time series from satellite data offer the possibility to infer aspects of resilience from real systems. The objective of this study was to investigate two aspects of ecological resilience, namely resistance to climate variability and recovery from drought, by applying a change detection method (Breaks For Additive Seasonal and Trend; BFAST) spatially on a 28-year Landsat NDVI time series in a dry rangeland in southern Cyprus. First, we used the number of breakpoints fitted by the BFAST model as an inverted proxy for long-term vegetation resistance to climate variability (the ability to withstand change during a disturbance reduces the likelihood to trigger a breakpoint in the time series). Second, we used the linear slope of the BFAST model after a known drought as a proxy of the recovery rate of the vegetation. This information was then used to analyse the spatial distribution of the total number of breakpoints and of the NDVI recovery trend in relation to grazing and environmental properties. Our results show that high NDVI and a northern orientation (i.e. favourable environmental conditions) were associated with a highly resilient system, due to high resistance to climate variability and fast recovery after drought. Intermediate conditions were associated with low resistance. Unfavourable conditions and high grazing intensities were associated with an unresponsive ecosystem state characterised by high resistance and slow recovery after a drought event. Low grazing intensities positively affected the NDVI recovery trend, but did not improve resistance. On northern slopes, terrain slope had a positive effect on the NDVI recovery trend, while on southern slopes it had a negative effect. Our satellite-driven approach has a strong potential for resilience monitoring, because it can be applied on broad spatial and temporal scales in areas with low availability of field data. Moreover, it allows to jointly extract two important components of resilience: resistance and recovery rate.

Beschreibung: The evolution of chemical bonding in ferropericlase, (Mg,Fe)O, with pressure may affect the physical and chemical properties of the Earth’s lower mantle. Here, we report highpressure optical absorption spectra of single-crystalline ferropericlase ((Mg0.87Fe0.13)O) up to 135 GPa. Combined with a re-evaluation of published partial fluorescence yield X-ray absorption spectroscopy data, we show that the covalency of the Fe−O bond increases with pressure, but the iron spin transition at 57−76.5 GPa reverses this trend. The qualitative crossover in chemical bonding suggests that the spin-pairing transition weakens the Fe−O bond in ferropericlase. We find, that the spin transition in ferropericlase is caused by both the increase of the ligand field-splitting energy and the decrease in the spin-pairing energy of high-spin Fe2+.

Beschreibung: he competition between the impact of inherited weaknesses and plate kinematics determines the location and style of deformation during rifting, yet the relative impacts of these ‘internal’ and ‘external’ factors remain poorly understood, especially in 3D. In this study, we used brittle-viscous analogue models to assess how multiphase rifting, that is changes in plate divergence rate or direction, and the presence and orientation of weaknesses in the competent mantle and crust, influences rift evolution. We find that the combined reactivation of mantle and crustal weaknesses without any kinematic changes already creates complex rift structures. Divergence rates affect the strength of the weak lower crustal layer and hence the degree of mantle-crustal coupling; slow rifting decreases coupling, so that crustal weaknesses can dominate deformation localisation and surface structures, whereas fast rifting increases coupling and deformation related to mantle weaknesses can have a dominant surface expression. Through a change from slow to fast rifting mantle-related deformation can overprint structures that previously formed along (differently oriented) crustal weaknesses. Conversely, a change from fast to slow rifting may shift deformation from mantle-controlled towards crust-controlled. When changing divergence directions, structures from the first rifting phase may control where subsequent deformation occurs, but only when they are sufficiently well developed. We furthermore place our results in a larger framework of brittle-viscous rift modelling results from previous experimental studies, showing the importance of general lithospheric layering, divergence rate, the type of deformation in the mantle, and finally upper crustal structural inheritance. The interaction between these parameters can produce a variety of deformation styles that may, however, lead to comparable end products. Therefore, careful investigation of the distribution of strain localisation, and to an equal extent of basin depocenter locations over time is required to properly determine the evolution of complex rift systems, providing an incentive to revisit various natural examples.

Beschreibung: The Kiaka orogenic gold deposit (Burkina Faso), located in the Paleoproterozoic domain of the West African Craton, is characterized by a two-stage gold mineralization hosted in volcano-sedimentary metamorphic rocks that was formed during the Eoeburnean (2.20–2.13 Ga) and Eburnean (2.13–2.05 Ga) orogenic cycles. These two stages include an early disseminated low-grade gold mineralization and a late vein-hosted high-grade gold mineralization. Paragenetic studies indicate that the first gold stage was coeval with the deposition of hydrothermal tourmaline. The aim of this paper is twofold: (i) to determine the processes responsible for deposition of the economic disseminated gold mineralization and (ii) to identify the source of the mineralizing fluids. For this purpose, we performed an in situ study on tourmaline by combining electron probe microanalysis and secondary ion mass spectrometry measurements of boron isotopes. Hydrothermal tourmaline hosted in metamafic volcanic rocks and metagreywackes has a dravite composition but shows different δ11B values falling within the two intervals of − 25.1 to − 22.0‰ and − 19.6 to − 15.1‰, respectively. Our results suggest that tourmaline formed from a distal, high-temperature (ca. 400 °C), reduced, and low-salinity hydrothermal fluid that interacted with the local host rocks. Based on the modeling of tourmaline–fluid boron isotope fractionation, we propose a metamorphic fluid origin derived from devolatilization of deeply buried muscovite schists during the regional prograde to peak metamorphism prior ca. 2.13 Ga. This metamorphic fluid–rock interaction model may possibly extend to other orogenic gold deposits in the West African Craton.

Beschreibung: Earthquakes are generally known to alter the stress field near seismogenic faults. Observations using YRY-four-gauge borehole strainmeters within Yushu (YSH) borehole near the Ganzi-Yushu fault in eastern Tibetan Plateau shows that the azimuth variation of maximum horizontal stress (SH) first decreased and then increased substantially when the earthquakes occurred during the measurement period from January 1, 2009 to December 31, 2018. In this period, 38 earthquakes (M ≥ 3) were detected near the fault and the SH orientation showed a drastic change after the 2010 Ms 7.3 Yushu mainshock. We present a discrete element modelling using Particle Flow Code 2D (PFC2D) to simulate a dynamic fault rupturing process and to use the modelling results for interpretation of the stress reorientation. The modelling reveals that dilatation and compression quadrants are formed around a fault rupturing in strike-slip model, resulting in different spatiotemporal changes of the orientation of maximum horizontal stress (Δθ). The value of Δθ in the compression quadrants shows a sharp drop at the time of coseismic slip, then approaches slowly to an asymptotic value. In the dilatation quadrants, Δθ drops by coseismic slip, then increases sharply and finally reaches a stable value. The modelled Δθ by coseismic fault slip agrees with in-situ observations at YSH borehole during 2010 Ms 7.3 Yushu mainshock. It is also found that, the value of Δθ decreases with increasing distance from the rupturing source. We modelled the effect of fault geometry and host rock properties on the Δθ, and found that structural complexity and off-fault damage by coseismic fault slip have significant impact on the stress field alteration near the rupturing source.

Beschreibung: The permeability characteristics of natural fracture systems are crucial to the production potential of shale gas wells. To investigate the permeability behavior of a regional fault that is located within the Wufeng Formation, China, the gas permeability of shale samples with natural micro-fractures was measured at different confining pressures and complemented with helium pycnometry for porosity, computed micro-tomographic (µCT) imaging, and a comparison with well testing data. The cores originated from a shale gas well (HD-1) drilled at the Huayingshan anticline in the eastern Sichuan Basin. The measured Klinkenberg permeabilities are in the range between 0.059 and 5.9 mD, which roughly agrees with the permeability of the regional fault (0.96 mD) as estimated from well HD-1 productivity data. An extrapolation of the measured permeability to reservoir pressures in combination with the µCT images shows that the stress sensitivity of the permeability is closely correlated to the micro-fracture distribution and orientation. Here, the permeability of the samples in which the micro-fractures are predominantly oriented along the flow direction is the least stress sensitive. This implies that tectonic zones with a large fluid potential gradient can define favorable areas for shale gas exploitation, potentially even without requirements for hydraulic fracture treatments.

Beschreibung: Clathrate hydrates—also known as gas hydrates—are ice-like compounds consisting of gas and water molecules. They occur wherever elevated pressures and low temperatures prevail; and where enough water and hydrate-forming gas molecules are available. Therefore, natural gas hydrates occur at all active and passive continental margins, in permafrost regions, in some deep lakes, and under unfavorable circumstances, also, in pipelines. This article provides an overview of the (thermodynamic) requirements and various models for the nucleation and growth of gas hydrates and the different gas hydrate structures that may occur and which have been detected in nature. Furthermore, this study also shows the influence of the properties of the enclosed gas molecules such as size and shape on the structure and thermodynamic properties of the resulting hydrate phase. Finally, the complexity of a natural environment with regard to the various influences of sediments, microbial activity, and salinity of the pore fluid on hydrate formation is also discussed.

Beschreibung: We believe that the transdisciplinary studies on water-related multi-hazards are innovative and critical research by the water community, thus answering the call of the recent Nature Sustainability Editorial ‘Too much and not enough’1 for water science ideas that are not derivative or stagnant. This domain of water studies focuses on the specific contexts where water-related hazardous events occur simultaneously, in cascade or cumulatively with other events. Characteristic of the field is the intensive collaboration of scientists and practitioners from different disciplines working together to better understand, assess and manage water-related multi-hazards. At the recent Asia Oceania Geosciences Society–European Geosciences Union Joint Conference on New Dimensions for Natural Hazards in Asia, we discussed the statement ‘Too much and not enough’1 and here suggest three reasons why transdisciplinary collaborations have led to many new ideas and notable advancements in the field of water-related multi-hazard research in recent years.

Beschreibung: Stibnite was mined until the end of the twentieth century in the Schlaining ore district, Austria, near the easternmost border of the Eastern Alps where windows of Penninic ophiolites and metasediments are exposed below Austroalpine tectonic units. In Early Miocene, structurally controlled small vein and metasomatic stibnite-quartz deposits were formed in Penninic Mesozoic calcareous marbles and calcite schists. Fluid inclusion studies identified two fluids involved in the mineralization: (i) a low-salinity, low-CO2 metamorphic fluid that precipitated quartz at approximately 240 °C and (ii) a stibnite-forming ore fluid that had a meteoric origin. There is no evidence of boiling or that the fluids mixed during mineralization. The ore components Sb and H2S were leached by fluid/rock interaction from buried rock units. Stibnite mineralization occurred by cooling the ore fluid to below 300 °C, at less than 2000 m depth. Quartz precipitated at slightly lower temperatures, approximately contemporaneous with stibnite. Fluid migration and ore deposition are probably related to high heat flow during the exhumation of the Rechnitz Window in response to Neogene extension and/or shallow Early Miocene andesitic magmatism. The study emphasizes that data obtained from the analyses of gangue minerals alone cannot routinely be used to infer the origin and depositional conditions of the associated ore minerals.

Beschreibung: The upper crustal structure of the individual tectono-magmatic segments within the Main Ethiopian Rift is poorly understood. The Gedemsa tectono-magmatic segment (GTMS) is considered to be within an intermediate stage of rift development with magma-assisted rifting contributing to overall extension. The interpretation of gravity data from Global Gravity Model plus2013 has provided an insight on the nature and structure of the shallow upper crust beneath the GTMS within the Central Main Ethiopian Rift. A detailed investigation of the structure of the shallow upper crust based on 2D gravity models, 3D geologic interpretation maps, and residual and upward continued gravity anomaly maps indicate the occurrence of dense, mafic intrusions at depths between 10-20 km and 5–10 km beneath the Tullu Moye-Gedemsa and Boku volcanic centers, respectively. The gravity data indicate at deeper depths, the GTMS (25 km wide and 60 km long) is elongated parallel to the Main Ethiopian Rift axis. Our 2D gravity modeling indicates that the Quaternary faults within the Wonji Fault Belt are closely associated with the dense magmatic intrusions as the faults aid in transporting melt from the deeper crustal melt areas to the shallow magma chambers.

Beschreibung: A Cretaceous paleo-accretionary wedge, the Ashin Complex, now exposed along the Zagros suture zone in southern Iran, exhibits mafic, metasedimentary and ultramafic lithologies. Field, geochemical and petrological observations point to an anomalous high-temperature event that gave rise to the formation of peritectic (trondhjemitic) melts associated with restitic garnet-bearing amphibolites. Lu-Hf isotopic dating of centimetre-sized garnet in amphibolite-facies metasediments yielded a crystallization age of 113.10 ± 0.36 Ma, possibly representing the age of prograde to near-peak metamorphic conditions. SHRIMP U-Th-Pb zircon dating from trondhjemitic leucosomes yields crystallization ages of 104 ± 1 Ma, interpreted as the age of the temperature peak, which occurred in the upper amphibolite-facies (c. 650–680 °C at 1.1–1.3 GPa), according to thermodynamic modelling and Ti-in-zircon thermometry. Rutile crystals from two leucosomes yield Zr-in-rutile temperatures in the range of 580–640 °C and a LA-ICP-MS U-Pb age range from 85 to 112 Ma, interpreted as a consequence of partial re-equilibration during incipient cooling. A late static recrystallization event is indicated by the presence of sodic-calcic clinopyroxene, sodic amphibole, Si-rich phengite, titanite overgrowths after rutile and lawsonite within former leucosomes and late fractures. This mineral assemblage is a typical blueschist-facies (high pressure-low temperature) paragenesis and is interpreted as reflecting long-term isobaric cooling that occurred until the end of the Cretaceous as a consequence of increasing slab thermal age. This first report of a melting event in the Zagros paleo-accretionary wedge reveals the presence of a transient, abnormally high thermal gradient of c. 18 °C/km that occurred at c. 105–113 Ma. We speculate that this could be explained by the subduction of a thermal anomaly such as a seamount chain, a transform fault system or, more likely, a spreading ridge under the southern Iranian margin. Indeed, paleogeographic reconstructions of the Tethyan realm suggest the entrance of the Northern Tethyan basin ridge into the subduction zone shortly after 120 Ma.

Beschreibung: We present a comprehensive Li and B isotope study of granites, aplites, and igneous enclaves from the multi-phase Eibenstock granite in the Western Erzgebirge-Vogtland metallogenic province of Germany. The studied samples cover the entire compositional range of the granites from moderately to highly evolved and include variably altered types as obtained by magmatic fractionation, post-magmatic high- to medium-temperature and near-surface low-temperature alteration. Fractionation and alteration processes are unequivocally documented by the chemical variability of the rocks. Despite the marked imprint of these processes on bulk-rock compositions, our granite samples show only little variation in δ7Li (−0.52 to 0.75‰) and δ11B (−17.46 to −14.78‰), with surface samples defining the lower end of the δ7Li range. The narrow range in δ7Li suggests that magmatic fractionation and high-temperature overprint have a very minor effect on δ7Li. The B budget of the samples is dominated by tourmaline, which makes δ11B values insensitive for later high- to medium-temperature overprint or surficial low-temperature alteration. Depending on whether tourmaline crystallized before or after exsolution and loss of magmatic fluids, whole-rock samples have higher or lower δ11B values. Granite enclaves have δ7Li and δ11B values ranging from −1.51 to −0.81‰ and − 14.55 to −13.89‰, respectively. Some samples have chemical and mineralogical evidence for wall-rock interaction during emplacement or later overprint by external fluids. These samples show broader ranges in δ7Li (−2.61 to 2.21‰) and δ11B (−21.58 to −9.85‰). These values show that wall-rock interaction via assimilation and external fluids may affect δ7Li and δ11B to a larger extent than intra-magmatic processes, such as fractional crystallization, fluid-mediated autometasomatic overprinting, or exsolution of fluids from the melt. The offset of δ7Li and δ11B values towards the compositions of the wall rocks reflects the contrasting composition of granite and country rock and the addition of country-rock material to the granite. The magnitude of the offset reflects both the relative contribution of wall-rock derived Li and B to the granite and the magnitude of the difference in the Li and B isotopic compositions between them.

Beschreibung: At the Decatur carbon capture and storage site (IL, USA) CO2 has been injected from 2011–2014 and from 2017 to present near the base of the Lower Mt. Simon Sandstone saline reservoir, resulting in microseismicity. Microseismicity is mainly located in the basement and distributed in distinct spatial clusters. The lack of significant impedance contrasts within the basement makes the interpretation of active-source seismic reflection data challenging, however, recent reprocessing allowed to resolve faults above and at the top of the basement. These faults generally do not coincide with the location of microseismic events and their continuation to the general depth of the seismic events cannot be assumed. This paper shows how the interpretation of the microseismicity can complement structural interpretations of active-source seismic reflection data. In particular, we analyze clusters and bursts (abrupt increases) of microseismicity, identify unresolved, smaller-scale weaknesses and extract statistical parameters. These parameters allow comparisons with the interpreted faults, and with fracture sets intercepted by boreholes. During injection at the Decatur site, the injection pressure was kept far below fracture pressure, nevertheless, seismic events were induced and spread far beyond the expected extent of the CO2 plume. We argue that local stress transfers related to the CO2 injection reactivated pre-existing fractures within the critically stressed basement. Finally, we conducted a slip tendency analysis for faults interpreted from active seismic, selected cluster, bursts and nodal planes from focal mechanisms to determine if the interpreted structures are optimally oriented with respect to the stress regime. Our results suggest that the orientation of fractures close to the injection well, generally shows slight deviations from the optimal orientation for slip. This might indicate either slight local deviations of the maximum horizontal stress azimuth from the average direction used in the analysis, or the lack of optimally oriented fractures at this location.

Beschreibung: The occurrence of felt earthquakes due to gas production in Groningen has initiated numerous studies and model attempts to understand and quantify induced seismicity in this region. The whole bandwidth of available models spans the range from fully deterministic models to purely empirical and stochastic models. In this article, we summarise the most important model approaches, describing their main achievements and limitations. In addition, we discuss remaining open questions and potential future directions of development.

Beschreibung: Py4HIP is an open-source software tool for Heat-In-Place calculations implemented as a self-explanatory Jupyter notebook written in Python (Py4HIP.ipynb) Calculating the Heat In Place (HIP) is a standard method for assessing the geothermal potential for a defined geological unit (e.g., Nathenson, 1975; Muffler and Cataldi, 1978; Garg and Combs, 2015). The respective implementation in Py4HIP is based on a volumetric quantification of contained energy after Muffler and Cataldi (1978), where the geological unit at hand is considered spatially variable in terms of its temperature, thickness, porosity, density and volumetric heat capacity of its solid and fluid (brine) components. The energy values

Beschreibung: This work presents a study of an unusual synthetic diamond sample, where calcite is the main phase in addition to diamond. Using Raman spectroscopy it can be shown that in addition to diamond carbon forms four different generations: highly ordered graphite-II, lower ordered graphite-I as spandrel between different mineral grains, weakly-ordered carbonaceous material as graphite-III and the hexagonal diamond lonsdaleite (here graphite-IV) as last formation in calcite. Graphite-III can be found mainly in the calcite body. According to the fine dispersion of this carbonaceous material and the arrangement on grain boundaries, we assume that carbon was dissolved in the calcite melt, and that by the activated state of carbon in the calcite melt the formation of diamond is favored near 1760°C, and 6.8 GPa. Lonsdaleite as minor phase may have been formed as an longlived intermediate state under standard conditions between graphite and diamond. Evidence shows that there are two different lonsdaleite phases (possibly hexagonal and monocline) present. The prevailing diamond is characterized by the first-order Raman line at 1333 cm-1. However, there are also present diamonds with the first-order Raman line down to 1310.6 cm-1, corresponding to 13C = 0.511. Significantly there is the strong decrease of the optical damage threshold with increase of the 13C content.

Beschreibung: Erosion by landslides is a common phenomenon in mountain regions around the globe, affecting all climatic zones. Landslides facilitate bedrock weathering, pedogenesis and ecological succession, being key drivers of biodiversity. Landslide chronosequences have long been used for studies of vegetation succession in initial ecosystems, but they further offer ideal model systems for studies of soil development and microbial community succession. In this review we synthesize the state of knowledge on the role of landslides in ecosystems, their influence on element cycles and interactions with biota. Further, we discuss feedback mechanisms between global warming, landslide activity and greenhouse gas emissions. In the view of increasing anthropogenic influence and climate change, soils are becoming a critical resource. Due to their ubiquity, landslide chronosequences have the potential to provide critical insights into soil development under different climates and thereby contribute to future soil restoration efforts.

Beschreibung: Llallagua is one of the world's biggest tin deposits and part of the metallogenic Bolivian Tin Belt, which occurs in the inner arc of the Central Andes. The Llallagua deposit formed through emplacement of a subvolcanic porphyry stock, of intermediate dacitic to rhyodacitic composition, metasomatism, and hydrothermal mineralization. The deposit is the subject of a well-established geochronologic controversy that includes contradictory ages (≈40 and 20 Ma) from multiple geochronometers. Geochronological characterization of phosphate minerals from the metasomatized igneous porphyry and hydrothermal vein assemblages, along with detailed petrography and chemical analyses, are used to reconcile the age controversy. The new interpretation is further supported by existing textural and geochemical data. The Usingle bondPb ages of unaltered igneous fluorapatite (e.g. 21.2 ± 2.9 Ma), monazite (e.g. 21.22 + 0.80/−0.66 Ma), and zircon (e.g. 21.15 ± 0.39 Ma) from the Llallagua porphyry are all equivalent within error at ≈20 Ma, and are interpreted to represent the age of porphyry stock emplacement. Ages determined from altered portions of these minerals (e.g. apatite, 18.8 ± 8.0 Ma) are within error the same as unaltered portions, suggesting alteration of the porphyry soon after emplacement. Usingle bondPb ages from unaltered fluorapatite (e.g. 21.4 ± 6.7 Ma), altered and unaltered monazite (e.g. 19.4 ± 1.4 Ma and 20.29 ± 0.30 Ma respectively), and unaltered xenotime (e.g. 19.32 ± 0.67 Ma) from the hydrothermal veins are also ≈20 Ma, indicating a very short hiatus or temporal continuity between stock emplacement and hydrothermal vein formation. The early Miocene age for tin mineralization at Llallagua is consistent with neighboring mineralization centers at the Morococala, Colquechaca, Japo, and Santa Fe mines, and with regional trends observed across the Bolivian Tin Belt, including stratigraphy, magmatic, and tectonic history of the Andean Eastern Cordillera.

Beschreibung: Scandium is important in modern technology and is regarded as a strategic metal in many countries. It is highly dispersed in Earth’s crust and rarely forms independent minerals. Clinopyroxene is the most important Sc-bearing mineral in some world-class deposits hosted in mafic–ultramafic intrusions, which are also the major source of laterite-hosted Sc deposits. However, the factors controlling Sc distribution in minerals have been little explored, impeding the understanding of the geochemical behavior of Sc and why it is common in some clinopyroxene grains. The newly discovered Mouding Sc deposit in SW China is hosted in a zoned intrusion composed, from core to rim, of monzogabbro, syenogabbro, gabbro, magnetite clinopyroxenite, and clinopyroxenite. Clinopyroxene in the intrusion is diopsidic in composition with high Sc contents (80–105 ppm). In-situ trace element mappings of diopside crystals reveal homogeneous, zoned, swallow-tailed, and hourglass internal Sc distribution patterns. These patterns can be produced through kinetically controlled incorporation of Sc on different crystal faces. The preferential substitution of Sc can take place on the {1 0 0}, {1 1 0} and {0 1 0} prism faces because of the high flexibility of the octahedral M1 protosites. The fast growth of diopside, which facilitates kinetically controlled crystallization, is dominated by textural coarsening and promoted by the hydrous parental magmas with low viscosities and active convection. The active flow and efficient interstitial communication of the magma can direct compatible elements from the magma into clinopyroxene, thus favoring formation of Sc-rich grains. Our study provides a feasible way to study intra-grain variations of Sc in minerals and emphasizes that kinetic effects may play a critical role in Sc distribution and enrichment in hydrous magmatic Sc deposits. We also show that disequilibrium crystallization may be more pervasive than previously thought, and the hourglass zoning of clinopyroxene can provide valuable information on this process.

Beschreibung: Iron-bearing carbonates play an important role in Earth’s carbon cycle. Owing to their stability at mantle conditions, recently discovered iron carbonates with tetrahedrally coordinated carbon atoms are candidates for carbon storage in the deep Earth. The carbonates’ iron oxidation and spin state at extreme pressure and temperature conditions contribute to the redox conditions and element partitioning in the deep mantle. By laser heating FeCO3 at pressures of about 83 GPa, Fe3+4C3O12 and Fe2+2Fe3+2C4O13 were synthesized and then investigated by x-ray emission spectroscopy to elucidate their spin state, both in situ and temperature quenched. Our experimental results show both phases in a high-spin state at all pressures and over the entire temperature range investigated, i.e., up to 3000 K. The spin state is conserved after temperature quenching. A formation path is favored where Fe3+4C3O12 forms first and then reacts to Fe 2+2Fe3+2C4O13, most likely accompanied by the formation of oxides. Density functional theory calculations of Fe2+2Fe3+2C4O13 at 80 GPa confirm the experimental findings with both ferric and ferrous iron in high-spin state with antiferromagnetic order at 80 GPa. As the intercrystalline cation partitioning between the Fe-bearing carbonates and the surrounding perovskite and ferropericlase depends on the spin state of the iron, an understanding of the redox conditions prevalent in subducted slab regions in the lower mantle has to take the latter into account. Especially, Fe2+2Fe3+2C4O13 may play a key role in subducted material in the lower mantle, potentially with a similar role as silicate perovskite.

Beschreibung: Transnational access (TNA) allows cross‑border, short‑term and frequently free‑of‑charge access to world-class research facilities, to foster collaborations and exchanges of experience. Specifically, TNA aims to encourage open science and innovation and to increase the efficient and effective use of scientific infrastructure. Within EPOS, the European Plate Observing System, the Volcano Observatories and Multi‑scale Laboratories communities have offered TNA to their high‑quality research facilities through national and European funding. This experience has allowed the definition, design, and testing of procedures and activities needed to provide transnational access in the EPOS context. In this paper, the EPOS community describes the main objectives for the provision of transnational access in the EPOS framework, based on previous experiences. It includes practical procedures for managing transnational access from a legal, governance, and financial perspective, and proposes logistical and technical solutions to effectively execute transnational access activities. In addition, it provides an outlook on the inclusion of new thematic communities within the TNA framework, and addresses the challenges of providing market‑driven access to industry.

Beschreibung: Semi-arid Mongolia is a highly sensitive region to climate changes, but the region’s Holocene paleoclimatic evolution and its underlying forcing mechanisms have been the subject of much recent debate. Here we present a continuous 7.4 ka sediment record from the high-altitude Shireet Naiman Nuur (Nuur = lake) in the central Mongolian Khangai Mountains. We extensively dated the sediments and analyzed elemental composition and bulk isotopes for lake sediment characterization. Our results show that 14C-dating of bulk organic carbon and terrestrial macrofossils provide a robust and precise chronology for the past 7.4 ± 0.3 cal ka BP at Shireet Naiman Nuur and 14C-ages are mostly in stratigraphic order. The 14C-based chronology is confirmed by paleomagnetic secular variations, which resemble the predictions of spherical harmonic geomagnetic field models. The very good chronological control makes paleomagnetic secular variation stratigraphy a powerful tool for evaluating and refining regional 14C-chronologies when compared to the record presented here. The lake sediment proxies TOC, N, log (Ca/Ti) and log (Si/Ti) reveal increased lake primary productivity and high growing season temperatures from 7.4 ± 0.3 to 4.3 ± 0.2 cal ka BP, which is likely the result of stronger summer insolation and pronounced warming. Reduced summer insolation thereafter results in decreased productivity and low growing season temperatures at Shireet Naiman Nuur from 4.3 ± 0.3 cal ka BP until present day. The globally acknowledged 4.2 ka event also appears as a pronounced cooling event at Shireet Naiman Nuur, and additional abrupt cooling events occurred during minima in total solar irradiance at ∼3.4, 2.8 and 2.4 ka BP. Low lake primary productivity and growing season temperatures are likely the result of longer ice cover periods at the high-altitude (2,429 m a.s.l.) Shireet Naiman Nuur. This leads to shorter mixing periods of the lake water which is supported by more positive δ13CTOC because of increased incorporation of dissolved HCO3− by aquatic producers during periods of longer ice cover.

Beschreibung: Soil moisture signatures provide a promising solution to overcome the difficulty of evaluating soil moisture dynamics in hydrologic models. Soil moisture signatures are metrics that quantify the dynamic aspects of soil moisture timeseries and enable process-based model evaluations. To date, soil moisture signatures have been tested only under limited land-use types. In this study, we explore soil moisture signatures' ability to discriminate different dynamics among contrasting land-uses. We applied a set of nine soil moisture signatures to datasets from six in-situ soil moisture networks worldwide. The dataset covered a range of land-use types, including forested and deforested areas, shallow groundwater areas, wetlands, urban areas, grazed areas, and cropland areas. Our set of signatures characterized soil moisture dynamics at three temporal scales: event, season, and a complete timeseries. Statistical assessment of extracted signatures showed that (1) event-based signatures can distinguish different dynamics for all the land-uses, (2) season-based signatures can distinguish different dynamics for some types of land-uses (deforested vs. forested, urban vs. greenspace, and cropped vs. grazed vs. grassland contrasts), (3) timeseries-based signatures can distinguish different dynamics for some types of land-uses (deforested vs. forested, urban vs. greenspace, shallow vs. deep groundwater, wetland vs. non-wetland, and cropped vs. grazed vs. grassland contrasts). Further, we compared signature-based process interpretations against literature knowledge; event-based and timeseries-based signatures generally matched well with previous process understandings from literature, but season-based signatures did not. This study will be a useful guideline for understanding how catchment-scale soil moisture dynamics in various land-uses can be described using a standardized set of hydrologically relevant metrics.

Beschreibung: Apatite is a ubiquitous accessory mineral in crustal rocks. The Sr-isotope record of apatite has a wide range of applications in earth science studies. However, apatite has been documented to be easily altered by certain fluids. Currently, the impact of metasomatic alteration on Sr-isotopic abundances in apatite is not well known. In order to better understand this issue, well-characterized fluorapatite grains have been metasomatized experimentally at a temperature of 800 or 600 °C and a pressure of 200 MPa. Fluids used included a CO2-H2O mixture and NaF-, CaCl2-, and HCl-bearing solutions, all of which were doped with a standard solution with a known 87Sr/86Sr ratio. In the fluorapatite + CO2-H2O experiments, the fluorapatite grains were not altered by the fluids; thus, their Sr isotope compositions were generally kept unchanged. However, the other fluids induced partial to complete alteration of fluorapatite. In experiments involving the NaF- or HCl-bearing solutions, the Sr content remained constant or was increased in altered areas of the fluorapatite, and the Sr isotopes underwent changes with partial isotopic signature of the reacting solutions. In experiments involving CaCl2, the Sr content was decreased in the altered fluorapatite because high Ca activity in the solution caused Ca to replace Sr on the Ca site. Notably, the Sr isotopic ratios are still changed, although to a relatively small extent. The efficient Sr isotopic exchange between fluid and apatite is attributed to the rapid transport rate of Sr between the reaction-interface fluid and bulk fluid surrounding the apatite. This experimental study demonstrates that the response of apatite Sr isotopes to metasomatic alteration is mainly controlled by the chemistry of fluids. Overall, Sr isotopes become susceptible to hydrothermal alteration once the apatite is chemically reactive with the fluids. Therefore, it is important to evaluate the fluid-rock history, especially the conjectured fluid composition before using Sr isotopes from apatite as a geochemical tracer.

Beschreibung: Thaw consolidation of degrading permafrost is a serious hazard to the safety and operation of infrastructure. Monitoring thermal changes in the active layer (AL), the proportion of the soil above permafrost that thaws and freezes periodically, is critical to understanding the conditions of the top layer above the permafrost and regulating the construction, operation, and maintenance of facilities. However, this is a very challenging task using ground-based methods such as ground-penetrating radar (GPR) or temperature sensors. This study explores the integration of interferometric measurements from high-resolution X-band Synthetic Aperture Radar (SAR) images and volumetric water content (VWC) data from SoilGrids to quantify detailed spatial variations in active layer thickness (ALT) in Iqaluit, the territorial capital of Nunavut in Canada. A total of 21 SAR images from COSMO Sky-Med (CSK) were first analyzed using the freely connected network interferometric synthetic aperture radar (FCNInSAR) method to map spatial and temporal variations in ground surface subsidence in the study area. Subsequently, we built an ALT retrieval model by introducing the thaw settlement coefficient, which takes soil properties and saturation state into account. The subsidence measurements from InSAR were then integrated with VWC extracted from the SoilGrids database to estimate changes in ALT. For validation, we conducted a comparison between estimated ALTs and in situ measurements in the airport sector. The InSAR survey identifies several sites of ground deformation at Iqaluit, subsiding at rates exceeding 80 mm/year. The subsidence rate changes along the runway coincide with frost cracks and ice-wedge furrows. The obtained ALTs, ranging from 0 to 5 m, vary significantly in different sediments. Maximum ALTs are found for rock areas, while shallow ALTs are distributed in the till blanket (Tb), the intertidal (Mi) sediments, and the alluvial flood plain (Afp) sediment units. The intersection of taxiway and runway has an AL thicker than other parts in the glaciomarine deltaic (GMd) sediments. Our study suggests that combining high-resolution SAR imagery with VWC data can provide more comprehensive ALT knowledge for hazard prevention and infrastructure operation in the permafrost zone.

Beschreibung: Cosmic-ray neutron sensing (CRNS) has become an effective method to measure soil moisture at a horizontal scale of hundreds of metres and a depth of decimetres. Recent studies proposed operating CRNS in a network with overlapping footprints in order to cover root-zone water dynamics at the small catchment scale and, at the same time, to represent spatial heterogeneity. In a joint field campaign from September to November 2020 (JFC-2020), five German research institutions deployed 15 CRNS sensors in the 0.4 km2 Wüstebach catchment (Eifel mountains, Germany). The catchment is dominantly forested (but includes a substantial fraction of open vegetation) and features a topographically distinct catchment boundary. In addition to the dense CRNS coverage, the campaign featured a unique combination of additional instruments and techniques: hydro-gravimetry (to detect water storage dynamics also below the root zone); ground-based and, for the first time, airborne CRNS roving; an extensive wireless soil sensor network, supplemented by manual measurements; and six weighable lysimeters. Together with comprehensive data from the long-term local research infrastructure, the published data set (available at https://doi.org/10.23728/b2share.756ca0485800474e9dc7f5949c63b872; Heistermann et al., 2022) will be a valuable asset in various research contexts: to advance the retrieval of landscape water storage from CRNS, wireless soil sensor networks, or hydrogravimetry; to identify scale-specific combinations of sensors and methods to represent soil moisture variability; to improve the understanding and simulation of land–atmosphere exchange as well as hydrological and hydrogeological processes at the hillslope and the catchment scale; and to support the retrieval of soil water content from airborne and spaceborne remote sensing platforms.

Beschreibung: Arctic permafrost landscapes have functioned as a global carbon sink for millennia. These landscapes are very heterogeneous, and the omnipresent water bodies within them act as a carbon source. Yet, few studies have focused on the impact of these water bodies on the landscape carbon budget. We deepen our understanding of carbon emissions from thermokarst ponds and constrain their impact by comparing carbon dioxide and methane fluxes from these ponds to fluxes from the surrounding tundra. We use eddy covariance measurements from a tower located at the border between a large pond and semi-terrestrial tundra. When we take the open-water areas of thermokarst ponds into account, our results show that the estimated summer carbon uptake of the polygonal tundra is 11 % lower. Further, the data show that open-water methane emissions are of a similar magnitude to polygonal tundra emissions. However, some parts of the pond's shoreline exhibit much higher emissions. This finding underlines the high spatial variability in methane emissions. We conclude that gas fluxes from thermokarst ponds can contribute significantly to the carbon budget of Arctic tundra landscapes. Consequently, changes in the water body distribution of tundra landscapes due to permafrost degradation may substantially impact the overall carbon budget of the Arctic.

Beschreibung: Major earthquakes, such as the Canterbury and Kaikoura events recorded in New Zealand in 2010-2011 and 2016 respectively, highlighted that floor systems can be heavily damaged. Quasi-static cyclic experimental tests of structural sub-assemblies can help to establish the seismic performance of structural systems. However, the experimental performance obtained with such tests is likely to be dependent on the loading protocol adopted. This paper provides an overview of the loading protocols which have been assumed in previous experimental activities, with emphasis on those adopted for testing floor systems. The paper also describes the procedure used to define the loading protocol applied in the testing of a large precast concrete floor diaphragm as part of the ReCast floor project jointly conducted by the University of Canterbury, the University of Auckland and BRANZ. Subsequently the limitations of current loading protocols for bi-directional testing are discussed. The relevance of local seismicity on bidirectional demand is demonstrated by examining a large dataset of records from the RESORCE database. It is concluded that bi-directional experimental testing be undertaken using at least two loading protocols that impose different ratios of demand in orthogonal directions.

Beschreibung: Barite (BaSO4) is a critical raw material with the major application as a weighing agent in drilling mud. Some of the largest barite deposits in the geological record are spatially associated with large clastic dominant (CD-type) Zn-Pb massive sulfide deposits. In the Selwyn Basin (Yukon, Canada), bedded barite deposits are hosted by black mudstones of the Late Devonian Canol Formation, where barite occurs in the matrix, laminations and nodules commonly in an assemblage with pyrite. The Late Devonian bedded barite also forms part of a larger province extending along western North America’s ancient continental margin. The bedded barite deposits have been variably interpreted, with genesis linked to either lower temperature diagenetic or high-temperature hydrothermal processes. In this study, in-situ isotopic microanalyses by secondary ion mass spectrometry (SIMS) have been used to determine the isotopic composition of sulfur (δ34S) and oxygen (δ18O) in barite and sulfur (δ34S) in pyrite. The δ34Sbarite (+37.1‰ to +67.9‰) and δ18Obarite (+8.8‰ and +23.9‰) values represent a substantial offset from Late Devonian seawater sulfate, consistent with precipitation from modified diagenetic pore fluids. The coexistence of the barite with the highly 34S-positive pyrite also demonstrates coprecipitation under progressively sulfate-limited conditions. A similar assemblage containing diagenetic pyrite and barite has been described in the pre-ore assemblage of Late Devonian CD-type deposits nearby in the Selwyn Basin, demonstrating this is the expression of a regional diagenetic process. The diagenetic assemblage formed at the sulfate methane transition zone (SMTZ) where opposing diffusional fluxes of methane (+ barium) and sulfate interact, leading to sulfate reduction coupled to anaerobic methane oxidation (SR-AOM). The replacement of diagenetic barite by other Ba-bearing minerals (e.g. witherite, cymrite and hyalophane) provides further evidence of severe sulfate depletion and conditions under which barite was soluble.

Beschreibung: Groundwater (GW) is the world’s largest distributed freshwater storage for mankind, ecosystems, and is a key resource for industrial and agricultural demands. Due to its fundamental role in the Earth's water and energy cycles, groundwater has been declared as an Essential Climate Variable (ECV) by GCOS, the Global Climate Observing System. However, within Copernicus - the European Earth Observation Programme - there is no service available yet to deliver data on this fundamental resource, nor is there any other data source worldwide that operationally provides information on changing groundwater resources in a consistent way, observation-based, and with global coverage. Therefore, the Global Gravity-based Groundwater Product (G3P) project aims at developing an operational global groundwater service as a cross-cutting extension of the existing Copernicus - the European Earth Observation Programme - portfolio. G3P capitalizes from the unique capability of GRACE and GRACE-FO satellite gravimetry as the only remote sensing technology to monitor subsurface mass variations, and from other satellite-based water storage products that are already part of the Copernicus portfolio, to provide a data set of groundwater storage change for large areas with global coverage. G3P is obtained by using a mass balance approach, i.e., by subtracting satellite-based water storage compartments (WSCs) such as snow water equivalent, root-zone soil moisture, glacier mass, and surface water storage from GRACE/GRACE-FO monthly terrestrial water storage anomalies (TWSA). Compatibility of the observation-based WSCs with TWSA is achieved by a filtering process, where optimal filter types were derived by analyses of spatial correlation patterns. G3P groundwater variations are provided for almost two decades (from 2002 to the present), with the monthly resolution, and at a 0.5-degree spatial resolution globally. In this contribution, we also illustrate preliminary results of the G3P data set and of its uncertainties, as well as its evaluation by independent groundwater data. This study has been run in the context of the European Union’s Horizon 2020 research project G3P (Global Gravity-based Groundwater Product, grant agreement nº 870353).

Beschreibung: While there has been a lot of work focusing on improving our understanding of divergent and convergent plate boundaries, the intricate nature of back-arc extension, where subduction and large-scale extension occur and interact in close, is yet to be explored properly. It has long been proposed that the strength of the subduction interface, which depends among others on the amount of subducted sediments, plays a pivotal role in subduction dynamics. Here, we investigate the role of back-arc rheology and subduction interface strength on the deformation style of the overriding plate. Using two-dimensional thermomechanical model experiments, we demonstrate, that the presence of a weak mantle–lithospheric domain in the overriding plate can result in back-arc breakup even during the subduction of narrow, land-locked oceanic basins such as those found in the Mediterranean region. The thinning of the back-arc mantle–lithosphere results in a weaker overriding plate, hence a lower slab-pull force is sufficient to initiate back-arc extension. Convective thinning at the subduction interface also reduces the length of the interface, reducing the portion of slab-pull lost as energy dissipation. A weak plate interface, can also reduce the energy dissipated along the subduction zone, leading to earlier extension. A detailed analysis of the forces shaping the overriding plate stress field shows that transmission of slab-pull force has a predominant role while viscous basal drag has a negligible effect in our experiments. Our results compare favorably with large-scale characteristics of land-locked Mediterranean back-arc basins such as the North Tyrrhenian basin and the Pannonian basin.

Beschreibung: Back-arc basins such as the ones behind the island-arcs of the Western Pacific Ocean or the ones in the Mediterranean Sea are ubiquitous structures of the Earth. They are extensional basins forming in the overriding plate behind subduction zones and similarly to continental rifts, they can exhibit different structural styles from narrow, localized rifting to wide-rift extension. While these different structural styles have been long recognized, the factors controlling the style of extension in these basins have not been explored properly. We use thermo-mechanical models to investigate how the relative rates of progressive build-up of slab-pull force and of convective thinning and thermal weakening of the overriding plate control the style of back-arc rifting. Following subduction-initiation, a high subducting plate velocity results in rapid build-up of the slab-pull force. The relatively low rate of convectively thinning and associated moderate weakening of the overriding plate require slab-pull to build up to close to its maximum value to overcome the high back-arc integrated strength resulting in a narrow back-arc rift. In turn a low subducting plate velocity in comparison with the timescale of convective thinning of the overriding plate allows for significant back-arc weakening before the slab-pull force becomes large enough to drive back-arc extension. In this case, the back-arc exhibits a wide rifting style as extension occurs at significantly reduced overriding plate integrated strength. Our model results provide an explanation why some subduction zones exhibit wide, distributed extension in the overriding plate such as for instance observed in the Pannonian basin.

Beschreibung: Mg-carbonate is one of the main carbonates subducting from the Earth’s surface into the mantle and possible reactions of this compound with the main mantle minerals are of paramount importance for our understanding of the deep carbon cycle. Recently, we have theoretically shown that at pressures and temperatures of the Earth’s lower mantle MgCO3 reacts with MgO, producing Mg2CO4. Here, using Raman spectroscopy and in situ X-ray powder and single crystal diffraction techniques we have confirmed that MgCO3 does react with MgO at pressures above 50 GPa. At 2000-3000 K, the phase similar to the earlier predicted Mg2CO4-P21/c was observed in X-ray powder diffraction and Raman spectra. Using transmission electron microscopy we have confirmed the preservation of the orthocarbonate in the fully decompressed sample consistent with the theoretical prediction. A single crystal of another new phase of trigonal or hexagonal symmetry with cell parameters a=2.693(4), c=6.517(12) Å was identified at 50-60 GPa in the reaction products. The crystal structure of this phase was not solved. At heating above 3000 K, peaks of Mg2CO4 disappeared from the experimental spectra and the phase similar to MgO2-I4/mcm was detected in the quenched products. Our discoveries indicate that carbon-bearing lower mantle mineralogy can be more complex than previously thought and phases such as Mg-orthocarbonates can be quenched to the ambient conditions.

Beschreibung: Van Allen Probes measurements revealed the presence of the most unusual structures in the ultra-relativistic radiation belts. Detailed modeling, analysis of pitch angle distributions, analysis of the difference between relativistic and ultra-realistic electron evolution, along with theoretical studies of the scattering and wave growth, all indicate that electromagnetic ion cyclotron (EMIC) waves can produce a very efficient loss of the ultra-relativistic electrons in the heart of the radiation belts. Moreover, a detailed analysis of the profiles of phase space densities provides direct evidence for localized loss by EMIC waves. The evolution of multi-MeV fluxes shows dramatic and very sudden enhancements of electrons for selected storms. Analysis of phase space density profiles reveals that growing peaks at different values of the first invariant are formed at approximately the same radial distance from the Earth and show the sequential formation of the peaks from lower to higher energies, indicating that local energy diffusion is the dominant source of the acceleration from MeV to multi-MeV energies. Further simultaneous analysis of the background density and ultra-relativistic electron fluxes shows that the acceleration to multi-MeV energies only occurs when plasma density is significantly depleted outside of the plasmasphere, which is consistent with the modeling of acceleration due to chorus waves.

Beschreibung: Located in the middle of Shanxi Province, northern China, Taiyuan basin is a dry and water-short region. This region is reaching alarming levels of aquifer depletion due to decades of groundwater overexploitation, which has caused severe land subsidence in the basin. The Wanjiazhai Water Diversion Project (WWDP) was designed to ease water scarcity by transporting water from the Yellow River to the Taiyuan basin through 452.4 km-long canals. By the end of 2018, the WWDP had supplied 2.87 billion m3 of water to Shanxi Province, which replenish the basin’s surface water body as well as the underground aquifer. The groundwater levels have continued to rise since 2003, with rising levels of more than 70 meters by 2018 in comparison with its low stand in 2000. In this study, we use 2007-2010 ENVISAT, ALOS-1 data, and 2017-2020 Sentinel-1 data to study the response of the basin’s aquifer to the groundwater rebound against the background of the water transfer project. We addressed the issue of tropospheric delay and its impact on the seasonal deformation by combing GACOS (Generic Atmospheric Correction Online Service) and a common-point stacking method. The accuracy improvement of deformation by this correction method was validated with measurements from seven continuous GPS stations in the basin. Groundwater rebound triggers ground uplift, which was identified in five areas by InSAR with a rate up to 25 mm/yr. The uplifting displacement time series are well correlated with groundwater level recovery. The land subsidence in the south of the basin continues but the rates decreased significantly in 2017-2020 detected from Sentinel-1 as compared to that in period 2007-2010 from ENVISAT and ALOS-1. All these uplifting signals and the decreasing rates of land subsidence found in Taiyuan city provide the indication that water management practices are successful in mitigating further subsidence. We found a significant seasonal displacement concentrated within the central region of the basin corresponding to the main irrigated areas in the Taiyuan basin. The maximum peak-to-peak amplitude is 43 mm observed from ENVISAT and decreases to 20 mm observed from Sentinel-1. The seasonal amplitudes change rapidly across faults, indicating that the fault is an effective barrier to across-fault fluid flow. To further quantify the causal relationships between water level and ground displacement, groundwater levels and ground displacement at three wells located near the area affected by significant seasonal land subsidence are analyzed by Cross Wavelet Transform (XWT) method. We found the time lags of about one month between land subsidence and the forcing groundwater level declines. Such a cross wavelet analysis with high spatial-temporal resolution therefore enables tracking the health of the aquifer system and highlights the system’s sustainability in aiding water resources allocation against the background of the water diversion project.

Beschreibung: We explore nonlinear site effects in the new Japanese ground-motion dataset compiled by Bahrampouri et al. (2020). Following the approach of Seyhan and Stewart (2014), we evaluate the decrease of soil amplification according to the increasing and corresponding ground motion on surface rock (VS30 760 m=s). To better predict the rock ground motion associated with each record, we take into account the between-event variability of the ground motion, and to better evaluate the impact of nonlinearity, we correct observed ground motion on soil by the site-specific linear amplification. Instead of grouping the stations by site-response proxy, we focus on individual stations with several strong-motion records. We develop a framework to test recently published nonlinear site amplification models against a linear site amplification model and compare the results with recent building codes that include nonlinearity. The results show that the site response varies greatly from site to site, indicating that conventional site proxies, such as VS30, are not sufficient to characterize nonlinear site response. Out of all of the Kiban–Kyoshin network stations, 20 stations are selected as having recorded sufficient data to be used in the test. Out of these 20 stations, five stations show signs of nonlinearity, that is, the nonlinear models performed better than the linear-amplification model for all periods T. For most sites, however, the linear site amplification models get the best score. This suggest that, for the range of predicted rock motion considered in this study (peak ground acceleration 〈0:2g), nonlinearity may not have a sufficiently large impact on soil ground motion to justify the use of nonlinear site terms in ground-motion functional forms and seismic building codes for such moderate-level shaking.

Beschreibung: Archaeomagnetic data are fundamental for our understanding of the evolution of Earth's magnetic field on centennial to millennial timescales. From the earliest studies of the Thelliers, Aitken, Nagata and others in the 1950s and 1960s, archaeomagnetic data have been vital for extending our knowledge of the field to times prior to observational measurements. Today, many thousands of archaeomagnetic data allow us to explore the geomagnetic field in more detail than ever before. Both regional time series of archaeomagnetic data and the inclusion of archaeomagnetic data in time-varying global spherical harmonic field models have revealed a range of newly discovered field behaviour. More sophisticated approaches to developing regional curves and global models have allowed us to resolve the field in certain regions more robustly and with greater resolution than previously possible. In this review we give an overview of the widely used global archaeomagnetic database GEOMAGIA50, discuss the methods used to obtain archaeomagnetic data, their challenges, and explore progress over the past twenty years in developing regional secular variation curves and global spherical harmonic models of the archaeomagnetic field. We end the review by covering what we see as the “grand challenges” in archaeomagnetism, including which regions of the world should be focussed on with regards to data acquisition.

Beschreibung: In both seismic and electromagnetic imaging the diffracted wavefield has gained importance in recent years. While seismic data is often acquired for a large range of different source-receiver offsets, ground-penetrating radar (GPR) acquisitions are mostly (near-) zero-offset. This characteristic inhibits the use of reflected waves for the estimation of depth velocities, which in turn increases the importance of a reliable imaging and characterization of the diffracted wavefield. In this study, we adapt a coherence-based workflow originally designed for seismic wavefields to ground-penetrating radar (GPR) data, which often exhibit similar wave propagation phenomena. The first step of the proposed workflow is the coherence-based imaging of the often predominant reflected wavefield, which in the second step is adaptively subtracted from the original data, resulting in an approximation of the diffracted wavefield. In the third step, we characterize the previously revealed diffracted wavefield by means of wavefront attributes, namely slopes and curvatures. In the fourth and final step, these wavefront attributes can be used for the estimation of depth velocities by means of wavefront tomography, an inversion scheme that provides both the localization of scatterers and a smooth velocity model of the subsurface. We demonstrate the wide applicability of the suggested workflow on two GPR field data examples provided by the USGS – one recorded in the aftermath of Hurricane Sandy on the shores of Long Beach Island, New Jersey, the other capturing the internal structure of Wolverine Glacier, Alaska.

Beschreibung: Seismic monitoring refers to the measurement of time-lapse changes of seismic wave velocities and is a frequently used technique to detect dynamic changes in the Earth‘s crust. Its applications include a broad range of topics, such as natural hazard assessment and structural health monitoring. To obtain reliable measurements, results are usually stacked over time. Thereby, temporal resolution is lost, which makes the measurement less sensitive to short-term environmental processes. Another problem is that conventional datasets often lack spatial density and velocity changes can only be attributed to large areas. Recently, distributed acoustic sensing (DAS) has gained a lot of attention as a way to achieve high spatial resolution at low cost. DAS is based on Rayleigh-scattering of photons within an optical fibre. Because measurements can be taken every few meters along the cable, the fibre is turned into a large seismic array that provides information about the Earth’s crust at unprecedented resolution. In our study, we explore the potential of DAS for monitoring studies. Specifically, we investigate how spatial stacking of DAS traces affects the measurements of velocity variations. We use data recorded by a 21-km-long dark fibre located on Reykjanes Pensinsula, Iceland. The cable is sampled with a channel spacing of 4 meters. We analyze the energy of the oceans microseism continuously recorded between March and September 2020. At first, we stack adjacent traces on the fibre in space. We then cross correlate the stacks to obtain approximations of the Green’s functions between different DAS-channels. By measuring changes in the coda waveform of the extracted seismograms, velocity variations can be inferred. Our analysis shows that spatial stacking improves the reliability of our measurements considerably. Because of that, less temporal stacking is required and the time resolution of our measurements can be increased. In addition, the enhancement of the data quality helps resolve velocity variations in space, allowing us to observe variations propagating along the cable over time. These velocity changes are likely linked to magmatic intrusions associated with a series of repeated uplifts on the Peninsula. Our results highlight the potential of DAS for improving the localization capabilities and accuracy of seismic monitoring studies.

Beschreibung: Volcanic explosions produce energy that propagates both in the subsurface as seismic waves and in the atmosphere as acoustic waves. We analyse thousands of explosions which occurred at different craters at Etna volcano (Italy) in 2018 and 2019. We recorded signals from infrasound sensors, geophones (GPH), broadband seismometers (BB) and Distributed Acoustic Sensing (DAS) with fibre optic cable. The instruments were deployed at Piano delle Concazze at about 2 to 2.5 km from the active craters, within (or onto) a ~300,000 m2 scoria layer deposited by recent volcanic eruptions. The DAS interrogator was setup inside the Pizzi Deneri Volcanic Observatory (~2800 m elevation). Infrasonic explosion records span over a large range of pressure amplitudes with the largest one reaching 130 Pa (peak to peak), with an energy of ca. 2.5x1011 J. In the DAS and the BB records, we find a 4-s long seismic “low frequency” signal (1-2 Hz) corresponding to the seismic waves, followed by a 2-s long “high-frequency” signal (16-21 Hz), induced by the infrasound pressure pulse. The infrasound sensors contain a 1-2 Hz infrasound pulse, but surprisingly no high frequency signal. At locations where the scoria layer is very thin or even non-existent, this high frequency signal is absent from both DAS strain-rate records and BB/GPH velocity seismograms. These observations suggest that the scoria layer is excited by the infrasound pressure pulse, leading to the resonance of lose material above more competent substratum. We relate the high frequency resonance to the layer thickness. Multichannel Analysis of Surface Wave from jumps performed along the fibre optic cable provide the structure of the subsurface, and confirm thicknesses derived from the explosion analysis. As not all captured explosions led to the observation of these high frequency resonance, we systematically analyze the amplitudes of the incident pressure wave versus the recorded strain and find a non-linear relationship between the two. This non-linear behaviour is likely to be found at other explosive volcanoes. Furthermore, our observations suggest it might also be triggered by other atmospheric pressure sources, like thunderstorms. This analysis can lead to a better understanding of acoustic-to-seismic ground coupling and near-surface rock response from natural, but also anthropogenic sources, such as fireworks and gas explosions.

Beschreibung: Proficiency testing (PT) is one of the few ways for an analytical laboratory to assess data quality under routine operating conditions. Here we report the results of Round 1 of the G-Chron PT programme, which is sponsored by the International Association of Geoanalysts. G-Chron is the first PT scheme devoted to the U-Th-Pb dating of mineral phases, primarily zircon, in geological materials. In this first round of G-Chron a total of 72 geochronology laboratories received the test material “Rak-17”, which previously had been characterized by seven well-established isotope dilution TIMS laboratories. A total of 63 of the PT participating laboratories reported data by the 15 December 2019 deadline. Here we both report and assess the measurement results submitted to this round. Our analysis provides a means for participating laboratories to assess their individual performance in relation to the isotope ages assigned, the experimental fitness-for-purpose criteria proposed by the scheme’s organisers and the results of similar laboratories participating in this round.

Beschreibung: The Earth is rotating around its rotation axis in an irregular manner. The Earth rotation axis and its orientation in space vary with respect to the reference system (both the terrestrial and the celestial) due to the wide range of processes that contribute to the rotation excitation. Therefore, the study of the Earth rotation can provide essential information concerning the Earth system. Spaceborne geodetic sensors can determine Earth orientation parameters (EOP), which fully describe Earth’s behaviour in space. The EOP are needed for several fields and applications such as fundamental astronomical and geodetic reference systems, precise satellite orbit determination, space navigation, and disaster prevention.\\ Over the past three decades, climate change has caused undesirable alterations in living organisms, human activities, and socio-economic aspects. Climate change is fluctuating and alters weather patterns such as precipitation patterns and sea and ocean levels. It also threatens the biodiversity of ecosystems, food security, and human health, and exacerbates natural disasters. The intensity and frequency of natural hazards are increasing with erratic distribution due to changes in the climate. Also, the level of vulnerability and zonation of risk are changed. Analysis of natural hazards, such as atmospheric and hydrological events, can help improve crisis management. Therefore, satellite observation data and simulated data derived from different atmospheric models are needed in order to model different types of hazards and risks, which can help early warning and prediction systems. Even though continuous sensor measurements and archive data (historical data/climate) are used for weather forecasting in developed countries, deadly flooding happened close to Stuttgart in southern Germany in May 2016, which might be avoided by a precise weather warning system. Therefore, real-time space geodetic technique data estimation is necessary to use as input data in weather prediction models. For the analysis of space geodetic techniques in (near) real-time, predictions of the EOP are required. EOP are made available by the International Earth Rotation and Reference Systems Service (IERS) Rapid Service Prediction Centre at USNO, Washington D.C., with a delay of hours to days. Accordingly, in the past, several methods were developed and applied for the EOP prediction. However, the accuracy of EOP prediction is still unsatisfactory, even for prediction of just a few days in the future.\\ To improve the EOP prediction accuracy, this study investigates the consistency between Earth rotation’s theories and observations. Moreover, the potentials of different geophysical phenomena are examined to better understand the interaction of different processes that affect the Earth rotation excitation with the time. Most of the Earth’s rotation theories and solutions are based on the location of the Earth’s principal axes of inertia (PAI). That location is defined by the second-degree Stokes coefficients of the geopotential, which are accurately observed by the Gravity Recovery and Climate Experiment (GRACE) and satellite laser ranging (SLR). In this study, the evolution of the Earth’s axes of inertia is analyzed for the first time. The presented results are remarkable, as the inertia axes do not move around a mean position fixed to a given terrestrial reference frame in the study period, but drift away from their initial location in a non-negligible manner.\\ Moreover, this study proposes a novel hybrid approach to predict EOP. There is a well-introduced stochastic method called copula-based analysis, and I combined it with singular spectrum analysis (SSA) for EOP prediction. I analyzed the potential of copula-based methods for predicting Earth rotation parameters that are derived from the combination of different satellite geodetic sensors and from other geophysical parameters like effective angular momentums. The copula is a statistical method that exploits linear and non-linear relationships between two or more variables by fitting a theoretical copula function into an empirical bivariate or multivariate distribution function. I introduced a hybrid prediction method that can be applied to other geophysical parameters is introduced in this thesis.\\ In this study, the interconnection between the celestial pole motion (CPM) and geomagnetic field (GMF) is investigated to improve the current CPM prediction methods. During the last decade, several investigations have been conducted in order to discuss a possible interconnection of polar motion and geomagnetic jerks, which are rapid changes in GMF secular variations. However, less attention has been paid to the impact of the GMF changes on the CPM, e.g., the interrelation of the geomagnetic jerks, geomagnetic dipole moment, geomagnetic field elements, and CPM variations. In this study, I use the CPM time series obtained from very long baseline interferometry (VLBI) observations and the latest GMF data to explore the correlation between CPM and the GMF. Our preliminary results revealed some impressive common features in the CPM and GMF variations, which show the potential to improve our understanding of the GMF’s contribution to the Earth’s rotation. All in all, the results mathematically illustrate the coherency between the GMF parameters and CPM, which helps improve EOP products.

Beschreibung: Die Erde dreht sich um ihre Rotationsachse auf unregelmäßige Art und Weise. Die Erdrotationsachse und ihre Orientierung im Raum variieren in Bezug auf das Referenzsystem (sowohl im terrestrischen wie auch im zälestischen System) aufgrund des breiten Spektrums von Prozessen, die zur Rotationsanregung beitragen. Daher kann die Untersuchung der Erdrotation wesentliche Informationen über das Erdsystem liefern. Weltraumgestützte geodätische Sensoren liefern Informationen über erdgebundene Orientierungsparameter (EOP), die das Verhalten der Erde im Weltraum vollständig beschreiben. EOP werden für verschiedene Bereiche und Anwendungen, wie beispielsweise für grundlegende astronomische und geodätische Referenzsysteme, benötigt. Des Weiteren sind sie auch für die präzise Bestimmung von Satellitenorbits und die weltraumgestützten Navigation, bis hin zu Anwendungen im Katastrophenschutz von Bedeutung.\\ In den letzten drei Jahrzehnten zeigen sich bereits die negativen Auswirkungen des Klimawandels auf Biosphäre, menschliche Aktivitäten und sozioökonomische Aspekte. Der Klimawandel unterliegt Fluktuationen und verändert die Wettermuster, wie z.B. die Niederschlagsverteilung, sowie Meeres- und Ozeanspiegel. Er bedroht aber auch die biologische Vielfalt der Ökosysteme, die Ernährungssicherheit, die menschliche Gesundheit und verschlimmert Naturkatastrophen. Die Intensität und Häufigkeit von Naturgefahren werden zunehmen, hierbei ist ihre Verteilung aufgrund von Klimaveränderungen allerdings unregelmäßig; auch der Grad von Schadensanfälligkeiten und die Einteilung von Risikozonen werden sich zukünftig ändern.\\ Die Analyse von Naturgefahren, wie atmosphärische und hydrologische Ereignisse, kann zur Verbesserung des Krisenmanagements beitragen. Daher werden Satellitenbeobachtungen und simulierte Daten, die von verschiedenen atmosphärischen Modellen abgeleitet werden, für die Gefahr- und Risikomodellierung benötigt; dies kann Frühwarn- und Vorhersagesysteme unterstützen. Obwohl kontinuierliche Sensormessungen und Archivdaten (historische Daten/Klimadaten) für die Wettervorhersage in entwickelten Ländern zur Verfügung stehen, kamen es im Mai 2016 in der Nähe von Stuttgart in Süddeutschland bei einer katastrophalen Überschwemmung zu Verlusten von Menschenleben, die möglicherweise durch ein präzises Wetterwarnsystem vermeidbar gewesen wären. \\ Eine Echtzeitschätzung der geodätischen Weltraumtechnik wäre notwendig, um sie als Eingangsdaten in Wettervorhersagemodellen zu verwenden. Für die Analyse von raumgeodätischen Techniken in (nahezu) Echtzeit sind Vorhersagen der EOP unerlässlich. EOP werden durch das Rapid Service Prediction Centre des International Earth Rotation and Reference Systems Service (IERS) am USNO, Washington D.C., mit einer Verzögerung von Stunden bis Tagen zur Verfügung gestellt. Demzufolge wurden in der Vergangenheit mehrere Methoden für die EOP-Vorhersage entwickelt und angewendet. Die Genauigkeit dieser EOP-Vorhersagen ist jedoch nach wie vor - selbst für einen Vorhersagezeitraum von nur wenigen Tagen - unbefriedigend.\\ Um die Genauigkeit der EOP-Vorhersage zu verbessern, beschäftigt sich die vorliegende Studie mit der Untersuchung von der Vereinbarkeit von Erdrotationstheorien mit Beobachtungen. Darüber hinaus wird das Potenzial verschiedener geophysikalischer Phänomene analysiert, um die Wechselwirkungen verschiedener Prozesse, die die zeitabhängige Anregung der Erdrotation beeinflussen, besser zu verstehen.\\ Die meisten Theorien und Lösungen zur Erdrotation basieren auf der Lage der Hauptträgheitsachsen der Erde (PAI). Diese Position wird durch die Stokes-Koeffizienten zweiten Grades des Geopotentials, das mit Hilfe der Satellitenmission Gravity Recovery and Climate Experiment (GRACE) und der Satelliten-Laser-Entfernungsmessung (SLR) genau beobachtet wird, definiert. In dieser Studie wird zum ersten Mal die Entwicklung der Trägheitsachsen der Erde analysiert. Bemerkenswert an den vorgestellten Ergebnissen ist, daß sich die Trägheitsachsen während der Untersuchungsperiode nicht um eine mittlere Position, die in einem bestimmten terrestrischen Bezugsrahmen festgelegt ist, bewegen, sondern sich in nicht zu vernachlässigender Weise von ihrem ursprünglichen Position entfernen.\\ Darüber hinaus schlägt diese Studie einen neuartigen hybriden Ansatz zur EOP-Vorhersage vor. Es gibt eine gut eingeführte stochastische Methode, die “kopula-basierte Analyse”, die wir mit der “Singulär-Spektrum-Analyse” (SSA) für die EOP-Vorhersage kombiniert haben . Wir analysierten das Potenzial kopula-basierter Methoden zur Vorhersage von Erdrotationsparametern, die aus der Kombination verschiedener geodätischer Satellitensensoren und aus anderen geophysikalischen Parametern, wie z.B effektiven Drehimpulsen, abgeleitet werden. Die Kopula ist eine statistische Methode, die lineare und nicht-lineare Beziehungen zwischen zwei oder mehreren Variablen nutzt, indem eine theoretische Kopula-Funktion an eine empirische, bivariate oder multivariate Verteilungsfunktion angepasst wird. Wir haben eine hybride Vorhersagemethode entwickelt, die auch auf andere geophysikalische Parameter angewendet werden kann.\\ In dieser Studie wird der Zusammenhang zwischen der Bewegung des Himmelspols (CPM) und dem geomagnetischen Feld (GMF) untersucht, um die derzeitigen CPM-Vorhersagemethoden zu verbessern. Während des letzten Jahrzehnts wurden mehrere Untersuchungen durchgeführt, um eine mögliche Verbindung zwischen polaren Bewegungen und geomagnetischen Ausbrüten - hierbei handelt es sich um rasche Veränderungen der säkularen Variationen des GMF - zu erörtern. Weniger Aufmerksamkeit wurde jedoch den Auswirkungen der GMF-Änderungen auf die CPM, z.B. der Wechselbeziehung der geomagnetischen Ausbrüte, des geomagnetischen Dipolmoments, der geomagnetischen Feldelemente und der CPM-Variationen, gewidmet. In dieser Studie verwenden wir CPM-Zeitreihen, die aus Beobachtungen der Very Long Baseline Interferometry (VLBI) gewonnen wurden und aktuelle GMF-Daten, um die Korrelation zwischen CPM und GMF zu untersuchen. Unsere vorläufigen Ergebnisse zeigen einige auffallente Gemeinsamkeiten in den CPM- und GMF-Variationen, die das Potenzial besitzen, unser Verständnis des GMF-Beitrags zur Erdrotation zu verbessern. Alles in allem veranschaulichen die Ergebnisse mathematisch die Kohärenz zwischen den GMF-Parametern und der CPM und weisen damit perspektivisch den Weg für eine Verbesserung der EOP-Produkte.

Beschreibung: Over the past decades, ground-penetrating radar (GPR) has become a fundamental tool in glaciological studies thanks to its tremendous capacity to provide high-resolution images in snow and ice. 3D acquisitions in particular can give detailed information on the internal structure, properties, and dynamics of glaciers. For imaging and highlighting important englacial and subglacial features such as meltwater tunnels and voids, an analysis of the spatial distribution of diffractions in the data holds great potential. However, the diffracted wavefield typically has low amplitude and is often masked by more prominent arrivals. Diffraction separation and imaging procedures have already become topics of significant interest in the field of exploration seismology, and may potentially open new possibilities for the analysis of glacier GPR data. Here, we explore the potential of recent advances in diffraction imaging for the analysis of alpine glacier GPR data. To this end, we consider a 3D data set acquired on the Haut Glacier d’Arolla (Valais, Switzerland) using a 70-MHz single-antenna real-time-sampling GPR system. The approach we use coherently approximates the dominant reflected wavefield and subtracts it from the data. The remaining diffracted wavefield is then enhanced using local coherent stacking. We find that this methodology is highly effective at isolating diffractions in glacier GPR data and provides clean images of the diffracting structures. Current work includes investigation of the correlation between these structures and the englacial and subglacial hydrological network.

Beschreibung: Cassiterite, the economically most important tin mineral, typically has moderate U and variable common Pb contents, making it amenable for U-Pb dating. Cassiterite has extremely low Th/U ratios (Th/U 〈 0.01) and its 208Pb is dominantly common Pb. This is particularly helpful as there is significant interference of tungsten oxides on 202Hg and 204Pb. The feasibility of the 208Pb correction procedure is discussed in detail. The 208Pb corrected LA-SF-ICP-MS data are in good agreement with intercept ages in the Tera-Wasserburg diagram and 207Pb corrected ages. Twelve cassiterite samples were investigated using the ID-TIMS and LA-SF-ICP-MS methods. The ID-TIMS results of Pit-AB, Rond-A, RG-114, BB#7 and 19GX cassiterite are reported for the first time in this study. RG-114, BB#7 and 19GX cassiterite have very low common Pb contents and are recommended for use as primary reference materials for in situ cassiterite. Pit-AB, Rond-A and Yankee cassiterite contain a small amount of common Pb, produce reliable and consistent ages and are suitable as primary reference materials. The remaining five cassiterite samples (Kard, Zinnwald, Els, XBD-W and Y724) were only investigated using the LA-SF-ICP-MS method and produce ages consistent with published age data from the host rocks associated with the tin deposits and with published U-Pb ages of cassiterite from the same deposits. We present an ID-TIMS Usingle bondPb of 154.3 ± 0.7 Ma for the commonly used cassiterite reference material AY-4. This age differs from previously reported ID-TIMS ages. This age discrepancy is caused by different initial common Pb compositions rather than age heterogeneity.

Beschreibung: Our planet is in crisis! The latest report of the United Nations Intergovernmental Panel on Climate Change (IPCC AR6) confirms that human influence is causing widespread, rapid, and intensifying changes in our weather and climate that are affecting every region on Earth in multiple ways. With every additional ton of carbon we emit, the frequency and intensity of storms, floods, droughts, and fires become greater and the effects on the environment and on human health and civilization become more severe. As geoscientists and journal editors, most of us have been accustomed to being on the leading edge of human knowledge and understanding of climate change, where we deal in objectivity, uncertainty, and debate, but now we find ourselves at the core of this climate crisis. It is no longer scientific discoveries at stake, but also humanity itself. This is an uncomfortable place for many of us. We are trained to be dispassionate observers and cautious thinkers, yet the alarming rate of recent climate change impels us to turn our attention directly toward mitigating this impending crisis. We are making a plea for collective action: we must make the switch to a green economy, put a just and effective price on carbon now, and consider a portfolio of other equitable public investments in climate solutions. These actions will ensure that the true costs and risks of burning fossil fuels are accounted for and global carbon emissions are rapidly reduced. Rich countries must lead the way in making drastic cuts to carbon emissions and in helping low- and middle-income countries to develop sustainably. We are running out of time. For decades, American Geophysical Union (AGU) journals have been at the forefront of documenting human-caused climate change and warning of a worsening climate crisis. Over 2,000 publications from AGU journals are cited in the new IPCC AR6 report. But we too can do more than just document and scientifically explain the ongoing crisis—our profession must help lead the way to solutions. Finding solutions and adapting to change have become not only necessary, but essential in ensuring safe, sustainable, and healthy human communities in the future. The geosciences have an essential role to play in these efforts by pivoting toward more cross-sector, solution-based science. To help lead this vision, the AGU is adding a new publication forum for community science in partnership with associations outside the geosciences. This forum will enhance interactions among AGU's existing, more disciplinary journals and give local communities a voice in leading solutions to global challenges. We are scientists, but we also have families and loved ones alongside our fellow citizens on this planet. The time to bridge the divide between scientist and citizen, head and heart, is now. The lead-up to the 2021 UN Climate Change Conference, COP26, being held in Glasgow in November, is our “last best chance” to urge world leaders to come together and commit to keeping climate change and its devastating impacts in check.