An Approximate Approach to Nonisothermal Emplacement of Kilometer‐Sized Kilometer‐Deep Sills at Calderas

Caldera unrest is often caused by kilometer‐sized kilometer‐deep sills. Still unanswered questions include the following: How do sills spread? Why can magma propagate for kilometers without solidifying? Why do ground deformation data rarely, if ever, detect sill propagation? We show that kilometer‐sized kilometer‐deep magmatic sills spread like hydraulic fractures in an infinite medium. How magma propagates depends on overburden pressure, magma viscosity, injection rate, and difference between magma and rock temperatures. A small lag, filled with vapors from the fluid and/or the rock, exists between the propagating magma and fracture fronts. If the sill spreads along an interface, the lag slightly affects isothermal sill spreading but takes a key role in the case of nonisothermal propagation: A sill would stop after few tens of meters without it, unless magma intrudes rocks that are as hot as the solidification temperature or has unrealistic overpressures, because spreading velocity decreases soon to the critical value at which the tip becomes blocked with solidified magma. The lag defers magma solidification as heat exchange between the magma and the rock is effective only behind the thermal‐insulating lag, where magma has some finite thickness and sill opening grows with distance from the tip faster than thickness of solidified magma. Thus, the critical velocity decreases, allowing greater maximum sill sizes. We also show that the ground deformation pattern does not change appreciably over time if the final sill radius is smaller than 2 to 3 km, explaining why deformation is usually attributed to the inflation of a stationary source.

[1]  M. Calò,et al.  Anatomy of the Campi Flegrei caldera using Enhanced Seismic Tomography Models , 2018, Scientific Reports.

[2]  E. Detournay,et al.  The Tip Region of a Near-Surface Hydraulic Fracture , 2018 .

[3]  L. Crescentini,et al.  Thermally-assisted Magma Emplacement Explains Restless Calderas , 2017, Scientific Reports.

[4]  Andrew P. Bunger,et al.  Laboratory measurement of tip and global behavior for zero-toughness hydraulic fractures with circular and blade-shaped (PKN) geometry , 2017 .

[5]  G. Wörner,et al.  Timescales of magmatic processes prior to the ∼4.7 ka Agnano-Monte Spina eruption (Campi Flegrei caldera, Southern Italy) based on diffusion chronometry from sanidine phenocrysts , 2017, Bulletin of Volcanology.

[6]  E. Dontsov An approximate solution for a penny-shaped hydraulic fracture that accounts for fracture toughness, fluid viscosity and leak-off , 2016, Royal Society Open Science.

[7]  C. Michaut,et al.  Elastic-plated gravity currents with a temperature-dependent viscosity , 2016, Journal of Fluid Mechanics.

[8]  Giuseppe Aiello,et al.  Magma transfer at Campi Flegrei caldera (Italy) before the 1538 AD eruption , 2016, Scientific Reports.

[9]  R. Moretti,et al.  Open-system magma evolution and fluid transfer at Campi Flegrei caldera (Southern Italy) during the past 5 ka as revealed by geochemical and isotopic data: The example of the Nisida eruption , 2016 .

[10]  E. Dontsov Propagation regimes of buoyancy-driven hydraulic fractures with solidification , 2016, Journal of Fluid Mechanics.

[11]  Emmanuel M Detournay,et al.  Mechanics of Hydraulic Fractures , 2016 .

[12]  A. Peirce,et al.  A non-singular integral equation formulation to analyse multiscale behaviour in semi-infinite hydraulic fractures , 2015, Journal of Fluid Mechanics.

[13]  Valerio Acocella,et al.  An overview of recent (1988 to 2014) caldera unrest: Knowledge and perspectives , 2015 .

[14]  Riccardo Lanari,et al.  Magma injection beneath the urban area of Naples: a new mechanism for the 2012–2013 volcanic unrest at Campi Flegrei caldera , 2015, Scientific Reports.

[15]  Luca D'Auria,et al.  Evidence of thermal-driven processes triggering the 2005–2014 unrest at Campi Flegrei caldera , 2015 .

[16]  R. Katz,et al.  A review of mechanical models of dike propagation: Schools of thought, results and future directions , 2015 .

[17]  Prospero De Martino,et al.  Clues to the cause of the 2011–2013 Campi Flegrei caldera unrest, Italy, from continuous GPS data , 2014 .

[18]  F. Giudicepietro,et al.  Sill intrusion as a source mechanism of unrest at volcanic calderas , 2014 .

[19]  L. Crescentini,et al.  Paired deformation sources of the Campi Flegrei caldera (Italy) required by recent (1980–2010) deformation history , 2014 .

[20]  J. Neufeld,et al.  Viscous control of peeling an elastic sheet by bending and pulling. , 2013, Physical review letters.

[21]  C. Freda,et al.  Clinopyroxene–liquid thermometers and barometers specific to alkaline differentiated magmas , 2013, Contributions to Mineralogy and Petrology.

[22]  Agust Gudmundsson Magma chambers: Formation, local stresses, excess pressures, and compartments , 2012 .

[23]  P. Scarlato,et al.  A general viscosity model of Campi Flegrei (Italy) melts , 2011 .

[24]  C. Michaut Dynamics of magmatic intrusions in the upper crust: Theory and applications to laccoliths on Earth and the Moon , 2011 .

[25]  I. Arienzo,et al.  The magmatic feeding system of the Campi Flegrei caldera: Architecture and temporal evolution , 2011 .

[26]  P. Papale,et al.  The feeding system of Agnano-Monte Spina eruption (Campi Flegrei, Italy): Dragging the past into present activity and future scenarios , 2010 .

[27]  Agust Gudmundsson Toughness and failure of volcanic edifices , 2009 .

[28]  Emmanuel M Detournay,et al.  Experimental validation of the tip asymptotics for a fluid-driven crack , 2008 .

[29]  L. Crescentini,et al.  Simultaneous inversion of deformation and gravity changes in a horizontally layered half-space: Evidences for magma intrusion during the 1982–1984 unrest at Campi Flegrei caldera (Italy) , 2008 .

[30]  A. Peirce,et al.  An implicit level set method for modeling hydraulically driven fractures , 2008 .

[31]  P. Segall,et al.  Magma compressibility and the missing source for some dike intrusions , 2008 .

[32]  L. Crescentini,et al.  Effects of crustal layering on the inversion of deformation and gravity data in volcanic areas: An application to the Campi Flegrei caldera, Italy , 2007 .

[33]  E. Detournay,et al.  Early-Time Solution for a Radial Hydraulic Fracture , 2007 .

[34]  W. Chadwick,et al.  Results from new gps and gravity monitoring networks at fernandina and sierra negra volcanoes , 2006 .

[35]  N. Balmforth,et al.  Elastic-plated gravity currents , 2005, European Journal of Applied Mathematics.

[36]  H. Behrens,et al.  Water solubility in trachytic melts , 2004 .

[37]  L. Crescentini,et al.  Effects of crustal layering on source parameter inversion from coseismic geodetic data , 2004 .

[38]  G. Orsi,et al.  Volcanic hazard assessment at the restless Campi Flegrei caldera , 2004 .

[39]  Emmanuel M Detournay,et al.  Propagation of a penny-shaped fluid-driven fracture in an impermeable rock: Asymptotic solutions , 2002 .

[40]  E. Ventsel,et al.  Thin Plates and Shells: Theory: Analysis, and Applications , 2001 .

[41]  Mark Simons,et al.  Deformation due to a pressurized horizontal circular crack in an elastic half-space, with applications to volcano geodesy , 2001 .

[42]  Emmanuel M Detournay,et al.  The Tip Region of a Fluid-Driven Fracture in an Elastic Medium , 2000 .

[43]  J. Lister,et al.  The effect of solidification on fluid–driven fracture, with application to bladed dykes , 1999, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[44]  A. Rubin On the thermal viability of dikes leaving magma chambers , 1993 .

[45]  R. C. Kerr,et al.  Fluid‐mechanical models of crack propagation and their application to magma transport in dykes , 1991 .

[46]  A. McBirney Volcanic processes. , 1982, Science.

[47]  Arvid M. Johnson,et al.  Mechanics of growth of some laccolithic intrusions in the Henry mountains, Utah, II: Bending and failure of overburden layers and sill formation , 1973 .

[48]  J. D. Eshelby The determination of the elastic field of an ellipsoidal inclusion, and related problems , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[49]  Valerio Acocella,et al.  Rock Fractures in Geological Processes: Index , 2011 .

[50]  Agust Gudmundsson,et al.  Rock Fractures in Geological Processes: Index , 2011 .

[51]  G. Rolandi,et al.  Ground deformation at Campi Flegrei, Italy: implications for hazard assessment , 2006, Geological Society, London, Special Publications.

[52]  F. Dobran Volcanic Processes: Mechanisms in Material Transport , 2001 .

[53]  John R. Rice,et al.  Mathematical analysis in the mechanics of fracture , 1968 .

[54]  Michael P. Poland,et al.  Chapter 5 Magma Supply , Storage , and Transport at Shield-Stage Hawaiian Volcanoes , 2022 .