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

Description: This data publication contains a seismic survey which was acquired in the Mont Terri Underground Rock Laboratory (URL) in January 2019. The aim of the SI-A experiment (Seismic Imaging Ahead of and around underground infrastructure) is to provide a seismic characterization at the meso scale and to investigate the feasibility of tomographic and reflection imaging in argillaceous environments. The survey covered the different facies types of Opalinus Clay: shaly facies, carbonate -rich sandy facies and sandy facies (Bossart et al. 2017). Three different seismic sources (impact, vibro, ELVIS) were used to acquire the seismic data. The impact and magnetostrictive vibro sources were particularly designed for seismic exploration in the underground (Giese et al. 2005, Richter et al. 2018). The ELVIS source was mainly designed for near-surface investigations on roads or in open terrain (Krawczyk et al. 2012). All data were recorded on 32 3-component geophones (GS-14-L3, 28 Hz) which were deployed in 2 m deep boreholes, fixed at the tip of rock anchors. The data publication covers raw and preprocessed data stored in SEG-Y format.

Description: BACKGROUND AND OBJECTIVES: Soil or rock types in a region are often interpreted qualitatively by visually comparing various geophysical properties such as seismic wave velocity and vulnerability, as well as gravity data. Better insight and less human-dependent interpretation of soil types can be obtained from a joint analysis of separated and independent geophysical parameters. This paper discusses the application of a neural network approach to derive rock properties and seismic vulnerability from horizontal-to-vertical seismic ratio and seismic wave velocity data recorded in Majalengka-West Java, Indonesia.METHODS: Seismic microtremors were recorded at 54 locations and additionally multichannel analyses of surface wave experiments were performed at 18 locations because the multichannel analyses of surface wave experiment needs more effort and space. From the two methods, the values of the average shear wave velocity for the upper 30 meters, peak amplitudes and the dominant frequency between the measurement points were obtained from the interpolation of those geophysical data. Neural network was then applied to adaptively cluster and map the geophysical parameters. Four learning model clusters were developed from the three input seismic parameters: shear wave velocity, peak amplitude, and dominant frequency. FINDINGS: Generally, the values of the horizontal to vertical spectral ratios in the west of the study area were low (less than 5) compared with those in the southeastern part. The dominant frequency values in the west were mostly low at around 0.1–3 Hertz, associated with thick sedimentary layer. The pattern of the shear wave velocity map correlates with that of the horizontal to vertical spectral ratio map as the amplification is related to the soil or rock rigidity represented by the shear wave velocity. The combination of the geophysical data showed new features which is not found on the geological map such as in the eastern part of the study area. CONCLUSION: The application of the neural network based clustering analysis to the geophysical data revealed four rock types which are difficult to observe visually. The four clusters classified based on the variation of the geophysical parameters show a good correlation to rock types obtained from previous geological surveys. The clustering classified safe and vulnerable regions although detailed investigation is still required for confirmation before further development. This study demonstrates that low-cost geophysical experiments combined with neural network-based clustering can provide additional information which is important for seismic hazard mitigation in densely populated areas.

Description: The Walvis Ridge (WR) is the most prominent hotspot track related to the opening in the South Atlantic Ocean. Several hypotheses have been developed to explain its origin and evolution. The presence of a massive magmatic structure at the landfall of the WR in Northwest Namibia raised speculation about the role of a hotspot during the opening of the South Atlantic ocean. To investigate its deeper velocity structure at the junction of the WR with the African continent was the focus of the amphibious seismological WALPASS experiment. In total 12 ocean-bottom seismometers and 28 broad-band land stations were installed between 2010 and 2012 to acquire seismological data. Here, we present the results of seismic ambient noise tomography to investigate to which extent the Tristan hotspot modified the crustal structure in the landward prolongation of the ridge and in the adjacent oceanic basins. For the tomography, vertical and hydrophone component cross correlations for 〉300 d for OBS stations and between 1 and 2 yr for land stations data were analysed. More than 49 000 velocity measurements (742 dispersion curves) were inverted for group velocity maps at 75 individual signal periods, which then had been inverted for a regional 3-D shear wave velocity model. The resulting 3-D model reveals structural features of the crust related to the continent–ocean transition and its disturbance caused by the initial formation of the WR ∼130 Ma. We found relatively thick continental crust below Northwest Namibia and below the near-shore part of the WR, a strong asymmetry offshore with typical, thin oceanic crust in the Namibe Basin (crossing over into the Angola Basin further offshore) to the North and a wide zone of transitional crust towards the Walvis Basin south of the WR.