Source Of Deformation Research Articles

Sensing expectation enables simultaneous proprioception and contact detection in an intelligent soft continuum robot

A high-level perceptual model found in the human brain is essential to guide robotic control when facing perception-intensive interactive tasks. Soft robots with inherent softness may benefit from such mechanisms when interacting with their surroundings. Here, we propose an expected-actual perception-action loop and demonstrate the model on a sensorized soft continuum robot. By sensing and matching expected and actual shape (1.4% estimation error on average), at each perception loop, our robot system rapidly (detection within 0.4 s) and robustly detects contact and distinguishes deformation sources, whether external and internal actions are applied separately or simultaneously. We also show that our soft arm can accurately perceive contact direction in both static and dynamic configurations (error below 10°), even in interactive environments without vision. The potential of our method are demonstrated in two experimental scenarios: learning to autonomously navigate by touching the walls, and teaching and repeating desired configurations of position and force through interaction with human operators.

Precursory crater contraction associated with the 2017 eruption of Shinmoe-dake volcano (Japan) detected by PALSAR-2 and Sentinel-1 InSAR

The time series of PALSAR-2 and Sentinel-1 images reveal inflation at the volcanic flank and contraction at the crater for approximately 5 months before the 2017 eruption of Shinmoe-dake volcano, Japan. While the observation of inflation at the volcano’s flank is ubiquitous, few studies have reported crater contraction preceding an eruption. The flank inflation stopped after the 2017 eruption, while the contraction at the crater continued until the 2018 eruption. We found that a pipe-shaped deformation source above sea level best fits the observation preceding the 2017 eruption. Suppose the flux of ejected materials constrains the conduit radius during the previous 2011 eruption. In that case, the amount of deformation of the pipe-shaped deformation source, whether open or closed at its top, is too large to be realistic. Although constraining the conduit radius from the eruption flux overestimates the pressure change of the pipe-shaped deformation source, water-saturated fractures along the volcanic conduit could extend the effective conduit radius of the pressure source. We propose one potential scenario for the mechanism of the crater contraction preceding volcanic eruptions based on the combination of compaction due to cooling by ambient groundwater and material withdrawal within the conduit. The groundwater inflows from the ambient aquifer through cracks in the porous conduit wall, which are generated by conduit expansion during the magma ascent. Decoupling from the conduit wall due to a decrease in volume of the material promotes material instability and crater contraction. The interaction between the groundwater and the magma triggers the 2017 eruption of Shinmoe-dake volcano, as previous studies have reported.Graphical

A Century of Deformation and Stress Change on Kīlauea's Décollement

AbstractKīlauea Volcano on Hawai'i Island is host to a complex volcanic and interwoven fault system. Over the last ∼120 years, a range of seismic events, including large earthquakes such as the 1975 7.7 Kalapana earthquake, creep, and slow slip events, have occurred along the décollement underlying Kilauea's south flank. We explore both the deformation and stress changes of Kīlauea from 1896 to 2018 by collating six geodetic data sets and creating an analytical model to determine the dominant deformation sources (i.e., fault planes, rifts, magma chambers) driving this system at different times. The 1975 Kalapana earthquake significantly altered the region's state of stress and deformation; we find the average slip along the décollement was reduced from 10 cm/yr prior, to 4 cm/yr after the rupture. Prior to 1975 no slip is resolved along the décollement where the earthquake nucleated, suggesting that this portion may have been locked leading up to the rupture. After 1975, décollement slip overall is smaller and more irregular, suggesting increased control by spatial variation of mechanical properties. We find increases in shear stress along the Kīlauea décollement and a decrease in normal compressive stress within the East Rift Zone prior to the Kalapana earthquake, creating favorable conditions for failure of the décollement and subsequent magmatic intrusion.

Crustal Deformation and Seismic Velocity Perturbations in the Alto Tiberina Fault Zone (Northern Apennines, Italy)

AbstractCrustal perturbations related to seismic activity can generally be observed with the occurrence of a large magnitude event. For less energetic seismic sequences though, the associated transient crustal variations are questionably measurable, and their observation gets easily obscured by relatively stronger perturbations such as the ones related to hydrological processes. In this study we reveal the significant role that terrestrial water‐storage variations play in governing temporal crustal changes in the tectonically active Northern Apennines of Italy, and discuss the potential of accounting for its correction in order to monitor the relatively weaker transient perturbations caused by local seismic swarms. This area is characterized by an extensive level of low‐energetic seismic activity, typically clustered in time and space, of which three main seismic swarms outstand during the 12 year period of study (2010–2021). Our analysis compares independent observations and processing methods of Global Navigation Satellite System measurements and ambient seismic noise recordings. We adopt a multivariate statistical approach to discriminate between independent sources of ground deformation, and seismic noise cross‐correlation analysis to monitor relative seismic‐velocity variations. The result shows how the perturbation effects produced by variations in total water content are dominant in both time series of ground deformations and seismic‐velocity variations. After correcting for the water‐related variation effects, our monitoring results reveal perturbations in the crustal properties whose activation time and depth range correlate with the occurrences of the seismic swarms.

Anisotropy in tensile properties of a high strength metastable β titanium alloy

High strength metastable β titanium alloys are widely employed in the aircraft industry due to their outstanding strength-to-weight ratio. While components can endure complex in-service mechanical loading, the anisotropy in tensile properties has been the subject of limited attention. In this study, its origin was investigated focusing on the role played by millimeter scale β grains as they were recently identified as a source of heterogeneous deformation. Tensile properties of Ti-10V-2Fe-3Al processed via different thermomechanical routes were assessed using multiple sampling directions. In particular, elongation values were observed to vary significantly depending on the testing direction. A combination of SEM, EBSD, µ-CT and in-situ DIC during tensile tests was employed to clarify the underlying causes of this behavior. Substantial differences in strain heterogeneity and localization were found related to features of β grains, including their crystallographic and morphologic orientations. Furthermore, multiple fracture mechanisms were observed to derive from the differences in deformation behavior, and eventually compete to trigger specimen failure. Elongation values are then determined by both the degree of strain heterogeneity and the operating fracture mechanisms. These findings provide a new understanding of the role of the microstructure in the tensile behavior of high strength metastable β titanium alloys.

Tracking the 2007–2023 magma-driven unrest at Campi Flegrei caldera (Italy)

Understanding and managing unrest at a volcano include i) ascertaining the magmatic distribution and migration, and ii) tracking the evolution of the shallow plumbing system. Here we use multi-technique geodetic data, mechanical models, and petrological simulations to define both aspects for the ongoing (2005-present) unrest at Campi Flegrei caldera, Italy. Results show a deformation source exhibiting progressive widening and shallowing, from 5.9 to 3.9 kilometres. Concurrently, a deeper tabular source at 8 km depth experiences limited but constant deflation. Petrological calculations explain inflation of the shallower source resulting from the rise of 0.06 to 0.22 cubic kilometres of magma from depth ≥8 kilometres. Our analysis provides strong evidence that magma ascent to depths shallower than 8 kilometres is the ultimate driver behind the ongoing unrest. This merging of geodetic and petrological approaches to track the evolution of a plumbing system better constrains magma ascent at volcanoes experiencing unrest.

Evaluation of recent tectonic movement in northeast Japan by using long-term GNSS and tide gauge measurements

The study investigates tectonic movements in northeast Japan by using long-term (2000–2022) global navigation satellite system (GNSS) and tide gauge measurements. The effect of the 2011 Tohoku-Oki earthquake including the other eight seismic events that occurred within this period is also discussed using GNSS time-series. The result showed the break in GNSS-time series because of occurred earthquakes and pointed out tectonic movements significantly. The GNSS sites located in the central and southern parts of northeast Japan showed that the velocity vectors have strong internal variation and suggest the existence source of alternative deformation because of geological terranes within the region. A least square approach was used, and the trend of sea-level measurements was fitted with the straight line. The obtained results from tide gauge measurements showed a rising trend at almost every site and indicated lithospheric uprising movement because of tectonic activities. This is possible because of the ongoing subduction of the Pacific and Philippine Sea plates beneath the Eurasian and North American plates.

THEPORE: A software package for modeling THErmo-PORo-elastic displacements

THEPORE (THErmo-POro-Elastic solutions) is an open source software to perform forward and inverse modeling of ground displacements induced by thermo-poro-elastic sources. The software, implemented in MATLAB, offers a library of analytical and semi-analytical solutions to compute ground displacements induced by thermo-poro-elastic deformation sources of different geometries, embedded in an elastic, homogeneous and isotropic half-space. The solutions have been verified against finite-element simulations. THEPORE includes also an inversion procedure of the deformation data to constrain the source parameters that better fit the observed signals.The software's functionality is showcased by inverting the GPS deformation data recorded on Vulcano Island at the onset of the 2021 unrest, in order to estimate the position and volume change of the source responsible for the observed deformations. The results encourage to consider THEPORE as a practical tool suitable for a fast preliminary estimation of the deformation source during a volcanic crisis.

Non‐isometric deformation in double curvature and symmetric composite laminates due to shrinkage

AbstractApplication of the global Gauss–Bonnet theorem for the prediction of changes in Gaussian curvature for double curvature laminate geometries due to shrinkage is developed for non‐isometric deformation states where in surface area change accompanies changes in curvature. The transition between isometric and non‐isometric deformation is strongly influenced by the ratio of the principal radii of curvature of the double curved surface and a relationship is proposed by the authors to accommodate the transition. Examples for the radii of curvature ratios of 2, 3, and 4 are compared to results from finite‐element simulations and the developed relationship is shown to be accurate within 0.55%. Furthermore, these results demonstrate that isometric deformation state can be taken as an adequate representation for double curvature toroidal laminate geometries of >10 and the class of materials studied in this work.Highlights Shrinkage induced deformation in multi‐curvature composite laminates is a recurring issue in polymer composite manufacturing. Thermal, crystallization and cure shrinkage are the dominant sources of shrinkage deformation in polymer composite systems. Isometric deformation of multi‐curvature, symmetric laminates composite is defined by the absence of membrane strain. The Gauss–Bonnet theorem provides a platform for determining radii of curvature changes in multi‐curvature symmetric laminates, Non‐isometric deformation of symmetric laminates can be treated by the Gauss–Bonnet theorem for radii of curvature ratios that correspond to geometric self‐constraint. The above highlights yield insight into the anticipated geometric distortions due to shrinkage in polymer composites.

Effect of Machine Structure on Cutting Force Considering Different Working Conditions in Micro-End-Milling

Machine structure determines the performance of machine tools. The study of machine deformation contributes to the stability of machining processes, which in turn improves machining accuracy and surface quality. Vertical structure is the weak part of the horizontal micro-milling machine, which is the main source of deformation. Considering that the micro-milling cutting parameters and machine deformations are at the micron scale, it remains to be investigated whether such deformations cause changes in the cutting volume and thus affect the cutting forces. For this purpose, this paper investigates the deformation of the machine structure caused by various operating conditions and its effect on the cutting force. The finite element method (FEM) was used to calculate the deformation [Formula: see text][Formula: see text] of the clamp-workpiece. To verify the effect of machine structure on cutting force under different working conditions, the Al 6061 full immersion slot cutting experiments were conducted by considering different cutting parameters. The results show that the deformation of the machine structure has an effect on the cutting forces, and the degree of effect varies depending on the working conditions.

Mapping geodetically inferred Antarctic ice surface height changes into thickness changes: a sensitivity study

Abstract. Determining recent Antarctic ice volume changes from satellite altimeter measurements of ice surface height requires a correction for contemporaneous vertical crustal deformation. This correction must consider two main sources of crustal deformation: (1) ongoing glacial isostatic adjustment (GIA), that is, the deformational, gravitational, and rotational response to Late Pleistocene and Holocene ice and ocean mass changes and (2) modern ice mass change. In this study, we seek to quantify the uncertainties associated with each of these corrections. Corrections of ice surface height changes for correction 1 have generally involved the adoption of global models of GIA defined by some preferred combination of ice history and mantle viscoelastic structure. We have computed the GIA correction generated from a coupled ice sheet–sea level model and a realistic Earth model incorporating three-dimensional viscoelastic structure. Integrating the difference between this correction and those from recent GIA analyses widely adopted in the literature yields an uncertainty in total present-day ice volume change equivalent to approximately 10 % of Antarctic ice mass loss inferred for the period 2010–2020. This reinforces earlier work indicating that ice histories characterized by relatively high excess ice volume at the Last Glacial Maximum may be introducing a significant error in estimates of modern melt rates. Regarding correction 2, a spatially invariant scaling has commonly been used to convert GIA-corrected ice surface height changes obtained from satellite altimetry to ice volume estimates. We adopt modeling results based on a projection of Antarctic ice mass change over the period 2015–2055 to demonstrate a spatial variability in the scaling of up to 10 % across the ice sheet. Furthermore, using these calculations, we find a systematic error of ∼ 3 % in the projected net ice volume change, with most of the difference arising in areas of West Antarctica above mantle zones of low viscosity.

Independent component analysis and finite element modelling of the 2004–2005 ground deformation in Tenerife (Canary Islands)

Historic volcanic activity in Tenerife was concentrated within two of the island’s three dorsals and on the Teide-Pico Viejo complex located inside Las Cañadas caldera. Eruptions on the island are primarily characterised by basaltic and trachybasaltic fissural eruptions. However, the Teide-Pico Viejo complex also hosted explosive and effusive phonolitic eruptions. Our study focused on the analysis of the 2004–2005 unrest in Tenerife, which was characterised by an onset of a ground deformation pattern, heightened on-land seismic activity, changes in the chemical composition of fumaroles of the Teide composite volcano, an increase in diffusive emissions of carbon dioxide along the NW rift, and, significant gravity changes. We used the Envisat-ASAR satellite images from 2003 to 2010 to generate the Line-Of-Sight SBAS-DInSAR deformation time series to investigate the source responsible for ground deformation. We applied the Independent Component Analysis (ICA) to separate distinct ground deformation patterns. Specifically, we selected four components for the Independent Component Analysis decomposition: the first one mainly affects the stratovolcano’s summit region of Teide and shows a circular symmetry; the second and third components are possibly related to the topography and atmospheric artifacts, while the fourth contains only a noisy signal. We employed a non-linear optimisation approach in a Finite Element modelling environment to determine the source geometry responsible for the first identified ICA pattern of ground deformation within Las Cañadas. Our results revealed the existence of an ellipsoidal ground deformation source oriented along the E-W axis, located beneath the Teide and Montaña Blanca volcanoes at 1,600 m a.s.l. This source became active during the seismic crisis of 2004–2005, which was associated with the degassing of a magma batch that intruded into the northwest rift of Tenerife island. We propose that the ground deformation observed at Teide volcano during the 2004–2005 crisis was related to hydrothermal activity within the volcano.

On the Initial Phase of the Ongoing Unrest at Campi Flegrei and Its Relation with Subsidence at Vesuvio (Italy)

The densely inhabited area of Naples (Italy), between the Campi Flegrei and Vesuvio volcanoes, is one of the most hazardous regions in the world. After two decades of sustained subsidence, Campi Flegrei has been experiencing an accelerating uplift since 2005. The uplift is currently associated with unusual seismicity and increased degassing. To try to identify the cause of the shift from subsidence to uplift and explore any connection between Campi Flegrei and Vesuvio, we analysed the ground displacement time series of the two volcanoes from 1993 to 2010, obtained from ERS/ENVISAT Synthetic Aperture Radar imagery. To distinguish between the various sources of deformation, we used simple scatter plots and a blind source separation technique called variational Bayesian independent component analysis (vbICA). We obtained consistent results using both approaches. Specifically, with vbICA, we identified two significant independent components (ICs). IC1 describes the subsidence that occurred at Campi Flegrei prior to 2000, including the mini-uplifts of 2000 and 2005, and part of the post-2005 uplift. The expansion and contraction of two volumes beneath Campi Flegrei satisfy IC1: a sill-shaped volume at a depth of approximately 3 km and a small volume at a depth of 1–2 km, respectively. The two sources of deformation reproduce the large-scale deformation in the Campi Flegrei area and the local deformation in the Solfatara area, respectively. In the Campi Flegrei area, IC2 exhibits primarily uplift, which is concentrated in the eastern part of the caldera. The deformation pattern is complex and difficult to interpret. If we model it using simple spheroidal deformation sources, the pattern suggests that two volumes at depths of approximately 9 and 8 km are experiencing opposite activity, namely contraction (beneath the southwestern part of the caldera) and expansion (beneath the central part of the caldera). In the Vesuvio area, IC2 is consistent with the deformation induced by the contraction of a volume at a depth of around 9 km. The contraction beneath Vesuvio is smaller in magnitude than the expansion/contraction beneath Campi Flegrei. The correlation observed after 2002 between uplift at Campi Flegrei and subsidence at Vesuvio suggests the transfer of magma and/or magmatic fluids between the two plumbing systems at 8–9 km depth. This implies that part of the ongoing unrest at Campi Flegrei may have been promoted by mass transfer from below Vesuvio.

Numerical experiment toward simultaneous estimation of the spatiotemporal evolution of interseismic fault slips and block motions in southwest Japan

Interseismic strain around southwest Japan can be attributed to slip phenomena on the plate interface and processes inside the continental plate, such as block motions and inland fault slips. Although many crustal deformation analyses have been carried out, the simultaneous estimation of these phenomena and comprehensive discussion regarding the effect on the total strain budget remain topics to be investigated. In this study, we conducted numerical experiment to evaluate the possibility of simultaneously monitoring the spatiotemporal evolution of interseismic fault slips and block motions based on state-space modeling. We aimed to estimate fault slips and block motions using the time series of dense continuous Global Navigation Satellite System sites covering southwest Japan, encompassing 25 years from 1996 onward. We calculated synthetic block motions of the forearc region, monotonically increasing back slips on block boundary faults, and long-term slow slip events in the Tokai region, Kii Channel, and Bungo Channel. Subsequently, we generated the expected synthetic displacement time-series at GNSS stations. Applying Kalman filtering, we successfully estimated the spatiotemporal evolution of block motions, cumulative back slips, multiple slow slip events occurring in different regions. Although the contributions of fault slips and block motions showed slight trade-off, main characteristics of the slip distributions and the direction of block motions were well-recovered. We recovered slow slip events with slips of 5–10 cm or larger. We investigated the estimation uncertainty and separation precision of the unknown parameters using the covariance matrix estimated by Kalman filtering. Focusing on the structure of the non-diagonal component of the covariance matrix, we evaluated the complex effects of site and subfault locations on the estimated slip spatial distribution bias. For instance, the extent of the correlation between subfaults suggested that the three slow slip regions have different tendency of the uncertainties of slip areas, extending toward landside, seaward, or both of them. Our framework enables the comprehensive evaluation of the contributions and uncertainties of various deformation sources covering the entire subduction zone.Graphical abstract

The effect of powder recycling on the mechanical performance of laser powder bed fused stainless steel 316L

Powder recycling refers to the reuse of unused powder feedstock in the laser powder bed fusion (PBF-LB/M) process. This approach is crucial for the economic viability and sustainability of PBF-LB/M, as powder accounts for a large proportion of the total production cost. However, through powder recycling, the physical and chemical properties of powder are liable to change. This variation in powder properties can subsequently lead to knock-on effects on the mechanical properties of a fully built component.This research has investigated the changes that occur to stainless steel 316L (SS316L) powder as a result of recycling. This includes changes to powder size distribution (PSD), flowability, chemistry and phase composition. Likewise, the impact that these changes have will also be assessed in PBF-LB/M SS316L components manufactured from powders after different levels of recycling and subjected to alternative post processing routes such as hot isostatic pressing (HIP). This comprehensive investigation involves a thorough examination of both macro- and microstructures, encompassing detailed analyses of chemical composition, microstructural features, and defects. The study aims to elucidate differences in mechanical behaviour through a series of experiments, including uniaxial tensile tests, Charpy impact assessments, and low cycle fatigue (LCF) experiments. Additionally, the investigation will be complemented by pitting potential tests, providing a holistic understanding of the material's performance and characteristics.Although moderate changes to powder were observed for both PSD and chemistry, this was found to be negligible and not enough to result in any adverse changes to part performance. In addition, the microstructure of SS316L remained stable across differing levels of powder recycling. Whereas the porosity content increased marginally as the fine particle content of powder was reduced, this was not found to be sufficient to affect the LCF performance of the material. After powder recycling, increases in ductility and Young’s modulus were attributed to a reduction in oxides present in the microstructure, which were sources of localised damage and deformation.

Strain Partitioning and Fault Kinematics in the Northern Qilian Shan (NE Tibet) Determined From Bayesian Inference of Geodetic Data

AbstractOblique convergence across the northern Qilian Shan is accommodated by sub‐parallel strike‐slip and thrust faults that ruptured simultaneously in the Mw 8 Gulang earthquake in 1927. We investigate the kinematics of fault loading in the northern Qilian Shan and provide insights into the conditions necessary for generating multi‐fault earthquakes. We perform Bayesian inversions for the geometry and creep rate on the fault network. We infer that all of the thrust faults are locked north of the Qilian‐Haiyuan strike‐slip fault and are accumulating elastic strain. Multi‐fault earthquakes may occurr in this fault system because the faults are simultaneously loaded by the same source of deformation and are linked together by locked fault segments. The interseismic velocity field alone can not contain the location or activity of individual faults visible in the geomorphology, therefore the short‐term geodetic measurements may not reliably indicate the long‐term behavior of the fault system.

Form of salt rises as an indicator of interaction of geodynamic systems

The forms of salt uplifts are formed under the influence of stresses developing in the earth’s crust and, therefore, can be considered as sensitive indicators of the interaction of geodynamic systems. In the North German basin, in the Cis-Ural trough on the border with the Caspian depression and in the Pripyat trough, under the influ- ence of various sources of stress and deformation, forms of salt uplifts arose, reflecting the superposition of these systems. The mechanisms of such a superposition were both the formation of forms associated with the resulting addition of deformation vectors and the interference of forms associated with each geodynamic system separately.

Versatile Modeling Of Deformation (VMOD) Inversion Framework: Application to 20 Years of Observations at Westdahl Volcano and Fisher Caldera, Alaska, US

AbstractWe developed an open source, extensible Python‐based framework, that we call the Versatile Modeling Of Deformation (VMOD), for forward and inverse modeling of crustal deformation sources. VMOD abstracts from specific source model implementations, data types and inversion methods. We implement the most common geodetic source models which can be combined to model and analyze multi‐source deformation. VMOD supports Global Navigation Satellite System (GNSS), InSAR, electronic distance measurement, Leveling and tilt data. To infer source characteristics from observations, VMOD implements non‐linear least squares and Markov Chain Monte‐Carlo Bayesian inversions, including joint inversions using different sources of data. VMOD's structure allows for easy integration of new geodetic models, data types, and inversion strategies. We benchmark the forward models against other published results and the inversion approaches against other implementations. We apply VMOD to analyze deformation at Unimak Island, Alaska, observed with continuous and campaign GNSS, and ascending and descending InSAR time series generated from Sentinel‐1 satellite radar acquisitions. These data show an inflation pattern at Westdahl volcano and subsidence at Fisher Caldera. We use VMOD to test a range of source models by jointly inverting the GNSS and InSAR data sets. Our final model simultaneously constrains the parameters of two sources. Our results reveal a depressurizing spheroid under Fisher Caldera ∼4–6 km deep, contracting at a rate of ∼2–3 Mm3/yr, and a pressurizing spherical source underneath Westdahl volcano ∼6–8 km deep, inflating at ∼5 Mm3/yr. This and past applications of VMOD to volcanic unrest benefit from an extensible framework which supports jointly inversions of data sets for parameters of easily composable multi‐source models.

Dynamics of the 2021 unrest at Changbaishan Tianchi volcano from ALOS-2/PALSAR-2 and seismic data

The Changbaishan Tianchi intraplate volcano is one of the most active and hazardous volcanoes of NE Asia, characterized by a summit caldera formed after the 946 CE ‘Millennium’ Plinian eruption. From December 2020 to June 2021, the frequency and magnitude of earthquakes at Tianchi were significantly higher than during background periods, with hundreds of earthquakes (46 events per month in average) and reaching a local magnitude of ML 3.1. This study reports a comprehensive deformation analysis and geophysical inversion scheme aimed to unveil the dynamics of this period. Multi-temporal InSAR analysis results of 32 ALOS-2 images from 2018 to 2022 show that the surface deformation is a combination of seasonal fluctuations (± 25 mm in average, with a maximum ± 45 mm) and a long-term positive component. The least squares linear regression of the deformation time series and temperature data, isolates the seasonal fluctuations, revealing a clear uplift-subsidence process from June 2020 to July 2021 in the caldera area. To constrain the Tianchi plumbing system dynamics, a combined inversion scheme consisting of three deformation sources is designed. The inversion results and the seismic records indicate that Tianchi volcano has experienced a low-level unrest episode from December 2020 to June 2021. The shallower plumbing system, located at about 5–9 km depth and modelled by pressurized spheroids, underwent a cumulative volume increase of 26 × 106 m3 from November 2018 to April 2021, followed by a volume decrease of 9 × 106 m3 from April to July 2021. This suggests magma uprising from the 14 km deep storage zone to the shallower plumbing system, followed by depressurization of the plumbing system due to the escape of fluids. This research provides a comprehensive understanding of the magma and fluid migration dynamics within the Tianchi multi-level plumbing system for the first time.

Geodetic imaging of magma ascent through a bent and twisted dike during the Tajogaite eruption of 2021 (La Palma, Canary Islands)

On Sept. 19th, 2021, the largest historical eruption on the island of La Palma began, which had a significant scientific, social, and economic impact. The 2021 Tajogaite eruption was characterised by short precursors, lasting only 8 days. The seismicity started on Sept. 11th with a westward and upward migration of hypocenters. Permanent GNSS stations started recording deformation on Sept. 12th on the island's western side, which reached more than 15 cm just before the eruption. After the eruption onset, the ground deformation increased, reaching a maximum on Sept. 22nd and showing a nearly steady deflation trend in the following months. To better understand the dynamics of the eruption, we exploited a joint dataset of GNSS and Sentinel-1 SBAS time series along both ascending and descending orbits. To obtain the geometry of the causative source of the ground deformation, we combined the result of a preliminary non-linear inversion and the precise location of the seismicity. The resulting geometry of the source is that of a twisted dike bending eastward. We performed inverse modelling to obtain the spatiotemporal kinematics of the opening function of the dike. The forward modelling has been realised using a 3D finite-element approach considering the island's topography. Our findings reveal a close correspondence between the magmatic intrusion and pre-eruptive seismicity. The ascent of the magma occurred along two branches, and the rheology of a previously identified ductile layer strongly affected the magma propagation process. Finally, we found evidence of an early shallow deformation, which we interpret as the effect of ascending hydrothermal fluids. Our findings highlight the need for advanced modelling to understand pre-eruptive processes in basaltic volcanoes.