Manuscript of Karen Strehlow is now online in the discussion forum of Solid Earth
http://www.solid-earth-discuss.net/7/1673/2015/sed-7-1673-2015.html
Abstract. Well water level changes associated with magmatic unrest can be interpreted as a result of pore pressure changes in the aquifer due to crustal deformation, and so could provide constraints on the subsurface processes causing this strain. We use Finite Element Analysis to demonstrate the response of aquifers to volumetric strain induced by pressurised magma reservoirs. Two different aquifers are invoked – an unconsolidated pyroclastic deposit and a vesicular lava flow – and embedded in an impermeable crust, overlying a magma chamber. The time-dependent, fully coupled models simulate crustal deformation accompanying chamber pressurisation and the resulting hydraulic head changes as well as porous flow in the aquifer. The simulated deformational strain leads to centimetres (pyroclastic aquifer) to meters (lava flow aquifer) of hydraulic head changes; both strain and hydraulic head change with time due to substantial porous flow in the hydrological system. Well level changes are particularly sensitive to chamber volume and shape, followed by chamber depth and the phase of the pore fluid. The Young's Modulus and permeability of the aquifer, as well as the strength of pressurisation also have significant influence on the hydraulic head signal. While source characteristics, the distance between chamber and aquifer and the elastic stratigraphy determine the strain field and its partitioning, flow and coupling parameters define how the aquifer responds to this strain and how signals change with time. We investigated a period of pre-eruptive head changes recorded at Usu volcano, Japan, where well data were interpreted using an analytical deformation model. We find that generic analytical models can fail to capture the complex pre-eruptive subsurface mechanics leading to well level changes, due to aquifer pressure changes being sensitive to chamber shape and lithological heterogeneities. In addition, the presence of a pore fluid and its flow have a significant influence on the strain signal in the aquifer and are commonly neglected in analytical models. These findings highlight the need for numerical models for the interpretation of observed well level signals. However, simulated water table changes do mirror volumetric strain and wells can therefore serve as comparatively cheap strain meters that could provide important insights into pre-eruptive dynamics.