3D inversion based on multi-grid approach of magnetotelluric data from Northern Scandinavia
Abstract
In this work we investigate the geoelectrical structure of the cratonic margin of Fennoscandian Shield by means of magnetotelluric (MT) measurements carried out in Northern Norway and Sweden during summer 2011-2012. The project Magnetotellurics in the Scandes (MaSca) focuses on the investigation of the crust, upper mantle and lithospheric structure in a transition zone from a stable Precambrian cratonic interior to a passive continental margin beneath the Caledonian Orogen and the Scandes Mountains in western Fennoscandia. Recent MT profiles in the central and southern Scandes indicated a large contrast in resistivity between Caledonides and Precambrian basement. The alum shales as a highly conductive layers between the resistive Precambrian basement and the overlying Caledonian nappes are revealed from this profiles. Additional measurements in the Northern Scandes were required. All together data from 60 synchronous long period (LMT) and about 200 broad band (BMT) sites were acquired. The array stretches from Lofoten and Bodo (Norway) in the west to Kiruna and Skeleftea (Sweden) in the east covering an area of 500x500 square kilometers. LMT sites were occupied for about two months, while most of the BMT sites were measured during one day. We have used new multi-grid approach for 3D electromagnetic (EM) inversion and modelling. Our approach is based on the OcTree discretization where the spatial domain is represented by rectangular cells, each of which might be subdivided (recursively) into eight sub-cells. In this simplified implementation the grid is refined only in the horizontal direction, uniformly in each vertical layer. Using multi-grid we manage to have a high grid resolution near the surface (for instance, to tackle with galvanic distortions) and lower resolution at greater depth as the EM fields decay in the Earth according to the diffusion equation. We also have a benefit in computational costs as number of unknowns decrease. The multi-grid forward solver is implemented within the framework of the modular system for EM inversion (ModEM by G. Egbert, A. Kelbert, N. Meqbel), using the ModEM 3D finite difference staggered grid forward solver (second order PDE in the electric field, with divergence correction) as a starting point for our development. The first 3D inversion model for the crust and upper mantle shows the highly conducting bodies in the crust which can be interpreted as alum shales. The eastern and central parts are presented by resistive Precambrian rocks of the Svecofennian and Archaean domains. The upper mantle is resistive and relates to the Baltica basement. We also compare 3D inversion model with the results of 2D inversion along several profiles. We are able to explain some of the features in the data (out of quadrant phase) with 3D model, thus providing more reliable results compared to routine 2D approach.
- Publication:
-
AGU Fall Meeting Abstracts
- Pub Date:
- December 2012
- Bibcode:
- 2012AGUFMGP33A1092C
- Keywords:
-
- 1515 GEOMAGNETISM AND PALEOMAGNETISM / Geomagnetic induction;
- 3260 MATHEMATICAL GEOPHYSICS / Inverse theory