Creep of CarbFix basalt: influence of rock–fluid interaction

Geological carbon sequestration provides permanent CO2 storage to mitigate the current high concentration of CO2 in the atmosphere. CO2 mineralization in basalts has been proven to be one of the most secure storage options. For successful implementation and future improvements of this technology, the time-dependent deformation behavior of reservoir rocks in the presence of reactive fluids needs to be studied in detail. We conducted load-stepping creep experiments on basalts from the CarbFix site (Iceland) under several pore fluid conditions (dry, H2O saturated and H2O + CO2 saturated) at temperature, T≈80 ∘C and effective pressure, Peff=50 MPa, during which we collected mechanical, acoustic and pore fluid chemistry data. We observed transient creep at stresses as low as 11 % of the failure strength. Acoustic emissions (AEs) correlated strongly with strain accumulation, indicating that the creep deformation was a brittle process in agreement with microstructural observations. The rate and magnitude of AEs were higher in fluid-saturated experiments than in dry conditions. We infer that the predominant mechanism governing creep deformation is time- and stress-dependent subcritical dilatant cracking. Our results suggest that the presence of aqueous fluids exerts first-order control on creep deformation of basaltic rocks, while the composition of the fluids plays only a secondary role under the studied conditions.

Forearc density structure of the overriding plate in the northern area of the giant 1960 Valdivia earthquake

The objective of this work is to analyse the density structure of the continental forearc in the northern segment of the 1960 Mw 9.6 Valdivia earthquake. Regional 2D and local 3D density models have been obtained from available gravity data in the area, complemented by new gravimetric stations. Models are constrained by independent geophysical and geological information and new TEM and MT soundings. The results show a segmentation of the continental wedge along and perpendicular to the margin, highlighting a high-density anomaly, below the onshore forearc basin, that limits the late Paleozoic–early Mesozoic metamorphic basement in the region where Chaitenia terrane has been proposed. A progressive landward shift of this anomaly correlates with the high slip patch of the giant 1960 Mw 9.6 Valdivia earthquake. Based on these results, we propose that the horizontal extension of the less rigid basement units conforming the marine wedge and Coastal Cordillera domain could modify the process of stress loading during the interseismic periods, and also that changes in position and extension of the late Paleozoic–early Mesozoic accretionary complex could be linked with the frictional properties of the interplate boundary. This analysis provides new evidence of the role of the overriding plate structure in the seismotectonic process in subduction zones.

Impact of Timanian thrust systems on the late Neoproterozoic–Phanerozoic tectonic evolution of the Barents Sea and Svalbard

The Svalbard Archipelago consists of three basement terranes that record a complex Neoproterozoic–Phanerozoic tectonic history, including four contractional events (Grenvillian, Caledonian, Ellesmerian, and Eurekan) and two episodes of collapse- to rift-related extension (Devonian–Carboniferous and late Cenozoic). Previous studies suggest that these three terranes likely accreted during the early to mid-Paleozoic Caledonian and Ellesmerian orogenies. Yet recent geochronological analyses show that the northwestern and southwestern terranes of Svalbard both record an episode of amphibolite (–eclogite) facies metamorphism in the latest Neoproterozoic, which may relate to the 650–550 Ma Timanian Orogeny identified in northwestern Russia, northern Norway, and the Russian Barents Sea. However, discrete Timanian structures have yet to be identified in Svalbard and the Norwegian Barents Sea. Through analysis of seismic reflection, as well as regional gravimetric and magnetic data, this study demonstrates the presence of continuous thrust systems that are several kilometers thick, NNE-dipping, deeply buried, and extend thousands of kilometers from northwestern Russia to northeastern Norway, the northern Norwegian Barents Sea, and the Svalbard Archipelago. The consistency in orientation and geometry, as well as apparent linkage between these thrust systems and those recognized as part of the Timanian Orogeny in northwestern Russia and Novaya Zemlya, suggests that the mapped structures are likely Timanian. If correct, these findings would imply that Svalbard's three basement terranes and the Barents Sea were accreted onto northern Norway during the Timanian Orogeny and should hence be attached to Baltica and northwestern Russia in future Neoproterozoic–early Paleozoic plate tectonics reconstructions. In the Phanerozoic, the study suggests that the interpreted Timanian thrust systems represent major preexisting zones of weakness that were reactivated, folded, and overprinted by (i.e., controlled the formation of new) brittle faults during later tectonic events. These faults are still active at present and can be linked to folding and offset of the seafloor.

A functional tool to explore the reliability of micro-earthquake focal mechanism solutions for seismotectonic purposes

Improving the knowledge of seismogenic faults requires the integration of geological, seismological, and geophysical information. Among several analyses, the definition of earthquake focal mechanisms plays an essential role in providing information about the geometry of individual faults and the stress regime acting in a region. Fault plane solutions can be retrieved by several techniques operating in specific magnitude ranges, both in the time and frequency domain and using different data. For earthquakes of low magnitude, the limited number of available data and their uncertainties can compromise the stability of fault plane solutions. In this work, we propose a useful methodology to evaluate how well a seismic network, used to monitor natural and/or induced micro-seismicity, estimates focal mechanisms as a function of magnitude, location, and kinematics of seismic source and consequently their reliability in defining seismotectonic models. To study the consistency of focal mechanism solutions, we use a Bayesian approach that jointly inverts the P/S long-period spectral-level ratios and the P polarities to infer the fault plane solutions. We applied this methodology, by computing synthetic data, to the local seismic network operating in the Campania–Lucania Apennines (southern Italy) aimed to monitor the complex normal fault system activated during the Ms 6.9, 1980 earthquake. We demonstrate that the method we propose is effective and can be adapted for other case studies with a double purpose. It can be a valid tool to design or to test the performance of local seismic networks, and more generally it can be used to assign an absolute uncertainty to focal mechanism solutions fundamental for seismotectonic studies.

De-risking the energy transition by quantifying the uncertainties in fault stability

The operations needed to decarbonize our energy systems increasingly involve faulted rocks in the subsurface. To manage the technical challenges presented by these rocks and the justifiable public concern over induced seismicity, we need to assess the risks. Widely used measures for fault stability, including slip and dilation tendency and fracture susceptibility, can be combined with response surface methodology from engineering and Monte Carlo simulations to produce statistically viable ensembles for the analysis of probability. In this paper, we describe the implementation of this approach using custom-built open-source Python code (pfs – probability of fault slip). The technique is then illustrated using two synthetic examples and two case studies drawn from active or potential sites for geothermal energy in the UK and discussed in the light of induced seismicity focal mechanisms. The analysis of probability highlights key gaps in our knowledge of the stress field, fluid pressures, and rock properties. Scope exists to develop, integrate, and exploit citizen science projects to generate more and better data and simultaneously include the public in the necessary discussions about hazard and risk.

Biotite supports long-range diffusive transport in dissolution–precipitation creep in halite through small porosity fluctuations

Phyllosilicates are generally regarded to have a reinforcing effect on chemical compaction by dissolution–precipitation creep (DPC) and thereby influence the evolution of hydraulic rock properties relevant to groundwater resources and geological repositories as well as fossil fuel reservoirs. We conducted oedometric compaction experiments on layered NaCl–biotite samples to test this assumption. In particular, we aim to analyse slow chemical compaction processes in the presence of biotite on the grain scale and determine the effects of chemical and mechanical feedbacks. We used time-resolved (4-D) microtomographic data to capture the dynamic evolution of the porosity in layered NaCl–NaCl/biotite samples over 1619 and 1932 h of compaction. Percolation analysis in combination with advanced digital volume correlation techniques showed that biotite grains influence the dynamic evolution of porosity in the sample by promoting a reduction of porosity in their vicinity. However, the lack of preferential strain localisation around phyllosilicates and a homogeneous distribution of axial shortening across the sample suggests that the porosity reduction is not achieved by pore collapse but by the precipitation of NaCl sourced from outside the NaCl–biotite layer. Our observations invite a renewed discussion of the effect of phyllosilicates on DPC, with a particular emphasis on the length scales of the processes involved. We propose that, in our experiments, the diffusive transport processes invoked in classical theoretical models of DPC are complemented by chemo-mechanical feedbacks that arise on longer length scales. These feedbacks drive NaCl diffusion from the marginal pure NaCl layers into the central NaCl–biotite mixture over distances of several hundred micrometres and several grain diameters. Such a mechanism was first postulated by Merino et al. (1983).

Virtual field trip to the Esla Nappe (Cantabrian Zone, NW Spain): delivering traditional geological mapping skills remotely using real data

The restrictions implemented to contain the spread of the COVID-19 pandemic during 2020 and 2021 have forced university-level educators from around the world to seek alternatives to the residential physical field trips that constitute a fundamental pillar of Geoscience programmes. The field-mapping course for second-year Geology BSc students from Cardiff University was replaced with a virtual mapping course set in the same area as previous years, the Esla Nappe (Cantabrian Zone, NW Spain). The course was designed with the aim of providing the students with the same methodology employed in physical mapping, including such skills as gathering discrete data at stops located along five daily itineraries. Data included bedding attitude, outcrop descriptions with a certain degree of ambiguity, photographs and/or sketches, panoramic photos, and fossil images. Data were provided to the students through georeferenced KMZ files in Google Earth. Students were asked to keep a field notebook, define lithological units of mappable scale, identify large structures such as thrust faults and folds with the aid of age estimations from fossils, construct a geological map on a hard-copy topographic map, draw a stratigraphic column and cross sections, and plot the data in a stereonet to perform structural analysis. The exercise allowed for successful training of diverse geological field skills. In light of the assessment of reports and student surveys, a series of improvements for the future is considered. Though incapable of replacing a physical field course, the virtual exercise could be used in preparation for the residential field trip.

The Mid Atlantic Appalachian Orogen Traverse: a comparison of virtual and on-location field-based capstone experiences

The Stratigraphy, Structure, Tectonics (SST) course at James Madison University incorporates a capstone project that traverses the Mid Atlantic region of the Appalachian Orogen and includes several all-day field trips. In the Fall 2020 semester, the SST field trips transitioned to a virtual format, due to restrictions from the COVID pandemic. The virtual field trip projects were developed in web-based Google Earth and incorporated other supplemental PowerPoint and PDF files. In order to evaluate the effectiveness of the virtual field experiences in comparison with traditional on-location field trips, an online survey was sent to SST students that took the course virtually in Fall 2020 and to students that took the course in person in previous years. Instructors and students alike recognized that some aspects of on-location field learning, especially those with a tactile component, were not possible or effective in virtual field experiences. However, students recognized the value of virtual field experiences for reviewing and revisiting outcrops as well as noting the improved access to virtual outcrops for students with disabilities and the generally more inclusive experience of virtual field trips. Students highlighted the potential benefits for hybrid field experiences that incorporate both on-location outcrop investigations and virtual field trips, which is the preferred model for SST field experiences in Fall 2021 and into the future.

Deep oceanic submarine fieldwork with undergraduate students: an immersive experience with the Minerve software

We present the content and scripting of an active tectonic lab session conceived for third-year undergraduate students studying Earth sciences at Observatoire des Sciences de l'Univers in Lyon. This session is based on a research project conducted on the submarine Roseau active fault in the Lesser Antilles. The fault morphology is particularly interesting to map as this structure in the deep ocean is preserved from weathering. Thus, high-resolution models computed from remotely operated vehicle (ROV) videos provide exceptional educational material to link fault morphology and coseismic displacement. This class includes mapping exercises on geographical information systems and virtual fieldwork to provide basic understanding of active tectonics and active fault morphology in particular. The work has been conducted either in a full remote configuration via 3D online models or in virtual reality (VR) in a dedicated room using the Minerve software. During the VR sessions, students were either alone in the VR environment or participated as a group that included the instructor (physically in the classroom or remotely from another location), which is to our knowledge one of the first attempts of this kind in France. We discuss the efficiency of virtual fieldwork using VR based on feedback from teachers and students. We conclude that VR is a promising tool to learn observational skills in Earth sciences, subject to certain improvements that should be possible in the years to come.

Sedimentary basins of the eastern Asia Arctic zone: new details on their structure revealed by decompensative gravity anomalies

In the present study, the structure of sedimentary basins in the eastern Asia Arctic zone is analysed by employing the approach based on decompensative gravity anomalies. Two obtained models, differing in their initial conditions, provide thickness and density of sediments in the study area. They demonstrate essentially new details on the structure, shape, and density of the sedimentary basins. Significant changes in the sedimentary thickness and the depo-centre location have been found for the Anadyr Basin in its continental part. Also, new details on the sedimentary thickness distribution have been revealed for the central part of the Penzhin and Pustorets basins; for the latter, the new location of the depo-centre has been identified. The new model agrees well with the seismic data on the sedimentary thickness for the offshore part of the Chauna Basin confirming that the method is robust. The most significant lateral redistribution of the thickness has been found for the Lower Cretaceous coal-bearing strata in the northern part of the Zyryanka Basin, where the connection of two coal-bearing zones, which was not previously mapped, has been identified. Also, the new details on the sedimentary thickness distribution have been discovered for the Primorsk Basin. Therefore, the new results substantially improve our knowledge about the region, since previous geological and geophysical studies were unsystematic, sparse, and limited in depth. Thus, the implementation of the decompensative gravity anomalies approach provides a better understanding of the evolution of the sedimentary basins and the obtained results can be used for planning future detailed studies in the area.