In-situ experiment reveals CO2 enriched fluid migration in faulted caprock.

The sealing characteristics of the geological formation located above a CO2 storage reservoir, the so-called caprock, are essential to ensure efficient geological carbon storage. If CO2 were to leak through the caprock, temporal changes in fluid geochemistry can reveal fundamental information on migration mechanisms and induced fluid-rock interactions. Here, we present the results from a unique in-situ injection experiment, where CO2-enriched fluid was continuously injected in a faulted caprock analogue. Our results show that the CO2 migration follows complex pathways within the fault structure. The joint analysis of noble gases, ion concentrations and carbon isotopes allow us to quantify mixing between injected CO2-enriched fluid and resident formation water and to describe the temporal evolution of water-rock interaction processes. The results presented here are a crucial complement to the geophysical monitoring at the fracture scale highlighting a unique migration of CO2 in fault zones.

Monitoring microseismicity of the Hengill Geothermal Field in Iceland.

Induced seismicity is one of the main factors that reduces societal acceptance of deep geothermal energy exploitation activities, and felt earthquakes are the main reason for closure of geothermal projects. Implementing innovative tools for real-time monitoring and forecasting of induced seismicity was one of the aims of the recently completed COSEISMIQ project. Within this project, a temporary seismic network was deployed in the Hengill geothermal region in Iceland, the location of the nation's two largest geothermal power plants. In this paper, we release raw continuous seismic waveforms and seismicity catalogues collected and prepared during this project. This dataset is particularly valuable since a very dense network was deployed in a seismically active region where thousand of earthquakes occur every year. For this reason, the collected dataset can be used across a broad range of research topics in seismology ranging from the development and testing of new data analysis methods to induced seismicity and seismotectonics studies.

Spectral boundary integral method for simulating static and dynamic fields from a fault rupture in a poroelastodynamic solid.

The spectral boundary integral method is popular for simulating fault, fracture, and frictional processes at a planar interface. However, the method is less commonly used to simulate off-fault dynamic fields. Here we develop a spectral boundary integral method for poroelastodynamic solid. The method has two steps: first, a numerical approximation of a convolution kernel and second, an efficient temporal convolution of slip speed and the appropriate kernel. The first step is computationally expensive but easily parallelizable and scalable such that the computational time is mostly restricted by computational resources. The kernel is independent of the slip history such that the same kernel can be used to explore a wide range of slip scenarios. We apply the method by exploring the short-time dynamic and static responses: first, with a simple source at intermediate and far-field distances and second, with a complex near-field source. We check if similar results can be attained with dynamic elasticity and undrained pore-pressure response and conclude that such an approach works well in the near-field but not necessarily at an intermediate and far-field distance. We analyze the dynamic pore-pressure response and find that the P-wave arrival carries a significant pore pressure peak that may be observed in high sampling rate pore-pressure measurements. We conclude that a spectral boundary integral method may offer a viable alternative to other approaches where the bulk is discretized, providing a better understanding of the near-field dynamics of the bulk in response to finite fault ruptures.

Anatomy of a fumarolic system inferred from a multiphysics approach.

Fumaroles are a common manifestation of volcanic activity that are associated with large emissions of gases into the atmosphere. These gases originate from the magma, and they can provide indirect and unique insights into magmatic processes. Therefore, they are extensively used to monitor and forecast eruptive activity. During their ascent, the magmatic gases interact with the rock and hydrothermal fluids, which modify their geochemical compositions. These interactions can complicate our understanding of the real volcanic dynamics and remain poorly considered. Here, we present the first complete imagery of a fumarolic plumbing system using three-dimensional electrical resistivity tomography and new acoustic noise localization. We delineate a gas reservoir that feeds the fumaroles through distinct channels. Based on this geometry, a thermodynamic model reveals that near-surface mixing between gas and condensed steam explains the distinct geochemical compositions of fumaroles that originate from the same source. Such modeling of fluid interactions will allow for the simulation of dynamic processes of magmatic degassing, which is crucial to the monitoring of volcanic unrest.