Piezomagnetic fields associated with a dislocation source in a layered elastic medium

Abstract

The piezomagnetic effect is defined as a change in magnetization with applied stress. Changes in the geomagnetic field caused by the piezomagnetic effect, referred to as the piezomagnetic field, have been theoretically estimated and compared by previous studies to interpret observed variations in the geomagnetic field. However, the piezomagnetic field estimated in previous studies may not provide an accurate estimation because they ignored spatial variations in elasticity, leading to only a rough approximation of the properties of Earth's crust. In this paper, a semi-analytical procedure for calculating the piezomagnetic field arising from a point dislocation source embedded in a layered elastic medium is derived. Following a well-established method of the vector surface harmonic expansion, all of the governing equations written in partial differential equations in a real domain, together with the linear constitutive law of the piezomagnetic effect, are converted to a set of ordinary differential equations in a wavenumber domain. Equations in the wavenumber domain are solved analytically, and each component of the piezomagnetic field in the real domain is obtained after applying the Hankel transform. By using the derived procedure, the piezomagnetic and displacement fields due to a finite fault with strike-slip, dip-slip, and tensile-opening mechanisms are estimated for media with layered elasticity structures. Results for a finite fault are obtained by integrating the point source solution over the fault plane. The results of the numerical analysis allow the effect of heterogeneities in rigidity on the piezomagnetic effect to be examined and implications for the findings of previous investigations to be drawn. In cases where the moment-release at the dislocation source is fixed, the effect of the rigidity differences between upper and lower layers on the piezomagnetic field is minor even in the case where the Curie point depth is near the source of dislocation. This result is in contrast to a previous study that assumed the Mogi model and suggested that heterogeneities in the horizontal direction may be of importance when combined with layered rigidity structures. A contrast is seen between the piezomagnetic and displacement fields corresponding to models with layered rigidity structures: the piezomagnetic field is roughly proportional to the moment-release on a source fault, whereas the displacement field is proportional to slip or opening of the fault. Provided that the rigidity of the crust increases with increasing depth, the calculated piezomagnetic field is likely to have been underestimated in many earlier studies, which assumed uniform rigidity and a geodetically inverted size of slip.