The magma plumbing systems of volcanoes in subduction and divergent tectonic settings are relatively well known, whereas those of intraplate volcanoes remain elusive; robust geophysical information on the magma pathways and storage zones is lacking. We inverted magnetotelluric data to image the magma plumbing system of an intraplate monogenetic volcanic field located above the stagnant Pacific slab in northeast China. We identified a complex, vertically aligned, low-resistivity anomaly system extending from the asthenosphere to the surface consisting of reservoirs with finger- to lens-like geometries. We show that magma forms as CO₂-rich melts in a 150-km-deep asthenospheric plume crossing the whole lithosphere as hydrated melt, inducing underplating at 50 km depth, evolving in crustal reservoirs, and erupting along dikes. Intraplate volcanoes are characterized by low degrees of melting and low magma supply rates. Their plumbing systems have a geometry not so different from that of volcanoes in subduction settings.

SUMMARY: Constraining the distribution of water in different regions of the mantle remains one of the significant challenges to comprehend the global deep water cycle. Geomagnetic depth soundings can provide such constraint through the electrical conductivity structure. Hence, this study aims to propose a regularization technique that can estimate previously unavailable C-response. In the method, the objective function comprised an L1-norm measured data prediction error and a spectral smoothness constraint term. We used the data error of C-response to weight the predicted error. The L-BFGS method was introduced to determine the minimum point of the objective function, and the regularization parameter decreased adaptively during inversion. Thus, the geomagnetic data processed yielded high-quality C-responses in 31 stations in Eastern China. In addition, we obtained 1-D electrical conductivity profiles in the mantle transition zone (MTZ) beneath Eastern China from C-responses using the L-BFGS method. Compared with the global 1-D model, the conductivity–depth profiles revealed that the MTZ beneath Eastern China is more conductive in the east but more resistive in the west. The conversion of these conductivities to water content based on the mineral physics suggested that the MTZ beneath Eastern China is characterized by a high water concentration, approximately 0.2 and 1 wt per cent in the upper and lower MTZ, respectively. Owing to the inclusion of more stations, the water-rich region could be constrained roughly to the east of the North–South Gravity Lineament (NSGL). Further considering seismic images in the same area, this water content distribution pattern suggested that the front of the stagnant Pacific Plate in the lower MTZ might have reached the NSGL. However, the dehydration reactions in the stagnant slab were more active in the eastern part. Perhaps, some of these fluids migrated into the upper MTZ and could be the source of the trapped water found in the xenoliths from the deep upper mantle beneath Eastern China.

A model of electrical conductivity in the mid–mantle transition zone was obtained with improved constraints. An L₁‐norm regularization inversion algorithm is proposed here that reduces the influence of noisy data on three‐dimensional geomagnetic depth sounding inversion from C‐responses. Here the regularization is implemented by using an L₁‐norm to measure the predicted data error, which is normalized by the C‐response covariance, but an L₂‐norm is used to measure the regularization term associated with model parameters. The limited‐memory quasi‐Newton method (L‐BFGS) is used to invert for the three‐dimensional electrical conductivity model. The model is discretized by curved rectangular prisms in spherical coordinates. Sensitivity tests show that for good‐quality data contaminated by Gaussian noise, L₁ inversion, which could perform as well as L₂ inversion, can adequately recover the main features of the electrical conductivity structure within the region of data coverage. When data errors are drawn from an exponential distribution, L₁ inversion obtains relatively reliable reconstruction of the electrical structure, even when the noise level is comparable to that of actual C‐responses. C‐responses from 129 low‐latitude and midlatitude geomagnetic observatories are inverted using L₁‐norm minimization of the data error. The resulting model reveals an electrically conductive feature in the lower mantle transition zone and upper lower mantle that is broadly coincident with that found in previous studies. The reduced influence of data with large variances on L₁‐norm misfits, along with inclusion of responses estimated from more observatories, makes L₁ inversion more clearly identify these deep conductive features while identifying previously obscured anoconductive zones. A feature of particular interest is the high electrical conductivity anomaly beneath the Bermuda‐Sargasso Sea region in the mid–mantle transition zone and the uppermost lower mantle. Rock physics analysis indicates that the anomaly is most possibly caused by the wet upwelling material with excessive ~650 K higher temperature, suggesting a narrow tail with a broad head.