P. Schmoldt, PhD - MTNet - DIAS

P. Schmoldt, PhD - MTNet - DIAS P. Schmoldt, PhD - MTNet - DIAS

mtnet.dias.ie
from mtnet.dias.ie More from this publisher
04.08.2013 Views

9. Data collection and processing 9.6. Compensating for distortion of the impedance tensor Removing distortion of the data, as described in Section 4, is a crucial, but difficult, step during the analysis for most MT datasets. Moreover, since previous geological and geophysical studies indicate a complex subsurface setting for the region under investigation by the PICASSO Phase I project (cf. Sec. 7), this step is carried out with great diligence. 9.6.1. Geoelectric strike estimation Determining whether a regional 2D approximation is justified, and (if that is found true) subsequent identification of the geoelectric 2D strike direction are major elements of MT impedance tensor analysis as it affects the decomposition of EM fields into TM and TE mode contribution. The approach used here for the PICASSO Phase I data utilises a script developed by Jan Vozar, displaying the RMS misfit calculated by the strike algorithm for different directions at all stations. The strike algorithm is based on the impedance tensor decomposition method by McNeice and Jones [2001], using the Groom and Bailey [1989] technique (Sec. 4.4.4), and is today commonly used my many in the MT community. The RMS misfit is therein calculated for a chosen depth range derived by the Niblett-Bostick depth approximation for the rotational invariant arithmetic average of the off-diagonal elements (also referred to as Berdichevsky average) [Berdichevsky and Dmitriev, 1976b] (cf. Sec. 6.3.1). A relative impedance error floor of 3.5% is used during the application of the strike algorithm, resulting in an error floor of 2.0 degrees for the phase and 7.12% for the apparent resistivity. The range of directions under investigation is limited to the interval 0 to 90 degrees, owing to the 90 degrees ambiguity of the geoelectric strike estimation. The ambiguity originates from the fact that during the calculation no assumption can be made about the orientation of the two modes, TE and TM, except that they are orthogonal. Hence the true geoelectric strike direction and its orthogonal fit the data equally well. For example, a calculated strike direction of N50E (50 degrees clockwise from North) indicates that the true strike of the geological features has a direction of either N50E or N140E. The decision about the direction has than to be made by considering additional sources of information, such as results of geological studies or the vertical field responses. For the strike analysis of the PICASSO Phase I data an increment of one degree is chosen in order to provide sufficient resolution. The geoelectric strike analysis using the above-described RMS misfit calculation is supplemented by a multisite, multifrequency (transformed into Niblett-Bostick depth) analysis with the same strike program. Therein, optimal strike direction and average RMS misfit are calculated for the chosen depth and stations ranges of; results are from here on referred to as multi-strike. These multi-strike directions are calculated separately for the whole profile and its northern region only, reasoned by the observed difference 210

0 km 10 km 30 km 100 km 300 km Depth range Depth: 0 – 300 km 90 75 60 45 30 15 0 pic041 pic040 pic037 pic035 pic033 pic031 Betic Cordillera 9.6. Compensating for distortion of the impedance tensor pic029 pic027 pic025 pic023 CIZ pic022 pic020 pic019 pic017 pic015 Multi-strike (Av. RMS-misfit in brackets) pic013 pic011 pic009 pic007 Tajo Basin pic006 All stations Tajo Basin 51.5 (1.6) 51.4 (1.5) Fig. 9.6.: RMS misfit for different geoelectric strike directions of the stations recorded during PICASSO Phase I, using data in the Niblett-Bostick depth range 0 – 300 km; shown in the inset. Also shown in the top right corner is the optimal common geoelectric strike direction clockwise from North calculated for all stations and for stations in the Tajo Basin only; indicated by the dashed white line in the main plot. Assignment of stations to the different geological regions, displayed on the bottom of the plot, is based on their location in respect to geological units of the USGS EnVision map for Europe (Fig. 9.1), CIZ: Central Iberian Zone, part of the Iberian Massif. See text for further details. in geoelectric strike characteristics described in the following paragraphs. For a 2D subsurface, satisfactory multi-strike directions should coincide with low RMS misfits for the same direction at a number of adjacent stations, bounded by higher RMS misfits for other directions. Cases where the RMS misfit is low for most directions, on the other hand, indicate a rather 1D nature of the subsurface. In the following paragraphs geoelectric strike direction is determined for different depth ranges in order to investigate whether the best-fitting direction varies for individual regions in crust and mantle. Complete depth range Initially, the RMS misfit for different directions is calculated for the whole depth range spanning from 0 to 300 km, revealing an unsystematic distribution of the RMS misfit (cf. Fig. 9.6). Low RMS misfits at stations pic013 – pic017 and pic023 – pic025 are due to the truncation of their long-period data affected by EM noise sources, as described in Section 9.4. The low RMS misfit of station pic009, on the other hand, observed throughout the whole frequency range, is due a higher degree of uncertainty of the responses for this station. Multi-strike analysis yields optimal strike direction and RMS misfit (in brackets) of N51.5E (1.6) and N51.4E (1.5) for the whole profile and the Tajo Basin, respectively; which however, is not significant given the observed unsystematic RMS misfit distribution along the profile. Crust and mantle depth As a second step, the depth range is divided into crust and mantle regions spanning from 0 km to 30 km and 35 km to 200 km, respectively (Fig. 9.7); therein, Moho depths of 32 km derived by seismic reflection data (cf. Sec. 7.3.2) are used as crustal-thickness proxy with a small error margin. The crustal depth range exhibits an unsystematic RMS misfit distribution and a multi-strike direction of N58.0E (1.3 av. RMS misfit) and N51.2E (1.5 av. RMS misfit) for the whole profile and the Tajo Basin, respectively. pic005 pic004 pic003 pic002 pic001 RMS misfit 3 2 1 0 211

0 km<br />

10 km<br />

30 km<br />

100 km<br />

300 km<br />

Depth range<br />

Depth: 0 – 300 km<br />

90<br />

75<br />

60<br />

45<br />

30<br />

15<br />

0<br />

pic041<br />

pic040<br />

pic037<br />

pic035<br />

pic033<br />

pic031<br />

Betic Cordillera<br />

9.6. Compensating for distortion of the impedance tensor<br />

pic029<br />

pic027<br />

pic025<br />

pic023<br />

CIZ<br />

pic022<br />

pic020<br />

pic019<br />

pic017<br />

pic015<br />

Multi-strike<br />

(Av. RMS-misfit in brackets)<br />

pic013<br />

pic011<br />

pic009<br />

pic007<br />

Tajo Basin<br />

pic006<br />

All stations Tajo Basin<br />

51.5 (1.6) 51.4 (1.5)<br />

Fig. 9.6.: RMS misfit for different geoelectric strike directions of the stations recorded during PICASSO Phase I, using data in the<br />

Niblett-Bostick depth range 0 – 300 km; shown in the inset. Also shown in the top right corner is the optimal common geoelectric<br />

strike direction clockwise from North calculated for all stations and for stations in the Tajo Basin only; indicated by the dashed white<br />

line in the main plot. Assignment of stations to the different geological regions, displayed on the bottom of the plot, is based on their<br />

location in respect to geological units of the USGS EnVision map for Europe (Fig. 9.1), CIZ: Central Iberian Zone, part of the Iberian<br />

Massif. See text for further details.<br />

in geoelectric strike characteristics described in the following paragraphs. For a 2D subsurface,<br />

satisfactory multi-strike directions should coincide with low RMS misfits for the<br />

same direction at a number of adjacent stations, bounded by higher RMS misfits for other<br />

directions. Cases where the RMS misfit is low for most directions, on the other hand,<br />

indicate a rather 1D nature of the subsurface. In the following paragraphs geoelectric<br />

strike direction is determined for different depth ranges in order to investigate whether the<br />

best-fitting direction varies for individual regions in crust and mantle.<br />

Complete depth range<br />

Initially, the RMS misfit for different directions is calculated for the whole depth range<br />

spanning from 0 to 300 km, revealing an unsystematic distribution of the RMS misfit (cf.<br />

Fig. 9.6). Low RMS misfits at stations pic013 – pic017 and pic023 – pic025 are due to the<br />

truncation of their long-period data affected by EM noise sources, as described in Section<br />

9.4. The low RMS misfit of station pic009, on the other hand, observed throughout the<br />

whole frequency range, is due a higher degree of uncertainty of the responses for this<br />

station. Multi-strike analysis yields optimal strike direction and RMS misfit (in brackets)<br />

of N51.5E (1.6) and N51.4E (1.5) for the whole profile and the Tajo Basin, respectively;<br />

which however, is not significant given the observed unsystematic RMS misfit distribution<br />

along the profile.<br />

Crust and mantle depth<br />

As a second step, the depth range is divided into crust and mantle regions spanning from<br />

0 km to 30 km and 35 km to 200 km, respectively (Fig. 9.7); therein, Moho depths of<br />

32 km derived by seismic reflection data (cf. Sec. 7.3.2) are used as crustal-thickness<br />

proxy with a small error margin. The crustal depth range exhibits an unsystematic RMS<br />

misfit distribution and a multi-strike direction of N58.0E (1.3 av. RMS misfit) and N51.2E<br />

(1.5 av. RMS misfit) for the whole profile and the Tajo Basin, respectively.<br />

pic005<br />

pic004<br />

pic003<br />

pic002<br />

pic001<br />

RMS<br />

misfit<br />

3<br />

2<br />

1<br />

0<br />

211

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!