Ground Uplift Research Articles

Study on Geological Deformation of Supercritical CO2 Sequestration in Oil Shale after In Situ Pyrolysis

After the completion of in situ pyrolysis, oil shale can be used as a natural place for CO2 sequestration. However, the effects of chemical action and formation stress-state changes on the deformation of oil shale should be considered when CO2 is injected into oil shale after pyrolysis. In this study, combined with statistical damage mechanics, a transverse isotropic model of oil shale with coupled damage mechanisms was established by considering the decreased mechanical properties and the chemical damage caused by CO2 injection. The process of injecting supercritical CO2 into oil shale after pyrolysis was simulated by COMSOL6.0. The volume distribution of CO2 and the stress evolution in oil shale were analyzed. It is found that CO2 injection into oil shale after pyrolysis will not produce new force damage, and the force damage caused by the decrease in the mechanical properties of oil shale after pyrolysis can offset the ground uplift caused by CO2 injection to a certain extent. Under the combined action of chemical damage and mechanical damage, the uplift of a formation with a thickness of 200 m is only 10 cm. The injection of supercritical CO2 is beneficial for maintaining the stability of oil shale after in situ pyrolysis.

Land surface response to groundwater drawdown and recovery in Taiyuan city, Northern China, analyzed with a long-term elevation change measurements from leveling and multi-sensor InSAR

In recent years, a comprehensive program of groundwater management was adopted in Taiyuan city, Northern China, with the aim to prevent the lowering of groundwater level and manage land subsidence due to excessive groundwater pumping in the past decades. Groundwater recovery has been observed in this city following the termination of groundwater pumping. While the link between land subsidence and groundwater pumping is well documented, how the land surface responds to the recovery of groundwater has not been studied in detail. In this study, we analyzed the leveling data from 1956 to 2000 to investigate changes in ground elevation in response to groundwater drawdown prior to the termination of groundwater pumping and combined them with InSAR (Interferometric Synthetic Aperture Radar) displacement data during 2003–2020 to evaluate the effect of groundwater recharge after the termination of groundwater pumping. A method was proposed to merge InSAR time series of displacement derived from multiple satellites (ENVISAT, COSMO-SkyMed, TerraSAR-X, and Sentinel-1) to achieve a consistent, long-term continuous deformation. We found that previous groundwater depression cones in Taiyuan (Xizhang, Wanbailin, Wujiabao, and Xiayuan) have turned into groundwater recovery and the corresponding subsidence centers have switched to land uplift. The subsidence center of Xizhang, north of Taiyuan city, experienced a total subsidence of 749 mm by Aug 2003; it then turned into uplift with an amount of rebound of 9.2 % of previous subsidence by the end of 2020. In central area of Taiyuan, the cumulated subsidence is 1127 mm in Wanbailin, 1817 mm in Xiayuan, and 3343 mm in Wujiabao. Beginning in Aug 2010, they all turned into land uplift, with a rebound of 7.9 %, 4.8 % and 3.6 % of the previous subsidence in each center, respectively. An analysis of the correlation between groundwater level and ground displacement suggests that the ground uplift is related to the poro-elastic rebound mechanism in the porous aquifer system, which is caused by the rise in pore pressure with groundwater recovery, and the amount of subsidence/uplift is controlled by the stratigraphic properties. The adverse impacts on the built environment due to the rising groundwater level were discussed. The findings improve our understanding of anthropogenic ground deformation in complex urban aquifers influenced by groundwater pumping and recharge and provides essential information that can contribute to future groundwater management and groundwater-related hazard mitigation over the Taiyuan city.

Delineation and Fine‐Scale Structure of Fault Zones Activated During the 2014–2024 Unrest at the Campi Flegrei Caldera (Southern Italy) From High‐Precision Earthquake Locations

AbstractIn the past two decades, the central portion of Campi Flegrei caldera has experienced ground uplift up to 15 mm/month, with an increase of rate, magnitude and extent of the seismicity. In this work, we perform multi‐scale precise earthquake relocation of the 2014–2024 seismicity, mapping in detail activated fault zones. We relate the geometry, extent, and depth of these zones with up‐to‐date structural reconstructions of the caldera. The current seismicity is mainly driven by ground‐uplift‐induced stress concentration on pre‐existing, weaker fault zones, some of which identified for the first time. These structures are not only related to the inner caldera and dome resurgence but also to volcano‐tectonic events of the last 10 ka. The extent of imaged fault segments suggests they can accommodate ruptures up to a moment magnitude 5.1, significantly increasing seismic hazard in the area.

Complex variable solution for ground deformation responses considering the grouting pressure of shallow shield

During shield tunneling in shallow strata, ground deformation is highly sensitive to grouting pressure response. Excessive grouting pressure can cause ground uplift, and an excessively low grouting pressure can cause excessive ground settlement, affecting the surrounding environment. Current studies seldom analyze the influence of shield excavation on ground deformation when the grouting pressure is considered. To ensure micro-disturbance of the surrounding environment during shield tunneling, this study presents a complex variable analysis method for analyzing the influence of grouting pressure on surface settlement. This solution maps the semi-infinite strata outside the segment into circular domains by introducing a complex function conformal mapping method. The nonuniform convergence deformation pattern BC-4 was introduced as the boundary condition for the displacement around the tunnel. The Laurent series expansion of the analytical function considering the grouting pressure in the circular domain was substituted into each boundary condition to solve the analytical function. Subsequently, the accuracy of the analytical method was verified by combining it with a finite element numerical method. Finally, the influence of grouting pressure on ground deformation under different working conditions was analyzed through parameterization. This analytical method can provide a theoretical basis for controlling grouting pressure parameters in shield tunnels.

Investigating Earth surface deformation with SAR interferometry and geomodeling in the transborder Meuse–Rhine region

Introduction: This study aims to differentiate local and regional ground uplift, as well as sub-regional subsidence induced by groundwater level drawdown, which are possibly enhanced across fault structures, as monitored by various synthetic aperture radar interferometry (InSAR) processing methods. A buoyant mantle plume under the Eifel may be responsible for the regional ground uplift, including the Weser–Geul (BE) and South Limburg regions (NL), which could negatively affect the area proposed for the future Einstein Telescope.Methods: Different InSAR processing techniques are compared to evaluate their limits in tracking fault structures on a time series of Copernicus Sentinel-1 images while detecting and measuring ground motion based on their phase signature. The results present an overall stable ground for the Euregio Meuse–Rhine region, especially at the Belgian–Dutch border, while tectonic activity is observed along the German side of the Rhine Graben.Results: As the current neotectonic activity in the target area was not well known, we performed a spatiotemporal analysis of ground deformation associated with the presence of NW–SE-trending normal faults where karst also develops, as well as along the Variscan NE–SW-trending thrust faults. This work demonstrates that the identification of deformation hazards using satellite remote sensing (and connected seismological) techniques is challenging mainly due to (very) small regional scale deformation, terrain conditions, and SAR properties.Discussion: Thus, the results mostly indicate ground stability over the area; however, also some agricultural activities were observed, as was deformation along some infrastructure such as railways. Displacements of millimetric order measured along the faults located near the Geul valley (BE) are probably related to old mining activities.

Ground Deformation of Shield Tunneling through Composite Strata in Coastal Areas

In order to mitigate ground deformation during shield construction in both upper soft and lower hard strata of coastal areas, a numerical simulation was executed. This simulation assessed surface deformation under varying stratum ratios, grouting pressures, and earth bin pressures. The evaluation was primarily based on the amount of ground deformation, which revealed that hard rock strata offer superior settlement control compared to soft rock strata. The excavation of the right tunnel line increased disturbance to the left line at higher stratum ratios. Surface deformation demonstrated a linear correlation with earth pressure, with 130 kPa identified as the optimal point. Higher pressures resulted in extrusion deformation and ground uplift. Grouting pressure had a minimal impact on stratum deformation over time. The stratum ratio exerted the most significant influence on settlement, followed by earth pressure, with grouting pressure having the least impact. In the context of coastal tunnel construction, hard rock excavation is favored. Earth pressure must be balanced to prevent subsidence or uplift, while excessive grouting pressure does not significantly reduce subsidence. Grouting pressure should ensure the complete filling of voids.

Assessing Mechanical Properties and Response of Expansive Soft Rock in Tunnel Excavation: A Numerical Simulation Study.

This study investigates the dynamics of moisture absorption and swelling in soft rock during tunnel excavation, emphasizing the response to support resistance. Utilizing COMSOL numerical simulations, we conduct a comparative analysis of various strength criteria and non-associated flow rules. The results demonstrate that the Mohr-Coulomb criterion combined with the Drucker-Prager model under compressive loads imposes stricter limitations on water absorption and expansion than when paired with the Drucker-Prager model under tensile loads. Restricted rock expansion leads to decreased horizontal displacement and ground uplift, increased displacement in the tunnel's bottom arch, and significantly reduced displacement in the top arch. The study also considers the effects of shear dilation, burial depth, and support resistance on the stress and displacement of the surrounding rock. Increased shear dilation angles correlate with greater rock expansion, resulting in increased horizontal displacement and ground uplift. The research study concludes that support resistance is critical in limiting the movement of the tunnel's bottom arch and impacting the stability of the surrounding rock. Additionally, the extent of rock damage during the excavation of expansive soft rock tunnels is found to be minimal. Overall, this study provides valuable insights into the processes of soft rock tunnel excavation and contributes to the development of more efficient support systems.

Sustainable construction and quality of improved columns with three types of water-cement ratios on deep mixing method in Saga Lowland, Kyushu, Japan

In this study, we investigated the application of the deep mixing method (DMM) to cohesive soil in the Saga Lowland of Kyushu, Japan. The study focused on examining three types of water-cement ratio (W/C) conditions, with a constant addition of cement-based binder (C). Due to the soft clay nature of the Saga Lowland, frequent ground settlement and deformation occur, necessitating measures to prevent adverse effects on the surrounding environment. The objective of this research is to provide a valuable approach to optimizing the quality of the improved columns while minimizing ground displacement in an environmentally considerate manner. In Saga Lowland, it is common to fix the W/C ratio at 1.0 and vary the cement content. Through experimental construction for improved columns on the field, the study confirmed that W/C values of 0.5, 1.0, and 1.5 influence the quality of the improved structure. A higher W/C value of 1.5 resulted in a more fluid cement slurry due to a higher injection rate (IR = 23.9%), as evidenced by statistical analysis revealing higher average unconfined compressive strength ( qu¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\overline{{q }_{u}}$$\\end{document}), and a lower coefficient of variation (CV). The defective rate of 10% (qudr) from the design standard strength shows that values are lowest for Case 2, followed by Case 3 and then Case 1. Comparing the values of Case 2 and Case 3, it is observed that in Case 3, with a higher W/C, the CV is lower.. Regarding horizontal ground displacement (Sh), Case 3 exhibited a Sh value of 2.0 to 6.5 mm, significantly lower than Saga prefecture standards (20 mm). This outcome is attributed to reduced viscosity during mixing, leading to improved fluidity and minimal lateral displacement of the soil–cement columns (which often results in lateral ground uplift). Even though with a higher W/C = 1.5, the implementation cost remains the same, but the constructed structure would be of higher quality and smaller displacement, with the overall structure corresponding to the standard quality. The study includes the specific geotechnical conditions of the Saga Lowland and the scope of experimentation. Nonetheless, in terms of the applicability and optimization of DMM in Saga Lowland, the findings provide practical guidance for engineers in selecting W/C ratio and IR during construction for future DMM implementations, thereby contributing to the development of long-lasting infrastructure and sustainable societal development.

Monitoring Analysis of Urban Subsidence in Northern Henan Province Based on TS-InSAR Technology

The protracted and pervasive incidence of land subsidence emerges as a pivotal factor exerting a substantial impact on the sustainable development of urban landscapes. A nuanced comprehension of the spatiotemporal evolution characteristics of land subsidence within the Northern Henan Plain assumes paramount significance in the context of mitigating potential urban geological disasters. This study endeavors to redress the deficiency in information concerning the temporal and spatial evolution characteristics of enduring deformation in cities within the northern plain of Henan Province. To this end, the authors leveraged Sentinel-1A radar data processed through persistent scatterer interferometric synthetic aperture radar (PS-InSAR) technology to elucidate the distribution patterns of ground deformation and temporal evolution characteristics within the expansive 24-scene coverage research area. Empirical findings illuminate conspicuous surface deformation in Anyang, Puyang, and Hebi throughout the monitoring period. Spatially, land subsidence in the study area predominantly clusters in the suburban peripheries of the cities, with Hebi and Puyang registering a maximum subsidence rate exceeding 25 mm per annum. Temporally, land subsidence manifests predominantly during autumn and winter, whereas spring and summer display relatively stable land subsidence interspersed with a slight ground uplift. In order to rectify the spatial disparities observed between leveling data and PS-InSAR monitoring data, this experiment employed an averaging procedure on the PS-InSAR monitoring data, subsequently subjecting it to comparative analysis with the leveling data. Additionally, through the integration of the singular spectrum analysis (SSA) method and the time series deformation model, this study aspires to attain a comprehensive understanding of the temporal dynamics manifested in the PS-InSAR monitoring outcomes, while concurrently elucidating the factors influencing the observed deformations. Ultimately, this analysis discloses that the monitoring outcomes derived via PS-InSAR technology exhibit a root mean square error of ±12.9 mm and a standard deviation of ±13.31 mm. These statistical metrics furnish valuable insights into the precision and consistency of the PS-InSAR monitoring data. Drawing upon a comparative scrutiny of on-site data and historical remote sensing imagery within the study area, it has been discerned that excessive groundwater extraction and expansive surface engineering initiatives stand as the principal instigators of land subsidence in the research domain. Consequently, this experiment assumes the role of a salient reference for the mitigation of urban ground subsidence within the study area.

Research On Soil Displacement Prediction Method Caused By High-pressure Jet Grouting Pile Group Construction Under Adjacent Pile Response

This paper analyzes the soil displacement effect during high-pressure jet grouting pile construction by simplifying the construction process as a series of pressure-controlled spherical cavity expansions in a semi-infinite soil medium. Based on the elastic-plastic solution of pressure-controlled spherical cavity expansion theory, a calculation method is proposed to estimate the soil displacement caused by high-pressure jet grouting pile construction. Non-linear contact theory between the pile and soil is introduced, and the finite difference method is employed to determine the internal forces and deformations of the pile. The proposed method is applied to field construction cases, and its validity is verified by comparing the results with on-site monitoring data. The research findings reveal that the soil displacement effect decreases with depth and primarily occurs in shallow soil layers during the construction of individual high-pressure jet grouting piles. The lateral displacement at the ground surface increases first and then decreases with the increase in horizontal distance, while the ground uplift exhibits an exponential decay trend.

Evolution in unrest processes at Campi Flegrei caldera as inferred from local seismicity

Changes in seismicity with time and location are diagnostic signals for understanding the dynamics of volcanic unrest. We used these signals at the Campi Flegrei caldera, in southern Italy, to investigate how structural changes have determined three styles of unrest since 1982, distinguished by a ground uplift (measured at Pozzuoli, near the centre of the caldera) of 178 cm in 1982–84; a subsidence of 93 cm in 1985–2005; and an uplift of 118 cm between 2005 and November 2023. Double-difference seismic locations and concentrations of seismic energy release have revealed impermeable horizons that correspond to the cap rock and self-sealed base of the geothermal system at depths of 1.5 and 3 km, respectively. Most earthquakes have been shallower than 3 km, consistent with the brittle upper crust being stretched over a zone of pressurization below the geothermal system. The 1982–84 uplift decayed after a major seismic swarm on April 1st, 1984, breached the lower impermeable horizon, which increased the flux of escaping gas and reduced the source pressure. Continued gas escape promoted subsidence until the lower horizon had resealed itself and initiated a new episode of uplift while gas from depth re-accumulated beneath. Compared with 1982–84, a greater proportion of recent seismicity has occurred at shallower depths below the actively degassing fumaroles of Solfatara-Pisciarelli, north-east of Pozzuoli. The associated high seismic b‐values, between 1 and 2, are consistent with high fluid pressures and, hence, with locations more favourable to fracturing. Continued uplift may thus persist until fracturing of the shallow crust allows faster rates of gas release and depressurization of the pressure source.

Transient ground bulge derived from a dense broadband seismic array during an aquifer step-drawdown pumping test

We deployed nine broadband seismometers within 32 m of a well where a four-hour step-drawdown pumping test was conducted at a site near Shuangxi, Taiwan in August 2018. The whole seismic array simultaneously recorded hour-long seismic signals that we interpret as ground tilt evolution during and after the pumping. Here we document data from the six seismic stations within 10 m of the boreholes to show the ground deformation closest to the pump site. Instead of ground subsidence due to groundwater over-draw as documented in many long-term and even short-term GPS and tiltmeter studies, four out of the six seismometers show ground bulge at the pumping site with radial directional tilts up to 1.0 × 10-5 rad for a drawdown of 10 m. We interpret the bulge as the dominant transient crustal response to the unloading of discharged groundwater, outpacing the minor subsidence due to reduced effective stress from reduced pore fluid head. After the pumping stopped, the ground tilt signals started to level back, though not simultaneously: some signals recorded delays up to an hour, possibly related to replenishment of groundwater from far-field or different aquifers and shallow depth hydrological heterogeneity including minor faults and fractures. We conducted seismic reflection surveys and found mostly layered structures without large fault offsets. Using an analytical solution that calculates the ground tilt due to surface load changes, we calculated the effects of unloading on ground deformation. The predicted tilt signals are in the range of 10-5 rad, the same order as the seismically derived tilts from this study. For 10 m of drawdown we estimate 10-5 to10-4 m of vertical ground uplift for this particular site. Our study demonstrates that dense broadband seismic arrays may offer a relatively less expensive and non-invasive tool to study temporal and spatial evolution of ground motions, and thus the ground water flow patterns, during pumping tests.

Ionospheric Total Electron Content (TEC) Anomalies as Earthquake Precursors: Unveiling the Geophysical Connection Leading to the 2023 Moroccan 6.8 Mw Earthquake

The study delves into the relationship between ionospheric total electron content (TEC) anomalies and seismic activity, with a focus on Morocco’s 6.8 Mw earthquake on 8 September 2023, lying within a tectonically active region at the convergence of the African and Eurasian Plates. To enhance the reliability of our findings, we incorporate space weather conditions, utilizing indices (Dst, Kp, and F10.7) to pinpoint periods of stable space weather. This minimizes the possibility of erroneously attributing natural ionospheric fluctuations to seismic events. Notably, our TEC analysis unveils positive and negative anomalies, with some occurring up to a week before the earthquake. These anomalies, exceeding predefined thresholds, provide compelling evidence of significant deviations from typical ionospheric conditions. Spatial mapping techniques employing both station-specific vTEC data and pseudorandom noise codes (PRNs) from multiple global navigation satellite system (GNSS) stations highlight a strong correlation between ionospheric anomalies and the earthquake’s epicenter. The integration of PRNs enhances coverage and sensitivity to subtle anomalies. Additionally, the analysis of satellite imagery and ground displacement data using Sentinel-1 confirms significant ground uplift of approximately 15 cm following the earthquake, shedding light on surface responses to seismic events. These findings underscore the potential of ionospheric science in advancing earthquake early warning systems and deepening our understanding of earthquake precursors, thus contributing to the mitigation of seismic event impacts and the protection of lives and infrastructure.

Potential for rupture before eruption at Campi Flegrei caldera, Southern Italy

Volcanoes reawakening after long repose must rupture the crust before magma can erupt. Rupture is preceded by repeatable variations in the rate of seismicity with ground movement, which trace the amount of applied stress that is released by local earthquakes. A rupturing sequence has been developing across four episodes of ground uplift at Italy’s Campi Flegrei caldera: in 1950-1952, 1969-1972, 1982-1984 and since 2004. We predicted in 2016 that the approach to rupture would continue after an additional uplift of 30-40 cm at the location of largest movement. We have updated our analysis with new data on changes in the numbers of local earthquakes with amounts of ground movement. Here we show that subsequent events have confirmed our prediction and that the unrest has been changing the structure of Campi Flegrei’s crust. The results provide new constraints for evaluating the volcano’s potential to erupt or to subside without eruption.

Shallow Geologic Storage of Carbon to Remove Atmospheric CO2 and Reduce Flood Risk.

Geologic carbon storage currently implies that CO2 is injected into reservoirs more than 1 km deep, but this concept of geologic storage can be expanded to include the injection of solid, carbon-bearing particles into geologic formations that are one to two orders of magnitude shallower than conventional storage reservoirs. Wood is half carbon, available in large quantities at a modest cost, and can be milled into particles and injected as a slurry. We demonstrate the feasibility of shallow geologic storage of carbon by a field experiment, and the injection process also raises the ground surface. The resulting CO2 storage and ground uplift rates upscale to a technique that could contribute to the mitigation of climate change by storing carbon as well as helping to adapt to flooding risks by elevating the ground surface above flood levels. A life-cycle assessment indicates that CO2 emissions caused by shallow geologic storage of carbon are a small fraction of the injected carbon.

Efficient SOM’s Application to Seismic Fumarolic Tremor for the Detection of Anomalous Hydrothermal Activity in Campi Flegrei Volcano (Italy)

In recent decades, the Campi Flegrei caldera (Italy) showed unrest characterized by increases in seismicity, ground uplift, and hydrothermal activity. Currently, the seismic and hydrothermal phenomena are mostly concentrated in the Solfatara–Pisciarelli area, which presents a wide fumarolic field and mud emissions. The main fumarole in Pisciarelli is associated with a boiling mud pool. Recently, episodes of a sudden increase in hydrothermal activity and expansion of mud and gas emissions occurred in this area. During these episodes, which occurred in December 2018 and September 2020, Short Duration Events (SDEs), related to the intensity of mud pool boiling, were recorded in the fumarolic seismic tremor. We applied a Self-Organizing Map (SOM) neural network to recognize the occurrence of SDEs in the fumarolic tremor of Campi Flegrei, which provides important information on the state of activity of the hydrothermal system and about the possible phreatic activity. Our method, based on an ad hoc feature extraction procedure, effectively clustered the seismic signals containing SDEs and separated them from those representing the normal fumarolic tremor. This result is useful for improving the monitoring of the Solfatara–Pisciarelli hydrothermal area which is a high-risk zone in Campi Flegrei.

Evolution of the magmatic-hydrothermal system at Lastarria volcano (Northern Chile) between 2006 and 2019: Insights from fluid geochemistry

One of the major problems in the volcanic surveillance is how data from several techniques can be correlated and used to discriminate between possible precursors of volcanic eruptions and changes related to non-eruptive processes. Gas chemical surveys and measurements of SO2 emission rates performed in the past (2006–2019) at Lastarria volcano in Northern Chile have revealed a persistent increment of magmatic sourced gas emissions since late November 2012, following a 13 years period of intense ground uplift. In this work, we provide new insights into the gas-chemical evolution of Lastarria’s fumarolic discharges obtained from direct sampling (2006–2019) and SO2 emission rates using UV camera and DOAS instruments (2018–2019) and link these to pre-existing information on ground deformation (1998–2016) in order to determine the origin of observed degassing and ground deformation processes. We revise the four mechanisms originally proposed as alternatives by Lopez et al. (Geosphere, 2018, 14 (3), 983–1007) to explain the changes observed in the fluid geochemistry and ground deformation between 2009 and 2012, in order to explain major changes in gas-geochemistry over an extended period between 1998 and 2019. We hypothesize that a continuous sequence of processes explains the evolution in the fluid geochemistry of fumarolic discharges. Two mechanisms are responsible of the changes in the gas composition during the studied period, corresponding to a 1) deep magma chamber (7–15 km depth) pressurized by volatile exsolution (1998–2020), which is responsible of the large-scale deformation; followed by 2) a crystallization-induced degassing (2001–2020) and pressurization of the hydrothermal system (2003-early November 2012), where the former process induced the changes in the gas composition from hydrothermal-dominated to magmatic-dominated, whereas the last produced the small-scale deformation at Lastarria volcano. The changes in the gas composition since late November 2012, which were strongly dominated by magmatic volatiles, produced two consecutive processes: 1) acidification (late November 2012–2020) and 2) depletion (2019–2020) of the hydrothermal system. In this work we have shown that a long-term surveillance of the chemistry of fluid discharges provides valuable insights into underlying magmatic/volcanic processes, and consequently, for forecasting future eruptions.

Modelling of geomechanical response for coal and ground induced by CO2-ECBM recovery

Carbon dioxide enhanced coalbed methane (CO2-ECBM) recovery is widely regarded as a premium choice for improving CH4 recovery and sequestrating CO2 in subterranean seams. Some geomechanical responses, such as the coal softening, the variety of reservoir stress state and the subsidence and uplift of ground, are induced by CO2-ECBM recovery. In this paper, a set of linear equations is established to characterize the impact of adsorptive gas and water in reservoir on mechanical property for coal. Then a fully coupled mathematical model comprised of the mechanical, hydraulic and thermal fields for coal reservoir as well as the mechanical and thermal fields for its overlying strata is constructed to estimate the geomechanical response of CO2-ECBM recovery. Based on the geometric model including coal reservoir and its overlying strata in Qinshui Basin, the changes of mechanical property for coal, effective stress for reservoir and vertical displacement for strata are evaluated with CO2 injection rates of 6, 8 and 10 t/d by numerical simulation. The simulation results indicate that: (1) Elastic modulus diminishes while Poisson's ratio raises because the softening effect induced by gas ad-desorption is stronger than that of water for the coal with large adsorption capacity and low water content, which generally weakens the ability of coal to resist shearing and compressing, (2) Inpouring CO2 into coal reservoir decelerates the compressive stress on the coal skeleton, and the risk of coal failure increases with the growth of CO2 injection rate, (3) Vertical displacement of coal reservoir top outstrips that of ground, which has little impact on ground buildings in the study area with CO2 injection rate not exceeding 10 t/d for 8000 d. This work analyzes the geomechanical response and safety for CO2-ECBM recovery and can provides guidance for the large-scale CO2-ECBM project.

Discriminating carbon dioxide sources during volcanic unrest: The case of Campi Flegrei caldera (Italy)

AbstractLarge calderas are among the main emitters of volcanic CO2, which is mainly supplied by the deep degassing of magmatic fluids. However, other sources of non-magmatic CO2 can also occur due to the intense interaction among magmatic fluids, wide hydrothermal systems, and their host rocks. In particular, massive amounts of CO2 are released by calderas during unrest phases and have been often detected before eruptions. An accurate assessment of CO2 sources is thus fundamental to properly understand gas monitoring signals during volcanic crises. We focused on the restless Campi Flegrei caldera, in southern Italy, where CO2 fluxes at the Solfatara-Pisciarelli hydrothermal site have been progressively increasing up to 4000–5000 t/d during the ongoing unrest that started in 2005. Theoretical models of magma degassing have been able to reproduce the CO2-N2-He variations at the Solfatara fumaroles. However, a time-dependent deviation between measured and modeled N2/CO2 and He/CO2, well correlated with the temporal evolution of ground uplift and temperature of the hydrothermal system, has been observed since 2005. We show that these variations are controlled by intense physical-chemical perturbation of the hydrothermal system, which is driving the decarbonation of hydrothermal calcite stored in reservoir rocks. This process is providing large volumes of non-magmatic CO2 during the current unrest, contributing up to 20%–40% of the total fumarolic CO2.

Space distribution of the seismic source energy between 1-1-2000 and 3-7-2022 at Campi Flegrei caldera

The spatial distribution of seismic sources is commonly visualized through location maps, where a symbol (commonly a dot or a cross) shows where the earthquake source is located inside the area under study. In the present note I propose a different visualization approach, based on the isoline representation of the Energy Space Density (ESD), which is conceptually similar to approaches already used in literature to represent the pattern of space distribution of strong earthquakes or to spatially identify earthquake swarms. Visualizing such a pattern jointly with the hypocentral distribution map, favors the conceptual association of the Seismic Energy Release to the known geological structures in the area. A test is carried out using the Campi Flegrei seismic catalogue between 2000 and, 2022, period in which the seismic activity in the area resumed after a long calm of several years. The catalogue reports >4000 seismic events with magnitude ranging from 0.5 to 3.5, all located inside a very small volume (5 × 5 × 5 km) and associated with the phenomenon called ”bradiseism”, a slow but intense ground uplift (>80 cm in the 22-year interval spanned by the catalogue). Results indicate that the ESD pattern enhances the visualization of the rock fracture density associated with the ground deformation, in agreement with other studies in the same area, producing a picture of the most fractured rock volumes, possibly associated with fluid storage under the active areas of Solfatara and Pisciarelli.