ALBERT

Description: 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.

Description: Due to the major role of greenhouse gas emissions in global climate change, the development of non-fossil energy technologies is essential. Deep geothermal energy represents such an alternative, which offers promising properties such as a high base load capability and a large untapped potential. The present work addresses barite precipitation within geothermal systems and the associated reduction in rock permeability, which is a major obstacle to maintaining high efficiency. In this context, hydro-geochemical models are essential to quantify and predict the effects of precipitation on the efficiency of a system. The objective of the present work is to quantify the induced injectivity loss using numerical and analytical reactive transport simulations. For the calculations, the fractured-porous reservoirs of the German geothermal regions North German Basin (NGB) and Upper Rhine Graben (URG) are considered. Similar depth-dependent precipitation potentials could be determined for both investigated regions (2.8-20.2 g/m3 fluid). However, the reservoir simulations indicate that the injectivity loss due to barite deposition in the NGB is significant (1.8%-6.4% per year) and the longevity of the system is affected as a result; this is especially true for deeper reservoirs (3000 m). In contrast, simulations of URG sites indicate a minor role of barite (〈 0.1%-1.2% injectivity loss per year). The key differences between the investigated regions are reservoir thicknesses and the presence of fractures in the rock, as well as the ionic strength of the fluids. The URG generally has fractured-porous reservoirs with much higher thicknesses, resulting in a greater distribution of precipitates in the subsurface. Furthermore, ionic strengths are higher in the NGB, which accelerates barite precipitation, causing it to occur more concentrated around the wellbore. The more concentrated the precipitates occur around the wellbore, the higher the injectivity loss. In this work, a workflow was developed within which numerical and analytical models can be used to estimate and quantify the risk of barite precipitation within the reservoir of geothermal systems. A key element is a newly developed analytical scaling score that provides a reliable estimate of induced injectivity loss. The key advantage of the presented approach compared to fully coupled reservoir simulations is its simplicity, which makes it more accessible to plant operators and decision makers. Thus, in particular, the scaling score can find wide application within geothermal energy, e.g., in the search for potential plant sites and the estimation of long-term efficiency.

Description: Aufgrund der tragenden Rolle der Treibhausgasemissionen für den globalen Klimawandel ist die Entwicklung von nicht-fossilen Energietechnologien essenziell. Die Tiefengeothermie stellt eine solche Alternative dar, welche vielversprechende Eigenschaften wie eine hohe Grundlastfähigkeit und ein großes ungenutztes Potenzial bietet. Die vorliegende Arbeit befasst sich mit Barytausfällungen inner- halb geothermaler Systeme und der damit einhergehenden Verringerung der Gesteinsdurchlässigkeit, welche ein Haupthindernis für die Aufrechterhaltung einer hohen Effizienz darstellen. Dabei sind hydro-geochemische Modelle unerlässlich, um die Auswirkungen von Ausfällungen auf die Effizienz eines Systems zu quantifizieren und vorherzusagen. Ziel der vorliegenden Arbeit ist es, mittels numerischer und analytischer reaktiver Transportsimulationen, den induzierten Injektivitätsverlust zu quantifizieren. Für die Berechnungen werden die klüftig-porösen Reservoire der deutschen Geothermieregionen Norddeutsches Becken (NDB) und Oberrheingraben (ORG) betrachtet. Für beide untersuchte Regionen konnte ein ähnliches, tiefenabhängiges Fällungspotenzial bestimmt werden (2,8–20,2 g/m3 Fluid). Die Reservoirsimulationen zeigen jedoch, dass der Injektivitätsverlust aufgrund von Barytablagerungen im NDB erheblich ist (1,8%–6,4% pro Jahr) und die Langlebigkeit der Anlage dadurch beeinträchtigt wird, dies gilt insbesondere für tiefere Reservoire (3000 m). Im Gegensatz dazu deuten die Simulationen der ORG-Standorte auf eine untergeordnete Rolle von Baryt hin (〈 0,1%–1,2% Injektivitätsverlust pro Jahr). Die entscheidenden Unterschiede zwischen den untersuchten Regionen sind die Reservoirmächtigkeiten und das Vorhandensein von Rissen im Gestein sowie die Ionenstärke der Fluide. Der ORG weist in der Regel klüftig-poröse Reservoire mit deutlich höheren Mächtigkeiten auf, was zu einer größeren Verteilung der Präzipitate im Untergrund führt. Weiterhin sind die Ionenstärken im NDB höher, was die Barytausfällung beschleunigt und diese dadurch konzentrierter um das Bohrloch herum entstehen lässt. Je konzentrierter die Präzipitate um die Bohrung herum auftreten, desto höher ist der Injektivitätsverlust. In dieser Arbeit wurde ein Workflow erarbeitet, innerhalb dessen mittels numerischer und analytischer Modelle das Risiko von Barytausfällungen innerhalb des Reservoirs geothermischer Systeme abgeschätzt und quantifiziert werden kann. Ein zentrales Element ist ein neu entwickelter, analytischer Scaling-Score, der eine zuverlässige Schätzung des induzierten Injektivitätsverlustes ermöglicht. Der entscheidende Vorteil des präsentierten Ansatzes im Vergleich zu voll-gekoppelten Reservoirsimulationen liegt in ihrer Einfachheit, die sie für Anlagenbetreiber und Entscheidungsträger zugänglicher macht. Somit kann insbesondere der Scaling-Score eine breite Anwendung innerhalb der Geothermie finden, z.B. bei der Suche nach potenziellen Anlagenstandorten und der Abschätzung der langfristigen Effizienz.

Description: The current western margin of the South American continent is an active subduction orogeny, which is famous for the ’Andean-type’ subduction. The Andes orogeny extends more than 6000 km along the margin. The Central Andes is governed by plateau-style deformation and the subduction of the Nazca plate below Chile and western Argentina, which has caused drastic crustal shortening and thickening, magmatism and periodical back-arc lithosphere and lower crust delamination. The deformation is also influenced by the periodical dip angle variations of the Nazca plate. However, the interaction zone between the overriding plate and the subducting Nazca slab is still enigmatic. In my dissertation, I use earthquake multi-scale full waveform inversion (FWI) to image the crust and upper mantle structure providing new robust images to resolve the interactions between the slab and the Central Andean plateau. First, I performed FWI using 117 earthquakes to investigate the seismic structure for the northern Central Andes which cover the normal-dip subduction zone from southern Peru to central Chile, proceeding from long period data (40–80 s) over several steps down to 12–60 s. In this study, the subducting Nazca slab is clearly imaged in the upper mantle, with dip-angle variations from the north to the south. Bands of low velocities in the crust and mantle wedge indicate intense crustal partial melting and hydration of the mantle wedge beneath the frontal volcanic arc, respectively. They are also linked to the vigorous dehydration from the subducting Nazca plate and intermediate-depth seismicity within the slab. These low-velocity bands are interrupted at 19.8 –21 S, both in the crust and uppermost mantle with an absence of active volcanoes, hinting at the lower extent of crustal partial melting and hydration of the mantle wedge, associated with the weak dehydration from the Nazca plate. The variation of lithospheric high-velocity anomalies below the backarc from north to south allows insight into the evolutionary foundering stages of the Central Andes. An extended high-velocity layer at lithospheric depths beneath the southern Altiplano suggests underthrusting of the leading edge of the Brazilian Shield following the removal of the autochthonous lithosphere. In contrast, a steeply westward dipping high-velocity block and low-velocity lithospheric uppermost mantle beneath the southern Puna plateau hints at the ongoing lithospheric delamination. Second, I took advantage of 134 events to perform FWI and started from even longer period data 60- 120 s to 12-100 s. In this study, a new seismic velocity model for the southern Central Andes is derived, covering the Pampean flat and adjacent Payenia steep subduction segments. Strong focused crustal low-velocity anomalies indicate partial melts in the Payenia segment along the volcanic arc, whereas weaker low-velocity anomalies covering a wide zone in Pampean possibly indicate remnant melts in the past. Thinning and tearing of the flat Nazca slab below the Pampean is inferred by gaps in the high-velocity slab along the inland projection of the Juan-Fernandez-Ridge. A high-velocity anomaly in the upper mantle below the flat slab is interpreted as a relic Nazca slab segment, which indicates an earlier slab break-off during the flattening process, triggered by the buoyancy of the Juan-Fernandez-Ridge. In Payenia, large-scale low-velocity anomalies atop and below the re-steepened Nazca slab are associated with the re-opening of the mantle wedge and sub-slab asthenospheric flow, respectively.

Description: The Scandes mountain range along the western rim of the Archean Baltic craton with elevation up to 2500 m forms an exceptional setting as the orogeny terminated 420 Ma ago and the Caledonides were deeply eroded afterwards. Since this region lacks recent compressional tectonic forces, a comprehensive explanation for the topography, which shows north-south and lateral variations along the Scandes, is missing. In my dissertation, I use earthquake surface waves and ambient noise to image the crustal and mantle structure aiming to provide new clues about the topography’s origin. The focus is also on exploring structural differences between the various tectonic domains. Here, I benefit from the seismic recordings by the ScanArray network supplemented by permanent and previous projects, distributed over entire Scandinavia. First, I performed a beamforming of Rayleigh surface waves which yielded average phase velocities for the study region and several of its sub-regions. An unusual 360° or sin(1θ) phase velocity variation with propagation azimuth is observed in northern Scandinavia and southern Norway/Sweden but not in the central area. For periods 〉35 s, a 5% variation between the maximum and minimum velocities was measured for opposite backazimuths of 120° and 300°, respectively. Such a variation is incompatible with the intrinsic azimuthal anisotropy and the path average approximation made in tomography. I assumed an eastward dipping lithosphere-asthenosphere boundary (LAB) to be the causing structure, inspired by some preliminary velocity models and observations made in previous studies. To test this hypothesis, I carried out 2D full-waveform modeling of the Rayleigh wave propagation. The models include a steep gradient at the LAB in combination with a pronounced reduction in the shear velocity below the LAB. The synthetic results are consistent with the observations: Faster phase velocities are obtained for propagation towards the thinning lithosphere, and lower ones for propagation in the direction of deepening LAB. The interference of reflected surface wave energy at the steep LAB with the forward propagating fundamental mode probably causes this peculiar effect. Second, the joint inversion of Rayleigh surface waves and ambient noise provides structural imaging down to 250 km depth. Resultant from my velocity model, I derive a new crustal model from which maps of the Moho depth as well as of the high-density lower crustal layer (LCL) are obtained. I observe crustal thickening from west to east below the Precambrian low-topography terranes, which is mainly a consequence of eastward thickening of the LCL. In contrast to the southern Scandes, with the overall highest topography (2,500 m), a crustal root below the northern Scandes (max. 2,100 m) is seen which diminish towards the central Scandes (max. 1,000 m). The LAB below the Scandes is deepening from west to east. The sharp steps in the LAB and strong velocity reductions both in the south (90–120 km LAB depth with 5.5% Vsv contrast) and the north (150 km LAB depth with 9% Vsv contrast) surprisingly correlate with the Caledonian mountain front. Whereas smoother laterally varying structures (150–170 km LAB depth with 4% Vsv contrast) are found below the central Scandes. The correlation of the lithosphere thickening with the Caledonian front might be related to metasomatism as result of the orogeny and/or the passive margin rifting. In Precambrian Scandinavia, low-velocity areas below 150 km depth are observed beneath the Archean Karelia craton in northern Finland. At mantle depth, the Paleoproterozoic Norrbotten craton can be separated from the Karelia craton, Caledonides and Paleoproterozic Svecofennian likely due to different degrees of metasomatism. Based on the structural differences, I conclude that different mechanisms are responsible for the compensation of the topography. The northern Scandes are likely compensated by a combined Airy-Pratt isostasy as implied by low-density rocks in the shallow crust (〈15 km depth), a high-density layer in the deep crust (〉10 km LCL thickness) and the mountain root. The strongly reduced velocities at sub-lithospheric depth additionally suggest an uplift contribution from the upper mantle. Since the southern Scandes lacks these crustal attributes, they experience mainly mantle-driven buoyancy. In both cases, however, I assume the influence of small-scale edge-driven convections (EDC) that can arise at sharp LAB gradients. EDC emplaces thereby low-density material at sub-lithospheric depths by the upwelling of hot asthenosphere which implies additional buoyancy of the lithosphere. Moreover, the lateral topography differences along the Scandes can be explained by varying EDC cell dimensions. Primarily, Pratt isostasy compensates the low topography central Scandes, but a contribution from dynamic support could act as well. Ultimately, I see the strong gradients at the LAB below the southern and northern Scandes as the cause of the observed 1θ phase velocity variation. While the smoother velocity structure in the central study area explains the absence of the 1θ effect.

Description: The occurrence and the style of volcanic eruptions are largely controlled by the ways in which magma is stored and transported from the mantle to the surface through the crust. Nevertheless, our understanding of the deep roots of volcano-magmatic systems remains very limited. Here, we use the sources of seismovolcanic tremor to delineate the active part of the magmatic system beneath the Klyuchevskoy Volcanic Group in Kamchatka, Russia. The tremor sources are distributed in a wide spatial region over the whole range of crustal depths connecting different volcanoes of the group. The tremor activity is characterized by rapid vertical and lateral migrations explained by fast pressure transients and dynamic permeability. Our results support the conceptual model of extended and highly dynamic trans-crustal magmatic systems.

Description: Seismic event detection and phase picking are the base of many seismological workflows. In recent years, several publications demonstrated that deep learning approaches significantly outperform classical approaches, achieving human-like performance under certain circumstances. However, as studies differ in the datasets and evaluation tasks, it is unclear how the different approaches compare to each other. Furthermore, there are no systematic studies about model performance in cross-domain scenarios, i.e., when applied to data with different characteristics. Here, we address these questions by conducting a large-scale benchmark. We compare six previously published deep learning models on eight datasets covering local to teleseismic distances and on three tasks: event detection, phase identification and onset time picking. Furthermore, we compare the results to a classical Baer-Kradolfer picker. Overall, we observe the best performance for EQTransformer, GPD and PhaseNet, with a small advantage for EQTransformer on teleseismic data. Furthermore, we conduct a cross-domain study, analyzing model performance on datasets they were not trained on. We show that trained models can be transferred between regions with only mild performance degradation, but models trained on regional data do not transfer well to teleseismic data. As deep learning for detection and picking is a rapidly evolving field, we ensured extensibility of our benchmark by building our code on standardized frameworks and making it openly accessible. This allows model developers to easily evaluate new models or performance on new datasets. Furthermore, we make all trained models available through the SeisBench framework, giving end-users an easy way to apply these models.

Description: The accurate knowledge of the Earth’s orientation and rotation in space is essential for a broad variety of scientific and societal applications. Among others, these include global positioning, near-Earth and deep-space navigation, the realisation of precise reference and time systems as well as studies of geodynamics and global change phenomena. In this paper, we present a refined strategy for processing and combining Very Long Baseline Interferometry (VLBI), Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) observations at the normal equation level and formulate recommendations for a consistent processing of the space-geodetic input data. Based on the developed strategy, we determine final and rapid Earth rotation parameter (ERP) solutions with low latency that also serve as the basis for a subsequent prediction of ERPs involving effective angular momentum data. Realising final ERPs on an accuracy level comparable to the final ERP benchmark solutions IERS 14C04 and JPL COMB2018, our strategy allows to enhance the consistency between final, rapid and predicted ERPs in terms of RMS differences by up to 50% compared to existing solutions. The findings of the study thus support the ambitious goals of the Global Geodetic Observing System (GGOS) in providing highly accurate and consistent time series of geodetic parameters for science and applications.

Description: This field campaign aimed at densifying the station coverage on the Armutlu Peninsula in the eastern Sea of Marmara. The Armutlu peninsula is directly crossed by the Armutlu fault, located roughly ~50 km away from the Istanbul metropolitan region. The main objective of this experiment is to characterize the seismic and aseismic deformation of this region. Waveform data are available from the GEOFON data centre, under network code 9P.

Description: This dataset contains subaquatic passive seismic recordings taken in September 2021 at 88 locations off Tuktoyaktuk Island as well as in a small lake (“Lake 3”) between the villages of Tuktoyaktuk and Inuvik, Northwest Territories, Canada. The measurements were part of the “Mackenzie Delta Permafrost Field Campaign” (mCan2021) within the “Modular Observation solutions for Earth Systems” (MOSES) program. Data is from a seismic intermediate-bandwidth seismic sensor lowered for few minutes to the bottom of the sea and lake, respectively, and from underwater short-period sensors deployed for a few days. The aim of the study was to determine the depth of the subaquatic permafrost (local lake and oceanic locations). Raw data is provided in proprietary “Cube” format and standard mseed format.