No abstract

Moving beyond the widely used kinematic models for the deformation sources, we present a new dynamic model to describe the process of injecting magma into an existing magma reservoir. To validate this model, we derive an analytical solution and compare its results to those calculated using the Finite Element Method. A Newtonian fluid characterized by its viscosity, density, and overpressure (relative to the lithostatic value) flows through a vertical conduit, intruding into a reservoir embedded in an elastic domain, leading to an increase in reservoir pressure and time‐dependent surface deformation. We apply our injection model to Interferometric Synthetic Aperture Radar (InSAR) data from the ongoing unrest episode at Laguna del Maule (Chile) volcanic field that started in 2007. Using a grid search optimization, we minimize the misfit to the InSAR displacement data and vary the three parameters governing the analytical solution: the characteristic timescale τP for magma propagation, the maximum injection pressure, and the inflection time when the acceleration switches from positive to negative. For a spheroid with semimajor axis a = 6200 m, semiminor axis c = 100 m, located at a depth of 4.5 km in a purely elastic half‐space, the best fit to the InSAR displacement data occurs for τP=9.5 years and an injection pressure rising up to 11.5 MPa for 2 years. The volume flow rate increased to 1.2 m3/s for 2 years and then decreased to 0.7 m3/s in 2014. In 7.3 years, at least 187 × 106 m3 of magma was injected.

Mid-ocean ridge morphology and crustal accretion are known to depend on the spreading rate of the ridge. Slow-spreading mid-ocean-ridge segments exhibit significant crustal thinning towards transform and non-transform offsets, which is thought to arise from a three-dimensional process of buoyant mantle upwelling and melt migration focused beneath the centres of ridge segments. In contrast, fast-spreading mid-ocean ridges are characterized by smaller, segment-scale variations in crustal thickness, which reflect more uniform mantle upwelling beneath the ridge axis. Here we present a systematic study of the residual mantle Bouguer gravity anomaly of 19 oceanic transform faults that reveals a strong correlation between gravity signature and spreading rate. Previous studies have shown that slow-slipping transform faults are marked by more positive gravity anomalies than their adjacent ridge segments, but our analysis reveals that intermediate and fast-slipping transform faults exhibit more negative gravity anomalies than their adjacent ridge segments. This finding indicates that there is a mass deficit at intermediate- and fast-slipping transform faults, which could reflect increased rock porosity, serpentinization of mantle peridotite, and/or crustal thickening. The most negative anomalies correspond to topographic highs flanking the transform faults, rather than to transform troughs (where deformation is probably focused and porosity and alteration are expected to be greatest), indicating that crustal thickening could be an important contributor to the negative gravity anomalies observed. This finding in turn suggests that three-dimensional magma accretion may occur near intermediate- and fast-slipping transform faults.

No abstract