Magnetic anisotropy produced by magma flow: theoretical model and experimental data from Ferrar dolerite sills (Antarctica)

Volcanic rocks forming sills, dykes or lava flows may display a magnetic anisotropy derived from the viscous flow during their emplacement. We model a sill as a steady-state flow of a Bingham fluid, driven by a pressure gradient in a horizontal conduit. The magma velocity as a function of depth is calculated from the motion and constitutive equations. Vorticity and strain rate are determined for a reference system moving with the fluid. The angular velocity and the orientation of an ellipsoidal magnetic grain immersed in the fluid are calculated as functions of time or strain. Magnetic susceptibility is then calculated for a large number of grains with a uniform distribution of initial orientations. It is shown that the magnetic lineation oscillates in the vertical plane through the magma flow direction, and that the magnetic foliation plane changes periodically from horizontal to vertical. The results are compared with the magnetic fabric of Ferrar dolerite sills (Victoria Land, East Antarctica) derived from low-field susceptibility measurements.     

Quasi-static deformation of a poroelastic half-space with anisotropic permeability by two-dimensional surface loads

An analytical solution is obtained of the fully coupled diffusion–deformation system of equations governing the quasi-static plane strain deformation of a poroelastic half-space with anisotropic permeability and compressible constituents. The stresses and the pore pressure are taken as the basic state variables. Displacements are obtained by integrating the coupled constitutive relations. The problem of surface loads is discussed in detail. Explicit analytical solutions are derived for normal line loading, shear line loading and normal strip loading. The permeability anisotropy is found to have a significant effect on the quasi-static deformation of the half-space. However, in the drained and undrained limits, the anisotropy has no effect. The stresses in the drained and undrained states are independent of the poroelastic parameters. Numerical computation of the pore pressure indicates that ignoring permeability anisotropy may lead to an overestimation of the pore pressure at points vertically below the point of normal loading. Further, anisotropy in permeability may lead to a dilution in the theoretical prediction of the Mandel–Cryer Effect.     

Ground deformation in a viscoelastic medium composed of a layer overlying a half-space: a comparison between point and extended sources

We obtain and compare analytical and numerical solutions for ground displacement caused by an overpressurized magma chamber placed in a linear viscoelastic medium composed of a layer over a half-space. Different parameters such as size, depth and shape of the chamber, crustal rheology and topography are considered and discussed. Numerical solutions for an axisymmetric extended source are computed using a finite element method (FEM). Analytical solutions for a point source are obtained using the dislocation theory and the propagator matrix technique. In both cases, the elastic solutions are used together with the correspondence principle of linear viscoelasticity to obtain the solution in the Laplace transform domain. Viscoelastic solutions in the time domain are derived inverting the Laplace transform using the Prony series method. The differences between the results allow us to constrain the applicability of the point source and the flat surface hypothesis, which are usually implicitly assumed when analytical solutions are derived. The effect of the topography is also considered. The results obtained show that neglecting the topographic effects may, in some cases, introduce an error greater than that implicit in the point-source hypothesis. Therefore, for an adequate modelling and interpretation of the time-dependent displacements, topography must be considered.