ALBERT

Description: Large igneous provinces (LIP) are vast (0.2 to 〉1 Mkm3) outpourings of basaltic lava and voluminous intrusions of magmas that have had important environmental consequences, in many cases leading to immense greenhouse gas release and mass extinctions. Magmatic oxygen fugacity (fO2) influences the chemistry of volcanic gases and is an important parameter for examining the links between LIP eruptions and environmental change. To constrain the fO2 of LIP magmas, we report olivine elemental chemistry of 399 crystals from a set of fifteen olivine-rich LIP samples, spanning in age from the Proterozoic (∼1270 Ma) to the Miocene (∼17 Ma). Concentrations of V in olivine are used to show that mafic LIP lavas erupted at +1.20 ± 0.95 ΔFMQ, on average more oxidized than mid ocean ridge basalts (MORB) at −0.28 ± 0.28 ΔFMQ. Mafic LIP magmas show a much larger range than MORB, however. Additionally, fO2 shows a negative correlation with parental magma MgO content, with high MgO lavas approaching the MORB range. This correlation is likely due to sampling of a heterogeneous mixture of oxidized and reduced lithologies, as also sampled by ocean island basalts (OIB). Correlation between fO2 and isotopic ratios such as 143Nd/144Nd demonstrates that the oxidized endmember is geochemically enriched, and may result from subduction recycling of oxidized surficial materials. The high fO2 of primitive LIP magmas demonstrate that they largely emitted oxidized gases during eruption, and furthermore, that LIP magmas associated with mass extinctions have similar magmatic fO2 to those that are not. Global plate tectonic position, magnitude and duration of LIP volcanic eruptions and magmatic degassing, as well as interaction with sedimentary basins in the crust - but not mantle source fO2 - are likely to be the critical factors for whether a LIP was associated with a mass extinction.

Description: Analysis and identifying the displacement characteristics play a key role in timely monitoring and detecting the physical responses of the bridge to ensure the safety of the human and structure. Many previous kinds of research used GNSS data to identify displacement and oscillation modelling of the bridge with different algorithms. This study uses GNSS time-series data to determine linear displacement and model oscillation of the bridge using a procedure including filtering outliers, linear regression, and sin function to identify amplitude in three directions, the plane displacement velocity, spatial displacement velocity, and vibration model of the bridge. The data in the research in the GNSS time-series data from three P5 GNSS receivers of the CHC brand on the Dachongyong bridge in Nanning, China with 1646 observations, at one-hour sample intervals in 68 consecutive days. The plane and spatial velocity of the three points DCQ01, DCQ02, and DCQ03 is 0.0181 mm/h, 0.0185 mm/h; 0.0114 mm/h, 0.0173 mm/h; and 0.0071 mm/h, 0.0082 mm/h respectively. The study results are significant in analyzing and identifying the bridge's displacement characteristics.

Description: Double seismic zones (DSZs) are a feature of some subducting slabs, where intermediate-depth earthquakes (~70–300 km) align along two separate planes. The upper seismic plane is generally attributed to dehydration embrittlement, whereas mechanisms forming the lower seismic plane are still debated. Thermal conductivity of slab minerals is expected to control the temperature evolution of subducting slabs, and therefore their seismicity. However, effects of the potential anisotropic thermal conductivity of layered serpentine minerals with crystal preferred orientation on slab’s thermal evolution remain poorly understood. Here we measure the lattice thermal conductivity of antigorite, a hydrous serpentine mineral, along its crystallographic b- and c-axis at relevant high pressure-temperature conditions of subduction. We find that antigorite’s thermal conductivity along the c-axis is ~3–4 folds smaller than the b-axis.Our numericalmodels further reveal thatwhen the lowthermal-conductivity c-axis is aligned normal to the slab dip, antigorite’s strongly anisotropic thermal conductivity enables heating at the top portion of the slab, facilitating dehydration embrittlement that causes the seismicity in the upper plane of DSZs. Potentially, the antigorite’s thermal insulating effect also hinders the dissipation of frictional heat inside shear zones, promoting thermal runaway along serpentinized faults that could trigger intermediatedepth earthquakes.

Description: Tourmaline is common in rare element pegmatites of the Nb-Y-F (NYF) type in the south-central part of the Proterozoic Sveconorwegian orogen in southern Norway. In the global context, however, tourmaline appears rare in this type of pegmatite. This study aims to explain the unusual tourmaline abundance in these pegmatites and the origin of boron (B) in the respective melts, and to raise awareness of tourmaline in NYF pegmatites generally. Tourmalines from six pegmatites in three Sveconorwegian lithotectonic units: Bamble, Kongsberg and Idefjorden, were investigated in terms of their mineral chemistry and δ11B values, in addition to bulk rock analyses of pegmatites and host rocks. Tourmalines in pegmatites from Bamble and Kongsberg record B isotopic compositions (δ11B = -1.0 to + 9.9 ‰) that are heavy relative to continental crust and mantle sources. In contrast, tourmaline in pegmatites and host rocks from Idefjorden have light B isotopic ratios (δ11B = -14.8 to −12.5 ‰) that are typical crustal values. We suggest that the latter melts were sourced from orthogneisses at depth. We relate the heavy B isotopic composition of Bamble and Kongsberg pegmatites to regional Na-metasomatism by fluids sourced from Mesoproterozoic shallow marine sediments. This is supported by previously published δ11B ratios from metasomatized Bamble host rocks. The spatial association of pegmatites with Na-metasomatism in the basement rocks suggests that metasomatism enhanced the fertility and B-concentration in the affected lithologies, favouring partial melting and the formation of tourmaline-bearing pegmatites. These findings contribute to understanding the petrogenesis of Sveconorwegian pegmatites but they also imply that B can play a greater role in the formation of NYF pegmatites than previously thought and that tourmaline has value as a petrogenetic tool in this type of pegmatites as well as in the Li-Cs-Ta (LCT) type to which is it is more commonly applied.

Description: Watershed management requires an understanding of key hydrochemical processes. The Pra Basin is one of the five major river basins in Ghana with a population of over 4.2 million people. Currently, water resources management faces challenges due to surface water pollution caused by the unregulated release of untreated household and industrial waste into aquatic ecosystems and illegal mining activities. This has increased the need for groundwater as the most reliable water supply. Our understanding of groundwater recharge mechanisms and chemical evolution in the basin has been inadequate, making effective management difficult. Therefore, the main objective of this work is to gain insight into the processes that determine the hydrogeochemical evolution of groundwater quality in the Pra Basin. The combined use of stable isotope, hydrochemistry, and water level data provides the basis for conceptualizing the chemical evolution of groundwater in the Pra Basin. For this purpose, the origin and evaporation rates of water infiltrating into the unsaturated zone were evaluated. In addition, Chloride Mass Balance (CMB) and Water Table Fluctuations (WTF) were considered to quantify groundwater recharge for the basin. Indices such as water quality index (WQI), sodium adsorption ratio (SAR), Wilcox diagram, and salinity (USSL) were used in this study to determine the quality of the resource for use as drinking water and for irrigation purposes. Due to the heterogeneity of the hydrochemical data, the statistical techniques of hierarchical cluster and factor analysis were applied to subdivide the data according to their spatial correlation. A conceptual hydrogeochemical model was developed and subsequently validated by applying combinatorial inverse and reaction pathway-based geochemical models to determine plausible mineral assemblages that control the chemical composition of the groundwater. The interactions between water and rock determine the groundwater quality in the Pra Basin. The results underline that the groundwater is of good quality and can be used for drinking water and irrigation purposes. It was demonstrated that there is a large groundwater potential to meet the entire Pra Basin’s current and future water demands. The main recharge area was identified as the northern zone, while the southern zone is the discharge area. The predominant influence of weathering of silicate minerals plays a key role in the chemical evolution of the groundwater. The work presented here provides fundamental insights into the hydrochemistry of the Pra Basin and provides data important to water managers for informed decision-making in planning and allocating water resources for various purposes. A novel inverse modelling approach was used in this study to identify different mineral compositions that determine the chemical evolution of groundwater in the Pra Basin. This modelling technique has the potential to simulate the composition of groundwater at the basin scale with large hydrochemical heterogeneity, using average water composition to represent established spatial groupings of water chemistry.

Description: Die Bewirtschaftung von Wassereinzugsgebieten erfordert ein Verständnis der wichtigsten hydrochemischen Prozesse. Das Pra-Becken ist eines der fünf großen Flusseinzugsgebiete Ghanas mit einer Bevölkerung von über 4,2 Millionen Menschen. Die Bewirtschaftung der Wasserressourcen wird derzeit durch die Verschmutzung der Oberflächengewässer erschwert, die durch die unkontrollierte Einleitung von unbehandelten Haushalts- und Industrieabfällen in die aquatischen Ökosysteme und durch illegale Bergbauaktivitäten entsteht. Dies hat den Bedarf an Grundwasser als zuverlässigste Wasserversorgung erhöht. Unser Verständnis der Mechanismen der Grundwasserneubildung und der chemischen Entwicklung im Einzugsgebiet ist bislang unzureichend, was eine wirksame Bewirtschaftung erschwert. Daher ist das Hauptziel dieser Arbeit Einblicke in die Prozesse zu bekommen, welche die hydrogeochemische Entwicklung der Grundwasserqualität im Pra-Becken bestimmen. Die kombinierte Verwendung von Daten stabiler Isotope, der Hydrochemie und von Wasserständen bildet die Grundlage für die Konzeption der chemischen Entwicklung des Grundwassers im Pra-Becken. Dafür wurden die Herkunft und die Verdunstungsraten des in die ungesättigte Zone infiltrierenden Wassers bewertet. Darüber hinaus wurden die Chlorid-Massenbilanz und die Wasserspiegelschwankungen betrachtet, um die Grundwasserneubildung für das Einzugsgebiet zu quantifizieren. Indizes wie der Wasserqualitätsindex (WQI), das Natriumadsorptionsverhältnis (SAR), das Wilcox-Diagramm und der Salzgehalt (USSL) wurden in dieser Studie verwendet, um die Qualität der Ressource für die Verwendung als Trinkwasser und zu Bewässerungszwecken zu bestimmen. Aufgrund der Heterogenität der hydrochemischen Daten wurden die statistischen Verfahren der hierarchischen Cluster- und Faktorenanalyse angewandt, um die Daten entsprechend ihrer räumlichen Korrelation zu unterteilen. Ein konzeptionelles hydrogeochemisches Modell wurde entwickelt und anschließend durch Anwendung kombinatorischer inverser und reaktionspfadbasierter geochemischer Modelle validiert, um plausible mineralische Assemblagen zu bestimmen, welche die chemische Zusammensetzung des Grundwassers kontrollieren. Die Wechselwirkungen zwischen Wasser und Gestein bestimmen die Grundwasserqualität im Pra-Becken. Die Ergebnisse unterstreichen, dass das Grundwasser eine gute Qualität aufweist und als Trinkwasser und für Bewässerungszwecke genutzt werden kann. Es wurde nachgewiesen, dass ein großes Grundwasserpotenzial vorhanden ist, um den derzeitigen und künftigen Wasserbedarf des gesamten Pra-Beckens zu decken. Als Hauptneubildungsgebiet wurde die nördliche Zone im Gebiet identifiziert, während die südliche Zone das Abflussgebiet ist. Der vorherrschende Einfluss der Verwitterung von Silikatmineralen spielt bei der chemischen Entwicklung des Grundwassers eine zentrale Rolle. Die hier vorgestellte Arbeit gibt grundlegende Einblicke in die Hydrochemie des Pra-Beckens und liefert für das Wassermanagement wichtige Daten für eine fundierte Entscheidungsfindung bei der Planung und Zuweisung von Wasserressourcen für verschiedene Zwecke. Ein neuartiger Ansatz zur inversen Modellierungwurde in dieser Studie eingesetzt, um unterschiedliche Mineralzusammensetzungen zu ermitteln, welche die chemische Entwicklung des Grundwassers im Pra-Becken bestimmen. Diese Modellierungstechnik hat das Potenzial, die Zusammensetzung eines Grundwassers auf der Skala eines Beckens mit großer hydrochemischer Heterogenität zu simulieren, wobei die durchschnittliche Wasserzusammensetzung zur Darstellung der etablierten räumlichen Gruppierungen der Wasserchemie verwendet wird.

Description: Watershed management requires an understanding of key hydrochemical processes. The Pra Basin is one of the five major river basins in Ghana with a population of over 4.2 million people. Currently, water resources management faces challenges due to surface water pollution caused by the unregulated release of untreated household and industrial waste into aquatic ecosystems and illegal mining activities. This has increased the need for groundwater as the most reliable water supply. Our understanding of groundwater recharge mechanisms and chemical evolution in the basin has been inadequate, making effective management difficult. Therefore, the main objective of this work is to gain insight into the processes that determine the hydrogeochemical evolution of groundwater quality in the Pra Basin. The combined use of stable isotope, hydrochemistry, and water level data provides the basis for conceptualizing the chemical evolution of groundwater in the Pra Basin. For this purpose, the origin and evaporation rates of water infiltrating into the unsaturated zone were evaluated. In addition, Chloride Mass Balance (CMB) and Water Table Fluctuations (WTF) were considered to quantify groundwater recharge for the basin. Indices such as water quality index (WQI), sodium adsorption ratio (SAR), Wilcox diagram, and salinity (USSL) were used in this study to determine the quality of the resource for use as drinking water and for irrigation purposes. Due to the heterogeneity of the hydrochemical data, the statistical techniques of hierarchical cluster and factor analysis were applied to subdivide the data according to their spatial correlation. A conceptual hydrogeochemical model was developed and subsequently validated by applying combinatorial inverse and reaction pathway-based geochemical models to determine plausible mineral assemblages that control the chemical composition of the groundwater. The interactions between water and rock determine the groundwater quality in the Pra Basin. The results underline that the groundwater is of good quality and can be used for drinking water and irrigation purposes. It was demonstrated that there is a large groundwater potential to meet the entire Pra Basin’s current and future water demands. The main recharge area was identified as the northern zone, while the southern zone is the discharge area. The predominant influence of weathering of silicate minerals plays a key role in the chemical evolution of the groundwater. The work presented here provides fundamental insights into the hydrochemistry of the Pra Basin and provides data important to water managers for informed decision-making in planning and allocating water resources for various purposes. A novel inverse modelling approach was used in this study to identify different mineral compositions that determine the chemical evolution of groundwater in the Pra Basin. This modelling technique has the potential to simulate the composition of groundwater at the basin scale with large hydrochemical heterogeneity, using average water composition to represent established spatial groupings of water chemistry.

Description: Rapidly growing seismic and macroseismic databases and simplified access to advanced machine learning methods have in recent years opened up vast opportunities to address challenges in engineering and strong motion seismology from novel, datacentric perspectives. In this thesis, I explore the opportunities of such perspectives for the tasks of ground motion modeling and rapid earthquake impact assessment, tasks with major implications for long-term earthquake disaster mitigation. In my first study, I utilize the rich strong motion database from the Kanto basin, Japan, and apply the U-Net artificial neural network architecture to develop a deep learning based ground motion model. The operational prototype provides statistical estimates of expected ground shaking, given descriptions of a specific earthquake source, wave propagation paths, and geophysical site conditions. The U-Net interprets ground motion data in its spatial context, potentially taking into account, for example, the geological properties in the vicinity of observation sites. Predictions of ground motion intensity are thereby calibrated to individual observation sites and earthquake locations. The second study addresses the explicit incorporation of rupture forward directivity into ground motion modeling. Incorporation of this phenomenon, causing strong, pulse like ground shaking in the vicinity of earthquake sources, is usually associated with an intolerable increase in computational demand during probabilistic seismic hazard analysis (PSHA) calculations. I suggest an approach in which I utilize an artificial neural network to efficiently approximate the average, directivity-related adjustment to ground motion predictions for earthquake ruptures from the 2022 New Zealand National Seismic Hazard Model. The practical implementation in an actual PSHA calculation demonstrates the efficiency and operational readiness of my model. In a follow-up study, I present a proof of concept for an alternative strategy in which I target the generalizing applicability to ruptures other than those from the New Zealand National Seismic Hazard Model. In the third study, I address the usability of pseudo-intensity reports obtained from macroseismic observations by non-expert citizens for rapid impact assessment. I demonstrate that the statistical properties of pseudo-intensity collections describing the intensity of shaking are correlated with the societal impact of earthquakes. In a second step, I develop a probabilistic model that, within minutes of an event, quantifies the probability of an earthquake to cause considerable societal impact. Under certain conditions, such a quick and preliminary method might be useful to support decision makers in their efforts to organize auxiliary measures for earthquake disaster response while results from more elaborate impact assessment frameworks are not yet available. The application of machine learning methods to datasets that only partially reveal characteristics of Big Data, qualify the majority of results obtained in this thesis as explorative insights rather than ready-to-use solutions to real world problems. The practical usefulness of this work will be better assessed in the future by applying the approaches developed to growing and increasingly complex data sets.

Description: Das rapide Wachstum seismischer und makroseismischer Datenbanken und der vereinfachte Zugang zu fortschrittlichen Methoden aus dem Bereich des maschinellen Lernens haben in den letzen Jahren die datenfokussierte Betrachtung von Fragestellungen in der Seismologie ermöglicht. In dieser Arbeit erforsche ich das Potenzial solcher Betrachtungsweisen im Hinblick auf die Modellierung erdbebenbedingter Bodenerschütterungen und der raschen Einschätzung von gesellschaftlichen Erdbebenauswirkungen, Disziplinen von erheblicher Bedeutung für den langfristigen Erdbebenkatastrophenschutz in seismisch aktiven Regionen. In meiner ersten Studie nutze ich die Vielzahl an Bodenbewegungsdaten aus der Kanto Region in Japan, sowie eine spezielle neuronale Netzwerkarchitektur (U-Net) um ein Bodenbewegungsmodell zu entwickeln. Der einsatzbereite Prototyp liefert auf Basis der Charakterisierung von Erdbebenherden, Wellenausbreitungspfaden und Bodenbeschaffenheiten statistische Schätzungen der zu erwartenden Bodenerschütterungen. Das U-Net interpretiert Bodenbewegungsdaten im räumlichen Kontext, sodass etwa die geologischen Beschaffenheiten in der Umgebung von Messstationen mit einbezogen werden können. Auch die absoluten Koordinaten von Erdbebenherden und Messstationen werden berücksichtigt. Die zweite Studie behandelt die explizite Berücksichtigung richtungsabhängiger Verstärkungseffekte in der Bodenbewegungsmodellierung. Obwohl solche Effekte starke, impulsartige Erschütterungen in der Nähe von Erdbebenherden erzeugen, die eine erhebliche seismische Beanspruchung von Gebäuden darstellen, wird deren explizite Modellierung in der seismischen Gefährdungsabschätzung aufgrund des nicht vertretbaren Rechenaufwandes ausgelassen. Mit meinem, auf einem neuronalen Netzwerk basierenden, Ansatz schlage ich eine Methode vor, umdieses Vorhaben effizient für Erdbebenszenarien aus dem neuseeländischen seismischen Gefährdungsmodell für 2022 (NSHM) umzusetzen. Die Implementierung in einer seismischen Gefährdungsrechnung unterstreicht die Praktikabilität meines Modells. In einer anschließenden Machbarkeitsstudie untersuche ich einen alternativen Ansatz der auf die Anwendbarkeit auf beliebige Erdbebeszenarien abzielt. Die abschließende dritte Studie befasst sich mit dem potenziellen Nutzen der von makroseismischen Beobachtungen abgeleiteten pseudo-Erschütterungsintensitäten für die rasche Abschätzung von gesellschaftlichen Erdbebenauswirkungen. Ich zeige, dass sich aus den Merkmalen solcher Daten Schlussfolgerungen über die gesellschaftlichen Folgen eines Erdbebens ableiten lassen. Basierend darauf formuliere ich ein statistisches Modell, welches innerhalb weniger Minuten nach einem Erdbeben die Wahrscheinlichkeit für das Auftreten beachtlicher gesellschaftlicher Auswirkungen liefert. Ich komme zu dem Schluss, dass ein solches Modell, unter bestimmten Bedingungen, hilfreich sein könnte, um EntscheidungsträgerInnen in ihren Bestrebungen Hilfsmaßnahmen zu organisieren zu unterstützen. Die Anwendung von Methoden des maschinellen Lernens auf Datensätze die sich nur begrenzt als Big Data charakterisieren lassen, qualifizieren die Mehrheit der Ergebnisse dieser Arbeit als explorative Einblicke und weniger als einsatzbereite Lösungen für praktische Fragestellungen. Der praktische Nutzen dieser Arbeit wird sich in erst in Zukunft an der Anwendung der erarbeiteten Ansätze auf wachsende und zunehmend komplexe Datensätze final abschätzen lassen.

Description: Large parts of the Earth’s interior are inaccessible to direct observation, yet global geodynamic processes are governed by the physical material properties under extreme pressure and temperature conditions. It is therefore essential to investigate the deep Earth’s physical properties through in-situ laboratory experiments. With this goal in mind, the optical properties of mantle minerals at high pressure offer a unique way to determine a variety of physical properties, in a straight-forward, reproducible, and time-effective manner, thus providing valuable insights into the physical processes of the deep Earth. This thesis focusses on the system Mg-Fe-O, specifically on the optical properties of periclase (MgO) and its iron-bearing variant ferropericlase ((Mg,Fe)O), forming a major planetary building block. The primary objective is to establish links between physical material properties and optical properties. In particular the spin transition in ferropericlase, the second-most abundant phase of the lower mantle, is known to change the physical material properties. Although the spin transition region likely extends down to the core-mantle boundary, the ef-fects of the mixed-spin state, where both high- and low-spin state are present, remains poorly constrained. In the studies presented herein, we show how optical properties are linked to physical properties such as electrical conductivity, radiative thermal conductivity and viscosity. We also show how the optical properties reveal changes in the chemical bonding. Furthermore, we unveil how the chemical bonding, the optical and other physical properties are affected by the iron spin transition. We find opposing trends in the pres-sure dependence of the refractive index of MgO and (Mg,Fe)O. From 1 atm to ~140 GPa, the refractive index of MgO decreases by ~2.4% from 1.737 to 1.696 (±0.017). In contrast, the refractive index of (Mg0.87Fe0.13)O (Fp13) and (Mg0.76Fe0.24)O (Fp24) ferropericlase increases with pressure, likely because Fe Fe interactions between adjacent iron sites hinder a strong decrease of polarizability, as it is observed with increasing density in the case of pure MgO. An analysis of the index dispersion in MgO (decreasing by ~23% from 1 atm to ~103 GPa) reflects a widening of the band gap from ~7.4 eV at 1 atm to ~8.5 (±0.6) eV at ~103 GPa. The index dispersion (between 550 and 870 nm) of Fp13 reveals a decrease by a factor of ~3 over the spin transition range (~44–100 GPa). We show that the electrical band gap of ferropericlase significantly widens up to ~4.7 eV in the mixed spin region, equivalent to an increase by a factor of ~1.7. We propose that this is due to a lower electron mobility between adjacent Fe2+ sites of opposite spin, explaining the previously observed low electrical conductivity in the mixed spin region. From the study of absorbance spectra in Fp13, we show an increasing covalency of the Fe-O bond with pressure for high-spin ferropericlase, whereas in the low-spin state a trend to a more ionic nature of the Fe-O bond is observed, indicating a bond weakening effect of the spin transition. We found that the spin transition is ultimately caused by both an increase of the ligand field-splitting energy and a decreasing spin-pairing energy of high-spin Fe2+.

Description: Geodynamische Prozesse werden von den physikalischen Materialeigenschaften unter den extremen Druck- und Temperaturbedingungen des Erdinneren gesteuert, gerade diese Areale sind aber faktisch nicht für direkte Beobachtungen zugänglich. Umso wichtiger ist es, die physikalischen Eigenschaften unter Bedingungen des Erdinneren zu untersuchen. Mit diesem Ziel vor Augen erlaubt das Studium der optischen Eigenschaften von Mineralen des Erdmantels, eine große Bandbreite an physikalischen Materialeigenschaften, in einer einfachen, reproduzierbaren und effizienten Art und Weise zu bestimmen. Dadurch bieten sich wichtige Einblicke in die physikalischen Prozessen des Erdinneren. Die vorliegende Arbeit konzentriert sich auf das System Mg-Fe-O, im Speziellen auf Periklas (MgO) und seine Eisen-haltige Variante Ferroperiklas ((Mg,Fe)O), ein wichtiger Baustein planetarer Körper. Das Hauptziel der Arbeit besteht darin Verbindungen zwischen optischen Eigenschaften und physikalischen Materialeigenschaften zu finden. Gerade der Spin-Übergang in Ferroperiklas, der zweithäufigsten Phase des unteren Erdmantels, ist dabei von Bedeutung, da damit Veränderungen in den physikalischen Materialeigenschaften einhergehen. Obwohl sich der Spinübergangsbereich vermutlich bis zur Kern-Mantel-Grenze erstreckt, sind die Auswirkungen des gemischten Spin-Zustandes, bei dem sowohl Hoch- als auch Tief-Spin präsent sind, nur unzureichend untersucht. Die hier vorgestellten Studien zeigen, wie optische Eigenschaften mit anderen wichtigen physikalischen Eigenschaften wie elektrischer und thermischer Leitfähigkeit, Viskosität oder auch mit der chemischen Bindung verbunden sind. Daraus lässt sich auch ableiten wie der Spin-Übergang in Ferroperiklas diese Eigenschaften beeinflusst. Von Raumbedingungen bis zu ~140 GPa sinkt der Brechungsindex von MgO um ~2.4 % von 1.737 auf 1.696 (±0.017). Im Gegensatz dazu steigt der Brechungsindex von (Mg0.87Fe0.13)O (Fp13) und (Mg0.76Fe0.24)O (Fp24) Ferroperiklas mit dem Druck an. Dies ist auf Fe-Fe Wechselwirkungen zwischen benachbarten Eisenpositionen zurückzuführen, die eine starke Verringerung der Polarisierbarkeit, wie im Falle von reinem MgO mit zunehmender Dichte, behindern. Eine Analyse der Dispersion des Brechungsindexes von MgO (Abnahme um ~23 % von 1 Atm zu ~103 GPa) offenbart eine Verbreiterung der Bandlücke von ~7.4 eV bei 1 Atm zu ~8.5 (±0.6) eV bei ~103 GPa. Die Messung der Dispersion (zwischen 550 und 870 nm) in Fp13 zeigt eine starke Abnahme über den Bereich des Spin-Überganges (~44–100 GPa) bis zu einem Faktor von ~3. Die Bandlücke nimmt in der Region des gemischten Spin-Zustandes signifikant auf bis zu ~4.7 eV zu (entspricht einer Zunahme um den Faktor ~1.7). Dies deutet auf eine Verringerung der Elektronen-Mobilität zwischen benachbarten Fe2+-Positionen mit unterschiedlichem Spin-Zustand hin, was die bereits in früheren Arbeiten beobachtete Abnahme der elektrischen Leitfähigkeit im Bereich des gemischten Spin-Zustandes erklärt. Absorptionsspektren an Fp13 zeigen eine Druck-bedingte Zunahme der Kovalenz der Fe-O Bindung für Ferroperiklas im Hoch-Spin Zustand, wohingegen Tief-Spin Ferroperiklas einen Trend zu einer mehr ionischen Fe-O Bindung auf-weist, was auf einen Bindungs-schwächenden Effekt des Spin-Wechsels hinweist. Der Übergang von Hoch- zu Tiefspin ist letztlich auf eine Zunahme der Ligandenfeldaufspaltungsenergie sowie eine abnehmende Spinpaarungsenergie von Hoch-Spin Fe2+ zurückzuführen.

Description: Ionospheric disturbances are associated with the propagation of seismic waves generated by earthquakes. We analyze total electron content (TEC) fluctuations from Global Navigation Satellite Systems (GNSS) data. This study investigates the global earthquakes from December 2004 to March 2020 with magnitude (Mw) 4.0 to 9.1. Data of 65 permanent GNSS stations are used to analyze the impacts of these earthquakes on the ionosphere. The experimental data based on the cluster sampling method ensure strict conditions such as accuracy, the distance from the monitoring stations to the epicenter, and the depth of the hypocenter. At least three GNSS stations near the epicenter participate in the analysis of each earthquake. Probability and statistics are applied to remove outliers and rough errors in the input data, to select datasets with similar quality, and to analyze TEC anomalies. The results show that when a strong earthquake occurs, the TEC values calculated at different GNSS stations surrounding the earthquake region tend to exhibit similar variations. Depending on the magnitude, epicenter, and hypocenter depth, these fluctuations can range from ± 3.2 TECU to ± 14.5 TECU for large earthquakes. These TEC fluctuations occur from 30 min to almost two hours before the mainshock of the earthquakes and last to the aftershock period. For earthquakes with a magnitude greater than 6.0 Mw, the TEC fluctuations are significant. The findings of this study contribute to GNSS applications in studying earthquakes in the future.

Description: In deformation analysis, irregularly spaced data, extreme values, and anomalies in time series can lead to misleading simulations for forecast models, such as overfitting and underfitting. Therefore, K-fold cross-validation is one of the hyperparameter tuning techniques used in machine learning (ML) to deal with these problems. In this study, we use data from 22 permanent GNSS stations to predict the motion trajectory of the Earth’s crust. Lag functions and sampling techniques are applied to generate 924-time series samples. Time series standardization techniques are also performed to improve the quality of data. To test the efficiency of the K-fold cross-validation method, we investigate 26 mathematical models based on six ML algorithms. The optimal K values are selected through trial methods. Root mean squared error (RMSE) of validation and test is the basis for determining the overfitting and underfitting models. The investigations show that the optimal intervals of K-fold range from five to ten folds for the GNSS time series with many anomalies, jumps, and significant variations, from three to ten for stable time series. The sensitivity of cross-validation is more effective on the time series of the Up component than those of the North and East components. In addition, cross-validation can also detect effectively overfitting and underfitting for forecast models in motion of permanent GNSS stations.

Description: Understanding the phase behavior and structural properties of salt water at high pressures is essential for understanding the dynamics and physical characteristics of icy planets. In this study, we employed high-pressure experimental and ab initio simulation techniques to investigate the impact of CaCl2 on the structure of ice VII. Our findings reveal that 1.8 mol% CaCl2 can be incorporated into the ice VII structure above 10 GPa. This CaCl2-bearing ice VII (Cb VII) exhibits a lower O-H stretching frequency in the Raman spectra as well as a reduced volume of the unit cell compared to pure ice VII. In contrast to doping ice VII with other salts such as LiCl and NaCl that leads to an increase of the ice VII to ice X transition pressure occurring at 100–150 GPa, CaCl2 doping stands out by reducing the transition pressure. It shifts the transition to a pressure of 52 GPa, which is significantly lower than the transition pressure of 80 GPa in the pure H2O ice system. This notable distinction highlights the unique influence of CaCl2 on the phase behavior of water under high pressure, and we attribute these effects to the phenomenon of chemical pressure induced by CaCl2 within the ice VII structure. Our study suggests that the presence of a modified ice VII phase, contaminated with salt and referred to as Cb VII, may influence the composition, structure, and evolution of planets.

Description: High-resolution thermospheric mass density low Earth orbit (LEO)-based measurements are valuable for accurately estimating short-term atmospheric abrupt disturbances triggered by solar flux forcing. To investigate the enhancing status of solar cycle 25 between August 29 and December 31, 2020, we processed Gravity Recovery and Climate Experiment Follow-on (GRACE-FO) 10-s accelerometer-derived thermospheric mass density (TMD) measurements normalized at 500 km altitude by the NRLMSIS-2.0 empirical model. These 4-month enhancing disturbance observations suggest a shift from relative quiescence to a much more active solar phase, revealing unexpected dependencies on temporal and spatial characteristics. The results indicated that the dominant driver is solar extreme ultraviolet radiation (EUV) during this ascending phase. Density enhancement was symmetric in both hemispheres around the autumn equinox. After the equinox, the neutral density enhancement intensity in the Southern Hemisphere surpasses that in the Northern Hemisphere. Density maxima occurred from high to low latitudes, accompanied by a 2–3 h time lag. The Wygant function was applied to model the response to solar wind geomagnetic field changes and quantify the impact of geomagnetic activities on upper atmosphere density, verifying the time lag of density disturbances. All these findings could potentially improve our understanding of the solar cycle and LEO orbital drag.

Description: Uraninite [UO2] is an increasingly recognized accessory mineral for geochronological studies of the mid to upper crust. Similar to what is seen for zircon and monazite, the U-Pb system of uraninite can become reset under relatively low temperatures in certain domains via the action of fluids through the process of coupled dissolution-reprecipitation. Whether or not the uraninite geochronometer is reset will be dependent on the chemistry of the fluid it interacts with as opposed to being purely dependent on P-T. This makes uraninite a mineral of interest for the dating of low- to mid-temperature, fluid-controlled geological processes. In order to better understand which factors cause the recrystallization and/or metasomatic alteration of uraninite, a set of 5 metasomatism experiments have been performed in cold seal autoclaves on a hydrothermal line involving a natural uraninite from Příbram, Czech Republic and a series of Na-, Ca-, OH-, and F-bearing fluids at 600°C and 200 MPa for 21 days. A second set of the same 5 experiments, to which elemental sulfur was added, were subsequently run at 450°C and 200 MPa for 66 days. Generally, little textural alteration of the starting material was observed in any of the experiments, which was independent of the fluid chemistry and temperature, except for an increase in the apparent porosity of the reacted grains. In the second set of experiments galena formed as small grains in four of the runs, indicating that Pb had migrated out from the uraninite into the solution and reacted with the sulfur to form galena. The excessive depletion of Pb in the metasomatized uraninite to negligible amounts in some of these fluids is especially evident if the solutions were NaF + H2O and 2M NaOH. This suggests that interaction of uraninite with F- or high pH Na-bearing fluids can metasomatically reset the uraninite geochronometer at 450 °C and mid to upper crustal pressures and by analogy to even lower temperatures given sufficient time.

Description: Rare earth elements (REE) include the lanthanides (La–Lu), Y, and Sc which are critical elements for the green energy transition. The REE show a decrease in ionic radii with increased atomic numbers, which results in a so-called lanthanide contraction systematically affecting crystal structures and mineral properties. Here we present a compilation of reference Raman spectra of ten REE sesquioxides (A-, B- and C-type), five REE hydroxides, eight xenotime-structured REE phosphate endmembers and two solid solutions, seven monazite-structured REE phosphate endmembers and two solid solutions and seven rhabdophane endmembers with up to five Ce1−xLREEx rhabdophane solid solutions (LREE = La–Gd). Raman mode assignment is based on a detailed literature review summarizing existing analytical work and theoretical calculations and systematic trends observed in this study by analyzing different REE-bearing solids. The wavenumbers of the main REE-O Raman band systematically increase with decreasing ionic radii forming discrete linear trends within isostructural mineral groups, that can be used to estimate the REE-O mode in other solids with known REE-O coordination numbers. Photoluminescence using 266 nm, 532 nm and 633 nm excitation laser wavelengths for REE-bearing oxides, hydroxides, anhydrous and hydrous phosphates is also presented providing a new framework for identifying REE-phases in phosphate-bearing natural mineral deposits.

Description: Iron oxide-apatite (IOA) deposits generally undergo extensive late-stage metasomatic overprinting that modifies the texture and geochemistry of the magnetite, apatite, and other ore-forming minerals. However, the onset of metasomatism in IOA deposits worldwide generally remains poorly constrained. The Heiyingshan IOA deposit is located in the Beishan area of the Central Asian Orogenic Belt, NW China. Fluorapatite from this IOA deposit has undergone extensive fluid-aided alteration resulting in the formation of numerous monazite inclusions via a coupled dissolution-reprecipitation process (CDRP) in altered domains of the fluorapatite. Petrographic studies show that these monazite grains are mainly subhedral to anhedral in shape, have a topotaxial relationship with the parent fluorapatite, and co-exist with magnetite, xenotime, and other mineral inclusions. The crystallization age of this monazite constrains the metasomatic activity following the original IOA mineralization. This study presents EMPA and trace element LA-ICP-MS analytical data for fluorapatite and monazite from the Heiyingshan IOA deposit. It also includes in-situ U–Pb isotopic data for unaltered, altered domains of fluorapatite, and monazite. U–Pb data obtained for unaltered fluorapatite, altered fluorapatite and monazite inclusions within altered fluorapatite grains have similar U–Pb ages (∼325 Ma), consistent with the regional Carboniferous volcanic rocks. Furthermore, this study also presents Sr–Nd isotopic data for fluorapatite, the unaltered and altered fluorapatite domains in the Heiyingshan IOA deposit have similar initial Sr isotope ratios, 143Nd/144Nd ratios, and εNd(t) values, with the Sr isotope ratios consistent with the regional Carboniferous volcanic rocks. This implies that magmatism, mineralization, and metasomatism in the Heiyingshan IOA deposit were essentially contemporaneous, the metasomatizing fluids were derived from the evolution of mineralizing fluids. A comparison of the results from Heiyingshan with the available geochronological data from other IOA ore deposits reveals a general consistency in the timing between the mineralization and the metasomatism of IOA deposits due to evolved later-stage fluids. In contrast, dating of monazite metasomatically derived from fluorapatite in some older (e.g., Kiirunavaara) IOA deposits indicates that these IOA deposits appear to have undergone additional episodes of metasomatic alteration during various later geological events after the original IOA mineralization. We suggest that the formation age of monazite inclusions within apatite in IOA deposits should be consistent with the timing of the IOA mineralization, as the original formation and subsequent fluid-aided modification of IOA deposit occur in the same magmatic-hydrothermal event/system. However, some IOA deposits may have undergone additional late episodes of metasomatic events.

Description: In high-precision space geodetic techniques data processing, the mapping function (MF) is a key factor in mapping the radio waves from the zenith direction down to the signal incoming direction. Existing MF products, either site-wise Vienna Mapping Function (VMF1 and VMF3) or grid-wise VMF1 and VMF3, are only available at the Earth surface. For overhead areas, height correction is always required, which is becoming increasingly important with growing airborne aircraft activity. In this contribution, we introduce a novel method aimed at providing a large number of MFs to the user in a simple and efficient manner, while minimizing the loss of precision. The approach effectively represents the vertical profile of the MFs from the Earth's surface up to altitudes of 14 km. In addition, the new model corrects for height in the assessment using the fifth generation of the European Centre for Medium-Range Weather Forecasts ReAnalysis (ERA5) ray tracing calculations for a global 5° × 5° grid with 54 layers in the vertical direction, a total of 8 azimuths in the plane, and 7 elevation angles, for each day in 2021. Specifically, for both polynomial and exponential model of order 2 and 3, the relative residuals are 〈 0.3% for the hydrostatic delay MF coefficient , and 〈 1% for the wet delay MF coefficient . The precision of the new model on the Earth’s surface is evaluated using site-wise VMF1 and VMF3 GNSS (Global Navigation Satellite System) products from Technische Universität Wien. The root mean square error of slant hydrostatic delay and slant wet delay at a 3° elevation angle is approximately 4–5 cm and 2–5 cm, respectively.

Description: Radio signals transmitted by Global Navigation Satellite System (GNSS) satellites experience tropospheric delays. While the hydrostatic part, referred to as zenith hydrostatic delay (ZHD) when mapped to the zenith direction, can be analytically modelled with sufficient accuracy, the wet part, referred to as zenith wet delay (ZWD), is much more difficult to determine and needs to be estimated. Thus, there exist several ZWD models which are used for various applications such as positioning and climate research. In this study, we present a data-driven, global model of the spatial ZWD field, based on the Extreme Gradient Boosting (XGBoost). The model takes the geographical location, the time, and a number of meteorological variables (in particular, specific humidity at several pressure levels) as input, and can predict ZWD anywhere on Earth as long as the input features are available. It was trained on ZWDs at 10718 GNSS stations and tested on ZWDs at 2684 GNSS stations for the year 2019. Across all test stations and all observations, the trained model achieved a mean absolute error of 6.1 mm, respectively, a root mean squared error of 8.1 mm. Comparisons of the XGBoost-based ZWD predictions with independently computed ZWDs and baseline models underline the good performance of the proposed model. Moreover, we analysed regional and monthly models, as well as the seasonal behaviour of the ZWD predictions in different climate zones, and found that the global model exhibits a high predictive skill in all regions and across all months of the year.

Description: Tin (Sn) and tungsten (W) behave incompatibly in reduced magmatic systems and may become enriched in late highly-evolved melts. Nonetheless, Sn and W rarely concentrate in the same deposit. In deposits formed by Sn- and W-bearing granites, this separation may be due to the contrasting behavior of Sn and W during exsolution of a magmatic fluid or the scavenging of Sn by silicate minerals. We illustrate the separation of Sn and W for the world-class Zhuxi W skarn deposit (South China). Although tin orebodies have not yet been identified within the Zhuxi deposit, tiny (commonly 〈 20 μm) cassiterite grains are widespread within the endoskarn and the retrogressed exoskarn. We analyzed the W and Sn contents of the magmatic minerals biotite and ilmenite in ore-forming granites and the prograde anhydrous skarn minerals garnet, pyroxene and vesuvianite. Our data show that (i) magmatic ilmenite (65.5–79.1 ppm Sn; 8.7–14.3 ppm W) and biotite (109–120 ppm Sn; 1.3–6.3 ppm W) from biotite monzogranite strongly enrich Sn relative to W, implying that W partitioned more strongly into the magmatic fluids than Sn, (ii) there is 100 Kt non-recoverable Sn within the Zhuxi deposit in addition to the certified 3.44 Mt WO3 reserves, and (iii) W is mainly hosted in scheelite, whereas Sn is dominantly sequestered in prograde skarn minerals, most importantly garnet (76–4086 ppm Sn, 〈 42 ppm W), pyroxene (3–103 ppm Sn, 〈 1 ppm W), and vesuvianite (43–361 ppm Sn, 〈 2 ppm W). The formation of secondary cassiterite requires the release of silicate-bound Sn by alteration of primary skarn minerals, which depends on the availability of magmatic or metamorphic fluids. Deep-seated granites such as those associated with the Zhuxi skarn deposit, which crystallized at 5 km to 12.6 km depth, do not release or mobilize copious amounts of fluid. Therefore, the Zhuxi deposit, like other deep-seated reduced skarn systems shows little alteration and most Sn remains in silicate minerals and is economically non-recoverable. Thus, reduced, deep-seated W skarn systems are unlikely to have associated Sn orebodies even if significant amounts of Sn are present.

Description: Carbonatites and their comagmatic silicate rocks related deposit provide significant resources of rare earth elements (REEs), niobium (Nb) and other elements such as U, Th, Mo, V, Ba, Sr, etc. However, the genesis of mineralization, especially for REEs and Nb, in carbonatite remains enigmatic. Previous liquid immiscibility experiments have demonstrated that both REEs and Nb are preferentially enriched in the silicate conjugate instead of carbonate melts under anhydrous conditions. Nevertheless, ligands other than carbonate ion appear to be abundant due to ubiquity of apatite, baryte, celestine, fluorite and sodalite in carbonate–silicate magmatic systems. Here, we experimentally investigate the trace element partitioning between natrocarbonate and silicate (nephelinite) melts in systems doped with varying amounts of additional F−, PO43−, Cl−, and SO42− ligands (0, 2, 4 and 6 wt%) to understand and constrain the role of ligands. The experiments were conducted at 850 °C and 0.1 GPa using rapid quench cold-seal pressure vessels (CSPVs). A comparison of experimental partition coefficients in this study reveals that the significant amounts F− and PO43− incorporated in the silicate melts can increase the D values for REE by influencing melt structure (DLaCM/SM = 0.85–7.42). In contrast, irrespective of the amount of added Cl− and SO42−, DCM/SM is not affected significantly by these species and the DREECM/SM values remain always lower than 1 (DLaCM/SM = 0.12–0.40). Notably, the DNbCM/SM values are all 〈1, with only one exception containing 6 wt% F. Besides, in all the investigated systems, Ba, Sr, Mo, V, Cs, Rb and Li preferentially partition into the conjugate carbonate melt. All the high field strength elements (Pb, Th, U, Zr, Hf, Nb, Ta), transition metals (Mn, Co, Cu, Zn) and common network formers (Ga, Ge) essentially partition into the silicate melt.

Description: This work studies variations of ionospheric total electron content (TEC) during four distinct solar eclipse events over the Ethiopia region. Dual‐frequency global positioning system (GPS) data obtained from UNAVCO over Addis Ababa (9.036°N, 38.76°E) and Bahir Dar (11.6°N, 37.34°E) stations are used to examine the ionospheric variability during two annular solar eclipses on 15 January 2010 and 1 September 2016, a partial solar eclipse on 4 January 2011, and a hybrid solar eclipse (the eclipse path starts out as annular but later changes to total) on 3 November 2013. The results show a significant decrease in TEC values during the occurrence of the solar eclipses. Specifically, the TEC values are reduced to

Description: An ideal target for geodetic very long baseline interferometry (VLBI) is a strong and point-like radio source. In reality, most celestial sources used in geodetic VLBI have spatial structure. This is as a major source of error in VLBI Global Observing System (VGOS) and also affects legacy S/X observations. Source structure causes a systematic delay, which can affect the geodetic estimates if not modelled or otherwise accounted for. In this work, we aim to mitigate its impact by extending the stochastic model used in the least-squares fitting of the VLBI group delays. We have developed a weighting scheme to re-weight the observations by parameterizing the source structure component in terms of closure delays and jet orientation relative to the observing baseline. It was implemented in the Vienna VLBI Software. To assess the performance of the extended stochastic model, we analysed the CONT17 legacy sessions and generated suitable reference solutions for comparison. The effects of re-weighting were evaluated with respect to the session fit statistics, source-wise residuals, and geodetic parameters. We find that this relatively simple noise model consistently improves the session fit by about 5% with moderate variation from session to session. The geodetic estimates are not affected to a significant level by this new weighting method. Source-wise we see improved post-fit residuals for 63 out of a total of 91 sources observed.

Description: Radio telescopes with dual linearly polarized feeds regularly participate in Very Long Baseline Interferometry. One example is the VLBI Global Observing System (VGOS), which is employed for high-precision geodesy and astrometry. In order to achieve the maximum signal-to-noise ratio, the visibilities of all four polarization products are combined to Stokes I before fringe-fitting. Our aim is to improve cross-polarization bandpass calibration, which is an essential processing step in this context. Here we investigate the shapes of these station-specific quantities as a function of frequency and time. We observed the extra-galactic source 4C 39.25 for 6 hours with a VGOS network. We correlated the data with the DiFX software and analyzed the visibilities with PolConvert to determine the complex cross-bandpasses with high accuracy. Their frequency-dependent shape is to first order characterized by a group delay between the two orthogonal polarizations, in the order of several hundred picoseconds. We find that this group delay shows systematic variability in the range of a few picoseconds, but can remain stable within this range for several years, as evident from earlier sessions. On top of the linear phase-frequency relationship there are systematic deviations of several tens of degrees, which in addition are subject to smooth temporal evolution. The antenna cross-bandpasses are variable on time scales of ∼1 hr, which defines the frequency of necessary calibrator scans. The source 4C 39.25 is confirmed as an excellent cross-bandpass calibrator. Dedicated surveys are highly encouraged to search for more calibrators of similar quality.

Description: Integer Ambiguity Resolution (IAR) can significantly improve the accuracy of GNSS Precise Orbit Determination (POD). Traditionally, the IAR in POD is achieved at the Double Differenced (DD) level. In this contribution, we develop an Un-Differenced (UD) IAR method for Global Positioning System (GPS)+ BeiDou Navigation Satellite System (BDS) + Galileo navigation satellite system (Galileo)+ Global'naya Navigatsionnaya Sputnikovaya Sistema (GLONASS) quad-system POD by calibrating UD ambiguities in the raw carrier phase and generating the so-called carrier range. Based on this method, we generate the UD ambiguity-fixed orbit and clock products for the Wuhan Innovation Application Center (IAC) of the International GNSS Monitoring and Assessment System (iGMAS). One-year observations in 2020 from 150 stations are employed to investigate performance of orbit and clock products. Notably, the UD Ambiguity Resolution (AR) yields more resolved integer ambiguities than the traditional DD AR, scaling up to 9%, attributable to its avoidance of station baseline formation. Benefiting from the removal of ambiguity parameters, the computational efficiency of parameter estimation undergoes a substantial 70% improvement. Compared with the float solution, the orbit consistencies of UD AR solution achieve the accuracy of 1.9, 5.2, 2.8, 2.1, and 2.7 cm for GPS, BeiDou-2 Navigation Satellite System (BDS-2), BeiDou-3 Navigation Satellite System (BDS-3), Galileo, and GLONASS satellites respectively, reflecting enhancements of 40%, 24%, 54%, 34%, and 42%. Moreover, the standard deviations of Satellite Laser Ranging (SLR) residuals are spanning 2.5–3.5 cm, underscoring a comparable accuracy to the DD AR solution, with discrepancies below 5%. A notable advantage of UD AR lies in its capability to produce the Integer Recovered Clock (IRC), facilitating Precise Point Positioning (PPP) AR without requiring additional Uncalibrated Phase Delay (UPD) products. To assess the performance of quad-system kinematic PPP based on IRC, a network comprising 120 stations is utilized. In comparison to the float solution, the IRC-based PPP AR accelerates convergence time by 31% and enhance positioning accuracy in the east component by 54%.

Description: The global navigation satellite system (GNSS)-interferometric reflectometry technique has been applied to retrieve snow depth, which has a high potential for application. The GNSS reflectometry classical algorithm retrieves the snow depth by extracting the frequency of the multipath signal and substituting it into an empirical formula. However, the retrieval errors of high and low snow depths are large due to the influence of factors such as surface vegetation and terrain environment. In this paper, we propose a snow depth retrieval algorithm based on a particle swarm optimized long short-term memory (PSO-LSTM) neural network. The algorithm extracted three characteristic parameters (frequency, amplitude, and phase) from the signal-to-noise ratio (SNR) data as inputs, and optimized the LSTM hyperparameters by the PSO algorithm to improve the retrieval accuracy for low snow depths and snow depths close to the antenna. The snow depth retrieval results of global positioning system L1 band SNR data collected from the P351 station in 2022 and AB33 station in 2017 were evaluated in this paper. The snow depth retrieval results of the PSO-LSTM algorithm for P351 station were in high agreement with the snow depth data provided by the snowpack telemetry network; the coefficient of determination () reached 0.986, and the root mean square error (RMSE) and mean absolute error (MAE) were 7.30 cm and 4.94 cm, respectively. Compared with the classical algorithm, the PSO-LSTM algorithm decreased the RMSE and MAE by 53.0% and 30.4% for the retrieval results of snow depths below 15 cm at the P351 station, and by 76.8% and 84.4% for the retrieval results of snow depths above 117 cm from the 1st day to the 137th day, respectively. Similarly, the RMSE, MAE, and for the 2017 retrieval results at AB33 station were 5.90 cm, 4.25 cm, and 0.965, respectively. Compared with the classical algorithm, the PSO-LSTM algorithm decreased the RMSE and MAE by 47.9% and 33.0% for the retrieval results of snow depths below 15 cm at the AB33 station, and by 75.4% and 82.3% for the retrieval results of snow depths above 56 cm from the day 46 to day 120. In addition, the snow depth retrieval algorithm was proposed in this paper does not require antenna height and empirical formulas to realize snow depth retrieval, and at the same time, the algorithm effectively improved the retrieval accuracy for both high and low snow depths with strong robustness.

Description: The Gram-positive, rod-shaped endophytic bacterium Cellulomonas sp. strain ATA003 was isolated from the endemic cactus Maihueniopsis domeykoensis seeds collected in the Coastal Atacama Desert, Chile. Here, we present a circular genome with a size of 4,084,881 bp and a GC content of 73.8% obtained by Nanopore sequencing.

Description: Geomechanics play an important role in any underground activity, such as carbon dioxide (CO2) and hydrogen (H2) geo-storage, owing to the considerable hazards linked to the injection and withdrawal of fluids into and from the subsurface. In order to quantify these risks, knowledge of full stress tensor is required. Yet, most of our stress information in the Australian target basins for geo-storage is limited to the stress orientations, while stress magnitude data is sparse. 3D geomechanical modelling has proved to be an invaluable tool for prediction of full stress tensor. Nevertheless, a model requires some stress magnitude data in order to tune the model to be representative of real stress state. In situations where stress magnitude data is lacking, this means that the model is susceptible to significant uncertainties. Herein, we present a novel strategy for stress modelling, which involves the utilisation of indirect data such as borehole breakouts, drilling-induced fractures, seismic activity records, and formation integrity tests to calibrate a 3D geomechanical model. We employ the northern Bowen Basin, an onshore basin in Queensland, Australia, as a case study for a comprehensive 3D geomechanical modelling approach. We assess all the indirect information in the model’s volume to narrow down the model predictions and find the most reliable stress state. This innovative approach is an important step forward in stress modelling of Australian basins, where lack of stress magnitudes is a great challenge for geomechanical assessment of geo-storage.

Description: Chromium (Cr) leached from iron (Fe) (oxyhydr)oxide-rich tropical laterites can substantially impact downstream groundwater, ecosystems, and human health. However, its partitioning into mineral hosts, its binding, oxidation state, and potential release are poorly defined. This is in part due to the current lack of well-designed and validated Cr-specific sequential extraction procedures (SEPs) for laterites. To fill this gap, we have (i) first optimized a Cr SEP for Fe (oxyhydr)oxide-rich laterites using synthetic and natural Cr-bearing minerals and laterite references, (ii) used a complementary suite of techniques and critically evaluated existing non-laterite and non-Cr-optimized SEPs, compared to our optimized SEP, and (iii) confirmed the efficiency of our new SEP through analyses of laterites from the Philippines. Our results show that other SEPs inadequately leach Cr host phases and underestimate the Cr fractions. Our SEP recovered up to seven times higher Cr contents because it (a) more efficiently dissolves metal-substituted Fe phases, (b) quantitatively extracts adsorbed Cr, and (c) prevents overestimation of organic Cr in laterites. With this new SEP, we can estimate the mineral-specific Cr fractionation in Fe-rich tropical soils more quantitatively and thus improve our knowledge of the potential environmental impacts of Cr from lateritic areas.

Description: A large dataset of detrital zircon U–Pb ages (N = 5940) of Aptian-Albian strata of the Colombian-Ecuadorian retroarc region suggests that these rocks were sourced from Proterozoic cratonic and Permian–Triassic to Cretaceous rocks of the Andean proto-Cordilleras. The nonconformity between Aptian-Albian strata and Proterozoic rocks of the Garzón Massif indicates the existence of positive relief in this region. Topographic highs could have caused local basin compartmentalization and the prevalence of a localized provenance of coeval strata. Areas of positive relief were seemingly exposed to intense chemical weathering as suggested by the high Chemical Index of Alteration values, and low Rb/Sr and SiO2/Al2O3 ratios of analysed strata. Our results highlight the value of provenance analysis to study the ancient topography of retroarc systems and open the avenue for further research on the role of extensional tectonics in the topographic and basin evolution in the Andes.

Description: Using the electron density (Ne) observations from the Defense Meteorological Satellite Program, and Constellation Observing Systems for Meteorology, Ionosphere, and Climate mission and simulations from the Thermosphere Ionosphere Electrodynamic General Circulation Model. we investigate the dynamic evolution of the polar tongue of ionization (TOI) from double to single structures at different altitudes during a geomagnetic storm. The modeled Ne depicted that double and single TOIs occurred at altitudes above 300 km, respectively. During the northward turning of interplanetary magnetic field (IMF) Bz, the afternoon TOI disappeared and the morning TOI was reduced. The plasma transport due to neutral winds and ambipolar diffusion facilitated (prevented) the depletion of plasma density of the morning TOI at 300 (500) km, with a relative contribution of 42.8% and 28.6% (−15.4% and −76.9%), respectively. Downward E × B drifts led to an enhancement/reduction of plasma density in the storm enhanced density region in the lower/upper ionosphere. During the duskward turning of IMF By, the morning TOI could be mostly attributed to the anti-sunward plasma drifts (75.8% at 300 km, 100% at 500 km), with a relatively stronger role of the zonal component than that of meridional E × B drifts. The upward E × B drifts were important/ignorable in the upper/lower ionosphere. Both the neutral winds and ambipolar diffusion resulted in an accumulation of plasma density of the morning TOI at 300 km indirectly (24.2%), however, their roles were minor at 500 km

Description: Natural alteration of zircon takes place in melts or fluids either via dissolution coupled with overgrowth or via a coupled dissolution-reprecipitation process. The latter results in the zircon being partially or totally replaced by new, compositionally re-equilibrated zircon or a new mineral phase or both. In this study, fragments (50–300 μm) from a large, inclusion-free, clear, 520–530 Ma euhedral zircon with light radiation damage from a nepheline syenite pegmatite, Seiland Igneous Province, northern Norway, were experimentally reacted in 20 mg batches with 5 mg of ThO2 + ThSiO2 + SiO2 and a series of alkali-bearing fluids in sealed Pt capsules at 900 °C and 1000 MPa for 6–11 days in the piston cylinder press using a CaF2 setup with a cylindrical graphite oven. ThO2 + ThSiO2 + SiO2 was present at the end of the experiment. In experiments involving H2O, H2O + NaCl, H2O + KCl, and 2 N KOH, no reaction textures formed other than a slight dissolution of the zircon grain fragments. Experiments involving 2 N NaOH, Na2Si2O5 + H2O, and NaF + H2O resulted in zircon reaction textures with varying degrees of intensity, which took the form of partial replacement by compositionally modified zircon via a coupled dissolution-reprecipitation process. In the NaF + H2O experiment some overgrowth also occurred. Altered zircon is separated by sharp compositional boundaries from unaltered zircon. Secondary ion mass spectrometry (SIMS) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) analysis indicates that, relative to the unaltered zircon, the altered zircon is strongly enriched in Th, and heavily to moderately depleted in U and (Y + REE). In all the experiments, 206Pb (3–5 ppm in unaltered zircon) is depleted in the altered zircon to below the SIMS detection limit and to at or below the LA-ICP-MS detection limit. Hafnium and Ti concentrations in the altered zircon retained the same approximate value (within error) as the original zircon. The results from these experiments demonstrate that zircon can be compositionally modified by alkali-bearing and alkali-F-bearing fluids via a coupled dissolution-reprecipitation process. Near to total loss of radiogenic Pb via such processes under high-grade conditions resets the internal zircon geochronometer. Although the end result is the same as with zircon overgrowth, i.e. the production of new generation zircon at the time of a metamorphic/metasomatic event, such replacement processes can explain incomplete isotopic ‘resetting’; inclusion production through unmixing of solid solutions in metastable zircon compositions; and ‘ghost’ textures that preserve initial growth features but with isotopic disturbance. Diagnostic replacement features produced in experiments, such as interface geometries between altered and unaltered zircon, provide markers of the mechanism and aid in zircon interpretation. A major implication from this study is that if zircon with low radiation damage can be metasomatically altered under high-grade conditions, this would have important consequences with respect to zircons presumed role as an impregnable container for mineral inclusions. Namely the mineral inclusions contained within zircon could also be altered, reset as a geochronometer, or even replaced by another mineral.

Description: To investigate the long-term stability of deep rocks, a three-dimensional (3D) time-dependent model that accounts for excavation-induced damage and complex stress state is developed. This model com- prises three main components: a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state, a quantitative relationship between critical irreversible deformation and complex stress state, and evolution characteristics of strength parameters. The proposed model is implemented in a self-developed numerical code, i.e. CASRock. The reliability of the model is validated through experiments. It is indicated that the time-dependent fracturing potential index (xTFPI) at a given time during the attenuation creep stage shows a negative correlation with the extent of excavation- induced damage. The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress, thereby highlighting the 3D stress-dependent characteristic of the model. Finally, the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed. Ó 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).

Description: The Lindero deposit is located in the Puna plateau, northwest Argentina, at the southern end of the Central Volcanic Zone of the Central Andes. The high-K calc-alkaline dioritic composition of the subvolcanic intrusions, the shallow emplacement depth (〈 1.5 km), and the gold-rich and copper-depleted mineralization style suggest that the Lindero deposit is a porphyry gold deposit. Porphyry gold deposits are scarce worldwide and the factors controlling their formation are still poorly known. Here we present a detailed study of fluid inclusions in order to characterize the mineralizing fluids that precipitated the Au mineralization at Lindero. Different types of fluid inclusions in quartz veins (A-type and banded quartz), which are associated with the K-silicate alteration, were analyzed using Raman spectroscopy, microthermometry, and LA-ICP-MS (laser ablation inductively coupled plasma mass spectrometry). Four inclusion types can be recognized in quartz veins: (i) Salt melt inclusions, which are characterized by a dense packing of daughter minerals (mainly Fe-chloride, sylvite, halite, anhydrite, and hematite), by a distorted vapor bubble, and by the lack of liquid phase; (ii) Halite-bearing inclusions which contain liquid, vapor, and halite; (iii) Two-phase aqueous inclusions that contain liquid and vapor; (iv) Vapor-rich inclusions containing only vapor. The inclusion types are related to different stages of hydrothermal evolution. Stage 1 is the main mineralization stage, characterized by vapor-rich inclusions coexisting with salt melt inclusions. Salt melt inclusions commonly show total homogenization temperature (ThL) 〉 1000 °C. This Na-K-Fe-Cl-rich highly saline brine (~ 90 wt% NaCl eq.) was of magmatic origin and responsible for the Au mineralization. Two later stages involving cooler fluids (ThL 〈 300 °C) and gradually lower salinities (from 36.1 to 0.2 wt% NaCl eq.) trapped by halite-bearing and two-phase aqueous inclusions during stages 2 and 3, respectively, correspond to a late magmatic-hydrothermal system, that is probably related to a deep supercritical fluid exsolution. Salt melt inclusions represent the most likely parental fluid of K-silicate alteration and associated Au mineralization at Lindero. This uncommon type of fluid must have played an important role in Au transport and precipitation in shallow porphyry gold deposits.

Description: Osmium isotope and highly siderophile element (HSE: Os, Ir, Ru, Pt, Pd, Re) abundance data are reported for picrites and basalts from the ∼132 Ma Etendeka large igneous province (LIP) and the ∼60 Ma North Atlantic Igneous Province (NAIP). Picrite dykes of the Etendeka LIP have HSE abundances and 187Os/188Os (0.1276 to 0.1323; γOsi = -0.5 to +3.1) consistent with derivation from high-degree partial melting (〉20 %) of a peridotite source with chondritic to modestly supra-chondritic long-term Re/Os. High-3He/4He NAIP picrites from West Greenland represent large-degree partial melts with similarly elevated HSE abundances and 187Os/188Os (0.1273 to 0.1332; γOsi = -0.2 to +3.9). Broadly chondritic Os isotope ratios have also been reported for the ∼132 Ma Paraná LIP and the ∼201 Ma Central Atlantic Magmatic Province (CAMP). Consequently, LIP associated with Atlantic Ocean opening derive, at least in part, from partial melting of peridotite mantle distinct from the depleted mantle associated with mid-ocean ridge basalt volcanism. Modern locations with high-3He/4He (〉25RA) include ocean island basalts (OIB) from Ofu (Samoa), Loihi (Hawaii) and Fernandina (Galapagos) in the Pacific Ocean, and from Iceland, which is considered the modern manifestation of NAIP magmatism. Unlike Etendeka and NAIP picrites, these modern OIB have Sr-Nd-Pb-Os isotopes consistent with contributions of recycled oceanic or continental crust. The lower degree of partial melting responsible for modern high-3He/4He OIB gives higher proportions of fusible recycled crustal components to the magmas, with radiogenic 187Os/188Os and low-3He/4He. The high-3He/4He, incompatible trace element-depleted mantle component in both LIP and OIB therefore also has long-term chondritic Re/Os, which is consistent with an early-formed reservoir that experienced late accretion. Atlantic LIP (CAMP; Paraná-Etendeka; NAIP) provide geochemical evidence for a prominent role for mantle plume contributions during continental break-up and formation of the Atlantic Ocean, a feature hitherto unrecognized in other ocean basin-forming events.

Description: Paleomagnetic records of middle Neoproterozoic (820 to 780 Ma) rocks display high amplitude directional variations that lead to large discrepancies in paleogeographic reconstructions. Hypotheses to explain these data include rapid true polar wander (TPW), a geomagnetic field geometry that deviates from a predominantly axial dipole field, a hyper-reversing field (〉10 reversals/Ma), and/or undiagnosed remagnetization. To test these hypotheses, we collected 1,057 oriented cores over a 85 m stratigraphic succession in the Laoshanya Formation (Yangjiaping, Hunan, China). High precision U-Pb dating of two intercalated tuff layers constrain the age of the sediments between 809 and 804 Ma. Thermal demagnetization isolates three magnetization components residing in hematite which are not time-progressive but conflated throughout the section. All samples possess a north and downward directed component in geographic coordinates at temperatures up to 660°C that is ascribed to a Cretaceous overprint. Two components isolated above 660°C reveal distinct directional clusters: one is interpreted as a depositional remanence, while the other appears to be the result of a mid-Paleozoic (460 to 420 Ma) remagnetization, which is likely widespread throughout South China. The high-temperature directions are subtly dependent on lithology; microscopic and rock magnetic analyses identify multiple generations of hematite that vary in concentration and distinguish the magnetization components. A comparison with other middle Neoproterozoic paleomagnetic studies in the region indicates that the sudden changes in paleomagnetic directions, used elsewhere to support the rapid TPW hypothesis (ca. 805 Ma), are better explained by mixtures of primary and remagnetized components, and/or vertical axis rotations.

Description: Volcanoes produce a variety of seismic signals and, therefore, continuous seismograms provide crucial information for monitoring the state of a volcano. According to their source mechanism and signal properties, seismo‐volcanic signals can be categorized into distinct classes, which works particularly well for short transients. Applying classification approaches to long‐duration continuous signals containing volcanic tremors, characterized by varying signal characteristics, proves challenging due to the complex nature of these signals. That makes it difficult to attribute them to a single volcanic process and questions the feasibility of classification. In the present study, we consider the whole seismic time series as valuable information about the plumbing system (the combination of plumbing structure and activity distribution). The considered data are year‐long seismograms recorded at individual stations near the Klyuchevskoy Volcanic Group (Kamchatka, Russia). With a scattering network and a Uniform Manifold Approximation and Projection (UMAP), we transform the continuous data into a two‐dimensional representation (a seismogram atlas), which helps us to identify sudden and continuous changes in the signal properties. We observe an ever‐changing seismic wavefield that we relate to a continuously evolving plumbing system. Through additional data, we can relate signal variations to various state changes of the volcano including transitions from deep to shallow activity, deep reactivation, weak signals during quiet times, and eruptive activity. The atlases serve as a visual tool for analyzing extensive seismic time series, allowing us to associate specific atlas areas, indicative of similar signal characteristics, with distinct volcanic activities and variations in the volcanic plumbing system.

Description: Saline water is a common fluid on the Earth‘s surface and in ice planets. Potassium chloride (KCl) is a common salt and is expected to be a ubiquitous solute in salt water in the Universe; however, few studies investigated the behavior of KCl-H2O system at high pressures and temperatures. In this study, powder and single-crystal X-ray diffraction (SC-XRD), Raman and Brillouin scattering combined with diamond anvil cells were used to investigate the phase relation in the KCl-H2O system for different KCl concentrations at 0–4 GPa and 298–405 K. The results of powder X-ray diffraction and Raman scattering demonstrate that a novel KCl hydrate is formed when KCl aqueous solutions transform to solid ice-VI and ice-VII at high pressure. Simultaneously, the single-crystal of KCl hydrate is synthesized from a supersaturated KCl solution at 298 K and 1.8 GPa. The structure is solved by SC-XRD, indicating a KCl monohydrate with the P21/n space group is formed. We have verified the phase stability of KCl monohydrate by using Raman spectroscopy and density functional theory. Our results indicate that KCl monohydrate is a stable phase under pressure and temperature conditions between 1.6 and 2.4 GPa and 298–359 K. By considering the thermal profile and composition of icy moons, we hypothesize that the formation and decomposition of KCl monohydrate might induce mantle convection in these moons.

Description: An experimental platform for dynamic diamond anvil cell (dDAC) research has been developed at the High Energy Density (HED) Instrument at the European X-ray Free Electron Laser (European XFEL). Advantage was taken of the high repetition rate of the European XFEL (up to 4.5 MHz) to collect pulse-resolved MHz X-ray diffraction data from samples as they are dynamically compressed at intermediate strain rates (≤103 s−1), where up to 352 diffraction images can be collected from a single pulse train. The set-up employs piezo-driven dDACs capable of compressing samples in ≥340 µs, compatible with the maximum length of the pulse train (550 µs). Results from rapid compression experiments on a wide range of sample systems with different X-ray scattering powers are presented. A maximum compression rate of 87 TPa s−1 was observed during the fast compression of Au, while a strain rate of ∼1100 s−1 was achieved during the rapid compression of N2 at 23 TPa s−1.

Description: Experiments accessing extreme conditions at x-ray free electron lasers (XFELs) involve rapidly evolving conditions of temperature. Here, we report time-resolved, direct measurements of temperature using spectral streaked optical pyrometry of x-ray and optical laser-heated states at the High Energy Density instrument of the European XFEL. This collection of typical experiments, coupled with numerical models, outlines the reliability, precision, and meaning of time dependent temperature measurements using optical emission at XFEL sources. Dynamic temperatures above 1500 K are measured continuously from spectrally- and temporally-resolved thermal emission at 450–850 nm, with time resolution down to 10–100 ns for 1–200 μs streak camera windows, using single shot and integrated modes. Targets include zero-pressure foils free-standing in air and in vacuo, and high-pressure samples compressed in diamond anvil cell multi-layer targets. Radiation sources used are 20-fs hard x-ray laser pulses at 17.8 keV, in single pulses or 2.26 MHz pulse trains of up to 30 pulses, and 250-ns infrared laser single pulses. A range of further possibilities for optical measurements of visible light in x-ray laser experiments using streak optical spectroscopy are also explored, including for the study of x-ray induced optical fluorescence, which often appears as background in thermal radiation measurements. We establish several scenarios where combined emissions from multiple sources are observed and discuss their interpretation. Challenges posed by using x-ray lasers as non-invasive probes of the sample state are addressed.

Description: The ambient-temperature compressibility and room-pressure thermal expansion of two Mg3(PO4)2 polymorphs (farringtonite=Mg3(PO4)2-I, with 5- and 6-fold coordinated Mg, and chopinite=“Mgsarcopside”=[6]Mg3(PO4)2-II), three Mg2PO4OH polymorphs (althausite, hydroxylwagnerite and ɛ- Mg2PO4OH, all with [5]Mg and [6]Mg) and phosphoellenbergerite ([6]Mg) were measured on synthetic powders using a synchrotron-based multi-anvil apparatus to 5.5 GPa and a laboratory high-temperature diffractometer, with whole-pattern fitting procedures. Bulk moduli range from 64.5 GPa for althausite to 88.4 GPa for hydroxylwagnerite, the high-pressure Mg2PO4OH polymorph. Chopinite, based on an olivine structure with ordered octahedral vacancies (K0=81.6 GPa), and phosphoellenbergerite, composed of chains of face-sharing octahedra (K0=86.4 GPa), are distinctly more compressible than their homeotypical silicate (127 and 133 GPa, respectively). The compressibility anisotropy is the highest for chopinite and the lowest for phosphoellenbergerite. First-order parameters of quadratic thermal expansions range from v1=2.19x10-5K-1 for ɛ-Mg2PO4OH to v1=3.58x10-5K-1 for althausite. Phosphates have higher thermal-expansion coefficients than the homeotypical silicates. Thermal anisotropy is the highest for farringtonite and the lowest for hydroxylwagnerite and chopinite. These results set the stage for a thermodynamic handling of phase-equilibrium data obtained up to 3 GPa and 1000°C in the MgO–P2O5–H2O and MgO–Al2O3–P2O5–H2O systems.

Description: Unicellular eukaryotic plankton communities (protists) are the major basis of the marine food web. The spring bloom is especially important, because of its high biomass. However, it is poorly described how the protist community composition in Arctic surface waters develops from winter to spring. We show that mixotrophic and parasitic organisms are prominent in the dark winter period. The transition period toward the spring bloom event was characterized by a high relative abundance of mixotrophic dinoflagellates, while centric diatoms and the haptophyte Phaeocystis pouchetii dominated the successive phototrophic spring bloom event during the study. The data shows a continuous community shift from winter to spring, and not just a dormant spring community waiting for the right environmental conditions. The spring bloom initiation commenced while sea ice was still scattering and absorbing the sunlight, inhibiting its penetration into the water column. The initial increase in fluorescence was detected relatively deep in the water column at ~55 m depth at the halocline, at which the photosynthetic cells accumulated, while a thick layer of snow and sea ice was still obstructing sunlight penetration of the surface water. This suggests that water column stratification and a complex interplay of abiotic factors eventually promote the spring bloom initiation.

Description: The advancement of the Global Geodetic Observing System (GGOS) has enabled monitoring of mass transport and solid-Earth deformation processes with unprecedented accuracy. Coseismic deformation is modelled as an elastic response of the solid Earth to an internal dislocation. Self-gravitating spherical Earth models can be employed in modelling regional to global scale deformations. Recent seismic tomography and high-pressure/high-temperature experiments have revealed finer-scale lateral heterogeneities in the elasticity and density structures within the Earth, which motivates us to quantify the effects of such finer structures on coseismic deformation. To achieve this, fully numerical approaches including the Finite Element Method (FEM) have often been used. In our previous study, we presented a spectral FEM, combined with an iterative perturbation method, to consider lateral heterogeneities in the bulk and shear moduli for surface loading. The distinct feature of this approach is that the deformation of the entire sphere is modelled in the spectral domain with finite elements dependent only on the radial coordinate. By this, self-gravitation can be treated without special treatments employed when using an ordinary FEM. In this study, we extend the formulation so that it can deal with lateral heterogeneities in density in the case of coseismic deformation. We apply this approach to a longer-wavelength vertical deformation due to a large earthquake. The result shows that the deformation for a laterally heterogeneous density distribution is suppressed mainly where the density is larger, which is consistent with the fact that self-gravitation reduces longer-wavelength deformations for 1-D models. The effect on the vertical displacement is relatively small, but the effect on the gravity change could amount to the same order of magnitude of a given heterogeneity if the horizontal scale of the heterogeneity is large enough.

Description: Strain energy from tectonic loading can be partly released through aseismic creep. Earthquake repeaters, repeatedly activated brittle fault patches surrounded by creep, indicate steady-state creep that affects the amount of seismic energy available for the next large earthquake along a plate contact. The offshore Main Marmara Fault (MMF) of the North Anatolian Fault Zone represents a seismic gap capable of generating a M 〉 7 earthquake in direct vicinity to the mega-city Istanbul. Based on a newly compiled seismicity catalog, we identify repeating earthquakes to resolve the spatial creep variability along the MMF during a 15-year period. We observe a maximum of seismic repeaters indicating creep along the central and western MMF segments tapering off toward the locked onshore Ganos fault in the west, and the locked offshore Princes Islands segment immediately south of Istanbul in the east. This indicates a high degree of spatial creep variability along the Istanbul-Marmara seismic gap.

Description: The Bakreswar geothermal province represents a medium enthalpy geothermal system with its Bakreswar and Tantloie hot springs. It lies within the Chotanagpur Granite Gneissic Complex in the eastern part of the Indian Peninsula. The province has a high heat flow and a high geothermal gradient of 90°C/km. Magnetotelluric data from 95 sites in a frequency range of 10 kHz–10 Hz were acquired over the Bakreswar geothermal province to obtain an electrical conductivity model and map the geothermal reservoir with its fluid pathways and related geological structures. Subsurface conductivity models obtained from three-dimensional inversions of the Magnetotelluric data exhibit several prominent anomalies, which are supplemented by gravity results. The conductivity model maps three features which act as a conduit (a) a northwest–southeast trending feature, (b) an east–west trending feature to the south of the northwest–southeast trending feature (which lies 1 km north of the Oil and Natural Gas Corporation fault marked by previous studies) and (c) shallow conducting features close to Bakreswar hot spring. The northwest–southeast trending feature coincides with the boundary of the high-density intrusive block. This northwest–southeast trending feature provides the pathway for the meteoric water to reach a maximum depth of 2.7 km, where it gets heated by interacting with deep-seated structures and then it rises towards the surface. The radiogenic process occurring within the granites of Chotanagpur Granite Gneissic Complex provides the heat responsible for heating the meteoric water. The northwest–southeast and east–west trending features are responsible for the transport of meteoric water to deeper depths and then towards the shallow regions of the Earth. The near surface features close to the Bakreswar hot spring are responsible for carrying the water further towards the hot spring. The resistivity of these structures plotted as a function of salinity and temperatures for saline crustal fluids suggests the involvement of meteoric water. Further, applying Archie's law to this resistivity suggests that the conduit path has a porosity greater than 10%. This study successfully maps the anomalous structures which might foster the migration of geothermal fluid in Bakreswar geothermal province.

Description: The twenty-eighth Conference of the Parties (COP28) to the United Nations Framework Convention on Climate Change (UNFCCC) in Dubai among many other items concluded the first Global Stocktake (GST) under the Paris Agreement. This article discusses the conference's outcomes in the areas of mitigation, loss and damage, adaptation, climate finance, and cooperation under Article 6 of the Paris Agreement. The conference arguably made history by for the first time ever recognising the need to "transition away" from fossil fuels, adopting specific targets for the scale-up of renewable energy and energy efficiency, and by operationalising a fund to support developing countries in dealing with loss and damage caused by climate change. However, the legal language in the call for an energy transition is relatively non-committal and the conference failed to underpin the new global objectives with adequate resources. Actual implementation of the Dubai outcomes will therefore to a large extent depend on whether COP29 in Baku, Azerbaijan - already billed as "finance COP" - will be able to cut the Gordian knot of finance.

Description: North Rhine-Westphalia (NRW) is the industrial center of Germany and one of the most important industrial locations in Europe. It is a key location for the energy-intensive basic materials industry like the production of steel and non-ferrous metals, (petro)chemicals, cement and lime, bricks, glass and ceramics, and paper. Around 20 % of NRW's total greenhouse emissions derive from industrial processes. By 2045, industry must achieve climate-neutrality, which requires a massive transformation effort. Technologically, this needs large-scale utilization of green hydrogen, carbon management, consequent circular economy, and climate-neutral production of process heat. Furthermore, various adjustments to the policy framework are essential.

Description: Over 175 million Nigerians rely on the use of traditional biomass for cooking, and it is estimated that more than 128,000 people died in Nigeria in 2019 from household air pollution related to these fuels. There is currently a gap in the study of possible pathways to meet Nigeria's goals in clean cooking and in understanding the health and climate impacts that different pathways can bring about. We explore clean cooking access scenarios for Nigeria until 2060 under a business-as-usual scenario, a moderate climate mitigation scenario, and an ambitious transformative scenario. We carry out a disaggregation at the state level for the period up to 2030 to better guide shorter-term policy development. Our analysis shows that under an ambitious scenario where 85 million households achieve access to clean cooking by 2060, annual premature deaths due to exposure to household air pollution would decrease by 7 % compared to 2018 levels. A baseline scenario, on the other hand, sees a dramatic 77 % increase, resulting in 209,000 people dying prematurely, of which 94,000 children under 5. Furthermore, we find that woodfuel removals from forestland would lead to a tripling of carbon dioxide emissions from land use change, reaching 602 Mt CO2 by 2060. Our findings stress the vital importance of a clean cooking transition in Nigeria and underline the urgent need for immediate acceleration in national efforts regarding access to clean cooking for all.

Description: The goal of this dissertation is to facilitate the assessment of impacts from sustainable measures and projects with an emphasis on impact reporting for Green, Social or Sustainability Bonds in the Sustainable Finance market. It does so by providing analysts with the means to develop, depict, formulate, and assess a causal hypothesis between an intervention and its subsequent effects in an impact-chain, represented by desired environmental (E), social (S) or governance (G) changes. This is achieved by developing a methodology for so-called ESG Logic Models or ESG-LM, that combine heuristic Theories-of-Change with propositional logic and Bayesian Reasoning. Three research questions are investigated and responded to. Research Question 1 asks how such Theories-of-Change can be developed for any type of ESG-related issue and how the different process steps in a causal chain can be classified, hierarchised, and prioritised regarding their efficacy towards overarching sustainability goals and their plausibility. Research Question 2 studies (a) the means by which the analyst or any other interested third party might be warranted in believing the causal claims from an ESG-LM, and (b) how an ESG-LM can be improved if this credence is low. Research Question 3 then looks at the reporting of impacts themselves regarding indicator selection, indicator assessment and indicator quantification as well as the provision of information on the contributions and attributions by different actors. The dissertation draws on a variety of theories and adapts existing methods to achieve that. It operationalises concepts from empirical Sustainable Finance research and already existing impact assessment methodologies. It adapts scholarly and practitioner approaches for theory-based evaluation and applies a qualitative social science perspective towards theory-building and evaluation, while some of the assessment tools in the dissertation are grounded in Logic, Set Theory and Bayesian Epistemology. Examples for such tools include rules for the Attribution by actors, heuristics for the abduction of plausible outcome pathways, or a four-stage Argument and Decision-Tree to assess the credibility of ESG-LM claims (based on Bayes Theorem). My assessment of the entire methodology is positive overall, as it provides solutions to each of the three research areas. Limitations of the approach, and thus opportunities for further research, are the additional expertise and time required by analysts compared to the existing, and somewhat more pragmatic, solutions in the current market. However, this is outweighed in my opinion by the ability of the framework to strongly mitigate impact washing by actors in the financial markets as well as biases by analysts. Its overall methodology also provides opportunities for new research angles in the area of sustainability indicators and assessments.

Description: Carrying out laboratory experiments is usually a time-consuming process. In addition, the options for varying parameter studies are limited and adjustments to the design of the measuring equipment are often not possible at all. In order to circumvent these limitations, we supplement our laboratory experiments with virtual experiments as best as possible. For this purpose, we have expanded our finite element library FEMALY [1] to include the so-called complete electrode model [2], which allows us to simulate electrodes of any shape for DC and IP applications and also provides us with explicit mathematical expressions for calculating sensitivities [3]. As a first case study, we consider an IP measurement on a measuring cylinder with embedded ring electrodes to virtually reproduce the time-varying change of the apparent resistivity for laboratory tracer experiments (Figure 1). We present the real and imaginary part of the sensitivity distribution of the underlying measurement configuration that confirms our initial assumption that the actual electrode surface shape has a relatively small influence on the observed measurement quantities.

Description: The Atacama Desert along the Chilean Coastal Cordillera is a unique landscape to understand the Earth's evolution in hyper-arid and arid environments. The Paranal clay pan has studied by the CRC 1211 project to recover a continuous climate record for paleoclimate research. The goal is to provide the sedimentary architecture and bedrock topography of the Paranal site by interpreting multidimensional inversion of loop source transient electromagnetic (TEM) data. A total of 133 TEM soundings were carried out using a central loop configuration, with a transmitter loop size of 40×40 m2 and a receiver of about 10×10 m2. The TEM data was processed and analyzed, exhibiting high-quality data, with an average of noise level of about ηnoi = 3·10−10V/Am2. The 1D Occam inversion results exhibits a clear three-layered resistivitydepth structure with a second conductive layer of roughly 20 Ωm. The clay pan's resistivity distribution is well-resolved with a global misfit of around 1.1. However, the study site showed 2D effects that were stronlgy visible at the edges of the clay pan, leading to misinterpretations of the TEM data. This was confirmed based on 2D forward modelling. In this manner, to better deal with the observed 2D distortions in the TEM data and to derive a more accurate geometry of the clay pan, the recently developed Julia Package (3DTEMinv) for time-domain 3D inversion and modeling data was performed. The resulting 3D inversion presents a high convergence rate, and acceptable solutions are obtained after ten iterations with a good misfit of about 1.6. The 3D model exhibits a well-resolved geometry of the clay pan, with a high resolution of the derived conductive body. The drill core results confirm the 1D and 3D TEM models at the center of the clay pan, which is in good agreement with the resulting lithology with a maximum thickness of about 171 m depth and a weathered granodioritic bedrock below. These results agree with the local and regional geological context, improving the understanding of sediment deposition and transportation in this hilly and arid environment.

Description: The transition towards renewable energies demands secure supply with critical raw material and requires efficient non-invasive methods for deep earth resources exploration. The novel DESMEX (Deep electromagnetic sounding for mineral exploration) semi-airborne electromagnetic (semi-AEM) exploration concept aims at efficient exploration of resources down to 1 km depth. Here we present a large-scale semi-AEM exploration study in a graphite mining district in eastern Bavaria, Germany. At the ground, several horizontal electrical dipole transmitters were deployed and helicopter-towed magnetic field sensors measure the EM fields along flight lines within several overlapping flight areas, providing a fast data acquisition and a high spatial coverage. Imaged shallow high conductivity structures can be correlated with graphite-rich zones and match well with existing helicopter-borne EM results. The presence of graphite leads to significant induced polarization (IP) effects with considerably high chargeabilities superposing electromagnetic induction. We include these effects in a realistic 3D inversion using a synthetic data study to analyse, if the IP effect alters the overall conductivity structure and demonstrate that the obtained 3D model is reliable.

Description: Transient electromagnetic (TEM) data can be significantly distorted by induced polarization (IP) effect, leading to a sign reversal feature and, if overlooked, false geological interpretation. The aim of this paper is to incorporate IP effects in the forward modelling and recover the distorted TEM data using an efficient inversion algorithm. To achieve this aim, we developed a 1D forward solver to incorporate the IP effects using various IP parameterizations including Cole-Cole, maximum phase angle (MPA), maximum imaginary conductivity (MIC) (Fiandaca et al., 2018) and the Jeffrey transform of Cole-Cole parameters (Ghorbani et al., 2007). For 1D inversion of distorted TEM data we used Levenberg-Marquardt and very fast simulated annealing algorithms. The result of 1D forward calculation and inversion of synthetic IPdistorted TEM data revealed that, for incorporation the IP effects into the TEM data, the Cole- Cole parametrization is more robust and reliable than MPA, MIC, and Jeffrey transform. Moreover, the result of inversion using Levenberg-Marquardt algorithm is strongly depends on the starting model. We successfully implemented these algorithms for 1D inversion of synthetic IP-affected TEM data (Fig. 1 ). For synthetic data generation, a 3-layered half space model with the thickness of the first and second layers of 5 m was considered. The resistivities of the layers from top to bottom are 10, 5 and 300 Ωm, respectively. To include the IP effect, second layer considered to be chargeable with Cole-Cole parameters of m = 0.5, τ = 0.01 s and c = 0.5. TEM central-loop configuration with a loop size of 50*50 m2 and step-off current of 1 A with a zero ramp time was used for data simulation. We evaluated the performance of our algorithm using field data, successfully.

Description: Multi-dimensional inversion of Transient electromagnetic data is a computationally expensive task. Only few developments and practical interpretation tools exist. Here, we present a multidimensional inversion framework for loop source time-domain electromagnetic data. The developed algorithm is a robust, efficient, and user-oriented tool for the multi-dimensional inversion of typical loop source time-domain electromagnetic configurations. A time-domain finite volume discretization and the direct solver MUMPS are utilized to solve the 3D TEM forward problem. An iterative Gauss-Newton optimization method is implemented for the inversion kernel. The code is parallelized for calculating multiple sources simultaneously to accelerate the inversion. Based on exploration tasks, different configurations exist for commonly used loop source TEM configurations and typical field scales. Synthetic examples are used to verify the effectiveness and benchmark the developed 3D algorithm. Considering that TEM data is often gathered along profiles, adjusting the model roughness along the different modeling domain directions, sufficiently constrains to allow for 2D imaging. In addition to the vertical signal components, we also included horizontal components for large scale fixed loop applications. Subsequent to synthetic validation, the inversion algorithm is further verified using ~120 dense TEM soundings collected over a clay pan site in the Atacama Desert, Chile, to provide bedrock geometry information and suitable coring sites. The 3D inversion result provided an excellent depth estimate of sedimentary infill as well as the bedrock topography and was later confirmed by deep coring. Another interesting site is the Roter Kamm impact crater in Namibia. Our preliminary results obtained from largescale multicomponent fixed loop TEM data reveal a sedimentary infill down to ~300 m depth. In conclusion, our presented 3D inversion code is capable to handle data from various exploration scenarios and provides a robust tool for advanced EM interpretation.

Description: In petrophysics, physical rock properties are typically established through laboratory measurements of individual samples. These measurements predominantly relate to the specific sample and can be challenging to associate with the rock as a whole since the physical attributes are heavily reliant on the microstructure, which can vary significantly in different areas. Thus, the obtained values have limited applicability to the entirety of the original rock mass. To examine the dependence of petrophysical measurements based on the variable microstructure, we generate sets of random microstructure representations for a sample, taking into account macroscopic parameters such as porosity and mean grain size. We show that the methodology can adequately mimic the physical behavior of real rocks, showing consistent emulation of the dependence of electrical conductivity on connected porosity according to Archie's law across different types of pore space (micro-fracture, inter-granular, and vuggy, oomoldic pore space). Furthermore, properties such as the internal surface area and its fractal dimension as well as the electrical tortuosity are accessible for the random microstructures and show reasonable behavior. Finally, the possibilities, challenges and meshing strategies for extending the methodology to 3D microstructures are discussed.

Description: This study investigates the optimization of Semi-Airborne Electromagnetic (SAEM) surveys for enhanced subsurface imaging in mineral exploration. It highlights the utility of multi-transmitter systems and explores real data utilization and the challenges of large-scale surveys. With emphasis on Data obtained from DESMEX project surveys. The use of multiple transmitters is crucial. Single transmitters can distort results and mask subsequent bodies. Employing two transmitters on both sides of the target enhances resolution and depth accuracy. results are based on finite element forward operator custEM and pyGIMLi’s inverse solver [1]. substantial advantages of combining single and multi-patch inversion data. This integration results in improved resolution, reduced artifacts, enhanced continuity of geological structures, superior anomaly detection, minimized edge effects, and improved depth penetration [2]. These findings open promising avenues for further exploration and research in geosciences, offering valuable insights into the Earth's subsurface and its intricate geological features. The next logical step involves expanding our methodology to large-scale inversion using more than three transmitters.

Description: We present a finite element software library written in Matlab for the numerical simulation and inversion of electromagnetic fields in two and three dimensions. It is designed in a modular way to easily plug together fundamental building blocks for various electromagnetic applications from DC to the inductive range in the frequency and even time domain. External modules comprise the mesh generator and the equation solver library. Through its homogeneous software concept the adoption to any field application is relatively simple and makes the code suitable to open source distribution. We introduce the key features of this library including higher-order Lagrange and Nédélec finite elements formulated on unstructured tetrahedral grids, a Gauss- Newton inversion approach using linear Raviart-Thomas elements for H1 regularization, and the ability to incorporate any geometric feature such as topography, bathymetry and internal voids like caves, tunnels and mine buildings. The library is currently being tested with large real data sets to confirm its usefulness as a tool for practical data interpretation. Therefore, case studies for the magnetotelluric, direct current resistivity, controlled source electromagnetic and induced polarization methods in the field and laboratory are briefly outlined as examples with challenging geometric features.

Description: play significant roles in coastal hydrologic systems. Despite the importance of these offshore groundwater systems and their interactions with onshore systems along global coastlines, a lack of understanding persists due to limitations in geophysical methodologies. Controlled-source electromagnetic (CSEM) techniques are one promising noninvasive avenue for identifying and characterizing OFG and SGD. However, the current availability of CSEM systems in academic research is limited, and applications are still restricted to specific regions. Existing CSEM systems are commonly associated with high deployment costs, logistical complexity, limited modification options and in case of seafloor-towed applications, slow data acquisition rates. To address these limitations, we introduce SWAN - a low-cost, modular, surface-towed hybrid time-frequency domain CSEM system capable of detecting OFG and SGD up to water depths of 100 m. A field test conducted in the central Adriatic Sea showcased the system's capabilities at water depths ranging from several tens to approximately 160 m. SWAN's ability to provide continuous measurements has proven effective in acquiring high-quality data while operating at towing speeds of 2.5 to 3 knots. The system's data coverage allows for the detection of subsurface resistivity variations to depths of approximately 150–200 m below the seafloor. With its user-friendly, modular design, SWAN offers a cost-efficient solution for investigating the hydrogeology of shallow offshore environments. The presentation shows the technological developments of SWAN, including illustrations of measured time series, processed data and first 2D inversion results from the Adriatic Sea.

Description: Semi-airborne EM (sAEM) transfer functions are estimated by regression over Fourier coefficients over a (small) number of time windows and binned frequencies. The estimation error is determined from the prediction error. However, Fourier coefficients themselves are afflicted with noise and thus erroneously estimated. Here, we make an attempt to incorporate the noise level as weighting coefficients in a weighted least squares regression (WLS). Reliable error estimates are crucial for inverting data; we show 2D inversion results with new error estimates for a test data set using a scalar magnetometer (OPM).

Description: As part of an investigation into applications of Neural Networks for EM problems, different approaches have been tested for DC resistivity modeling and inversion. The first approach consists of using Convolutional Neural Networks (CNN) for DC resistivity inversion. For this purpose, dipole-dipole pseudosections were simulated using the in-house FEMALY toolbox and used as input data for a CNN, which was trained to output underground resistivity. Training results showed qualitatively good match with the ground truth. However, the predictions are characterized by lack of extrapolation to unseen types of data (e.g. homogeneous half-spaces) and coarse grid enforced by the approach. The second approach laid in the use of Physics Informed Neural Networks (PINN). In this approach, the relevant partial differential equation is included as a regularization term in the loss function, leading to a network whose outputs are guided by physics. Derivatives for the PDE termare obtained via automatic differentiation, removing the need for discretization. This also necessitated a move to solving the forward problem. While this approach has the benefits of being mesh-free and incorporating physics into the training process, in practice it failed at even elementary modeling cases, particularly involving resistivity anomalies. A third approach aimed at combining the previous two, by creating a physics-informed Convolutional Neural Network. This was achieved by replacing the previous loss approaches by a convolution with a Laplace-operator Kernel. This approach produces results that look promising qualitatively for homogeneous half-spaces, however full Dirichlet boundary conditions are required and resistivity anomalies can again not be easily incorporated.

Description: Although physical models are improving our understanding of the crustal processes that lead to large earthquakes, observing their preparatory phases is still challenging. We show that the spatio-temporal evolution of the ground motion of small magnitude earthquakes can shed light on the preparatory phase of three main earthquakes that occurred in central Italy between 2016 and 2017. We analyze systematic deviations of peak ground accelerations generated by each earthquake from the values predicted by a reference ground motion model calibrated for background seismicity and refer to such deviations as event-specific ground motion anomalies (eGMAs). The eGMA temporal behavior indicates that during the activation phase of the main earthquakes, the ground shaking level deviates, positively or negatively, from the values expected for the background seismicity. eGMA can be exploited as beacons of stress change and help to monitor the mechanical state of the crust and the nucleation of large earthquakes.

Description: All datasets provided in the operational dataset (Heubeck et al., 2024) of the ICDP project BASE (ICDP 5069) consist of metadata, data and/or images. Here, a summary of explanations of the tables, data and images exported from the database of the project (mDIS BASE) are given and are complimented by additional information on data from measurements done in the laboratory prior to the sampling party. Finally, the sampling data from the first two sam-pling parties are added. Some basic definitions of identifiers used in ICDP, depths corrections and measurements are also introduced.

Description: To gain new insights into ground-motion phenomena in New Zealand (NZ), we apply the non-parametric generalized inversion technique (GIT) in the Fourier domain to isolate the systematic source, path, and site effects from 20 813 seismograms, recorded by 693 sensors at 439 unique locations, from 1200 shallow crustal events (Mw 〉 3) during the period 2000–2021. From the inverted source spectra, we derive Brune's stress parameter, ∆σ, which is found to follow a lognormal distribution with a log10 standard deviation of 0.36 or equivalently 0.83 in natural log unit. ∆σ slightly increases with focal depth and is practically independent of earthquake size (i.e. self-similar), but displays a statistically significant spatial clustering. Based on the inverted attenuation, a trilinear geometric-spreading function, and a distance-dependent quality-factor Q(f) model are found to well describe the attenuation in NZ; though a single model is also obtained for the whole distance range: ⁠. Using the site response decomposed from GIT, we find that the soil classification scheme specified in NZ seismic code, NZS1170.5, has a limited capability in discerning the site-specific frequency-dependent amplification functions in comparison to a non-parametric clustering with the same number of discrete classes. The potential use of the spatial variation in source parameters from this GIT analysis in region-specific physics-based simulations is discussed.

Description: Surface heat flow is a geophysical variable that is affected by a complex combination of various heat generation and transport processes. The processes act on different lengths scales, from tens of meters to hundreds of kilometers. In general, it is not possible to resolve all processes due to a lack of data or modeling resources, and hence the heat flow data within a region is subject to residual fluctuations. We introduce the REgional HEAT-Flow Uncertainty and aNomaly Quantification (REHEATFUNQ) model, version 2.0.1. At its core, REHEATFUNQ uses a stochastic model for heat flow within a region, considering the aggregate heat flow to be generated by a gamma-distributed random variable. Based on this assumption, REHEATFUNQ uses Bayesian inference to (i) quantify the regional aggregate heat flow distribution (RAHFD) and (ii) estimate the strength of a given heat flow anomaly, for instance as generated by a tectonically active fault. The inference uses a prior distribution conjugate to the gamma distribution for the RAHFDs, and we compute parameters for a uninformed prior distribution from the global heat flow database by Lucazeau (2019). Through the Bayesian inference, our model is the first of its kind to consistently account for the variability in regional heat flow in the inference of spatial signals in heat flow data. Interpretation of these spatial signals and in particular their interpretation in terms of fault characteristics (particularly fault strength) form a long-standing debate within the geophysical community. We describe the components of REHEATFUNQ and perform a series of goodness-of-fit tests and synthetic resilience analyses of the model. While our analysis reveals to some degree a misfit of our idealized empirical model with real-world heat flow, it simultaneously confirms the robustness of REHEATFUNQ to these model simplifications. We conclude with an application of REHEATFUNQ to the San Andreas fault in California. Our analysis finds heat flow data in the Mojave section to be sufficient for an analysis and concludes that stochastic variability can allow for a surprisingly large fault-generated heat flow anomaly to be compatible with the data. This indicates that heat flow alone may not be a suitable quantity to address fault strength of the San Andreas fault.

Description: The BASE (Barberton Archean Surface Environments) scientific drilling project focused on recovering unweathered continuous core through strata of the Paleoarchean Moodies Group (ca. 3.2 Ga), central Barberton Greenstone Belt (BGB), South Africa. They comprise some of the oldest well-preserved sedimentary strata on Earth, deposited within only a few million years in alluvial, fluvial, coastal-deltaic, tidal, and prodeltaic settings and represent a very-high-resolution record of Paleoarchean surface conditions and processes. Moodies Group strata consist of polymict conglomerates, widespread quartzose, lithic and arkosic sand-stones, siltstones, shales, and rare BIFs and jaspilites, interbedded with tuffs and several thin lavas. This report describes operations from preparations to the sampling workshop and complements the related scientific report. Eight inclined boreholes between 280 and 495 m length, drilled during November 2021 through July 2022, obtained a total of 2903 m of curated core of variable quality through steeply to subvertically dipping, in part overturned stratigraphic sections. All drilling objec-tives were reached. Boreholes encountered a variety of conglomerates, diverse and abun-dant, mostly tuffaceous sandstones, rhythmically laminated shale-siltstone and banded-iron formations, and several horizons of early-diagenetic sulfate concretions. Oxidative weather-ing reached far deeper than expected; fracturing was more intense, and BIFs and jaspilites were thicker than anticipated. Two km-long mine adits and a water tunnel, traversing four thick stratigraphic sections within the upper Moodies Group in the central BGB, were also sampled. All boreholes were logged by geophysical instruments. Core was processed (ori-ented, slabbed, photographed, described, and archived) in a large, publicly accessible hall in downtown Barberton. An exhibition provided background explanations for visitors and relat-ed the drilling objectives to the recently established Barberton-Makhonjwa Mountains World Heritage Site. A substantial education, outreach and publicity program addressed the information needs of the local population and of local and regional stakeholders.

Description: Highly fragmented rocks (i.e., pulverized rocks) in the fault damage zone presumably develop during co-seismic deformation processes. These pulverized rocks close to the fault core are generally thought to originate from high strain rates, whereas the genesis of pulverized rocks that can be found several hundred meters away from the fault core – where quasi-static conditions prevail – remains unclear. We thus conducted uniaxial cyclic loading experiments with axial strain rate of ∼10−3 s−1 on Leiyang marble in a stress-controlled manner in order to produce crushed rocks for analysis. We found that cyclic loading between 0.8 σc and 1.3 σc can simultaneously compact pre-existing cracks and generated new cracks in marble, which strengthened and stiffened the rock. The stiffened marble developed a higher crack density and energy density before rupture, thereby facilitating rock fragmentation compared with the reference sample, which was fractured monotonically in one cycle. Our results provide a plausible explanation for the genesis of pulverized marble at quasi-static strain rate in the field.